# 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)
pca.obj = getPCAFunc(es = es.sim,
labelVariable = "subjID",
hybName = "memSubj",
requireLog2 = FALSE,
corFlag = FALSE
)
pca2DPlot(pcaObj = pca.obj,
plot.dim = c(1,2),
labelVariable = "subjID",
hybName = "memSubj",
plotOutPutFlag = FALSE,
cex.legend = 0.5,
legendPosition = "topright")
Run the code above in your browser using DataLab