# qqstuff: Quantile-quantile information about Haar-NN and Anscombe... In binhf: Haar-Fisz functions for binomial data

## 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 ([email protected])

`qqnormy`
 ```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. ```