plot.change.points {seq2R} | R Documentation |
Useful for drawing the estimation and first derivative of the skew profile.
## S3 method for class 'change.points'
plot(x = model, y = NULL, base.pairs = NULL, der = NULL,
xlab = "x", ylab = "y", col = "black", CIcol = "black", main = NULL, type = "l",
CItype = "l", critical = FALSE, CIcritical = FALSE,ylim=NULL,...)
x |
|
y |
|
base.pairs |
Character string about the skew profile for A vs. T or C vs. G. |
der |
Number which determines inference process to be drawing into the plot. By default |
xlab |
Title for x axis. |
ylab |
Title for y axis. |
col |
A specification for the default plotting color. |
CIcol |
A specification for the default confidence intervals plotting color. |
main |
An overall title for the plot. |
type |
Type of plot should be drawn. Possible types are, |
CItype |
Type of plot should be drawn for confidence intervals. Possible types are, |
critical |
A logical value. If |
CIcritical |
A logical value. If |
ylim |
The y limits of the plot. |
... |
Other options. |
Simply produce a plot.
Nora M. Villanueva and Marta Sestelo.
library(seq2R)
#mtDNAhum <- read.genbank("NC_012920")
data(mtDNAhum)
DNA <- transform(mtDNAhum)
seq1 <- find.points(DNA)
plot(seq1,der=0,base.pairs="CG",CIcritical=TRUE,ylim=c(0.08,0.67))
plot(seq1,der=1,base.pairs="CG",CIcritical=TRUE,ylim=c(-0.0005,0.00045))
abline(h=0)
plot(seq1,critical=TRUE, CIcritical=TRUE)