# NOT RUN {
# generate simulated data set from t distribution
set.seed(1234567)
es.sim = genSimData.tDistr(nCpGs = 100, nCases = 20, nControls = 20,
df0 = 10, ncp0 = 0, df1 = 6, ncp1 = 2.393, testPara = "var",
eps = 1.0e-3, applier = lapply)
print(es.sim)
print(exprs(es.sim)[1:2,1:3])
# do AW score test for the first probe
dat = exprs(es.sim)
pDat = pData(es.sim)
print(pDat[1:2,])
res = FTest(value = dat[1,], group = pDat$memSubj)
print(names(res))
print(res)
# }
Run the code above in your browser using DataLab