Quantile-quantile information about Haar-NN and Anscombe samples

Description

A function to generate Q-Q plots (from simulations) for the Anscombe and (binomial) Haar-Fisz transforms.

Usage

1
2
qqstuff(intensity, binsize = 4, paths = 100, respaths = 1000, plot.q = FALSE, 
plot.sq = FALSE)

Arguments

intensity

an Bernoulli intensity vector, e.g. pintens.

binsize

a binomial size to generate a binomial mean vector.

paths

the number of paths sampled from the mean vector to use in Q-Q calculations.

respaths

the number of residual paths to use in squared residual calculations.

plot.q

A boolean variable, indicating whether simulation Q-Q plots should be outputted or not.

plot.sq

A boolean variable, indicating whether simulation squared residual plots should be outputted or not.

Details

respaths paths are sampled from the mean intensity vector. From these, the first paths are used to generate Q-Q data, which are then averaged for the Q-Q plots. The original paths are used to calculate a squared residual vector corresponding to the mean intensity vector.

Value

qqinfo. A 8 component list of quantile and residual plot information.

vmat

A matrix of dimensions respathsxlength(intensity), each row being a path from the intensity vector.

Av

A matrix of dimensions respathsxlength(intensity), each row an Anscombe-transformed path.

bfv

A matrix of dimensions respathsxlength(intensity), each row a binomial Haar-Fisz-transformed path.

vminusl

A matrix of the difference between the paths and the mean intensity.

vminusl

A matrix of the difference between the Anscombe-transformed paths and the mean intensity.

vminusl

A matrix of the difference between the binomial Haar-Fisz-transformed paths and the mean intensity.

Asqres

vector of squared residuals of Anscombe-transformed paths.

bfsqres

vector of squared residuals of binomial Haar-Fisz-transformed paths.

Note

This function requires the package wavethresh. N.B. Since this function returns a lot of information, assign the output to a variable, to avoid printing endless information in the console.

Author(s)

Matt Nunes (m.nunes@ucl.ac.uk)

See Also

qqnormy

Examples

1
2
3
4
5
data(pintens)

a<-qqstuff(intensity=pintens,binsize=4,paths=100,respaths=100,plot.q=TRUE,plot.sq=TRUE)

#plots some interesting graphs.