This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error in running DESeqDataSetFromTximport

Hello,

I am facing an error while running DESeqDatasetFromTximort after StringTie quantification step

> library(DESeq2)
> sampleTable <- data.frame(condition = factor(rep(c("mock", "mock"), each = 1)))
> rownames(sampleTable) <- colnames(txi$counts)
> dds <- DESeqDataSetFromTximport(txi = txi, sampleTable ,design = ~ Condition + Time)
Error in DESeqDataSet(se, design = design, ignoreRank) : 
  all variables in design formula must be columns in colData

My sample contains:

Sample_ids  Condition   Time
t_data_9.ctab   Mock    6hpi
t_data_10.ctab  Mock    6hpi

Am I supposed to add any argument?

tximport r rna-seq stringtie deseq2
• 3,990 views • link
updated by Ram Ram • written by ankur • ankur

1 answer

You need to redo how sampleTable is made. Firstly, it needs to have Condition rather than condition (or you need to change that in your design) and secondly it needs to contain Time, which it currently doesn't. Fixing that will remove the error.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this
  • DESEQ2 Code
    written by ananta.kapoor • ananta.kapoor
    # Load the DESeq2 library library(DESeq2) # Read your count data counts <- read.csv("C:/Users/anant/OneDrive/Documents/RNASeq/count_data.csv", row.names = 1) # Read the metadata metadata <- read.csv("C:/Users/anant/OneDrive/Documents/RNASeq/Metadata.csv") # …
  • Deseq2 dataframematrix error
    written by Mewz1k • Mewz1k
    Error in DESeqDataSet(se, design = design, ignoreRank) : some values in assay are negative I have raw counts of rna seq data with approxx 68 …
  • TCGA single dataset normalization
    written by ruby4bioinfo • ruby4bioinfo
    how can i resolve it dds <- DESeqDataSetFromMatrix(countData = PAADgenomicMatrix, colData = Samples, design = ~ submitter_id.samples) Error in DESeqDataSet(se, design = design, ignoreRank) : …
  • DESeq2 for Differential expression Ensembl gene ID
    written by Rob • Rob
    Hi friends I use DESeq2 for differential expression analysis for HTSeq read count data with about 20000 genes. my genes have ENSEMBL ID. I have …
  • tximport output should be integer values, IS NOT.
    written by c_u c_u
    Hi, I am analyzing RNA-seq data and I used Stringtie to assemble transcripts and quantify them. Now to perfom DE analysis on them (using DESeq2), …
  • DESeqDataSetFromMatrix giving an error
    written by natheshka • natheshka
    Hi, I am trying to follow the example in MLSeq documentation. But I keep on getting errors. The latest one I am getting is > …
  • including batch in design~
    written by Morris_Chair Morris_Chair
    Hello, I want to add the batch option in DESeqDataSet but I get error dds=DESeqDataSetFromTximport(txi,colData = samples,design=~ batch + condition) **Error in DESeqDataSet(se, design = …
  • adding batch variable in colData(dds)
    written by Mozart Mozart
    Dear all, I am quite struggling with some sort of coding here. Essentially, what I am trying to do is to simply add the batch …
  • DESeq2 and integers values
    written by Morris_Chair Morris_Chair
    Dear all, I have an experiment with two different conditions, two treated and two untreated samples from where I performed RNA sequencing. here is the …
  • Tximport > DESeq2 (sample names)
    written by Barry Digby Barry Digby
    Using R and tximport, I have read in 6 t_data.ctab files to create a DESeqDataSet by using DESeqDataSetFromTximport. When plotting heatmaps etc, the sample names …