適合在R 3參數中修改的weibull

Question

我想將以下數據擬合為Weibull分布乘以a。

datos： 在此處輸入鏈接描述

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

但是，我無法使用R中的nls函數來找到解決方案。有人可以指導我找到解決方法嗎？

使用的代碼如下：

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 。 它是為實現參數生存回歸而設計的（當然包括Weibull模型）。 這是代碼：

library(survival)

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

您看到的而不是y的Surv()對象是R“生存對象”，並且像生存回歸中的因變量一樣工作。 time和event變量必須分別表示持續時間和事件發生（0或1）。 請用您適當的變量集替換explanatory_variables變量。