Question: Model selection for differential expression
gravatar for rmf
5 months ago by
rmf700 wrote:

I am using limma to test for differential expression. And I am trying to find the best model to use. My variable of interest is condition. I have several batch variables plate, well and gender which seems to be highly associated with top three principal components in a PCA.

I have found this function in limma, selectModel() which should help me to pick the best fitting model.

d1 <- model.matrix(~0+condition,data=md1)
d2 <- model.matrix(~0+condition+gender,data=md1)
d3 <- model.matrix(~0+condition+plate,data=md1)
d4 <- model.matrix(~0+condition+well,data=md1)
d5 <- model.matrix(~0+condition+plate+well,data=md1)
d6 <- model.matrix(~0+condition+gender+plate,data=md1)
d7 <- model.matrix(~0+condition+gender+well,data=md1)
d8 <- model.matrix(~0+condition+gender+plate+well,data=md1)
dlist <- list(d1,d2,d3,d4,d5,d6,d7,d8)

sm1 <- selectModel(dfr,dlist,criterion="aic")
sm2 <- selectModel(dfr,dlist,criterion="bic")

Here are barplots summarising the best models as picked by selectModel() using AIC and BIC criterion respectively.


enter image description here

X axis shows the model number and y axis shows the number of probes for which a model was picked as the best. Both the AIC and BIC metrics seem to be differing in opinion. What would you say is the best model here? or am I doing this wrong?

ADD COMMENTlink modified 5 months ago • written 5 months ago by rmf700
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: 654 users visited in the last hour