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

时间:2016-11-23 19:25:10

标签: r

我想使用包drc将4参数逻辑回归模型拟合到一些毒性数据中。这是一个示例数据集:

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,hessian = TRUE,method = optMethod,   control = list(maxit = maxIt,:非有限的有限差分值   [4] drmOpt错误(opfct,opdfct1,startVecSc,optMethod,   约束,warnVal,:收敛失败

如果我运行3参数logistic regresison,它可以工作。

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

我可以先使用try或tryCatch测试4参数logistic,然后使用3参数logistic回归吗?还是有另一种解决方案吗?

1 个答案:

答案 0 :(得分:0)

您的cluster(id)期限是多少?在我添加了假的后,您的数据已同时使用LL.3()LL.4()解决。