2017-06-21 42 views
1

嗨我想在R的drc軟件包中擬合一個相當陡峭的劑量反應曲線。Moles是我的x函數,pct_control是我的函數。我嘗試將它安裝在Log-logistic 3參數函數上,並且配合非常糟糕。我嘗試了其他模型,我能得到的最接近的數據是logistic(5參數)。任何人有任何其他模型擬合的想法?擬合陡峭的劑量反應曲線

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

錯誤的Optim(startVec,opfct,粗麻布= TRUE,方法= optMethod, 控制=列表(麥克斯特= 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 

回答

0

我試圖用您的DRM代碼,以適應所提供的數據,並沒有得到你所得到的錯誤我能。繪製結果和我包括下面的圖片。我認爲你遇到的困難可能是由於「痣」點數很少 - 你有3個。根據我的經驗,如果你想要擬合一個合理的曲線,你將需要〜5點,跨越較低的漸近線,上漸近線和過渡區域。原因在於,如此小的覆蓋範圍,有許多方法來擬合曲線。例如,給定數據,尚不清楚7「pct_control」處的平臺是否代表較低的漸近線(例如,全劑量反應將跨越100%降至7%),或者如果它是上部漸近線(例如全劑量響應將跨越7%降至0%)或一些中間結果。

enter image description here