Dear all,

 I'm using the latest version of the package lme4 and R version 2.7.2
(2008-08-25).

After I run the model, I get the results of the model (cf. below). Then, I
run an ANOVA using the "anova" function and I get the following message
"Error in anova(lmer1) : Calculated PWRSS for a LMM is negative".

 I went trough the R-mailing list and a similar error message was reported
but only when using "weights". As I did not use "weights" in my model, I was
wondering what could cause this error message!

Thanks in advance for any advice,

 Sebastien.

R-Console

> (lmer1=lmer(ano$qf~ano$p3+(ano$p3|ano$cxx)+(ano$cxx|ano$id)))
Linear mixed model fit by REML
Formula: ano$qf ~ ano$p3 + (ano$p3 | ano$cxx) + (ano$cxx | ano$id)
   AIC   BIC logLik deviance REMLdev
 34167 34670 -16995    34029   33991
Random effects:
 Groups   Name         Variance Std.Dev.
Corr

 ano$id   (Intercept)  605250
777.98

          ano$cxxcII   674734   821.42
-0.648

          ano$cxxcIII  603555   776.89   -0.279
0.181
          ano$cxxcIV   484725   696.22    0.288 -0.187
-0.080
          ano$cxxcIX   699289   836.24    0.401 -0.260 -0.113
0.115
          ano$cxxcV    642216   801.38   -0.387  0.251  0.108 -0.111
-0.155
          ano$cxxcVI   602526   776.23   -0.433  0.281  0.121 -0.125 -0.173
0.167
          ano$cxxcVII  352721   593.90   -0.491  0.318  0.137 -0.141 -0.196
0.190  0.212
          ano$cxxcVIII 851755   922.91   -0.746  0.484  0.208 -0.215 -0.299
0.289  0.323  0.366
          ano$cxxcX    499538   706.78   -0.189  0.122  0.053 -0.054 -0.075
0.073  0.082  0.093  0.141
          ano$cxxcXI   513851   716.83   -0.125  0.081  0.035 -0.036 -0.050
0.048  0.054  0.061  0.093  0.023
          ano$cxxcXII  802332   895.73   -0.658  0.427  0.184 -0.190 -0.264
0.254  0.285  0.323  0.491  0.124  0.082
 ano$cxx  (Intercept)  251374
501.37

          ano$p32       12140   110.18
-0.866

          ano$p33       48679   220.63   -0.874
0.757
 Residual              150767
388.29

Number of obs: 2238, groups: ano$id, 254; ano$cxx, 12

Fixed effects:
            Estimate Std. Error t value
(Intercept)  76748.8      312.9  245.28
ano$p32      -1343.0      409.3   -3.28
ano$p33      -3054.5      349.6   -8.74

Correlation of Fixed Effects:
        (Intr) an$p32
ano$p32 -0.764
ano$p33 -0.895  0.684
> anova(lmer1)
Error in anova(lmer1) : Calculated PWRSS for a LMM is negative


Sebastien.

	[[alternative HTML version deleted]]

