Plot the number of jumps
# S3 method for njump
hist(x, breaks = NULL, ...)Arguments
- x
output of
compute_number_jumpsfunction- 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.
See also
Other Descriptive statistics:
boxplot.timeSpent(),
compute_duration(),
compute_number_jumps(),
compute_time_spent(),
estimate_pt(),
hist.duration(),
plot.pt(),
plotData(),
statetable(),
summary_cfd()
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)
nJump <- compute_number_jumps(d_JK)
hist(nJump)
# modify the plot using ggplot2
library(ggplot2)
hist(nJump, fill = "#984EA3") +
labs(title = "Distribution of the number of jumps")
