Anteil der Armen

Question

Ich versuche, einen Anteil von Armen mit dem Umfragepaket zu berechnen. Also erstelle ich eine Teilmenge mit der Variablen der Gesamtbevölkerung (eins) und n genannt (ich bin mir nicht sicher, ob es der richtige Weg ist, dies zu tun). So , ich habe:

str(base2015$q) 
num [1:117939] 0 0 0 0 0 0 0 0 0 0 ... 
str(base2015$one) 
int [1:117939] 1 1 1 1 1 1 1 1 1 1 ... 

und eine kategorische Variable Region. (y ist mein Svydesign) Ich habe versucht, die SvyRatio-Funktion so zu verwenden, aber er gibt mir das Ergebnis für alle Datensätze, und ich will nach Untergruppe (Region). In der Tat mag ich dies für jede Region (das Bild sehen) proportion of poor

svyratio(numerator=~q, 
     ~one, 
     design = y, separate=FALSE, na.rm=TRUE, 
     formula, covmat=FALSE, deff=FALSE) 

Ich habe auch versucht, eine Teilmenge erstellen mit der Bevölkerung, und dann die svyby verwenden:

n <- subset(y , one == 1) 
View(n) 
svyby(~ q , ~ region , n , vartype=c("se","ci","cv"), svytotal) 
svyby(~ q , ~ region , n , vartype=c("se","ci","cv"), svymean) 
    region q se ci_l ci_u cv 
1  1 NA NaN NaN NaN NaN 
2  2 NA NaN NaN NaN NaN 
3  3 NA NaN NaN NaN NaN 
4  4 NA NaN NaN NaN NaN 
5  5 NA NaN NaN NaN NaN 

Er gib mir diesen. Ich weiß nicht, was ich falsch mache. Kann mir jemand dabei helfen?

Answer 1

data(scd) 

scddes<-svydesign(data=scd, prob=~1, id=~ambulance, strata=~ESA,nest=TRUE, fpc=rep(5,6)) 

svyby(~ alive , denominator = ~ arrests , by = ~ ambulance , scddes , svyratio) 

# which missings do you want to remove? read `?svyby` 
svyby(~ alive , denominator = ~ arrests , by = ~ ambulance , scddes , svyratio , na.rm = TRUE , na.rm.by = TRUE , na.rm.all = TRUE)