cuffdiff outputs empty "run.info" file - and readCufflinks() complains
0
0
Entering edit mode
3.7 years ago
murraym • 0

Hi, We are analysing an RNAseq experiment using cuffquant and then cuffdiff on the resultant .cxb files, and then trying to load the data in R using the "readCufflinks" command in the cummeRbund package. readCufflinks() fails and writes an error message saying the run.info file has no lines. run.info is indeed empty, as are several other files in the cuffdiff output directory, but there are also other files that look fine (e.g. "gene_exp.diff"), so we are not sure if those empty files are normal. We suspect not, since cuffdiff eventually crashes with a segmentation fault 11 (just after printing "Writing run info"). Any suggestions as to what might be wrong would be appreciated. The commands used are below. thanks, Mike

here's the R command.

cuff<-readCufflinks() Creating database /Users/fred/Desktop/cuffdiff_output_dir/run.info Error in (function (file, header = FALSE, sep = "", quote = "\"'", dec = ".", : no lines available in input In addition: There were 50 or more warnings (use warnings() to see the first 50)

and here are the Unix commands.

cuffquant -o ./Ctl5A_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Ctl5A/accepted_hits.bam cuffquant -o ./Ctl5B_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Ctl5B/accepted_hits.bam cuffquant -o ./Ctl5C_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Ctl5C/accepted_hits.bam cuffquant -o ./Exp5A_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Exp5A/accepted_hits.bam cuffquant -o ./Exp5B_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Exp5B/accepted_hits.bam cuffquant -o ./Exp5C_quant -p4 ~/Desktop/Genome/gencode.v21.annotation.gtf ./Exp5C/accepted_hits.bam

cuffdiff ~/Desktop/Genome/gencode.v21.annotation.gtf -o ./cuffdiff_output_dir Ctl5A_quant/abundances.cxb,Ctl5B_quant/abundances.cxb,Ctl5C_quant/abundances.cxb Exp5A_quant/abundances.cxb,Exp5B_quant/abundances.cxb,Exp5C_quant/abundances.cxb

RNA-Seq cuffdiff cummeRbund readCufflinks • 1.2k views
ADD COMMENT
0
Entering edit mode

Directory name missing in cuff<-readCufflinks( dir=?).

ADD REPLY

Login before adding your answer.

Traffic: 2863 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6