[R-sig-ME] Problem_Downdated VtV is not positive definite

C. AMAL D. GLELE @|te@@ed@c2 @end|ng |rom gm@||@com
Sun Jan 19 14:29:54 CET 2020


Hello, everyone.
In the models below, mydata has 32 rows
village and Year are grouping factor variables with respectively 8 and 4
levels.
1)
model1.1 fits well
model1.1<-glmer(cbind(F,L)~village+Year+(1|Year)+(1|village),
family="binomial",data=mydata)
But, when I replace village+Year by village*Year, I get the warning and
error message below
model1.2<-glmer(cbind(F,L)~village*Year+(1|Year)+(1|village),
family="binomial",data=mydata)
Warning:
fixed-effect model matrix is rank deficient so dropping 2 columns /
coefficients
Error message:
Error in (function (fr, X, reTrms, family, nAGQ = 1L, verbose = 0L, maxit =
100L,  :
  Downdated VtV is not positive definite
I noticed a similar behavior when trying with glmmadmb:
model2.1 fits well
model2.1<-glmmadmb(cbind(F,L)~village+Year+(1|Year)+(1|village),family="binomial",data=na.omit

(mydata))
But, with village*Year, Iget the error message below:
model2.2<-glmmadmb(cbind(F,L)~village*Year+(1|Year)+(1|village),family="binomial",data=na.omit

(mydata))
Error:
Error in glmmadmb(cbind(F, L) ~ village * Year + (1 | Year) + (1 |
village),  :
  rank of X = 30 < ncol(X) = 32
##
in advance, thanks for helping solve these issues(I'm very interested in
interaction terms).
Kinds regards,

	[[alternative HTML version deleted]]



More information about the R-sig-mixed-models mailing list