[BioC] differential expression- edgeR

Ashutosh [guest] guest at bioconductor.org
Mon Dec 23 16:50:48 CET 2013


I am new to R and edgeR. I am trying to follow to example  
9. Case Study: deep-sequenced short tags from the edgeR manual. 

I download the data from http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM272105
by clicking the full table and then saving webpage as .txt file containing 

First 4 rows as below folled by fifth which contains the tag seq and count. 

#SEQUENCE = 
#COUNT = 
#TPM = tags per million
SEQUENCE	COUNT	TPM
CATCGCCAGCGGGCACC	1	0.37

Now we I follow the steps of 9.3 reading data and creating the DGElist: After running 
< d<- readDGE(targets, skip = 5, comment.char = "#"), I get 
Error in taglist[[i]] : subscript out of bounds

Can anyone please help, how I can solve this issue. 

Best regards,
Ashutosh


 -- output of sessionInfo(): 

R version 3.0.2 (2013-09-25)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] edgeR_3.4.2  limma_3.18.7


--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list