### R code from vignette source 'ChAMP.Rnw' ################################################### ### code chunk number 1: InstallLibraries (eval = FALSE) ################################################### ## source("http://bioconductor.org/biocLite.R") ## biocLite(c('minfi', 'DNAcopy', 'impute', 'marray', 'limma', 'preprocessCore', ## 'RPMM', 'sva', 'IlluminaHumanMethylation450kmanifest','wateRmelon')) ################################################### ### code chunk number 2: loadChAMPLibrary ################################################### library(ChAMP) ################################################### ### code chunk number 3: loadTest ################################################### testDir=system.file("extdata",package="ChAMPdata") ################################################### ### code chunk number 4: loadTest2 ################################################### data(testDataSet) myLoad=testDataSet ################################################### ### code chunk number 5: processFunctionIDAT (eval = FALSE) ################################################### ## champ.process(directory = testDir) ################################################### ### code chunk number 6: processSAVE (eval = FALSE) ################################################### ## save(myLoad,file="currentStudyloadedData.RData") ## load("currentStudyloadedData.RData") ## champ.process(fromIDAT=FALSE) ################################################### ### code chunk number 7: processSTEPbySTEP (eval = FALSE) ################################################### ## myLoad=champ.load(directory = testDir) ## myNorm=champ.norm() ## champ.SVD() ## batchNorm=champ.runCombat() ## limma=champ.MVP() ## lasso=champ.lasso(fromFile =TRUE, limma=limma) ## champ.CNA() ################################################### ### code chunk number 8: SampleSheet ################################################### myLoad$pd ################################################### ### code chunk number 9: combatFunction (eval = FALSE) ################################################### ## myLoad=champ.load(directory = testDir, filterBeads=TRUE) ################################################### ### code chunk number 10: loadFunction ################################################### myLoad = champ.load(directory=testDir) ################################################### ### code chunk number 11: normFunction ################################################### myNorm=champ.norm() ################################################### ### code chunk number 12: svdFunction ################################################### champ.SVD() ################################################### ### code chunk number 13: combatFunction (eval = FALSE) ################################################### ## batchNorm=champ.runCombat() ################################################### ### code chunk number 14: mvpFunction (eval = FALSE) ################################################### ## limma=champ.MVP() ## head(limma) ################################################### ### code chunk number 15: lassoFunction (eval = FALSE) ################################################### ## lasso=champ.lasso(fromFile=TRUE, limma=limma, image=FALSE,bedFile=FALSE) ## if(!is.null(lasso)) ## { ## head(lasso) ## } ################################################### ### code chunk number 16: cnaFunction (eval = FALSE) ################################################### ## CNA=champ.CNA()