如何处理R中的重复权重等于零？

Question

我对某些复制和采样权重的值为零。因此，当我从“调查”包中使用Svycoxph时，我收到一条错误消息“包无效权重，必须> 0”。我认为一种方法可能是排除这些观察结果。我想知道是否有办法保留Cox比例风险模型的观察结果？

谢谢！ 朱莉娅

Answer 1

零复制权重是完全正常的 - 随着折刀重复，它们将被PSU的观察结果遗漏;通过bootstrap复制，它们将来自在特定重采样中出现零次的PSU。另一方面，零采样权重并没有多大意义，可能表明样本中没有的观察结果或者不在采样框架中（个人而言，我不是这样）将这些代码编码为NA）

coxph()包中的survival功能无法处理负值或零权重，以及svycoxph()调用的内容。对于重复权重，svycoxph()会向权重添加一个小值（1e-10）以停止零问题。

但是，svycoxph()无法处理零采样权重。您可能希望在构造设计对象之前删除它们。