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

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

0 голосов
/ 06 декабря 2018

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

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

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

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

Это код, который я пробовал, но я продолжаю получать уведомление "Время события не может превышать ctimes для правильной цензуры" , потому что дляПациенты с оценкой полезности> 0 и <1, оценка выше, чем переменная 0/1 yc, которую я создал. </p> 

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)

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

1 Ответ

0 голосов
/ 13 декабря 2018

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

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

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

...