在库QUANTREG中使用CRQ审查了最小绝对偏差（CLAD）回归

Question

我想使用以下方法执行CLAD回归

• y = EQ-5D-5L实用程序得分（从上到下限为1.0）
• x =各种患者特征

我已经发现我需要在库 QUANTREG 中使用 CRQ ，但是到目前为止我还不清楚。我的问题是：

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

这是我尝试的代码，但是我不断收到通知“事件时间不能超过正确审查的c次” ，因为对于效用得分> 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)

提前谢谢！

Answer 1

如果有人遇到相同的障碍：对于每个y，yc必须是y被检查的值，而不是0/1检查指标。

在我的情况下（y = EQ-5D-5L的实用分数），yc必须为1。

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