merMod
object create an Rfunction that
takes a single argument, which is the new parameter value, and
returns the deviance.mkdevfun(rho, nAGQ = 1L, maxit = 100, verbose = 0, control = list())
pp
, a
prediction module, typically of class
merPredD
and resp
, a response
module, e.g., of class lmerResp
.lmerControl
(tolPwrss
and compDev
for GLMMs,
tolPwrss
for NLMMs)function
of one numeric argument.mkdevfun
evaluates the
deviance of the model represented by the predictor
module, pp
, and the response module, resp
. For lmer
model objects the argument of the
resulting function is the variance component parameter,
theta
, with lower bound. For glmer
or
nlmer
model objects with nAGQ = 0
the
argument is also theta
. However, when nAGQ > 0,
the argument is c(theta, beta)
.
lmer
, glmer
and
nlmer
(dd <- lmer(Yield ~ 1|Batch, Dyestuff, devFunOnly=TRUE))
dd(0.8)
minqa::bobyqa(1, dd, 0)
Run the code above in your browser using DataLab