मैं प्रिवेंशन स्कोर वेटिंग (विशेष रूप से आईपीटीडब्ल्यू) का उपयोग करके अवलोकन डेटा से औसत उपचार प्रभाव का अनुमान लगाने की कोशिश कर रहा हूं। मुझे लगता है कि मैं सही ढंग से ATE की गणना कर रहा हूं, लेकिन मुझे पता नहीं है कि व्युत्क्रम स्कोर के वज़न को ध्यान में रखते हुए ATE के आत्मविश्वास अंतराल की गणना कैसे करें।
यहां वह समीकरण है जिसका उपयोग मैं औसत उपचार प्रभाव (संदर्भ स्टैट मेड। 10 सितंबर, 2010; 29 (20): 2137-2148।): गणना के लिए कर रहा हूं। जहां विषयों की कुल संख्या, Z_i = उपचार की स्थिति, Y_i = परिणाम स्थिति, और p_i = प्रवृत्ति स्कोर।
क्या कोई ऐसे आर पैकेज के बारे में जानता है जो वजन को ध्यान में रखते हुए, औसत उपचार प्रभाव के विश्वास अंतराल की गणना करेगा? क्या survey
पैकेज यहां मदद कर सकता है? मैं सोच रहा था कि क्या यह काम करेगा:
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
मुझे नहीं पता कि अनुपात के अंतर के अंतर अंतराल (यानी औसत उपचार प्रभाव) को खोजने के लिए यहां से कहां जाना है।