Question: FeatureCounts from Galaxy to DESeq2 in RStudio
gravatar for Newbieish
2.7 years ago by
Newbieish0 wrote:

I'm having difficulty performing DESeq2 analysis on FeatureCounts data from Galaxy. I've tried both these websites, but I can't seem to properly format my FeatureCounts data for DESeq2. My FeatureCounts output file from Galaxy lists only the GeneId and Counts columns and it doesn't have the headers. Does anyone have the code for that? I've basically been trying:

# Convert to matrix
countdata <- as.matrix(countdata)

# Assign condition (first four are controls, second four contain the expansion)
(condition <- factor(c(rep("ctl", 4), rep("exp", 4))))

# Analysis with DESeq2 ----------------------------------------------------


# Create a coldata frame and instantiate the DESeqDataSet. See ?DESeqDataSetFromMatrix
(coldata <- data.frame(row.names=colnames(countdata), condition))
dds <- DESeqDataSetFromMatrix(countData=countdata, colData=coldata, design=~condition)
rna-seq • 1.9k views
ADD COMMENTlink modified 2.7 years ago by WouterDeCoster44k • written 2.7 years ago by Newbieish0
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.3.0
Traffic: 1600 users visited in the last hour