如何使用try (或tryCatch)来拟合4参数逻辑模型或3参数逻辑模型

Question

我想用drc软件包对一些毒性数据拟合一个4参数的logistic回归模型。这是一个示例数据集：

exp.df <- as.data.frame(matrix(nrow = 14, ncol = 2))

exp.df[,1] <- c(3200, 9600, 25600, 25600, 25600, 1600, 800,
                6400, 19200, 0, 12800, 1200, 400, 2400)

exp.df[,2] <- c(0.5855615, 0.9625668, 0.4064171, 0.4973262,
                0.4732620, 1.0000000, 0.6764706, 0.4652406,
                0.5106952, 0.7566845, 0.5427807, 0.5106952, 
                0.5935829, 0.4759358)

names(exp.df) <- c("Dose", "Response")

Dose_Response <- drm(Response ~ Dose, data = exp.df,
                     fct = LL.4(), type = "binomial")

本警告的结果如下：

optim中的误差(startVec，opfct，h hessian = TRUE，method = optMethod，control =list)(maxIt= maxIt，：drmOpt中的非有限差分值4的误差)(opfct，opdfct1，startVecSc，optMethod，约束，warnVal，：收敛失败)

如果我运行一个3参数的逻辑重新调整，它可以工作。

Dose_Response <- drm(Response ~ Dose, data = exp.df,
                     fct = LL.3(), type = "binomial")

我可以先使用try或tryCatch测试4个参数的逻辑，然后再使用3个参数的logistic回归吗？还是有别的解决办法？

Answer 1

你的cluster(id)术语是什么？在我添加了一个假的数据之后，您的数据就用LL.3()和LL.4()解决了。