Displaying 1 result from an estimated 1 matches for "setomega".
Did you mean:
geomega
2007 Jan 26
0
R crash with modified lmer code
...Zt <- if (length(Ztl) == 1)
Ztl[[1]]
else do.call("rbind", Ztl)
fl <- FL$fl
if (fltype < 0) {
mer <- .Call(mer_create, fl, Zt, X, as.double(Y), method ==
"REML", nc, cnames)
if (!is.null(start))
mer <- setOmega(mer, start)
.Call(mer_ECMEsteps, mer, cv$niterEM, cv$EMverbose)
LMEoptimize(mer) <- cv
return(new("lmer", mer, frame = if (model) fr$mf else
data.frame(),
terms = mt, call = mc))
}
if (method %in% c("ML", "REML"))...