Search
Question: how to read htseq-count
0
gravatar for Bioinformatics
4 months ago by
Bioinformatics30 wrote:

I am trying to read many files 

I do this and I get the error 

file.list <- list.files(pattern='*.txt')​
dge <- readDGE(file.list, header=FALSE)
Meta tags detected: __no_feature, __ambiguous, __too_low_aQual, __not_aligned, __alignment_not_unique
Error in colSums(x$counts) : 'x' must be numeric
ADD COMMENTlink modified 4 months ago by Gordon Smyth33k • written 4 months ago by Bioinformatics30
0
gravatar for Gordon Smyth
4 months ago by
Gordon Smyth33k
Walter and Eliza Hall Institute of Medical Research, Melbourne, Australia
Gordon Smyth33k wrote:

The error is telling you that at least one of the files has a non-numeric character in the second column.

Using list.files(pattern='*.txt') is convenient but dangerous. It is dangerous because it can easily include files that are not htseq output files and because, even if it selects the right files, you don't have any control over the order of files.

In this case, you might have one or more .txt files in the directory that are not a htseq-count output files. Have you checked that?

ADD COMMENTlink modified 4 months ago • written 4 months ago by Gordon Smyth33k

I see, now I get this warning

Meta tags detected: __no_feature, __ambiguous, __too_low_aQual, __not_aligned, __alignment_not_unique

I have a question, how can i change the annotation to gene name by your package? 

ADD REPLYlink modified 4 months ago by Gordon Smyth33k • written 4 months ago by Bioinformatics30

No, the message about meta tags is not a "warning". It's just an informative message. If it was a warning, then it would include the word: "Warning:".

The message alerts you to the fact that you will need to remove the meta tag rows from the DGEList object before you do any analysis.

To see how to add gene symbols to your annotation, follow the documentation in any example workflow. This will of course depend on what annotation you are using now, which is something you haven't told us.

 

ADD REPLYlink modified 4 months ago • written 4 months ago by Gordon Smyth33k
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.2.0
Traffic: 275 users visited in the last hour