chisq.stat: SAM Analysis for Categorical Data

Description

Generates the required statistics for a Significance Analysis of Microarrays of categorical data such as SNP data. Should not be called directly, but via sam(..., method = chisq.stat). Replaces cat.stat

Usage

chisq.stat(data, cl, approx = NULL, B = 100, n.split = 1,  check.for.NN = FALSE, lev = NULL, B.more = 0.1,  B.max = 50000, n.subset = 10, rand = NA)

Arguments

data

a matrix, data frame, or list. If a matrix or data frame, then each row must correspond to a variable (e.g., a SNP), and each column to a sample (i.e.\ an observation). If the number of observations is huge it is better to specify data as a list consisting of matrices, where each matrix represents one group and summarizes how many observations in this group show which level at which variable. These matrices can be generated using the function rowTables from the package scrime. For details on how to specify this list, see the examples section on this man page, and the help for rowChisqMultiClass in the package scrime.

a numeric vector of length ncol(data) indicating to which class a sample belongs. Must consist of the integers between 1 and $c$, where $c$ is the number of different groups. Needs only to be specified if data is a matrix or a data frame.

approx

should the null distribution be approximated by a $ChiSquare$-distribution? Currently only available if data is a matrix or data frame. If not specified, approx = FALSE is used, and the null distribution is estimated by employing a permutation method.

the number of permutations used in the estimation of the null distribution, and hence, in the computation of the expected $d$-values.

n.split

number of chunks in which the variables are splitted in the computation of the values of the test statistic. Currently, only available if approx = TRUE and data is a matrix or data frame. By default, the test scores of all variables are calculated simultaneously. If the number of variables or observations is large, setting n.split to a larger value than 1 can help to avoid memory problems.

check.for.NN

if TRUE, it will be checked if any of the genotypes is equal to "NN". Can be very time-consuming when the data set is high-dimensional.

lev

numeric or character vector specifying the codings of the levels of the variables/SNPs. Can only be specified if data is a matrix or a data frame. Must only be specified if the variables are not coded by the integers between 1 and the number of levels. Can also be a list. In this case, each element of this list must be a numeric or character vector specifying the codings, where all elements must have the same length.

B.more

a numeric value. If the number of all possible permutations is smaller than or equal to (1+B.more)*B, full permutation will be done. Otherwise, B permutations are used.

B.max

a numeric value. If the number of all possible permutations is smaller than or equal to B.max, B randomly selected permutations will be used in the computation of the null distribution. Otherwise, B random draws of the group labels are used.

n.subset

a numeric value indicating how many permutations are considered simultaneously when computing the expected $d$-values.

rand

numeric value. If specified, i.e. not NA, the random number generator will be set into a reproducible state.

Value

A list containing statistics required by sam.

Warning

This procedure will only work correctly if all SNPs/variables have the same number of levels/categories. Therefore, it is stopped when the number of levels differ between the variables.

Details

For each SNP (or more general, categorical variable), Pearson's Chi-Square statistic is computed to test if the distribution of the SNP differs between several groups. Since only one null distribution is estimated for all SNPs as proposed in the original SAM procedure of Tusher et al. (2001) all SNPs must have the same number of levels/categories.

References

Schwender, H. (2005). Modifying Microarray Analysis Methods for Categorical Data -- SAM and PAM for SNPs. In Weihs, C. and Gaul, W. (eds.), Classification -- The Ubiquitous Challenge. Springer, Heidelberg, 370-377.

Tusher, V.G., Tibshirani, R., and Chu, G. (2001). Significance analysis of microarrays applied to the ionizing radiation response. PNAS, 98, 5116-5121.