0
私は世帯調査のミクロデータから多次元貧困指数(MPI)を作成しています。そのために私は人員数比(H)と貧困(A)の強度を必要とするので、私は何ができる: MPI = Hの* A.svyを使って多次元貧困指数(MPI)を構築する方法
私のデータベースには、そのようなものです:すでに
str(base2015)
Classes ‘tbl_df’, ‘tbl’ and 'data.frame': 117939 obs. of 33 variables:
$ X1 : int 1 2 3 4 5 6 7 8 9 10 ...
$ domicilio : num 1.1e+10 1.1e+10 1.1e+10 1.1e+10 1.1e+10 ...
$ agua : int 0 0 6 2 0 0 0 0 0 0 ...
$ ind_agua : int 0 0 1 1 0 0 0 0 0 0 ...
$ esgoto : int 1 1 6 2 1 2 3 5 3 1 ...
$ ind_cond_sanitaria: int 1 1 1 1 1 1 1 1 1 1 ...
...
$ pre_wgt : int 200 200 200 200 200 200 200 200 200 200 ...
$ one : int 1 1 1 1 1 1 1 1 1 1 ...
$ region : int 1 1 1 1 1 1 1 1 1 1 ...
$ c : num 0.2 0.2 0.2 0.3 0.3 0.3 0.2 0.2 0.1 0.1 ...
$ q : num 0 0 0 0 0 0 0 0 0 0 ...
私は調査パッケージを使用してHとAを行いました。
#proporção de pobres para região (H)
svyby(~ q , ~ region , n , vartype=c("se","cv","ci"), svymean)
と、この:
# A by Região
# 1-norte 2-nordeste 3-sudeste 4-sul 5-centro-oeste
svyby(~ c , ~ region , subset_q , vartype=c("se","ci","cv"), svymean)
しかし、私は、アカウントへのサンプルデザインを取るMPI自体(MPI = Hの*のA)を計算する方法がわかりません。 誰かがこれを手伝ってくれますか?
PS:もう一度試してみましたが、わかりません。
#MPI
svyby(~ q , ~ region , n , vartype=c("se","cv","ci"), svymean)*svyby(~ c , ~ region , subset_q , vartype=c("se","ci","cv"), svymean)