[BioC] negative correlation?
Naomi Altman
naomi at stat.psu.edu
Wed Jun 29 19:24:54 CEST 2011
I am running limma 2.8 on a technical replicate dye-swap
design. Here is my design matrix:
(Intercept) as.factor(targets$time)24h as.factor(targets$time)4h
as.factor(targets$time)8h
1 -1 0 0
0
2 1 0 0
0
3 -1 0 0
0
4 1 0 0
0
5 -1 0 1
0
6 1 0 1
0
7 -1 0 1
0
8 1 0 1
0
9 -1 0 0
1
10 1 0 0
1
11 -1 0 0
1
12 1 0 0
1
13 -1 1 0
0
14 1 1 0
0
15 -1 1 0
0
16 1 1 0
0
corfit <- duplicateCorrelation(MA, design=designTimeR,ndups = 1,
block = targets$biorep)
targets$biorep
[1] 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8
> corfit$con
[1] -0.3414428
Why am I getting a negative correlation between bio reps? Is this
expected due to the dye swap?
--Naomi
More information about the Bioconductor
mailing list