PoweR (version 1.0.7)

plot.sizepower: size-power curves.

Description

This function produces a size-power curves plot.

Usage

# S3 method for sizepower
plot(x, xnull,legend.pos=NULL,...)

Arguments

x

Fx object as returned by function calcFx.

xnull

Fx object as returned by function calcFx, but computed under the null.

legend.pos

If NULL, position of the legend will be computed automatically. Otherwise, it should be either a character vector in "bottomright", "bottom", "bottomleft", "left", "topleft", "top", "topright", "right" and "center". Or a numeric vector of length 2 giving the x-y coordinates of the legend.

...

further arguments passed to the plot or points functions.

Details

See Section 2.3 in Lafaye de Micheaux, P. and Tran, V. A. (2014).

References

Pierre Lafaye de Micheaux, Viet Anh Tran (2016). PoweR: A Reproducible Research Tool to Ease Monte Carlo Power Simulation Studies for Studies for Goodness-of-fit Tests in R. Journal of Statistical Software, 69(3), 1--42. doi:10.18637/jss.v069.i03

See Also

See plot.pvalue, plot.discrepancy, graph.

Examples

# NOT RUN {
## You can increase M for better results:
stind <- c(43,44,42)   # Indices of test statistics.
alter <-list(stat43=3,stat44=3,stat42=3) # Type for each test.
# Several p-values computed under the null.
pnull <- many.pval(stat.indices=stind,law.index=1,
                   n=100,M=100,N=10,alter=alter,null.dist=2,
                   method="MC")$pvals
Fxnull <- calcFx(pnull)
p <- many.pval(stat.indices=stind,law.index=4,n=100,
               M=100,N=10,alter=alter,null.dist=2,
               method="MC")$pvals
Fx <- calcFx(p)
plot.sizepower(Fx,Fxnull)
# }