plot.scaleboot {scaleboot} | R Documentation |
plot
method for class "scaleboot"
.
## S3 method for class 'scaleboot': plot(x, models = names(x$fi), xval = c("inverse","sigma","square"), yval = c("zvalue", "pvalue", "psi"), xlab = NULL, ylab = NULL,log.xy = "", xlim = NULL, ylim = NULL, add = F, length.x = 300, col =1:6, lty = 1:5, lwd = par("lwd"), pch = 1, cex = 1, pt.col = col[1],pt.lwd = lwd[1], ...) ## S3 method for class 'scalebootv': plot(x,models=attr(x,"models"),...) ## S3 method for class 'scaleboot': lines(x,z,models=names(x$fi), length.x=z$length.x, col=z$col,lty=z$lty,lwd=z$lwd,... ) sblegend(x="topright",y=NULL,z,inset=0.1,...)
x |
an object used to select a method.
For sblegend ,
x is a numeric or character such as "lefttop" or "righttop", which
is passed to legend .
|
models |
character vector of model names. If numeric,
names(object$fi)[models] is used for each "scaleboot"
object. |
xval |
specifies x-axis. "inverse" for 1/σ or
1/sqrt(sa[i]) , "sigma" for σ or sqrt(sa[i]) . |
yval |
specifies y-axis. "zvalue" for
psi(σ^2|β)/σ or qnorm(1-bp[i]) ,
"pvalue" for 1-Phi(psi(σ^2|β)/σ) or bp[i] ,
"psi" for psi(σ^2|β) or sqrt(sa[i])*qnorm(1-bp[i]) . |
xlab |
label for x-axis. |
ylab |
label for y-axis. |
log.xy |
character to specify log-scale. "", "x", "y", or "xy". |
xlim |
range for x-axis. |
ylim |
range for y-axis. |
add |
logical for adding another plot. |
length.x |
the number of segments to draw curves. |
col |
color for model curves. |
lty |
lty for model curves. |
lwd |
lwd for model curves. |
pch |
pch for bp points. |
cex |
cex for bp points. |
pt.col |
col for bp points. |
pt.lwd |
lwd for bp points. |
... |
further arguments passed to or from other methods. |
z |
output from previous plot.scaleboot . |
y |
numeric passed to legend . |
inset |
inset distance from the margins, which is passed to
legend . |
plot
method plots bootstrap probabilities and call lines
method, which draws fitted curves for models.
Hidetoshi Shimodaira
data(mam15) ## a single plot a <- mam15.relltest[["t4"]] # an object of class "scaleboot" sblegend("topleft",z=plot(a)) # x=1/sigma, y=zvalue sblegend("topright",z=plot(a,xval="sigma",log="x")) # x=log(sigma) sblegend("bottomright",z=plot(a,xval="sigma", yval="pval",log="xy")) # x=log(sigma), y=log(pvalue) ## multiple plots b <- mam15.relltest[1:15] # an object of class "scalebootv" plot(b) # x=1/sigma, y=zvalue plot(b,xval="sigma",log="x") # x=log(sigma)