Displaying 2 results from an estimated 2 matches for "mer_coefgets".
2006 Mar 16
1
lme4/Matrix: Call to .Call("mer_update_y"...) and LMEoptimize gives unexpected side effect...
...e <- function(x, value)
{
if (value$msMaxIter < 1) return(x)
nc <- x at nc
constr <- unlist(lapply(nc, function(k) 1:((k*(k+1))/2) <= k))
fn <- function(pars)
deviance(.Call("mer_coefGets", x, pars, 2, PACKAGE = "Matrix"))
gr <- if (value$gradient)
function(pars) {
if (!isTRUE(all.equal(pars,
.Call("mer_coef", x,...
2006 Jun 14
1
lmer and mixed effects logistic regression
...tions are not
paired). For example:
id response age
1 0 30
1 0 30
2 1 55
2 1 55
3 0 37
4 1 52
5 0 39
5 0 39
etc.
I get the following error:
> (lmer(response~(1|id)+age,data=gdx,family=binomial))
Error in deviance(.Call("mer_coefGets", x, pars, 2, PACKAGE =
"Matrix")) :
Leading minor of order 2 in downdated X'X is not positive
definite
Similar problem if I use quasibinomial.
If I use glm, of course it thinks I have roughly twice the number of
subjects so the standard errors would be smaller than the...