seq2R (version 2.0.1)

critical: Critical points (maxima and minima)

Description

Function that maximizes or minimizes the first derivative of the model obtained with find.points function. Also, it is included their 95% confidence intervals.

Usage

critical(model, base.pairs = NULL)

Value

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.

Arguments

model

change.points object.

base.pairs

Character string for A vs. T or C vs G.

Author

Nora M. Villanueva and Marta Sestelo.

References

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.

Examples

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