眼镜

CRAN_状态_标签 CRAN_下载_徽章

光学重量包含用于估计权重的函数将处理平衡到给定的平衡阈值。它解出了一个二次方程最小化权重目标函数的规划问题使用solve_osqp()在中开放式服务质量计划包裹。这个是Zubizarreta(2015)中描述的方法。眼镜将该方法扩展到多项式、连续和纵向并提供简单的用户界面和与这个包裹。

下面是使用眼镜并评估协变量的平衡.

开发工具时间:安装github(“ngreifer/optweight”)#开发版本
图书馆(“optweight”)
图书馆(“钴”)
数据(“拉隆德”)

#估计权重
现在<- 眼镜(治疗~年龄+教育+比赛+节点度+已婚的+
第74页+re75标准+ (第74版== 0)+ (参考75== 0),
                数据=拉隆德,估计= “收件人”,托尔斯=.01)
现在
optweight对象-obs.数量:614-采样权重:无-治疗:2类-估计需求:ATT(焦点:1)-协变量:年龄、教育程度、种族、nodegree、已婚、re74、re75、I(re74==0)、I(re 75==0
总结(现在)
重量汇总:-重量范围:最小值最大值处理1||1.0000控件0 |-------------------------|7.4319-按组列出的最大重量为5的单位:2      3      4      5      6处理过的1 1 1 1608    574    559    573    303控制7.2344 7.3161 7.4058 7.4058 4.4319Var平均系数Abs Dev处理0.0000 0.0000控制1.9019 1.3719总体1.5897 0.9585-有效样本量:对照治疗组未加权429.000 185重量92.917 185
平衡标签(现在)
呼叫optweight(公式=治疗~年龄+教育+种族+无障碍+已婚+re74+re75+I(re74==0)+I(re 75==0,tols=0.01,estimand=“ATT”)平衡措施类型差异调整年龄控制。0.01导出Contin。0.01race_black二进制0.01race_hispan二进制0.00race_white二进制-0.01nodegree二进制0.01已婚二进制-0.01re74康定。0.01re75容器。0.01I(re74==0)二进制0.01I(re75==0)二进制0.01有效样本量对照治疗组未调整429.000 185调整92.917 185
#估计治疗效果
图书馆(“jtools”)
总和(勒姆(参考78~治疗,数据=拉隆德,重量=现在$重量),
     限制(confint)= 真的,坚固耐用的= 真的,模型.拟合= 错误的,
     model.info(型号.info)= 错误的)

美国东部时间。

2.5%

97.5%

t值。

第页

(截距)

5342.94

4058.75

6627.13

8.17

0

对待

1006.20

-710.34

2722.74

1.15

0.25

标准误差:稳健,类型=HC3

低级函数optwight.合身在上操作协变量和治疗变量。

除了估计用于估计治疗的平衡权重影响,光学重量可以估计采样权重将估计推广到由协变量定义的新目标人群使用函数的力矩女装.svy.