cfda (version 0.12.1)

hist.njump: Plot the number of jumps

Description

Plot the number of jumps

Usage

# S3 method for njump
hist(x, breaks = NULL, ...)

Value

a ggplot object that can be modified using ggplot2 package.

Arguments

x

output of compute_number_jumps function

breaks

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

...

parameters for geom_histogram

Author

Quentin Grimonprez

See Also

compute_number_jumps

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")