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)