[R-sig-ME] false positives?
John Kingston
jkingston at linguist.umass.edu
Mon May 26 05:34:52 CEST 2014
I received the warning messages below in response to running the
second model displayed below. The data file is attached to this message.
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.212882 (tol = 1e-06)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: very large eigenvalue
- Rescale variables?
sfHIndR <- glmer(fRspRN ~ cStep +
h1Context +
h2Context +
cQRT +
(1 + cStep +
h1Context +
h2Context +
cQRT | participant),
family = "binomial",
control = glmerControl(optCtrl = list(maxfun=1000000)),
data = sfCompiled)
summary(sfHIndR)
# interactions with RT
sfHSxRtCxRtR <- glmer(fRspRN ~ cStep * cQRT +
h1Context * cQRT +
h2Context * cQRT +
(1 + cStep +
h1Context +
h2Context +
cQRT | participant),
family = "binomial",
control = glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun=1000000)),
data = sfCompiled)
The test suggested in a earlier post produced the output below. These
values look small.
gg <- sfHSxRtCxRtR at optinfo$derivs$grad
hh <- sfHSxRtCxRtR at optinfo$derivs$Hessian
vv <- sqrt(diag(solve(hh/2)))
summary(abs(gg*vv))
# Min. 1st Qu. Median Mean 3rd Qu. Max.
# 6.720e-07 3.108e-05 8.866e-05 2.542e-04 1.745e-04 3.345e-03
Any advice about how to proceed? Thanks.
Best,
John
John Kingston
--------------------------------
Professor and Head
Linguistics Department
University of Massachusetts
150 Hicks Way, 226 South College
Amherst, MA 01003-9274
1-413-545-6837, fax -2792
jkingston at linguist.umass.edu
https://blogs.umass.edu/jkingstn/
More information about the R-sig-mixed-models
mailing list