plot.scaleboot {scaleboot}R Documentation

Plot Diagnostics for Multiscale Bootstrap

Description

plot method for class "scaleboot".

Usage


## 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,...)

Arguments

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.

Details

plot method plots bootstrap probabilities and call lines method, which draws fitted curves for models.

Author(s)

Hidetoshi Shimodaira

See Also

sbfit,legend.

Examples

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)


[Package scaleboot version 0.1-2 Index]