Question: DESeq2 RNA analysis work flow with NCBI GEO data
gravatar for maedakus
15 months 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 15 months ago by andrew.j.skelton73280 • written 15 months ago by maedakus10
gravatar for andrew.j.skelton73
15 months ago by
United Kingdom
andrew.j.skelton73280 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 15 months ago • written 15 months ago by andrew.j.skelton73280
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: 277 users visited in the last hour