[R-sig-ME] False Convergence warning
Robert Ackerman
r@ckerm@n27 @end|ng |rom gm@||@com
Sat May 4 21:38:12 CEST 2019
Hi everyone,
I have data from a speed-dating study, where groups of men and women
interacted with each other for five minutes and subsequently provided
ratings of the interaction and their partner.
These are the first lines of my data set:
GROUPID MaleID FemaleID MF FM Agender AATTRACT
1 1 1 2 1 0 1 3.00
2 1 1 4 1 0 1 3.25
3 1 1 6 1 0 1 6.00
4 1 1 8 1 0 1 3.50
5 1 1 10 1 0 1 3.50
6 1 1 2 0 1 -1 5.00
7 1 3 2 0 1 -1 3.50
8 1 5 2 0 1 -1 5.00
9 1 7 2 0 1 -1 4.75
10 1 9 2 0 1 -1 2.25
11 1 3 2 1 0 1 1.25
12 1 3 4 1 0 1 3.25
Altogether, I have 904 observations (33 speed-dating groups comprised of
3-5 men and women each).
I’m trying to use glmmTMB to get a model with crossed random effects and an
unstructured covariance matrix for the residuals to run. Please note that I
was able to get this model to run in SPSS without problems (estimates of
fixed effects and random effects are also virtually identical). However, I
wanted to use R for its ease of checking multilevel modeling assumptions.
Here is the syntax for the model I tried to run in glmmTMB:
glmmTMB(AATTRACT ~ 1 + Agender+ (0 + MF + FM|GROUPID:MaleID) +
(0 + FM + MF|GROUPID:FemaleID) + us(MF + FM + 0|
GROUPID:MaleID:FemaleID),
data = SD_data, family = gaussian(link = "identity"),
dispformula=~0, REML = FALSE,
verbose = TRUE)
I received the following output and error message after running this model:
Formula: AATTRACT ~ 1 + Agender + (0 + MF + FM | GROUPID:MaleID) + (0 +
FM + MF | GROUPID:FemaleID) + us(0 + MF + FM + 0 | GROUPID:MaleID:FemaleID)
Dispersion: ~0
Data: SD_data
AIC BIC logLik df.resid
2891.436 2944.311 -1434.718 893
Random-effects (co)variances:
Conditional model:
Groups Name Std.Dev. Corr
GROUPID:MaleID MF 0.5572656
FM 0.7916470 -0.02
GROUPID:FemaleID FM 0.6732786
MF 0.7529185 -0.16
GROUPID:MaleID:FemaleID MF 0.8523621
FM 0.9419247 0.26
Residual 0.0001221
Number of obs: 904 / Conditional model: GROUPID:MaleID, 120;
GROUPID:FemaleID, 120; GROUPID:MaleID:FemaleID, 452
Fixed Effects:
Conditional model:
(Intercept) Agender
3.9048 0.1112
Warning message:
In fitTMB(TMBStruc) :
Model convergence problem; false convergence (8). See
vignette('troubleshooting')
Because I used verbose = TRUE, I also received the following mgc values
(I’m only including the last few lines here).
mgc: 5.077952e-08
outer mgc: 0.544794
iter: 1 value: -5887.447 mgc: 0.003153025 ustep: 1
iter: 2 value: -5887.447 mgc: 6.546202e-08 ustep: 1
mgc: 9.447453e-08
outer mgc: 0.2579731
iter: 1 value: -5887.346 mgc: 0.003152465 ustep: 1
iter: 2 value: -5887.346 mgc: 5.733777e-08 ustep: 1
mgc: 7.849359e-08
outer mgc: 0.2594169
Does anyone have any insights into what might be going wrong?
Thanks!
Rob
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list