| hot.comp {hot} | R Documentation |
Computes Hotelling T^2 statistics
hot.comp(expr.dat, status, set.size, nset)
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) |
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.
Gilles Guillot
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
## 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)