1

尝试在 R 中使用 Weibull 分布运行生存分析时,我遇到了上述错误消息。

我的数据有点棘手,因为它包含左右删失的观察结果。

我按照https://stat.ethz.ch/R-manual/R-devel/library/survival/html/Surv.html上的说明进行操作:

“间隔删失数据可以用两种方式表示。对于第一种使用 type = "interval" 和上面显示的代码。在这种用法中,除非 event=3,否则 time2 参数的值将被忽略。第二种方法是考虑每个观察作为时间间隔,左删失为 (-infinity, t),右删失为 (t, infinity),精确为 (t,t),间隔为 (t1, t2)。这是用于类型 =区间2。无穷大值可以用实际无穷大 (Inf) 或 NA 来表示。已证明第二种形式更有用。

我决定采用第二种形式,并使用以下标题对我的数据进行编码:

t1, t2 (开始和结束时间), status (0,1,2), factor1, factor2

为了运行我的模型,我用 R 编写:

    model1 <- survreg(Surv(t1, t2, status, type='interval2')~factor(factor1) + factor(factor2), dist='weibull', data=data)

但是,它一直给我错误消息:

    Error in Surv(paduration, status, type = "interval") : Wrong number of args for this type of survival data

我不知道我做错了什么,因为 surv() 对象只有 5 个可能的参数,而 origin 参数与我的情况无关。当我尝试使用较少的参数时,它仍然给我同样的错误信息。我尝试改用“type=interval”,但仍然收到相同的错误消息。

请帮忙。

4

0 回答 0