拟合陡峭的剂量反应曲线

时间:2017-06-21 19:02:28

标签: r drc

嗨我想在R中的drc包中拟合一个相当陡峭的剂量 - 反应曲线.Moles是我的x函数,pct_control是我的y函数。我尝试将它安装在Log-logistic 3参数函数上,并且拟合非常糟糕。我尝试了其他模型拟合,我能得到的最接近的是logistic(5参数)。任何人都有其他关于模型拟合的想法吗?

 dm2005.test1<-drm(pct_control~moles, data=DONMON,fct=LL.3())
  

optim中的错误(startVec,opfct,hessian = TRUE,method = optMethod,   control = list(maxit = maxIt,:非有限的有限差分值   [3] drmOpt错误(opfct,opdfct1,startVecSc,optMethod,   约束,warnVal,:收敛失败

moles   pct_control
0.00344408  7.124090958
0.00344408  8.491257546
0.00344408  7.091147185
0.00344408  5.625149276
0.00344408  10.12197432
0.00344408  7.955921231
0.00344408  7.3794052
0.00344408  9.463098856
0.00344408  3.747354203
0.00344408  4.61212825
0.00344408  5.180408338
0.013748484 6.061654272
0.013748484 9.784300645
0.013748484 9.430155083
0.013748484 5.394542864
0.013748484 7.824146139
0.013748484 7.947685288
0.013748484 5.691036823
0.013748484 8.079460381
0.034387728 1.523649511
0.034387728 2.001334223
0.034387728 2.264884409
0.034387728 2.388423558
0.034387728 2.124873372
0.034387728 1.573065171
0.034387728 1.99309828
0.034387728 0.831830274

1 个答案:

答案 0 :(得分:0)

我尝试使用您的drm代码来适应所提供的数据,但没有得到您获得的错误。我能够绘制结果,并在下面包含该图像。我认为你遇到的困难可能是因为“鼹鼠”点的数量很少 - 你有3个。根据我的经验,如果你想要拟合一个合理的曲线,你将需要~5个点,跨越较低的渐近线,上渐近线和过渡区域。这样做的原因是,如此小的覆盖范围,有很多方法可以适应曲线。例如,鉴于数据,尚不清楚7“pct_control”的平台是否代表较低的渐近线(例如,全剂量反应将跨越100%降至7%),或者是否是上渐近线(例如全剂量)响应将延伸7%至0%) - 或某些中间结果。

enter image description here