Tag: dgelist

rna seq – R – [DESeq2] – How use TMM normalized counts (from EdgeR) in inputs for DESeq2?

I have several RNAseq samples, from different experimental conditions. After sequencing, and alignment to reference genome, I merged the raw counts to get a dataframe that looks like this: > df_merge T0 DJ21 DJ24 DJ29 DJ32 Rec2 Rec6 Rec9 G10 421 200 350 288 284 198 314 165 G1000 17208…

Continue Reading rna seq – R – [DESeq2] – How use TMM normalized counts (from EdgeR) in inputs for DESeq2?

Unsuccessful DE analysis using limma

This might be a bit long, please bare with me. I’m conducting a differential expression analysis using limma – voom. My comparison is regarding response vs non-response to a cancer drug. However, I’m not getting any DE genes, absolute zeros. Someone here once recommended not to use contrast matrix for…

Continue Reading Unsuccessful DE analysis using limma

Removing replicate not clustering and group with replicate Vs without -edgeR rnaseq analysis

Removing replicate not clustering and group with replicate Vs without -edgeR rnaseq analysis 0 I am working with bacteria samples – in 3 groups that include the control, Treatment A, and Treatment B. From the PCA I find that the replicates are far apart. So I have removed the treatment…

Continue Reading Removing replicate not clustering and group with replicate Vs without -edgeR rnaseq analysis

r – Contrast for Limma – Voom

I’m doing a differential expression analysis for RNA-seq data with limma – voom. My data is about a cancer drug, 49 samples in total, some of them are responders some of them are not. I need some help building the contrast. I’m dealing with only one factor here, so two…

Continue Reading r – Contrast for Limma – Voom

Limma Differential Analysis on Proteomics data

Limma Differential Analysis on Proteomics data 1 @3c9b3fdc Last seen 7 hours ago United States Hi, I have a proteomics data set and I am doing the differential analysis on that. I used the Limma package to do that. I first removed the negative counts and did the analysis but…

Continue Reading Limma Differential Analysis on Proteomics data

DGE from tumor-adjacent normal pair RNA-seq data. For an individual, no replicate

Single sample analysis can be done in edgeR, while it is deprecated in DNASeq2 in 2018 probably (if you can use old DESeq version then it will work in single sample also). For egdeR you just do library(edgeR) setwd() rawdata <- read.delim(“filename”, check.names=FALSE, stringsAsFactors=FALSE) ngenes <- 10000 #no. of genes…

Continue Reading DGE from tumor-adjacent normal pair RNA-seq data. For an individual, no replicate

keep only genes expressed in sample at the same time as a particular gene of interest

Filter DGElist object: keep only genes expressed in sample at the same time as a particular gene of interest 0 Hello, need some help here as I’m kind of stuck with the edgeR DGElist format. I have a DGE list named x with the following dimensions: > dim(x ) [1]…

Continue Reading keep only genes expressed in sample at the same time as a particular gene of interest

Calculate fold change in edgeR with one sample per condition

Hi, We have run a pilot RNA-Seq study with one sample per condition, this is just a test run. I understand there is no valid statistical test in this case, however just curious to obtain differential expression through edgeR package in R assuming dispersion = 0.4 for the human data….

Continue Reading Calculate fold change in edgeR with one sample per condition