CLAD回帰を実行したい
- y = EQ-5D-5L ユーティリティ スコア (上から 1.0 で制限)
- x = さまざまな患者の特徴
ライブラリQUANTREGでCRQを使用する必要があることはすでにわかりましたが、これまでのところ詳細を把握できませんでした。私の質問は次のとおりです。
- パウエル法を使用する必要がありますか?
- もしそうなら、時間変数ではなく 0/1 打ち切り変数を持っている場合、どのように "yc" (打ち切り時間) を指定すればよいですか?
これは私が試したコードですが、ユーティリティ スコアが 0 を超えて 1 未満の患者の場合、作成した 0/1 yc 変数よりもスコアが高いため、「正しい打ち切りではイベント時間は ctimes を超えることはできません」という通知が引き続き表示されます。 .
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)
前もって感謝します！