Konfidenzintervall für den durchschnittlichen Behandlungseffekt aus der Neigungsbewertung?

9

Ich versuche, den durchschnittlichen Behandlungseffekt anhand von Beobachtungsdaten mithilfe der Neigungsbewertung (speziell IPTW) abzuschätzen. Ich denke, ich berechne die ATE korrekt, aber ich weiß nicht, wie ich das Konfidenzintervall der ATE unter Berücksichtigung der inversen Neigungsbewertungsgewichte berechnen soll.

Hier ist die Gleichung, die ich zur Berechnung des durchschnittlichen Behandlungseffekts verwende (Referenz Stat Med. 10. September 2010; 29 (20): 2137–2148.): Wobei Gesamtzahl der Probanden, Behandlungsstatus, Ergebnisstatus und .

EINT.E.=1N.1N.Z.ichY.ichpich- -1N.1N.(1- -Z.ich)Y.ich1- -pich
N.=Z.ich=Y.ich=pich=

Kennt jemand ein R-Paket, das das Konfidenzintervall des durchschnittlichen Behandlungseffekts unter Berücksichtigung der Gewichte berechnet? Könnte das surveyPaket hier helfen? Ich habe mich gefragt, ob das funktionieren würde:

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

Ich weiß nicht, wohin ich von hier aus gehen soll, um das Konfidenzintervall der Differenz zwischen den Anteilen (dh den durchschnittlichen Behandlungseffekt) zu ermitteln.

JJM
quelle
Ich kann nicht spezifisch antworten, aber das Buch "Komplexe Umfragen: Ein Leitfaden zur Analyse mit R" des Autors des Umfragepakets behandelt IPTW und kann hilfreich sein. books.google.com/…
kaz_yos

Antworten:

11

Sie brauchen das surveyPaket oder etwas Kompliziertes nicht. Wooldridge (ab 2010, S. 920) "Ökonometrische Analyse von Querschnitts- und Paneldaten" verfügt über ein einfaches Verfahren, mit dem Sie die Standardfehler erhalten können, um die Konfidenzintervalle zu erstellen.

p(xich,γ)

dich=γp(xich,γ)'[Z.ich- -p(xich,γ)]]p(xich,γ)[1- -p(xich,γ)]]
ASSich=[Z.ich- -p(xich,γ)]]Y.ichp(xich,γ)[1- -p(xich,γ)]]
ASS^ichd^icheichN.(ASS^- -ASS)Var(eich)
[1N.ich=1N.eich2]]12N.

Sie können dann das Konfidenzintervall auf die übliche Weise berechnen (siehe beispielsweise die Kommentare zur Antwort hier für ein Codebeispiel). Sie müssen das Konfidenzintervall für die inversen Neigungsbewertungsgewichte nicht erneut anpassen, da dieser Schritt bereits in die Berechnung der Standardfehler einbezogen wurde.

Leider bin ich kein R-Typ, daher kann ich Ihnen den spezifischen Code nicht zur Verfügung stellen, aber das oben beschriebene Verfahren sollte einfach zu befolgen sein. Nebenbei bemerkt, dies ist auch die Art und Weise, wie der treatrewBefehl in Stata funktioniert. Dieser Befehl wurde von Cerulli (2014) im Stata Journal geschrieben und eingeführt . Wenn Sie keinen Zugriff auf den Artikel haben, können Sie seine Folien überprüfen, in denen auch die Berechnung der Standardfehler aus der Gewichtung der inversen Neigungsbewertung beschrieben wird. Dort erörtert er auch einige geringfügige konzeptionelle Unterschiede zwischen der Schätzung des Neigungsscores über Logit oder Probit, aber für diese Antwort war es nicht übermäßig wichtig, so dass ich diesen Teil weggelassen habe.

Andy
quelle