Entering edit mode
Luda Druskin
▴
10
@luda-druskin-757
Last seen 10.2 years ago
WE are trying to compare genes expression at different time points to
common reference (Ctl) and between each other.
This is a targets file:
slides Cy3 Cy5
367 96h Ctl
368 48h Ctl
381 48h Ctl
382 96h Ctl
383 48h Ctl
384 96h Ctl
561 12h Ctl
562 12h Ctl
563 12h Ctl
564 24h Ctl
565 24h Ctl
566 24h Ctl
This is a code I used:
"targets<- readTargets(file="targets.txt")
design<- designMatrix(targets, ref="Ctl")
design
96h 48h 12h 24h
1 -1 0 0 0
2 0 -1 0 0
3 0 -1 0 0
4 -1 0 0 0
5 0 -1 0 0
6 -1 0 0 0
7 0 0 -1 0
8 0 0 -1 0
9 0 0 -1 0
10 0 0 0 -1
11 0 0 0 -1
12 0 0 0 -1
fit <- lmFit(MA, design, weights = RG$weights)
colnames(design) <- c("96h", "48h", "12h", "24h")
contrast.matrix <- makeContrasts(96h, 48h,
24h,12h,96h-48h,48h-24h,24h-12h,levels=design)
Error: syntax error
contrast.matrix <- makeContrasts(96h-48h, 48h-24h, 24h-12h,
levels=design)
Error: syntax error"
What do I do wrong?
Thanks.
Luda