Question: Limma contrasts matrix drops columns and considers them as reference column
0
gravatar for fr
9 weeks ago by
fr150
fr150 wrote:

In the limma vignette page 44, a simple 2 factorial example is given. My case follows something similar, where I have treatment and time.

TS <- paste(targets$Strain, targets$Treatment, sep=".")
TS <- factor(TS, levels=c("WT.U","WT.S","Mu.U","Mu.S"))
design <- model.matrix(~0+TS)
colnames(design) <- levels(TS)

The resulting design contains all possible columns.

In my case, I'm trying to do the same, but specifying different columns (strain and treatment), and then specifying design = model.matrix(~0 + strain + treatment, data=targets), but this results in 1 less column in strain and in treatment. I already know that a workaround is to merge everything as shown above and in the vignette, but why can't I do it in the same way when I specify data=targets?

Note that I have assigned the intercept as 0 just like in the design above.

limma design • 148 views
ADD COMMENTlink modified 9 weeks ago by Luiz30 • written 9 weeks ago by fr150
0
gravatar for Luiz
9 weeks ago by
Luiz30
Brazil / São Paulo / UNIFESP
Luiz30 wrote:

You can explicitly tell your comparison design through contrast.matrix() function.

ADD COMMENTlink written 9 weeks ago by Luiz30

From which package does that come from? And taking the example of the vignette for instance, how would you do that?

ADD REPLYlink modified 8 weeks ago • written 8 weeks ago by fr150
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1103 users visited in the last hour