Removing replicate not clustering and group with replicate Vs without -edgeR rnaseq analysis
0
0
Entering edit mode
24 months ago
veronica • 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 A replicate r2 and treatment B r2 replicate and used edgeR (qCML) n=2 Vs n=1 analysis. I know rnaseq without sufficient replicates is not of any value But I want a suggestion if I can still use the data only to select those genes with log2fold change >2 and <-2?

or use GFOLD ( or gfold diff -s1 sample1,sample2 -s2 sample3 -suf .read_cnt -o 12VS3.diff) and present only the generalized log fold for selected genes. Any advice? enter image description here

group<-c(1,1,2,3)
y <- DGEList(counts=x, group=group)
y <- estimateDisp(y)
et <- exactTest(y,pair=c(1:2)
 topTags(et)
RNASeq edgeR • 367 views
ADD COMMENT

Login before adding your answer.

Traffic: 1650 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6