Hello!
I am new in the mailing list for R help and I hope to be able to formulate a
good question easy to understand.
After trying to do many different things to solve an error I do not find the
solution.
I am modeling my data set with hurdle negative binomial mixed effects, to find
the correlation of some bird counts with environmental (categorical and
continuous) variables.
When I run different models I have always an error. For instance:
-For the truncated modeling of the non-zero counts:
> HURgvgsw <- glmmadmb(count~ GVG*sward + (1|week),
> data=subset(SW_GVG,count>0), + family= "truncnbinom")
Error en glmmadmb(count ~ GVG * sward + (1 | week), data = subset(SW_GVG, :
rank of X = 6 < ncol(X) = 10
-Or the binomial part where the zeros are modeled with the non-zeros:
> HURgvgsw <- glmmadmb(count~ sward*GVG + (1|week) + (1|cluster), data=SW_GVG,
> family= "binomial")
Error en glmmadmb(count ~ sward * GVG + (1 | week) + (1 | cluster), data =
SW_GVG, :
rank of X = 13 < ncol(X) = 15
Would you have the solution to this?
When I run the model with many variables I also have that error. So what I have
been done so far is to split the model in to two differentiated by smaller set
of variables as follows so I do not get that error. I do not think this is the
solution though:
> hnb6 <- glmmadmb(count~ Year * (dist + sward + inject + fresh + veg) +
> (1|week), data=subset(newdata_,count>0), family= "truncnbinom")
Error en glmmadmb(count ~ Year * (dist + sward + inject + fresh + veg) + :
rank of X = 22 < ncol(X) = 24
However this model works (taking off one variable):
> hnb6 <- glmmadmb(count~ Year * (dist + sward + inject + veg) + (1|week),
> data=subset(newdata_,count>0),
+ family= "truncnbinom")
The same with no variable interactions (controling for week, nested in Year):
> hnb6 <- glmmadmb(count~ dist + sward + inject + fresh + veg + (1|Year/week),
> data=subset(newdata_,count>0),
+ family= "truncnbinom")
Error en glmmadmb(count ~ dist + sward + inject + fresh + veg + (1 |
Year/week), :
rank of X = 11 < ncol(X) = 12
When I do the binomial part of the model above, I get another error very
frequent as well:
> hnb6 <- glmmadmb(count~ dist + sward + inject + fresh + veg + (1|Year) +
> (1|week), data=newdata_,
+ family= "binomial")
Parameters were estimated, but not standard errors were not: the most likely
problem is that the curvature at MLE was zero or negative
Error en glmmadmb(count ~ dist + sward + inject + fresh + veg + (1 | Year) + :
The function maximizer failed (couldn't find STD file) Troubleshooting steps
include (1) run with 'save.dir' set and inspect output files; (2) change run
parameters: see '?admbControl'
Además: Mensajes de aviso perdidos
comando ejecutado 'C:\Windows\system32\cmd.exe /c
"C:/Users/Marta/Documents/R/win-library/2.15/glmmADMB/bin/windows32/glmmadmb.exe"
-maxfn 500 -maxph 5 -noinit -shess' tiene estatus 1
I hope it is easy to understand and you could help me with it.
Thanks a lot!!
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
