hot.comp {hot}R Documentation

Hotelling T2

Description

Computes Hotelling T^2 statistics

Usage

hot.comp(expr.dat, status, set.size, nset)

Arguments

expr.dat Expression data: a matrix with genes in rows and individuals in columns
status Label of individuals(0 or 1) with length the number of column in expr.dat
set.size Number of genes used to discriminate between groups: an integer (between 2 and 6).
nset Number of set of genes for which the T^2 statistic is stored in the returned matrix: an integer (should be > 1)

Value

A matrix with set.size+2 columns and nset lines. The index of genes are stored in columns 1 to set.size. Column set.size+1 contains the T^2 statistics. Column set.size+2 contains the corresponding nominal p-value.

Author(s)

Gilles Guillot

References

Discrimination and scoring using small sets of genes for two-sample microarray data, G. Guillot, M. Olsson, M. Bensson, M. Rudemo. To appear in Mathematical Biosciences

Examples


## Not run: 
## Assuming expression data are stored in a matrix named expr.dat
## with 10 controls + 10 patients
## the following instruction
res <- hot.comp(expr.dat=expr.dat,
                status=c(rep(0,10),rep(1,10)),
                set.size=4,
                nset=100)
## will return the 100 sets of 4 genes with highest T^2 statistic
## End(Not run)


[Package hot version 0.3 Index]