貧困層の割合

Question

調査パッケージを使用して貧困層の割合を計算しようとしています。 したがって、人口（1）とn（私はそれが正しい方法であるかどうかはわかりません）の合計の変数でサブセットを作成します。 だから、私は持っている：

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 ... 

とカテゴリ変数領域。 （yは私のsvydesignです） 私はそのようなsvyratio関数を使用しようとしましたが、彼は私にすべてのデータセットの結果を返し、私はサブグループ（region）で欲しいです。

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 

彼は私にこれを与える：行為で私は各地域のためにこれをしたい（画像を参照） proportion of poor

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

私はまた、集団を用いてサブセットを作成し、svybyを使用してみました。私は何が間違っているのか分からない。 誰かがこれを手伝ってくれますか？

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)