Entering edit mode
9.6 years ago
ifudontmind_plzz
▴
200
I was trying to run CrossMap for a set of bam files using this script. But its showing error.
bam_list <- list.files(path = "home_hg19_Bam/", pattern = ".bam$", full.names=T)
for (i in 1:length(bam_list)){
system(paste0("CrossMap.py bam ./hg18ToHg19.over.chain.gz", bam_list[i], folder, "home/SCLC_hg19_Bam/" ,gsub(".bam", "", bam_list[i]), "_hg19.bam"))}
why using R ?
set the full patch to CrossMap.py, add the python interpreter :" /usr/bin/python /path/to/CrossMap.py ...."
@Pierre, what would be the best option?
@Pierre i could not get you
a simple bash script...
Do you mind providing a example?
http://tldp.org/HOWTO/Bash-Prog-Intro-HOWTO-7.html