[BioC] Using limma with contrast matrix ,replicate spots,
Pita
pwilkinson_m at xbioinformatics.org
Tue Jan 25 00:30:13 CET 2005
What is the significance of this error? What does this tell me about these
two samples?. I added the coeficients for the donor then re-calculated the fit.
> donor.fit = lmFit(MA, design.III)
Coefficients not estimable: HISS0039 HISS0050
The design matrix I used was as follows, you can see that there is a column
of 1's that represent that the first 16 samples are for donor HISS0033, the
rest of the matrix for each donor looks the same but with the 1's in the
rescpective donor column:
X A_0h A_24h A_6h A_72h B_0h B_24h B_6h B_72h C_0h C_24h C_6h C_72h
N_0h N_24h N_6h N_72h HISS0033 HISS0036 HISS0037 HISS0038 HISS0039 HISS0050
1 1 0 0 0 0 0 0 0 0 0 0 0 0
1 0 0 0 1 0 0 0 0 0
2 2 0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0
3 3 0 0 0 0 0 0 0 0 1 0 0 0
0 0 0 0 1 0 0 0 0 0
4 4 1 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0
5 5 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 0 1 0 0 0 0 0
6 6 0 0 0 0 0 0 1 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0
7 7 0 0 0 0 0 0 0 0 0 0 1 0
0 0 0 0 1 0 0 0 0 0
8 8 0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0 0 0
9 9 0 0 0 0 0 0 0 0 0 0 0 0
0 1 0 0 1 0 0 0 0 0
10
10 0 0 0 0 0 1 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0
11
11 0 0 0 0 0 0 0 0 0 1 0 0 0
0 0 0 1 0 0 0 0 0
12
12 0 1 0 0 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0
13
13 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 1 0 0 0 0 0
14
14 0 0 0 0 0 0 0 1 0 0 0 0 0
0 0 0 1 0 0 0 0 0
15
15 0 0 0 0 0 0 0 0 0 0 0 1 0
0 0 0 1 0 0 0 0 0
16
16 0 0 0 1 0 0 0 0 0 0 0 0 0
0 0 0 1 0 0 0 0 0
More information about the Bioconductor
mailing list