Hi!
I am trying to do a glmer.nb but get this error:
Error in eval(expr, envir, enclos) :
..2 used in an incorrect context, no ... to look in
My data structure is data.frame
Here my data set: d1_2
Year ID Age Reproductive_status Rank_Residuals Asso.Y1 TotalY1
1994 109 8 Lactating 0.23947902 9 48
1994 110 6 Barren -0.030977951 56 62
1995 1000 9 Summer_yeld 0.08624 31 55
1995 103 4 Barren -0.06861 29 52
1996 23 8 Barren -0.015726749 64 64
1997 1000 11 Lactating 0.148665859 5 35
1997 103 6 Barren -0.116539886 51 51
2001 105 9 Lactating -0.036222842 27 28
2001 106 9 Barren 0.073777158 40 40
2001 116 8 Barren -0.015710059 5 31
library(lme4); library(MASS)
nb.glmer <- glmer.nb(Asso.Y1 ~ Reproductive_status * Age +
Reproductive_status * Rank_Residuals + Reproductive_status * I(Age^2) +
(1|ID) + (1|Year), offset = log(d1_2$TotalY1), data=d1_2)
I have also tried because I found that glmer.nb is not well understood:
library(glmmADMB)
random <- MCMCglmm(cbind(Asso.Y1,TotalY1-Asso.Y1) ~ Reproductive_status
* Age + Reproductive_status * Rank_Residuals + Reproductive_status *
I(Age^2), random = ~ ID+Year, data = d1_2)
BUT: Error in MCMCglmm(cbind(Asso.Y1, TotalY1 - Asso.Y1) ~
Reproductive_status * :
family must have the same length as the number of responses
I am not able to solve those problems and I have looked on forums too.
Tkx in advance!
--
*Karina Charest Castro*
Candidate à la maitrise en biologie - Masters student
Laboratoire de Steeve Côté
Département de biologie, Centre d'études nordiques
Université Laval
Pavillon Alexandre-Vachon
1045, av. de la Médecine, local 3076
Québec (Qc) G1V 0A6
Tel: 418-656-2131 poste 8152
[[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.
