Limma design matrix
2
0
Entering edit mode
R ▴ 40
@r-5604
Last seen 3.7 years ago
Germany

I have data with Treated and Control and three time-groups within each.
I want to compare across the time-groups for Treated vs Control.
How do I set up my desing matrix?

    groupTC = c("C","C","C","C","C","C","T","T","T")

    groupTC_Time = c("Early","Early","Mid,"Mid","Late,"Late","Early","Early","Mid,"Mid","Late,"Late",)

I want to compare Early vs Late and Early vs Mid for Treated vs Control?.

 

 

limma design matrix • 862 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 days ago
United States

Have you read the limma User's Guide? Section 9.6 specifically deals with timecourse experiments.

ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen just now
WEHI, Melbourne, Australia

I usually advise people to setup the design matrix in the standard way, as in Section 9.5.2 of the User's Guide:

Group <- paste(GroupTC, GroupTC_Time, sep=".")
Group <- factor(Group)
design <- model.matrix(~0+Group)
colnames(design) <- levels(Group)

Then use makeContrasts() to make whatever comparisons you like.

ADD COMMENT

Login before adding your answer.

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