[R] weights error in svyglm

Milan Bouchet-Valat nalimilan at club.fr
Fri Apr 11 17:36:52 CEST 2014


Le vendredi 11 avril 2014 à 13:38 +0000, Lee, Laura a écrit :
> Hi!
> 
> I am just learning the 'svyglm' package and have run into an error for
> which I have not found a solution. My data have been collected from a
> stratified random survey. Here is the code:
> 
> NB3<-glm.nb(COLNUM~Year+Depth+MESH+offset(LogEffort),data=data)
> 
> dstrat <- svydesign(id=~1, strata=~STRATA, weights = ~weight,
> data=data)
> 
> SNB3 <- svyglm(NB3, design=dstrat)
> 
> The error that is returned is:
> 
> Error in model.frame.glm(formula = list(coefficients =
> c(-6.96858807641624,  :
>   object '.survey.prob.weights' not found
> 
> I would appreciate any assistance in solving this problem.
The error is somewhat strange, but you should not pass the result of
glm.nb() to svyglm (where have you found this idea?). Just do
SNB3 <- svyglm(COLNUM ~ Year + Depth + MESH + offset(LogEffort), design=dstrat)


Regards




More information about the R-help mailing list