[BioC] maSigPro error in getting a design matrix

lohith madireddy lohithreddym at gmail.com
Fri Oct 30 22:15:20 CET 2009


Dear Bioconductors,

I have the following design for my experiment. I am trying to use
maSigPro to assess temporal differential expression across 3
experimental groups I want to consider (control, vehicle and Var). The
design matrix is as shown below in the maSigPro format.


design
                  X Time Replicate Control Vehicle  Var
1   Veh Day0 A1.CEL    0         1       0       1   0
2   Veh Day0 A2.CEL    0         1       0       1   0
3   Veh Day0 A4.CEL    0         1       0       1   0
4   Veh Day0 A5.CEL    0         1       0       1   0
5  Veh Day12 A1.CEL   12         2       0       1   0
6  Veh Day12 A2.CEL   12         2       0       1   0
7  Veh Day12 A3.CEL   12         2       0       1   0
8   Veh Day0 B1.CEL    0         3       0       1   0
9   Veh Day0 B3.CEL    0         3       0       1   0
10  Veh Day0 B4.CEL    0         3       0       1   0
11 Veh Day12 B1.CEL   12         4       0       1   0
12 Veh Day12 B2.CEL   12         4       0       1   0
13 Veh Day12 B3.CEL   12         4       0       1   0
14  Veh Day0 C1.CEL    0         5       0       1   0
15  Veh Day0 C2.CEL    0         5       0       1   0
16  Veh Day0 C3.CEL    0         5       0       1   0
17  Veh Day0 C4.CEL    0         5       0       1   0
18  Veh Day0 C5.CEL    0         5       0       1   0
19 Veh Day12 C1.CEL   12         6       0       1   0
20 Veh Day12 C2.CEL   12         6       0       1   0
21 Veh Day12 C3.CEL   12         6       0       1   0
22     Day 0 A1.CEL    0         7       0       0   1
23     Day 0 A2.CEL    0         7       0       0   1
24     Day 0 A3.CEL    0         7       0       0   1
25     Day 5 A1.CEL    5         8       0       0   1
26     Day 5 A2.CEL    5         8       0       0   1
27     Day 5 A3.CEL    5         8       0       0   1
28    Day 12 A1.CEL   12         9       0       0   1
29    Day 12 A2.CEL   12         9       0       0   1
30    Day 12 A3.CEL   12         9       0       0   1
31    Day 31 A1.CEL   31        10       0       0   1
32    Day 31 A2.CEL   31        10       0       0   1
33    Day 31 A3.CEL   31        10       0       0   1
34     Day 0 B1.CEL    0        11       0       0   1
35     Day 0 B2.CEL    0        11       0       0   1
36     Day 0 B3.CEL    0        11       0       0   1
37     Day 0 B4.CEL    0        11       0       0   1
38     Day 5 B1.CEL    5        12       0       0   1
39     Day 5 B2.CEL    5        12       0       0   1
40     Day 5 B3.CEL    5        12       0       0   1
41    Day 12 B1.CEL   12        13       0       0   1
42    Day 12 B2.CEL   12        13       0       0   1
43    Day 12 B3.CEL   12        13       0       0   1
44    Day 31 B1.CEL   31        14       0       0   1
45    Day 31 B2.CEL   31        14       0       0   1
46    Day 31 B3.CEL   31        14       0       0   1
47     Day 0 C1.CEL    0        15       0       0   1
48     Day 0 C2.CEL    0        15       0       0   1
49     Day 0 C3.CEL    0        15       0       0   1
50     Day 0 C4.CEL    0        15       0       0   1
51     Day 5 C1.CEL    5        16       0       0   1
52     Day 5 C2.CEL    5        16       0       0   1
53     Day 5 C3.CEL    5        16       0       0   1
54    Day 12 C1.CEL   12        17       0       0   1
55    Day 12 C2.CEL   12        17       0       0   1
56    Day 12 C3.CEL   12        17       0       0   1
57    Day 31 C1.CEL   31        18       0       0   1
58    Day 31 C2.CEL   31        18       0       0   1
59    Day 31 C3.CEL   31        18       0       0   1
60     Naive A1.CEL    0        19       1       0   0
61     Naive A2.CEL    0        19       1       0   0
62     Naive A3.CEL    0        19       1       0   0
63     Naive B1.CEL    0        20       1       0   0
64     Naive B2.CEL    0        20       1       0   0
65     Naive C1.CEL    0        21       1       0   0
66     Naive C2.CEL    0        21       1       0   0
67     Naive C3.CEL    0        21       1       0   0

I am getting an error when I am going to create a maSigPro design
matrix and the error states
> design<-make.design.matrix(
design.init, degree=3)
Error in Summary.factor(c(47L, 48L, 49L, 50L, 59L, 60L, 61L, 51L, 52L,  :
  min not meaningful for factors

> design<-make.design.matrix(design.init)
Error in Summary.factor(c(47L, 48L, 49L, 50L, 59L, 60L, 61L, 51L, 52L,  :
  min not meaningful for factors

Is the error due to not having arrays at all 4 time points for the
vehicle and Naive? In other words, can maSigPro analyze when the
number of replicates for each experimental group is different and the
arrays are not available at all time points (0, 5, 12, 31) for every
experimental group (Control, Veh, Var)?
I have 4 time points overall and I am giving 4-1 = 3 as degree for the model.

Please let me know what can be wrong in my analysis.

sessionInfo()
R version 2.10.0 (2009-10-26)
i386-pc-mingw32

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United
States.1252    LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                           LC_TIME=English_United States.1252

attached base packages:
[1] tcltk     stats     graphics  grDevices utils     datasets
methods   base

other attached packages:
[1] maSigPro_1.17.5    DynDoc_1.23.1      widgetTools_1.23.1
genefilter_1.26.4  rgu34acdf_2.5.0    affy_1.23.12       Biobase_2.5.8

loaded via a namespace (and not attached):
 [1] affyio_1.13.5        annotate_1.23.4      AnnotationDbi_1.7.20
DBI_0.2-4            limma_3.0.3          Mfuzz_2.3.4
preprocessCore_1.7.9
 [8] RSQLite_0.7-3        splines_2.10.0       survival_2.35-7
tkWidgets_1.23.2     tools_2.10.0         xtable_1.5-5


Thanks in advance!

With Best Regards,
Lohith


More information about the Bioconductor mailing list