[BioC] edgeR: error in readDGE
Wolfgang Huber
whuber at embl.de
Wed Nov 10 10:33:38 CET 2010
村田 卓也,
It looks like the second column ("COUNT") in your files is not always
numeric. Can you send the output of the following:
for (f in targets$files)
str(read.delim(f, stringsAsFactors=FALSE))
Best wishes
Wolfgang
Il Nov/10/10 9:36 AM, 村田 卓也 ha scritto:
> I am a beginner of everything of R and edgeR; just like "copy and paste of the command" of the manual.
>
> I tried to make DGEList object like descriptions on page 9.
> I could make object, "target", as follows.
>> targets
> files group description
> 1 GAT.txt MT FT9
> 2 ACA.txt WT FT9
> 3 ACT.txt MT FT10
> 4 ACG.txt WT FT10
>
> On the next command, I failed.
> d<- readDGE(targets)
>
> (I deleted "skip = 5, comment.char = "#"" from original command, because on my understanding,
> there were no rows to be skipped, nor comments.)
>
> error colSums(x$counts) : 'x' must be numeric
>
> This error message might be inaccurate, because it is localized.
>
> GAT.txt, and other txt files were made as a tab-delimited, plain-text, as follows, for example.
> SEQUENCE COUNT
> AAACAGCAGAGAGAGAGAGAGAGAGAG 0
> AAACAGCAGGTACAAGATCTACCCGGG 0
> AAACAGCAGTGTTCTGAAGCCAAACTC 0
> ....
>
> I use R for Mac OSX GUI 1.35 and newest edgeR.
>
> I am very happy, if you have any comments.
> Thank you.
>
> Takuya
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list