miRNA differential expression analysis
1
0
Entering edit mode
@kamyasingh90-11058
Last seen 7.0 years ago

Hi,

I am working on elucidating pathologically relevant miRNAs responsible for disease progression and chemoresistance in chronic lymphocytic leukemia. For differential expression of miRNAs I am using the edgeR package. I am new to statistics, edger and any help would be highly appreciated. The samples are classified as 1) Indolent which has 7 patient samples 2) Chemo-responsive, 9 patient samples 3) Non-responsive, 7 patient samples 4) Control, 2 samples. I want to check the differential expression of miRNAs as follows: (a) Control vs Indolent+Chemo-responsive+ Non-responsive (b) Indolent vs Chemo-responsive+ Non-responsive  (c)  Indolent vs Chemo-responsive (d) Indolent vs chemo-responsive. I went through the user guide of edgeR and in section 3.2.3 GLM approach by using contrast argument of glmLRT contrast(-1,1,0) will compare B to A sample. This works well when comparing two groups. How do I work when I have to compare A vs B+C and A vs B+C+D.

edger glm.fit • 1.2k views
ADD COMMENT
2
Entering edit mode
@gordon-smyth
Last seen 2 hours ago
WEHI, Melbourne, Australia

Do you mean  A vs (B+C)/2? That would compare A to the average of B & C, and would be done by contrast = c(1,-1/2,-1/2,0) if there are four groups in total.

A vs (B+C+D)/3 would be done by contrast = c(1,-1/3,-1/3,-1/3).

ADD COMMENT
0
Entering edit mode

Thank you so much for your help. In my analysis I have my samples divided into four groups Indolent (7 samples), chemo-responsive(9 samples), non-responsive(7 samples) and control(2 samples) . I wanted to find the differentially expressed miRNA in the following sub-groups :

1. Ctrl vs indolent,chemo-responsive and non- responsive will this contrast matrix  A vs (B+C+D)/3 would be done by contrast = c(1,-1/3,-1/3,-1/3) work for me ?

Similarly with Indolent vs chemo-responsive and non-responsive. A vs (B+C)/2  contrast = c(1,-1/2,-1/2,0) 

ADD REPLY

Login before adding your answer.

Traffic: 422 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