Как проверить предположение PH с остатками Шенфельда с моделируемыми путями в присутствии конкурирующих событий - PullRequest
Новая
       62

Как проверить предположение PH с остатками Шенфельда с моделируемыми путями в присутствии конкурирующих событий

0 голосов
/ 18 мая 2018

Я получил остатки Шенфельда для тестирования PH для регрессионной модели Fine & Gray's Cox при наличии конкурирующих событий с cmprsk.

Вот код: 

 fg<-crr( fu_m, event, diabetes,  failcode=1, cencode=0,
 na.action=na.omit, gtol=1e-06, maxiter=10, variance=TRUE)

fg$res
fg$uft
par(mfrow = c(1,1), mar = c(4.5,4,2,1))
for(j in 1:ncol(fg$res))
  fg$u
scatter.smooth(fg$uftime, fg$res[,j],
           main = names(fg$diabetes)[j],
           xlab = "Failure time",
           ylab = "Schoenfeld residuals")

Вот вывод:

> fg$res
           [,1]
[1,]  0.5498603
[2,] -0.3957394
[3,] -0.4024953
[4,]  0.5905142
[5,] -0.3421397

> fg$uft
[1]  4.238193  6.275154 16.131417 27.498973 46.817248

Вот график:

r1

Я сделал то же самое, используя SAS, и получил p-значение 0,55 на 1000 смоделированных путей и этот график для наблюдаемого пути и первых 20 смоделированных путей:

sas1

Как проверить предположение PH следующим образомЯ делаю с SAS?Как я могу представить имитированные пути?

Спасибо всем !!

...