hist.duration {cfda}R Documentation

Plot the duration

Description

Plot the duration

Usage

## S3 method for class 'duration'
hist(x, breaks = NULL, ...)

Arguments

x

output of compute_duration function

breaks

number of breaks. If not given, use the Sturges rule

...

parameters for geom_histogram

Value

a ggplot object that can be modified using ggplot2 package.

Author(s)

Quentin Grimonprez

Examples

# Simulate the Jukes-Cantor model of nucleotide replacement 
K <- 4
PJK <- matrix(1/3, nrow = K, ncol = K) - diag(rep(1/3, K))
lambda_PJK <- c(1, 1, 1, 1)
d_JK <- generate_Markov(n = 10, K = K, P = PJK, lambda = lambda_PJK, Tmax = 10)


# compute duration of each individual
duration <- compute_duration(d_JK)

hist(duration)

# modify the plot using ggplot2
library(ggplot2)
hist(duration) + 
   labs(title = "Distribution of the duration")


[Package cfda version 0.9.9 Index]