Displaying 1 result from an estimated 1 matches for "base_nps".
Did you mean:
base_npc
2008 Oct 22
1
Package survey
...yby for total statistics, for example:
svyby(~p_igov,~div_a,desenho_nps,svytotal,drop.empty.groups=TRUE,vartype
=c("se","var","cvpct"))
In the numerical variable p_igov (and others) I have many non responses
but if I maintain the NA it doesn’t work.
summary(base_nps$p_igov)
Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
1.0 9.0 51.0 208.7 185.0 19470.0 1687.0
desenho_nps <-
svydesign(id=~npc,strata=~estr_final,weights=~pond_nps,fpc=~nps_univ,nes
t=T,data=base_nps)
svyby(~p_igov,~div_a,desenho_nps,svytotal,drop.empty.grou...