Interval kepercayaan untuk efek pengobatan rata-rata dari bobot skor kecenderungan?

9

Saya mencoba memperkirakan efek pengobatan rata-rata dari data pengamatan menggunakan bobot skor kecenderungan (khususnya IPTW). Saya pikir saya menghitung ATE dengan benar, tetapi saya tidak tahu bagaimana menghitung interval kepercayaan ATE sambil memperhitungkan bobot skor kecenderungan terbalik.

Inilah persamaan yang saya gunakan untuk menghitung efek pengobatan rata-rata (referensi Stat Med. 10 Sep 2010, 29 (20): 2137–2148.): Di mana jumlah total subjek, status perawatan, status hasil, dan skor kecenderungan. N=Zi=Yi=pi=

ATE=1N1NZiYipi1N1N(1Zi)Yi1pi
N=Zi=Yi=pi=

Adakah yang tahu tentang paket R yang akan menghitung interval kepercayaan dari efek perawatan rata-rata, dengan mempertimbangkan bobotnya? Bisakah surveypaket membantu di sini? Saya bertanya-tanya apakah ini akan berhasil:

library(survey)
sampsvy=svydesign(id=~1,weights=~iptw,data=df)
svyby(~surgery=='lump',~treatment,design=sampsvy,svyciprop,vartype='ci',method='beta')

#which produces this result:
  treatment surgery == "lump"      ci_l      ci_u
   No         0.1644043 0.1480568 0.1817876
   Yes         0.2433215 0.2262039 0.2610724

Saya tidak tahu harus ke mana dari sini untuk menemukan interval kepercayaan dari perbedaan antara proporsi (yaitu efek pengobatan rata-rata).

JJM
sumber
Saya tidak bisa menjawab secara spesifik, tetapi buku "Survei Kompleks: Panduan untuk Analisis Menggunakan R" oleh penulis paket survei ini mencakup IPTW, dan mungkin bisa membantu. books.google.com/…
kaz_yos

Jawaban:

11

Anda tidak perlu surveypaket atau apapun yang rumit. Wooldridge (2010, hal. 920 dan seterusnya) "Analisis Ekonometrik dari Potongan Melintang dan Data Panel" memiliki prosedur sederhana yang darinya Anda dapat memperoleh kesalahan standar untuk membuat interval kepercayaan.

Di bawah asumsi bahwa Anda telah menentukan skor kecenderungan dengan benar yang kami nyatakan sebagai , tentukan skor dari estimasi skor kecenderungan (mis. Logit atau regresi probit pertama Anda) ) sebagai dan biarkan seperti yang Anda miliki dalam ekspresi Anda di atas. Kemudian ambil analog sampel dari dua ekspresi ini dan mundur onp(xi,γ)

di=γp(xi,γ)[Zip(xi,γ)]p(xi,γ)[1p(xi,γ)]
^ ATE i d iei
ATEi=[Zip(xi,γ)]Yip(xi,γ)[1p(xi,γ)]
ATE^id^i. Pastikan Anda menyertakan intersep dalam regresi ini. Biarkan menjadi residu dari regresi itu, maka varians asimptotik dari hanyalah . Jadi kesalahan standar asimptotik ATE Anda adalah eiVar(ei) [ 1N(ATE^ATE)Var(ei)
[1Ni=1Nei2]12N

Anda kemudian dapat menghitung interval kepercayaan dengan cara biasa (lihat misalnya komentar untuk jawaban di sini untuk contoh kode). Anda tidak perlu menyesuaikan interval kepercayaan lagi untuk bobot skor kecenderungan terbalik karena langkah ini sudah termasuk dalam perhitungan kesalahan standar.

Sayangnya saya bukan seorang cowok R jadi saya tidak bisa memberi Anda kode khusus tetapi prosedur yang diuraikan di atas harus langsung diikuti. Sebagai catatan tambahan, ini juga cara di mana treatrewperintah di Stata bekerja. Perintah ini ditulis dan diperkenalkan di Stata Journal oleh Cerulli (2014) . Jika Anda tidak memiliki akses ke artikel Anda dapat memeriksa slide- nya yang juga menguraikan prosedur menghitung kesalahan standar dari bobot skor kecenderungan terbalik. Di sana ia juga membahas beberapa perbedaan konseptual kecil antara memperkirakan skor kecenderungan melalui logit atau probit tetapi demi jawaban ini itu tidak terlalu penting sehingga saya menghilangkan bagian ini.

Andy
sumber