Question: DESeq2 RNA analysis work flow with NCBI GEO data
gravatar for maedakus
2.3 years ago by
maedakus10 wrote:


i am very new to Bioconductor, and am preparing for the analysis in DESeq2 packages for RNA data.

when data handling with GEO data(GSE22219), at first i got  GSE22219 data with ExpressionSet object,

so i need to convert it into DESeqDataSet class in order to use the DESeq2 package.

when i handled data as below, the error message appeared.

dds <- DESeqDataSetFromMatrix(countData = data,colData=pdata)

Error in DESeqDataSet(se, design = design, ignoreRank) :
  some values in assay are not integers

i am sure all data are integers,  if you know it , would you mind teaching me about it ?

thanks in advance !


ADD COMMENTlink modified 2.3 years ago by andrew.j.skelton73310 • written 2.3 years ago by maedakus10
gravatar for andrew.j.skelton73
2.3 years ago by
United Kingdom
andrew.j.skelton73310 wrote:

This GEO dataset is from an expression array, not suitable for DESeq2. DESeq2 requires integer counts, typically found through programs like htseq_count, which takes an RNA seq alignment. For the data you have, it's worth checking out Limma

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by andrew.j.skelton73310
Please log in to add an answer.


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