critical {seq2R} | R Documentation |
Function that maximizes or minimizes the first derivative of the model obtained with find.points
function. Also, it is included their 95% confidence intervals.
critical(model, base.pairs = NULL)
model |
|
base.pairs |
Character string for A vs. T or C vs G. |
The returned list has two component ($AT
, $CG
). Both of them containing a matrix with values about their critical points (maxima and minima), lower and upper 95% confidence intervals.
AT |
Critical points for AT. |
CG |
Critical points for CG. |
Nora M. Villanueva and Marta Sestelo.
N. M. Villanueva, M. Sestelo, M. M. Fonseca and J. Roca-Pardinas (2023). seq2R: An R package to detect change points in DNA sequences. Mathematics, 11 (10), 2299.
library(seq2R)
#mtDNAhum <- read.genbank("NC_012920")
data(mtDNAhum)
DNA <- transform(mtDNAhum)
seq1 <- find.points(DNA, nboot = 10)
critical(seq1,base.pairs="CG")
critical(seq1,base.pairs="AT")