适合在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变量。