Fit Weibull модифицируется в параметре R 3

Question

Я хочу подогнать следующие данные к распределению Вейбулла, умноженному на.

данные: введите описание ссылки здесь

                          y=b1*(1-exp(-(x/b2)^b3)

Однако я не смогнайти решение с помощью функции nls в R. Может ли кто-нибудь подсказать мне путь, по которому нужно идти, чтобы найти решение?

Используется следующий код:

ajuste_cg<-nls(y~b1*(-exp(-((x/b2)^b3))),data=d,start=list(b1=1000,b2=140,b3=20), trace=T,control = list(maxiter=10000000))

Спасибо!

Answer 1

Предлагаю вам воспользоваться пакетом survival.Это сделано для реализации параметрических регрессий выживания (конечно, в модели Вейбулла).Вот код:

library(survival)

weibull_model = survreg(Surv(time, event) ~ expalatory_variables, dist="weibull")

Объект Surv(), который вы видите вместо y, является R "объектом выживания" и работает как ваша зависимая переменная в регрессии выживания.Переменные time и event должны представлять продолжительность и возникновение события (0 или 1) соответственно.Пожалуйста, замените explanatory_variables на соответствующий набор переменных.