# generate simulated data set from conditional normal distribution
set.seed(1234567)
es.sim = genSimData.BayesNormal(nCpGs = 100,
nCases = 20, nControls = 20,
mu.n = -2, mu.c = 2,
d0 = 20, s02 = 0.64, s02.c = 1.5, testPara = "var",
outlierFlag = FALSE,
eps = 1.0e-3, applier = lapply)
print(es.sim)
# plot trajectories of the first 5 genes
plotCurves(
dat = exprs(es.sim)[1:5,],
curveNames = featureNames(es.sim)[1:5],
plotOutPutFlag=FALSE,
cex = 0.5,
requireLog2 = FALSE)
Run the code above in your browser using DataLab