[BioC] limma makeContrasts
Lorenzo Bomba
lory.bomb at gmail.com
Wed Oct 6 18:56:13 CEST 2010
Dear all,
I have a conceptual problem with the "makeContrast" function in Limma:
doing a factorial design and the follow are the targetfile and the design matrix
targetfile
FileName Scrofa Treatment
S10M S127 DR
S11 S52 DR
S12 S304 DR
S13 S9 NODR
S14 S52 DR
S15 S304 DR
S17rt S304 DR
S18 S66 DR
S19 S66 NODR
S20 S304 NODR
S21 S93 NODR
S22 S93 NODR
> design
S127.DR S304.DR S304.NODR S52.DR S66.DR S66.NODR S9.NODR S93.NODR
1 1 0 0 0 0 0 0 0
2 0 0 0 1 0 0 0 0
3 0 1 0 0 0 0 0 0
4 0 0 0 0 0 0 1 0
5 0 0 0 1 0 0 0 0
6 0 1 0 0 0 0 0 0
7 0 1 0 0 0 0 0 0
8 0 0 0 0 1 0 0 0
9 0 0 0 0 0 1 0 0
10 0 0 1 0 0 0 0 0
11 0 0 0 0 0 0 0 1
12 0 0 0 0 0 0 0 1
attr(,"assign")
[1] 1 1 1 1 1 1 1 1
attr(,"contrasts")
attr(,"contrasts")$TS
[1] "contr.treatment"
I want see the difference between the NODR and the DR group taking in account that they belong to different mothers ( S127 S304 .....etc ). How can I set the makeContrast....
I tried in this two ways.... but I don't know the difference between the first and the second
1)cont.matrix <- makeContrasts( DRvsNODR=(S127.DR+S304.DR+S52.DR)-(S304.NODR+S66.NODR+S9.NODR+S93.NODR), levels=design)
2) cont.matrix <- makeContrasts( DRvsNODR=(S127.DR-S304.DR-S52.DR)-(S304.NODR-S66.NODR-S9.NODR-S93.NODR), levels=design)
Could you help me please???
Thanks in advance!
Lorenzo
More information about the Bioconductor
mailing list