CNCDriver combined mutation recurrence and functional impact to identify coding and non-coding cancer drivers
User will need to install devtools in R for running CNCDriver package
library("remotes")
remotes::install_github("khuranalab/CNCDriver", ref="master", build_vignette=TRUE)
library(CNCDriver)
#####
# global parameters setup
#####
funseq2OutFile<-"path/to/funseq2/annotatedfile"
replicationTimingCutOff<-0.2
filterOutBlacklistMutations<-TRUE
seedNum<-42
reSampleIterations<-10000
reRunPvalueCutOff<-0.1
minPoints<-2
dRadius<-50
useCores<-4
debugMode<-FALSE
taskNum<-0
unitSize<-100
#####
mutationType<-"cds_cluster"
elementKeyWord<-"CDS"
cdsOutputDf<-getCDSPvalueWithPreFilter2(inputFileDir,outputFileDir,
codingRegionBedFile,elementKeyWord,
proteinDomainFile,proteinLengthFile,
minPoints,dRadius,
triNucleotideDistributionFile,
filterOutBlacklistMutations,
mutationBlacklistFile,
replicationTimingGenomeBinnedFile,
replicationTimingElementBinnedFileCDS,
tumorType,mutationType,cellType,
replicationTimingCutOff,
seedNum,reSampleIterations,
reRunPvalueCutOff,
useCores,taskNum,unitSize,debugMode)
makeQQplot(cdsOutputDf)
#####
mutationType<-"promoter_cluster"
elementKeyWord<-"Promoter"
promoterOutputDf<-getPromoterPvalueWithPreFilter2(inputFileDir,outputFileDir,
promoterRegionBedFile,elementKeyWord,
minPoints,dRadius,
triNucleotideDistributionFile,
filterOutBlacklistMutations,
mutationBlacklistFile,
replicationTimingGenomeBinnedFile,
replicationTimingElementBinnedFilePromoter,
tumorType,mutationType,cellType,
replicationTimingCutOff,
seedNum,reSampleIterations,
reRunPvalueCutOff,
useCores,taskNum,unitSize,debugMode)
makeQQplot(promoterOutputDf)
#####
mutationType<-"lincRNA_cluster"
elementKeyWord<-"lincRNA"
lincRNAOutputDf<-getLincRNAPvalueWithPreFilter2(inputFileDir,outputFileDir,
lincRNARegionBedFile,elementKeyWord,
minPoints,dRadius,
triNucleotideDistributionFile,
filterOutBlacklistMutations,
mutationBlacklistFile,
replicationTimingGenomeBinnedFile,
replicationTimingElementBinnedFileLincRNA,
tumorType,mutationType,cellType,
replicationTimingCutOff,
seedNum,reSampleIterations,
reRunPvalueCutOff,
useCores,taskNum,unitSize,debugMode)
makeQQplot(lincRNAOutputDf)
#####
mutationType<-"enhancerUnit_cluster"
elementKeyWord<-"Distal"
enhancerUnitOutputDf<-getEnhancerUnitPvalueWithPreFilter2(inputFileDir,
outputFileDir,
enhancerRegionBedFile,
elementKeyWord,
minPoints,dRadius,
triNucleotideDistributionFile,
filterOutBlacklistMutations,
mutationBlacklistFile,
replicationTimingGenomeBinnedFile,
replicationTimingElementBinnedFileEnhancer,
tumorType,mutationType,cellType,
replicationTimingCutOff,
seedNum,reSampleIterations,
reRunPvalueCutOff,
useCores,taskNum,unitSize,
debugMode)
makeQQplot(enhancerUnitOutputDf)
#####
For any questions, comments and suggestions, please email
Copyright © 2016-2019 Ekta Khurana Lab, WCMC
This project is licensed under the
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.