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. Thanks, Laura [[alternative HTML version deleted]]
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
Thanks! That worked. However, what error distribution is assumed? I was hoping to use a negative binomial. -----Original Message----- From: Milan Bouchet-Valat [mailto:nalimilan@club.fr] Sent: Friday, April 11, 2014 11:37 AM To: Lee, Laura Cc: r-help@r-project.org Subject: Re: [R] weights error in svyglm 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 ______________________________________________ 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.