Intervalle de confiance pour l'effet moyen du traitement à partir de la pondération du score de propension?

J'essaie d'estimer l'effet moyen du traitement à partir de données d'observation en utilisant la pondération du score de propension (spécifiquement IPTW). Je pense que je calcule correctement l'ATE, mais je ne sais pas comment calculer l'intervalle de confiance de l'ATE tout en tenant compte des poids de score de propension inverse.

Voici l'équation que j'utilise pour calculer l'effet moyen du traitement (référence Stat Med. 10 sept. 2010; 29 (20): 2137–2148.): Oùnombre total de sujets,état du traitement,état du résultat etscore de propension.

A T E = \frac{1}{N} \sum_{1}^{N} \frac{Z_{i} Y_{i}}{p_{i}} - \frac{1}{N} \sum_{1}^{N} \frac{(1 - Z_{i}) Y_{i}}{1 - p_{i}}

$ATE=\frac1N\sum_1^N\frac{Z_iY_i}{p_i}-\frac1N\sum_1^N\frac{(1-Z_i)Y_i}{1-p_i}$

N =

$N=$

Z_{i} =

$Z_i=$

Y_{i} =

$Y_i=$

p_{i} =

$p_i=$

Quelqu'un connaît-il un ensemble R qui calculerait l'intervalle de confiance de l'effet moyen du traitement, en tenant compte des poids? Le surveypackage pourrait-il aider ici? Je me demandais si cela fonctionnerait:

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

Je ne sais pas où aller à partir d'ici pour trouver l'intervalle de confiance de la différence entre les proportions (c'est-à-dire l'effet moyen du traitement).

r survey causality propensity-scores observational-study JJM
la source

Je ne peux pas répondre spécifiquement, mais le livre "Complex Surveys: A Guide to Analysis Using R" de l'auteur du module d'enquête couvre IPTW et peut être utile. books.google.com/…

kaz_yos

Réponses:

Vous n'avez pas besoin du surveypaquet ou de quelque chose de compliqué. Wooldridge (à partir de 2010, p. 920) «Analyse économétrique des données des sections transversales et des panneaux» a une procédure simple à partir de laquelle vous pouvez obtenir les erreurs standard afin de construire les intervalles de confiance.

$p(\textbf{x}_i,\textbf{$\gamma$})$

d_{i} = \frac{\nabla_{γ} p (x_{i}, γ)^{'} [Z_{i} - p (x_{i}, γ)]}{p (x_{i}, γ) [1 - p (x_{i}, γ)]}

$\textbf{d}_i = \frac{\nabla_\gamma p(\textbf{x}_i,\textbf{$\gamma$})'[Z_i-p(\textbf{x}_i,\textbf{$\gamma$})]}{p(\textbf{x}_i,\textbf{$\gamma$}){[1-p(\textbf{x}_i,\textbf{$\gamma$})]}}$

{ATE}_{i} = \frac{[Z_{i} - p (x_{i}, γ)] Y_{i}}{p (x_{i}, γ) [1 - p (x_{i}, γ)]}

$\text{ATE}_i = \frac{[Z_i-p(\textbf{x}_i,\textbf{$\gamma$})]Y_i}{p(\textbf{x}_i,\textbf{$\gamma$}){[1-p(\textbf{x}_i,\textbf{$\gamma$})]}}$

{\hat{ATE}}_{i}

$\widehat{\text{ATE}}_i$

{\hat{d}}_{i}

$\widehat{\textbf{d}}_i$

e_{i}

$e_i$

\sqrt{N} (\hat{ATE} - ATE)

$\sqrt{N}(\widehat{\text{ATE}} - \text{ATE})$

Var (e_{i})

$\text{Var}(e_i)$

\frac{{[\frac{1}{N} \sum_{i = 1}^{N} e_{i}^{2}]}^{\frac{1}{2}}}{\sqrt{N}}

$\frac{\left[ \frac{1}{N}\sum^N_{i=1}e_i^2 \right]^{\frac{1}{2}}}{\sqrt{N}}$

Vous pouvez ensuite calculer l'intervalle de confiance de la manière habituelle (voir par exemple les commentaires de la réponse ici pour un exemple de code). Vous n'avez pas besoin d'ajuster à nouveau l'intervalle de confiance pour les pondérations du score de propension inverse, car cette étape était déjà incluse dans le calcul des erreurs standard.

Malheureusement, je ne suis pas un gars R donc je ne peux pas vous fournir le code spécifique mais la procédure décrite ci-dessus devrait être simple à suivre. En remarque, c'est aussi la manière dont la treatrewcommande dans Stata fonctionne. Cette commande a été écrite et introduite dans le Stata Journal par Cerulli (2014) . Si vous n'avez pas accès à l'article, vous pouvez consulter ses diapositives qui décrivent également la procédure de calcul des erreurs standard à partir de la pondération du score de propension inverse. Là, il discute également de légères différences conceptuelles entre l'estimation du score de propension via logit ou probit, mais pour cette réponse, elle n'était pas trop importante, j'ai donc omis cette partie.

Andy
la source