Merge pull request #37 from chapmanb/nointernet-seqinfo

hg38: non-internet access and plot function export
gavinha · Oct 19, 2018 · b39a30d · b39a30d
2 parents 8396703 + 840c977
commit b39a30d
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 2 deletions.
diff --git a/NAMESPACE b/NAMESPACE
@@ -4,7 +4,7 @@ export(loadAlleleCounts,loadDefaultParameters,
 		filterData,getPositionOverlap,wigToRangedData,
 		outputTitanResults,outputModelParameters,correctReadDepth,
 		runEMclonalCN, viterbiClonalCN,computeSDbwIndex,
-		plotAllelicRatio,plotClonalFrequency,plotCNlogRByChr,
+		plotAllelicRatio,plotClonalFrequency,plotCNlogRByChr, plotIdiogram.hg38,
 		plotSubcloneProfiles, outputTitanSegments,
 		setGenomeStyle,plotSegmentMedians,
 		loadHaplotypeAlleleCounts,plotHaplotypeFraction, getHaplotypesFromVCF,

diff --git a/scripts/R_scripts/titanCNA.R b/scripts/R_scripts/titanCNA.R
@@ -181,7 +181,12 @@ if (is.null(outplot)){
 dir.create(outplot, showWarnings=verbose)
 outImage <- gsub(".titan.txt", ".RData", outfile)
 
-seqinfo <- Seqinfo(genome=genomeBuild)
+bsg <- paste0("BSgenome.Hsapiens.UCSC.", genomeBuild)
+if (!require(bsg, character.only=TRUE, quietly=TRUE, warn.conflicts=FALSE)) {
+	seqinfo <- Seqinfo(genome=genomeBuild)
+} else {
+	seqinfo <- seqinfo(get(bsg))
+}
 seqlevelsStyle(chrs) <- genomeStyle
 ## exclude chrX if gender==male ##
 if (gender == "male" || gender == "Male" || gender == "MALE"){