2017-06-21 12 views
1

こんにちは私は、Rのdrcパッケージにかなり急な用量反応曲線を収めようとしています。Molesは私のx関数で、pct_controlはmy y関数です。 Log-Logistic 3パラメータ関数にフィットさせてみましたが、フィット感はひどいものでした。私は他のモデルのフィットを試して、私が得ることができる最も近いフィットはロジスティック(5パラメータ)でした。誰にもモデル適合の他のアイデアはありますか?急勾配用量応答曲線をフィッティング

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

OPTIMでのエラー(startVec、opfct、ヘッセ= TRUE、方法= optMethod、 コントロール=リスト(maxit回= maxit回、非有限有限差分値 drmOpt [3]エラー( 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

関連する問題