我要回到几个月前这个问题解决的话题。
我有一个 coxph 模型,其中包括一个具有三个水平的因子协变量(分类变量)。
当我使用 cox.zph 测试违反比例风险假设的模型时,它表明因子的一个水平违反了假设。
通常,我会将违规协变量与时间函数进行交互,但在这种情况下 - 因为我正在处理一个因素 - 我不确定这是否有意义。
在任何情况下,尝试将因子与时间函数交互都会返回错误。
你能建议做什么吗?有没有办法通过与时间的相互作用来“瞄准”违规因素水平?谢谢。
我要回到几个月前这个问题解决的话题。
我有一个 coxph 模型,其中包括一个具有三个水平的因子协变量(分类变量)。
当我使用 cox.zph 测试违反比例风险假设的模型时,它表明因子的一个水平违反了假设。
通常,我会将违规协变量与时间函数进行交互,但在这种情况下 - 因为我正在处理一个因素 - 我不确定这是否有意义。
在任何情况下,尝试将因子与时间函数交互都会返回错误。
你能建议做什么吗?有没有办法通过与时间的相互作用来“瞄准”违规因素水平?谢谢。
嗯,有趣 - 时间(或其函数)与一个因素的交互应该没问题恕我直言,例如参见http://www.ats.ucla.edu/stat/examples/asa/test_proportionality.htm的以 R 为中心的位
...但我假设您正在创建一个新的“真实年龄”变量以与之交互,即通过将 DX 年龄添加到时间。我不确定出现在生存对象和公式右侧的“时间”变量会发生什么。
您可以发布一些代码(和错误消息)以更好地帮助您.....