

Hi Wolfgang,
Could you please elaborate the earlier email .
> teststat <- mt.teststat(compSet, label, test="t", nonpara="n")
> rawp <- 2 * (1-pt(abs(teststat)))
Do i need to include a df .
For my study I have 12 pre and 12 post samples. (n=24)
file 1 0
file 2 1
file 3 0
file 4 1
file 5 0
file 6 1
I used this for classlabel and the data frame .
I am using the function :
Results <- mt.teststat(X,classlabel,test="wilcoxon")
my script
Results <- mt.teststat(X,classlabel,test="wilcoxon")
rawp.test <- 2 * (1-pt(abs(Results),df=23))
procedures = c("Bonferroni", "Hochberg", "BH")
adjusted = mt.rawp2adjp(rawp.pnorm, procedures)
adjusted$adjp[1:10, ]
Could you please advice if the syntax is correct. I am not sure where I can put the factors you mentioned in the earlier mail above.
ctest::t.test
var.equal=TRUE.
Thank you
smeeta