CummeRbund error

[steindorff]

Hi,
I'm stuck in this cummeRbund error:

library(cummeRbund)
cuff_data <- readCufflinks('diff_out')
Creating database diff_out/cuffData.db
Reading diff_out/genes.fpkm_tracking
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 1 did not have 33 elements

If someone can help...
Thanks

[mikeyt]

Hi,

I am having the same error:

cuff_data <- readCufflinks ('diff_out')
Creating database diff_out/cuffData.db
Reading diff_out/genes.fpkm_tracking
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 1 did not have 25 elements


Were you able to resolve this problem?

Thanks.

[steindorff]

I look all the files searching the problem...the problem was in the GFF3 (or GTF) file. The '#' character appeared in some gene names.
Probably your annotation file have some character that the R don't recognize.

[mikeyt]

Thanks. I deleted all the '#' characters and received some progress but another error.

> cuff_data <- readCufflinks ('diff_out')
Creating database diff_out/cuffData.db
Reading diff_out/genes.fpkm_tracking
Checking samples table...
Populating samples table...
Writing genes table
Reshaping geneData table
Recasting
Writing geneData table
Reading diff_out/gene_exp.diff
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 1 did not have 14 elements

It may just be the same problem but with the gene_exp.diff file this time.
Thanks again and I appreciate your help!

[steindorff]

I'm glad that I can help.
I ran all the TopHat and cufflinks pipeline with the GFF3 corrected to avoid errors and it worked well for me.

[s869023]

Hi
I have similar problem while loading database in commeRbund

> library("cummeRbund")
> cuff_data<-readCufflinks('cuffdiff_output_sample-ATvsBC')
Creating database /cuffdiff_output_sample-ATvsBC/cuffData.db
Reading Run Info File cuffdiff_output_sample-ATvsBC/run.info
Writing runInfo Table
Reading Read Group Info cuffdiff_output_sample-ATvsBC/read_groups.info
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 1 did not have 7 elements

Then I try

> cuff_data
CuffSet instance with:
0 samples
0 genes
0 isoforms
0 TSS
0 CDS
0 promoters
0 splicing
0 relCDS

The cuffData.db file is only 70kb, so I assume there are some error in loading db.
Any help will be highly appreciated.

[jp.]

Hi there
I got the same problem... any solution ?

Hi
I have similar problem while loading database in commeRbund

> library("cummeRbund")
> cuff_data<-readCufflinks('cuffdiff_output_sample-ATvsBC')
Creating database /cuffdiff_output_sample-ATvsBC/cuffData.db
Reading Run Info File cuffdiff_output_sample-ATvsBC/run.info
Writing runInfo Table
Reading Read Group Info cuffdiff_output_sample-ATvsBC/read_groups.info
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, :
line 1 did not have 7 elements

Then I try

> cuff_data
CuffSet instance with:
0 samples
0 genes
0 isoforms
0 TSS
0 CDS
0 promoters
0 splicing
0 relCDS

The cuffData.db file is only 70kb, so I assume there are some error in loading db.
Any help will be highly appreciated.