Выживание :: определенное пользователем распределение бета-версий - PullRequest
Новая
       61

Выживание :: определенное пользователем распределение бета-версий

0 голосов

У меня есть набор данных с зависимой цензурированной переменной, распределенной по бета-распределению. Пожалуйста, не предлагайте flexsurv пакет. Необходимо описать функцию моего бета-распределения для использования в модели выживания регрессии AFT Survreg . В примере из книги М. Чжу я подставил известные мне функции, но не могу описать соотношение первой и второй производных функции плотности бета-распределения: 

# Beta Distribution - User-defined code for 'survreg'
a = 1.2; b = 0.2

myb <- list (name = "beta",

init = function (x, weights, ...)
{
mean <- sum (x * weights) / sum (weights)
var <- sum (weights * (x - mean) ^ 2) / sum (weights)
c (mean, var)
},

density = function (x, parms)
{
  cbind (
  pb (x, shape1 = a, shape2 = b, ncp = 0), # F - Cumulative Distribution Function
  1 - pb (x, shape1 = a, shape2 = b, ncp = 0), # S = 1 - F - Survival Function
  db (x, shape1 = a, shape2 = b, ncp = 0), # f - Density Function
  ?, # f'/ f - the first derivative of the Density function
  ? # f''/ f - the second derivative of the Density function
)
},

quantile = function (p, parms)
qb(x, shape1 = a, shape2 = b, ncp = 0),

deviance = function (...)
stop("deviance residuals not defined")
)

Я считаю, что это происходит из-за производных функции beta, функции gamma, digamma & trigamma в базовом пакете r , но мой уровень математики не позволять их применять.

...