hi, could you create a reproducible example starting from http://asdfree.com/pesquisa-nacional-de-saude-pns.html ? thanks
On Mon, Dec 4, 2017 at 9:56 AM, Luciane Maria Pilotto via R-help < r-help@r-project.org> wrote: > Hi, > I am trying to run analyzes incorporating sample weight, strata and > cluster (three-stage sample) with PNS data (national health survey) and is > giving error. I describe below the commands used. I could not make the code > reproducible properly. > Thanks, > ################################################# > library(survey)####change to 0 and 1 variable outcomedent2<-ifelse( > consdentcat2==2,0,1)table(dent2) > dent2<-as.factor(dent2)str(dent2)reg<-cbind(reg, > dent2) > #tchange to factor str(sexo)reg$sexo <- as.factor((reg$sexo)) > ################# > pns2013design<-svydesign(id=~upa, nest=TRUE, strata = estrato, weight = > peso, data = reg) > > PNS<-svyglm(dent2~sexo,design=pns2013design, method="logistic", data = > reg) > > Error in logistic(x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > : unused arguments (x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, > 1, 1, 1, 1, 1, 1,############################ > ############################################################ > ######################################################################## > > ###NOT running SVYGLM command#data for R-helpbteste <- reg[1:10, c(1, 4, > 26, 27, 28)]#dput(bteste)###omitting most of the lines "9967.223281", > "9967.870849", "9968.207979", "997.0045537", "997.1451224", > "997.1574275", "997.1782368", "997.1812338", "997.2480231", > "997.2531051", "997.3803145", "997.4345619", "997.4363662", > "997.5344599", "997.5964572", "997.7439813", "997.8360773", > "997.8770935", "997.8811374", "997.9147096", "997.9563562", > "9974.30392", "9974.344482", "9975.656981", "9977.382263", > "9979.999691", "998.0053953", "998.1069038", "998.2140192", > "998.2655421", "998.308316", "998.3090242", "998.3579509", > "998.3656231", "998.3766007", "998.6844831", "998.7030027", > "998.7112321", "998.8021132", "998.8839799", "998.9225688", > "998.9270228", "998.9337225", "9983.555066", "9985.353117", > "9989.517638", "999.0713699", "999.0771916", "999.1021413", > "999.1779133", "999.2539765", "999.3435971", "999.3809978", > "999.6348707", "999.7597985", "999.8002746", "999.8819267", > "999.8821907", "999.8921074", "999.9211427", "9991.102816", > "9991.440035", "9994.626337", "9994.723654", "9996.637923", > "9998.491819"), class = "factor")), .Names = c("consdentcat2", "sexo", > "estrato", "upa", "pesomorcc"), row.names = c(NA, 10L), class = > "data.frame")################################################## > bteste1 <-bteste[1:10, ]#bteste1 consdentcat2 sexo estrato upa > pesomorcc1 2 1 1110011 1100002 418.76819022 2 > 1 1110011 1100002 317.13175793 2 1 1110011 1100002 > 467.09452884 1 1 1110011 1100002 209.38409515 2 > 1 1110011 1100002 209.38409516 2 1 1110011 1100002 > 418.76819027 2 1 1110011 1100002 233.54726448 2 > 1 1110011 1100002 628.15228539 2 1 1110011 1100002 > 317.131757910 2 2 1110011 1100002 321.5014524> > _________________________________________ > Luciane Maria Pilotto > > > > | | Livre de vĂrus. www.avast.com. | > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.