[R-sig-ME] Warning message with lmer function
FMH
kagba2006 at yahoo.com
Wed Sep 16 12:17:28 CEST 2009
Dear All,
I have a set of data which consist of 1575 groups with 11 temperature values in each group. This temperature is recorded across 11 different depths of the sea. The R script are shown below:
############################################################
#Temp : Temperature
#dp1, dp2, dp3 : covariate with respect to linear, quadratic and cubic terms
#group : 1575 groups in which there are 11 observations in each group
#sub3 : Data set 1
#set3 : Data set 2
dp1 <- rep(rev(seq(1,51, by = 5)),1575)
dp2 <- dp1^2
dp3 <- dp1^3
group <- rep(1:1575, each = 11)
set3 <- data.frame(sub3, dp1, dp2, dp3)
(lm.lme3 <- lmer(Temp ~ dp1 + dp2 + dp3 + (dp1 + dp2 + dp3|group), data = set3))
############################################################
I tried to fit a linear mixed model via lmer function, with all the fixed and random effects are included, but there is a 'Warning' message given after the output, as shown below.
##############################################
Linear mixed model fit by REML
Formula: Temp ~ dp1 + dp2 + dp3 + (dp1 + dp2 + dp3 | group)
Data: set3
AIC BIC logLik deviance REMLdev
65627 65743 -32799 65539 65597
Random effects:
Groups Name Variance Std.Dev. Corr
group (Intercept) 4.8336e-01 6.9524e-01
dp1 5.2199e-04 2.2847e-02 0.000
dp2 3.0528e-07 5.5252e-04 0.000 0.000
dp3 7.8888e-11 8.8819e-06 -0.966 0.000 0.000
Residual 1.9939e+00 1.4120e+00
Number of obs: 17325, groups: group, 1575
Fixed effects:
Estimate Std. Error t value
(Intercept) 1.363e+01 4.043e-02 337.1
dp1 -2.930e-01 6.328e-03 -46.3
dp2 3.924e-03 2.870e-04 13.7
dp3 -1.900e-05 3.628e-06 -5.2
Correlation of Fixed Effects:
(Intr) dp1 dp2
dp1 -0.739
dp2 0.616 -0.959
dp3 -0.563 0.904 -0.982
Warning message:
In mer_finalize(ans) : false convergence (8)
#############################################
Does the output given is valid? Could someone please advice on this message.
Thank you
Fir
More information about the R-sig-mixed-models
mailing list