Displaying 1 result from an estimated 1 matches for "chaidata".
Did you mean:
caudata
2010 Apr 06
0
Strange error
Someone just sent me a data set that causes the lmekin function, part of
the kinship package, to fail. In chasing it down I get an error I have
never seen before.
fit <- lmekin(icam1 ~ factor(center) + age + factor(sex),
random= ~1|iid, data=chaidata, varlist=kmat)
Error in Y - fitted : non-numeric argument to binary operator
Add the recover option, and the offending lines are
fitted <- c(X %*% lfit$coef) #fitted, on the original scale
residuals <- Y - fitted
> options(error=recover)
> fit <- lmekin(icam1 ~ factor(cent...