[R-sig-ME] Regarding mcmcsamp in the most recent alpha release
Dave Warren
dwarren1 at cyrus.psych.uiuc.edu
Sat Mar 29 05:39:51 CET 2008
Hi all,
I've been trying out the most recent alpha release and keep running
into the following exception when I try using mcmcsamp.
> path.lmer = lmer( log( PathLength ) ~ Trial * Status + (Trial |
Patient), data = path_agg )
> summary( path.lmer )
Linear mixed model fit by REML
Formula: log(PathLength) ~ Trial * Status + (Trial | Patient)
Data: path_agg
AIC BIC logLik deviance REMLdev
640.4 674.1 -312.2 594.9 624.4
Random effects:
Groups Name Variance Std.Dev. Corr
Patient (Intercept) 1.3567e-03 0.0368331
Trial 2.1394e-05 0.0046254 -1.000
Residual 1.9524e-01 0.4418560
Number of obs: 502, groups: Patient, 10
Fixed effects:
Estimate Std. Error t value
(Intercept) 2.995689 0.060844 49.24
Trial -0.001977 0.002406 -0.82
StatusNormal 0.245692 0.084061 2.92
Trial:StatusNormal 0.008709 0.003382 2.58
Correlation of Fixed Effects:
(Intr) Trial SttsNr
Trial -0.661
StatusNorml -0.724 0.479
Trl:SttsNrm 0.470 -0.711 -0.660
> a = mcmcsamp( path.lmer, 10000 )
Error in .local(object, n, verbose, ...) :
crossproduct matrix 1 is not positive definite
Have I broken something? Is my model revealed to be invalid?
Thanks for any thoughts,
Dave
More information about the R-sig-mixed-models
mailing list