库QUANTREG中CRQ的截尾最小绝对偏差回归

Question

我想用

• Y= EQ-5D-5L实用分数(从高于1.0开始)
• X=各种病人特征

我已经发现我需要在库QUANTREG中使用CRQ，但是到目前为止我还没有弄清楚具体的内容。我的问题是：

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

这是我尝试过的代码，但我一直得到通知：“事件时间不能超过正确检查的时间”，因为对于效用分数>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需要被审查的值，而不是0/1的审查指标。

在我的例子中(EQ-5D-5L的y=utility分数) yc需要为1。

下面的命令非常有用: daten$d <- rep(1.000377)(因为我有377个观察)