Question: cuffdiff outputs empty "run.info" file - and readCufflinks() complains
0
gravatar for murraym
2.1 years ago by
murraym0
murraym0 wrote:

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

ADD COMMENTlink written 2.1 years ago by murraym0

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

ADD REPLYlink written 2.1 years ago by Arup Ghosh2.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2086 users visited in the last hour