Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
Refined support for bubble plots
  • Loading branch information
jenzopr committed Jul 9, 2015
1 parent 03bc4c5 commit e5f7991
Show file tree
Hide file tree
Showing 2 changed files with 10 additions and 10 deletions.
4 changes: 2 additions & 2 deletions src/admire
Expand Up @@ -289,7 +289,7 @@ base = file.path("./")
results = file.path("results/")
# WORKING ENVIRONMENT BASICS
sampleInfo = read.450k.sheet(base = base)
sampleInfo = read.450k.sheet(base = base, pattern="SampleSheet.csv")
# PARAMETERS
normalization="$N"
Expand Down Expand Up @@ -417,7 +417,7 @@ EXCEL=excel
mkdir -p $EXCEL
set -- ${groups[@]}
cat > additional_images.txt << EOL
file group region chr
file group region chr
EOL
for a; do
shift
Expand Down
16 changes: 8 additions & 8 deletions src/renderAdditionalImages.R
@@ -1,7 +1,7 @@
library("RColorBrewer")
#source(".Rprofile")
#source("config.R")
#source("meth_plot.R")
source(".Rprofile")
source("config.R")
source("meth_plot.R")

driver = read.table("additional_images.txt",header=T,sep="\t", stringsAsFactors = F)
intervals = data.frame(c(0,.1,.2,.3,.4,.5,.6,.7,.8,.9),c(.1,.2,.3,.4,.5,.6,.7,.8,.9,1))
Expand All @@ -21,12 +21,12 @@ for(i in 1:nrow(driver)) {
betaValues = t(betaValues)
rownames(betaValues) = sampleInfo[groupSamples,"Sample_Group"]

png(file=paste(dirname(driver[i,"file"]),"/",driver[i,"chr"],"-",min(probePositions),"-",max(probePositions),".png",sep=""))
pdf(file=paste(dirname(driver[i,"file"]),"/",driver[i,"chr"],"-",min(probePositions),"-",max(probePositions),".pdf",sep=""))
meth_plot(betaValues,intervals,"by-meth","proportional",probePositions,row.names(betaValues),probeNames)
dev.off()
png(file=paste(dirname(driver[i,"file"]),"/hm-",driver[i,"chr"],"-",min(probePositions),"-",max(probePositions),".png",sep=""))
meth_plot(betaValues,intervals,"by-meth","grid",probePositions,row.names(betaValues),probeNames)
dev.off()
#unlink(paste(driver[i,"file"],"txt",sep="."))
#pdf(file=paste(dirname(driver[i,"file"]),"/hm-",driver[i,"chr"],"-",min(probePositions),"-",max(probePositions),".pdf",sep=""))
#meth_plot(betaValues,intervals,"by-meth","grid",probePositions,row.names(betaValues),probeNames)
#dev.off()
unlink(paste(driver[i,"file"],"txt",sep="."))
}

0 comments on commit e5f7991

Please sign in to comment.