Moje dane są opisane tutaj Co może powodować, że „Model błędu () jest pojedynczym błędem” w aov przy dopasowywaniu ANOVA z powtarzanymi pomiarami?
Próbuję zobaczyć efekt interakcji przy użyciu, lmer
więc mój podstawowy przypadek to:
my_null.model <- lmer(value ~ Condition+Scenario+
(1|Player)+(1|Trial), data = my, REML=FALSE)
my.model <- lmer(value ~ Condition*Scenario+
(1|Player)+(1|Trial), data = my, REML=FALSE)
Uruchomienie anova
daje mi znaczące wyniki, ale kiedy próbuję uwzględnić przypadkowe nachylenie ( (1+Scenario|Player)
), model nie działa z tym błędem:
Warning messages:
1: In commonArgs(par, fn, control, environment()) :
maxfun < 10 * length(par)^2 is not recommended.
2: In optwrap(optimizer, devfun, getStart(start, rho$lower, rho$pp), :
convergence code 1 from bobyqa: bobyqa -- maximum number of function evaluations exceeded
3: In commonArgs(par, fn, control, environment()) :
maxfun < 10 * length(par)^2 is not recommended.
4: In optwrap(optimizer, devfun, opt$par, lower = rho$lower, control = control, :
convergence code 1 from bobyqa: bobyqa -- maximum number of function evaluations exceeded
5: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 36.9306 (tol = 0.002)
6: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge: degenerate Hessian with 1 negative eigenvalues
Alternatywnie, jeśli nie zbiegnie się po wielu iteracjach (ustawiłem na 100 000
), a następnie otrzymuję te same wyniki 50k
i 100k
oznacza to, że jest bardzo zbliżony do rzeczywistej wartości, po prostu jej nie osiąga. Czy mogę w ten sposób zgłosić moje wyniki?
Zauważ, że gdy ustawię tak wysokie iteracje, otrzymuję tylko te ostrzeżenia:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 43.4951 (tol = 0.002)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge: degenerate Hessian with 1 negative eigenvalues