Нахождение эмпирического отношения правдоподобия - PullRequest
Новая
       5

Нахождение эмпирического отношения правдоподобия

0 голосов
/ 21 апреля 2019

У меня есть домашняя работа о статье, и результат должен быть в интервале (0.10024, 1.0917).

В качестве конкретного примера мы взяли данные о времени ремиссии для пациентов с солидной опухолью (n = 10), которые представляют собой слегка модифицированную (прерывистую) версию Статистические методы анализа данных выживания, Elisa T. Lee, 1992 , Пример 4.2 :

3, 6,5, 6,51, 10, 12, 15, 8,4+, 4 +, 5,7 + и 10 +.

Предположим, что мы заинтересованы в получении 95% доверительного интервала для кумулятивная опасность в момент времени t = 9,8, ∆o (9,8). Следовательно, θo = ∆o (9.8). В этом Если функция g является индикаторной функцией: g (t) = I [t9.8].

95% доверительный интервал с использованием эмпирического отношения правдоподобия, -2logALR, для ∆o (9,8) составляет (0,10024, 1,0917)
Пожалуйста, помогите мне получить результат выше. Спасибо. 

remissiontime<-(3,4,5.7,6.5,6.51,8.4,10,10,12,15) 
status <- c(1,0,0,1,1,0,1,0,1,1)

и мой код (на самом деле я не уверен в этом коде) 

library(survival)
library(emplik)
x1 = c(3,4,5.7,6.5,6.51,8.4,10,10,12,15) 
d1 = c(1,0,0,1,1,0,1,0,1,1)

KM0 <- survfit(Surv(x1,d1) ~ 1,  type="kaplan-meier", conf.type="log")
summary(KM0)

myfun <-function(t){as.numeric(t <=9.8)}

emplikH1.test(x=x1,d=d1,theta=-log(0.643),fun=myfun)


myULfun <-function(theta,x,d){

  emplikH1.test(x=x1,d=d1,theta=theta,fun=function(t){as.numeric(t <= 9.8)})}

findUL(fun=myULfun,MLE =-log(0.643),x=x1,d=d1)

1 Ответ

0 голосов
/ 29 апреля 2019

Последняя строка представленного вами кода 

findUL(fun=myULfun,MLE =-log(0.643),x=x1,d=d1)

выдает ошибку:

Ошибка в emplikH1.test (x = x1, d = d1, theta = theta, fun = function (t) {: заданная тета слишком далеко от тета0

Изменение MLE аргумента findUL с -log(0.643) на 1 устраняет проблему. Пожалуйста, смотрите ниже: 

findUL(fun = myULfun, MLE = 1, x = x1, d = d1)

Вывод, который совпадает с требуемым: 

$`Low`
[1] 0.1002516

$Up
[1] 1.09165

$FstepL
[1] 6.103516e-05

$FstepU
[1] 6.103516e-05

$Lvalue
[1] 3.839313

$Uvalue
[1] 3.839971
...