critical {seq2R}R Documentation

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)

Arguments

model

change.points object.

base.pairs

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

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.

Author(s)

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

[Package seq2R version 2.0.1 Index]