Estoy tratando de estimar el efecto promedio del tratamiento a partir de datos de observación utilizando la ponderación de puntaje de propensión (específicamente IPTW). Creo que estoy calculando el ATE correctamente, pero no sé cómo calcular el intervalo de confianza del ATE teniendo en cuenta los pesos del puntaje de propensión inversa.
Aquí está la ecuación que estoy usando para calcular el efecto promedio del tratamiento (referencia Stat Med. 10 de septiembre de 2010; 29 (20): 2137–2148.): Donde número total de sujetos, estado del tratamiento, estado del resultado y puntaje de propensión.
¿Alguien sabe de un paquete R que calcule el intervalo de confianza del efecto promedio del tratamiento, teniendo en cuenta los pesos? ¿Podría survey
ayudar el paquete aquí? Me preguntaba si esto funcionaría:
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
No sé a dónde ir desde aquí para encontrar el intervalo de confianza de la diferencia entre las proporciones (es decir, el efecto promedio del tratamiento).