[BioC] p-values for multtest

Marcus Davy MDavy at hortresearch.co.nz
Wed Oct 1 15:12:42 MEST 2003

heres a function to extract a vector of dfs using the input objects X,
and classlabel 
inputs for mt.teststat.

welchdf <- function(X, classlabel)
    varx1 <- apply(X[,as.logical(classlabel)], 1, var)
    varx2 <- apply(X[,!as.logical(classlabel)], 1, var)
    n <- table(classlabel)
    df <-

It gives similar figures to
t.test(x,y, var.equal=F)$parameter
for any row in X.


>>> Ann Hess <hess at stat.colostate.edu> 1/10/2003 8:41:39 AM >>>
I am using mt.teststat (multtest) to obtain the teststatistics for
t-tests.  I was wondering if there is a way to get the Satterthwaite
df for the tests using multtest.

Any help would be appreciated.

Bioconductor mailing list
Bioconductor at stat.math.ethz.ch 

The contents of this e-mail are privileged and/or confidential to the
named recipient and are not to be used by any other person and/or
organisation. If you have received this e-mail in error, please notify 
the sender and delete all material pertaining to this e-mail.

More information about the Bioconductor mailing list