Интерпретация параметров Вейбулла из суррег

Question

Я пытаюсь сгенерировать обратное распределение Вейбулла, используя параметры, оцененные по коэффициенту выживания в R. Под этим я подразумеваю, что хотел бы, для данной вероятности (которая будет случайным числом в небольшой имитационной модели, реализованной в MS Excel) , верните ожидаемое время до отказа, используя мои параметры. Я понимаю, что общая форма обратного распределения Вейбулла:

X=b[-ln(1-rand())]^(1/a)

где a и b - параметры формы и масштаба соответственно, а X - время до отказа, которое я хочу. Моя проблема заключается в интерпретации параметров перехвата и ковариации из surreg. У меня есть эти параметры, единица времени - дни:

             Value   Std. Error    z    p
(Intercept)     7.79    0.2288  34.051  0.000
Group 2        -0.139   0.2335  -0.596  0.551
Log(scale)     0.415    0.0279  14.88   0.000
Scale= 1.51 

Weibull distribution
Loglik(model)= -8356.7   Loglik(intercept only)= -8356.9 
Chisq = 0.37 on 1 degrees of freedom, p= 0.55 
Number of Newton-Raphson Iterations: 4 
n=1682 (3 observations deleted due to missing values)

Я прочитал в файлах справки, что коэффициенты от R взяты из "распределения экстремальных значений", но я не уверен, что это на самом деле означает и как я "возвращаюсь" к стандартному параметру масштаба, используемому непосредственно в формулах. Использование b = 7,79 и a = 1,51 дает бессмысленные ответы. Я действительно хочу быть в состоянии генерировать время как для основной группы, так и для «Группы 2». Я также должен отметить, что я сам не проводил анализ и не могу больше запрашивать данные.

Answer 1

Это объясняется на странице руководства, ?survreg (в разделе "примеры").

library(survival)
y <- rweibull(1000, shape=2, scale=5)
r <- survreg(Surv(y)~1, dist="weibull")
a <- 1/r$scale      # Approximately 2
b <- exp( coef(r) ) # Approximately 5
y2 <- b * ( -ln( 1-runif(1000) ) ) ^(1/a)
y3 <- rweibull(1000, shape=a, scale=5)
# Check graphically that the distributions are the same
plot(sort(y), sort(y2))
abline(0,1)

Answer 2

Ключ в том, что параметр формы, сгенерированный rweibull, является обратным параметру формы, который вводит Surregreg