statgen: Statistics generator

Description

This function generates useful simulation statistics for NN and Anscombe transforms.

Usage

statgen(valuelist, xgrid = seq(0, 1, length = 21), ygrid = seq(0, 1, length = 21), 
binsize = 32, plot.m = FALSE, plot.v = FALSE, plot.ks = FALSE, ptype = "persp")

Arguments

valuelist

a two component list as produced by afgen.

xgrid

a vector of x coordinate binomial proportions.

ygrid

a vector of x coordinate binomial proportions.

binsize

binomial size to use in simulations.

plot.m

A boolean variable, indicating whether mean simulation plots should be outputted.

plot.v

A boolean variable, indicating whether variance simulation plots should be outputted.

plot.ks

A boolean variable, indicating whether Kolmogorov-Smirnov simulation plots should be outputted.

ptype

where appropriate, the type of plots to be produced. Possible values are "persp" for 3D persective plots or "contour" for corresponding contour plots.

Value

afm

matrix of sample mean values for binomial Haar-Fisz samples.

anm

matrix of sample mean values for Anscombe samples.

afv

matrix of sample variance values for binomial Haar-Fisz samples.

anv

matrix of sample variance values for Anscombe samples.

afk

matrix of Kolmogorov-Smirnof statistics for binomial Haar-Fisz samples.

ank

matrix of Kolmogorov-Smirnof statistics for Anscombe samples.

Details

The function does several sample variance plots, Kolmogorov-Smirnov and mean plots for the data in the variable valuelist (for both Anscombe and binomial Haar-Fisz transforms).

Examples

Run this code

# NOT RUN {
a<-afgen(xgrid = seq(0, 1, length = 21), ygrid = seq(0, 1, length = 21), 
samples = 1000, binsize = 32)

b<-statgen(a,xgrid=seq(0,1,length=21),ygrid=seq(0,1,length=21),binsize=32,plot.m=FALSE,
plot.v=TRUE,plot.ks=FALSE,ptype="persp")

# }