[R] Bug in lmer?

From: Mark Lyman <mlyman_at_byu.edu>
Date: Fri 30 Sep 2005 - 06:44:38 EST


I am relatively new to R so I am not confident enough in what I am doing to be certain this is a bug. I am running R 2.1.1 on a Windows XP machine and the lme4 package version 0.98-1. The following code fits the model I want using the nlme package version 3.1-60.

mltloc$loc <- factor(mltloc$loc)
mltloc$block <- factor(mltloc$block)
mltloc$trt <- factor(mltloc$trt)

Mltloc <- groupedData(adg~trt|loc,mltloc) plot(Mltloc)
plot(Mltloc,outer=~trt)

lme(adg~trt,
random=pdBlocked(list(pdIdent(~1),pdIdent(~block-1),pdIdent(~trt-1))),Mltloc)

The problem is that when I try fitting the model using the lmer function with the following code:

lmer(adg~trt+(1|loc)+(1|block:loc)+(1|loc:trt),mltloc)

I get this message from Windows and R closes.

R for Windows GUI front-end has encountered a problem and needs to close. We are sorry for the inconvenience.

This same code works on a Macintosh. So it doesn't seem that I have made an error in my code. Also if anyone of the random effect terms is removed there is no problem. Is this something that is being looked at? Or I have I made a mistake somewhere? I have included the data that I am using below.

      X obs loc block trt     adg      fe
1     1   3   A     1   3 3.16454  7.1041
2     2   4   A     1   4 3.12500  6.6847
3     3   6   A     1   2 3.15944  6.8338
4     4   7   A     1   1 3.25000  6.5254
5     5   9   A     2   2 2.71301  8.2505
6     6  10   A     2   1 3.20281  7.5922
7     7  12   A     2   3 3.02423  7.3894
8     8  16   A     2   4 2.87245  7.4604
9     9  19   A     3   1 2.68878  8.2785
10   10  20   A     3   2 2.86862  7.9470
11   11  21   A     3   3 2.89923  7.9739
12   12  22   A     3   4 3.02806  8.4331
13   13  25   B     1   3 2.18131  6.6691
14   14  27   B     1   4 2.51914  5.6281
15   15  28   B     1   2 1.88739  7.0723
16   16  31   B     1   1 2.34685  6.0295
17   17  33   B     2   4 2.45608  5.6195
18   18  35   B     2   1 2.25225  6.3978
19   19  36   B     2   3 2.23649  6.1799
20   20  40   B     2   2 2.47523  5.9985
21   21  41   B     3   2 1.94200  7.2975
22   22  44   B     3   3 2.43243  6.4350
23   23  47   B     3   4 2.30180  6.3339
24   24  48   B     3   1 2.53378  6.1564
25   25  50   C     1   4 2.96014  7.5110
26   26  51   C     1   2 3.23551  7.4762
27   27  54   C     1   3 3.24638  7.2063
28   28  56   C     1   1 3.04710  7.6389
29   29  58   C     2   3 3.26449  7.5466
30   30  59   C     2   2 2.71377  9.0895
31   31  61   C     2   1 3.06522  7.8723
32   32  62   C     2   4 2.71739  8.2318
33   33  66   C     3   4 3.03623  7.9426
34   34  67   C     3   3 3.10507  8.4608
35   35  69   C     3   1 3.16304  8.5549
36   36  70   C     3   2 3.02899  8.5038
37   37  74   D     1   1 2.49164  9.5758
38   38  77   D     1   3 2.51833  9.5121
39   39  79   D     1   2 2.35631 10.3264
40   40  80   D     1   4 2.30331  9.7715
41   41  81   D     2   3 2.72688  9.5628
42   42  83   D     2   2 2.59512  9.9414
43   43  85   D     2   1 2.56516  9.3887
44   44  88   D     2   4 2.91523  8.3158
45   45  89   D     3   3 2.57943 10.4416
46   46  90   D     3   4 2.98159  8.7710
47   47  93   D     3   2 2.35370 11.0148
48   48  94   D     3   1 2.21953 11.2417
49   49  99   E     1   3 2.84158  8.7886
50   50 100   E     1   4 2.65264  8.6946
51   51 102   E     1   2 2.47112  9.7143
52   52 103   E     1   1 2.89769  9.2401
53   53 105   E     2   2 2.57343  9.5353
54   54 106   E     2   1 2.99752  8.7538
55   55 110   E     2   4 2.95380  8.8210
56   56 112   E     2   3 3.08663  8.9427
57   57 114   E     3   1 2.72525  9.4308
58   58 115   E     3   2 2.75825  9.7721
59   59 116   E     3   3 3.08333  8.9010
60   60 117   E     3   4 3.12129  8.4852
61   61 122   F     1   1 3.20600  6.3983
62   62 123   F     1   2 2.89500  6.6569
63   63 125   F     1   4 3.36900  6.0821
64   64 126   F     1   3 3.12000  6.5349
65   65 130   F     2   2 3.19300  6.6729
66   66 131   F     2   1 3.29800  6.5488
67   67 133   F     2   4 3.09700  6.6598
68   68 135   F     2   3 3.38500  6.2998
69   69 139   F     3   3 3.44900  6.2849
70   70 140   F     3   2 3.05000  6.9957
71   71 141   F     3   1 3.43500  6.7302
72   72 143   F     3   4 3.60600  6.3827
73   73 145   G     1   2 2.58669  8.1394
74   74 147   G     1   1 3.17892  7.0972
75   75 148   G     1   4 2.95284  7.3140
76   76 151   G     1   3 3.17924  6.9430
77   77 154   G     2   4 2.62344  7.5150
78   78 155   G     2   3 2.64286  8.0237
79   79 157   G     2   1 3.12760  7.3169
80   80 160   G     2   2 2.54993  8.1957
81   81 163   G     3   4 2.58322  7.9687
82   82 164   G     3   3 2.84813  7.9284
83   83 166   G     3   2 2.69279  8.5303
84   84 167   G     3   1 3.14424  7.3564
85   85 169   H     1   3 3.39974  6.5945
86   86 173   H     1   2 3.12370  6.7530
87   87 175   H     1   4 3.17969  6.4279
88   88 176   H     1   1 3.70052  6.4830
89   89 177   H     2   2 2.95192  7.3809
90   90 179   H     2   3 3.44661  6.7929
91   91 182   H     2   4 3.28906  6.4807
92   92 184   H     2   1 3.37500  6.8139
93   93 188   H     3   4 3.65104  6.3068
94   94 190   H     3   1 3.27734  7.4789
95   95 191   H     3   3 3.42708  6.9327
96   96 192   H     3   2 3.04818  7.7264
97   97 193   I     1   2 2.22105  8.4243
98   98 196   I     1   1 3.15526  6.7119
99   99 197   I     1   4 2.40263  7.7486
100 100 198   I     1   3 3.00000  6.9215
101 101 203   I     2   1 2.29079  8.8861
102 102 204   I     2   4 2.25395  8.6850
103 103 205   I     2   3 2.60526  8.4150
104 104 208   I     2   2 2.34737  8.7866
105 105 209   I     3   2 2.50505  8.5975
106 106 212   I     3   1 2.31316  8.9267
107 107 213   I     3   4 2.42105  8.7750
108 108 214   I     3   3 2.74211  8.1599

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Sep 30 06:47:08 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:40:34 EST