Цензура регрессии наименьшего абсолютного отклонения (CLAD) с использованием CRQ в библиотеке QUANTREG

Question

Я хотел бы выполнить регрессию CLAD с

• y = баллы полезности EQ-5D-5L (ограничены сверху 1,0)
• х = различные характеристики пациента

Я уже узнал, что мне нужно использовать CRQ в библиотеке QUANTREG, но я до сих пор не мог выяснить специфику. Мои вопросы:

1. Нужно ли использовать метод Пауэлла?
2. Если да, то как мне указать "yc" (время цензуры), если у меня нет переменной времени, но есть переменная цензуры 0/1?

Это код, который я пробовал, но я продолжаю получать уведомление "Время событий не может превышать 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)

Заранее спасибо!

Answer 1

Если кто-то когда-либо сталкивался с одним и тем же препятствием: для каждого y, yc должен быть значением, при котором y подвергается цензуре, а не показателем цензуры 0/1.

В моем случае (y= оценка полезности EQ-5D-5L) yc должно быть 1.

Следующая команда делает трюк: daten$d <- rep (1.000,377) (потому что у меня 377 наблюдений)