使用ipw计算权重的无限值错误

Question

我正在尝试使用R中的 ipw 包估计倾向评分模型的一组稳定的反概率权重。我有一个包含34517个观测值的数据框，其中155个在我的治疗中组。

mfg_stabilized_full <- ipwpoint(exposure = pmd_dummy, family = "binomial", link = "logit",
                            numerator = ~1,
                            denominator = ~ mfgshare + owner_per + dist_km + network_density,
                            data = city_lehd_acs04)

运行时，我收到一条错误消息，说明权重中有NA。当我尝试通过从 ipwpoint（）中删除分子项来计算不稳定权重时，返回的权重是无限的。

在两次调用中，基础逻辑模型收敛并返回每个协变量的估计值。考虑到基础逻辑模型的完成，为什么我会继续获得不稳定权重和NA的无限值以求稳定？

Answer 1

检查曝光代码。根据包装手册，二进制曝光应编码为0和1。将曝光编码为一个因素时，我遇到了同样的错误。