0

我想执行 CLAD 回归

  • y = EQ-5D-5L 效用分数（以 1.0 为界）
  • x = 各种患者特征

我已经发现我需要在库QUANTREG中使用CRQ，但到目前为止我无法弄清楚具体细节。我的问题是：

  1. 我需要使用鲍威尔方法吗？
  2. 如果是这样，如果我没有时间变量而是 0/1 审查变量，我该如何指定“yc”（审查时间）？

这是我尝试过的代码，但我不断收到通知“事件时间不能超过正确审查的 ctimes”，因为对于效用得分 >0 和 <1 的患者，得分高于我创建的 0/1 yc 变量.

daten <- read.table ("P:/XXX.csv", header=TRUE, sep=";")

attach(daten)

x=cbind(factor(qlq) , AGE , SEX)

daten$c <- 1

daten$d <- ifelse (daten$UTILITY<1,0,1)

yc <- daten$d

y <- daten$UTILITY

clad <- crq (Curv(UTILITY, d, "right") ~ x, tau=0.5, method="Powell", data=daten)

先感谢您！

4

1 回答 1

0

如果有人遇到同样的障碍：对于每个 y，yc 需要是 y 被审查的值，而不是 0/1 审查指标。

在我的情况下（y=EQ-5D-5L 的效用分数）yc 需要为 1。

以下命令可以解决问题： daten$d <- rep(1.000,377) （因为我有 377 个观察值）

于 2018-12-13T16:13:45.653 回答