私はここに来たばかりで、あなたの助けが必要です。van genuchten モデルによる保水曲線のフィッティングに取り組んでおり、R のパッケージのnls
:nlsLM
メソッドの修正版を使用しようとしていますminpack.lm
が、その方法について質問があります。 「開始」オプションで作業します。つまり、このモデルの初期パラメーターを決定する方法は?、現在、ランダムな開始値を設定しており、うまく機能します (警告やエラー コードはありません)。私のジレンマは、van Genuchten モデルのパラメータの初期値を計算する方法ですか?, ここに私のコード:
sample <- data.frame(h = c(61,82,173,337,683,2356,7283,15300),
theta_h = c(0.481,0.41,0.374,0.235,0.24,0.155,0.084,0.096))
mod <- nlsLM(theta_h ~ theta_r + (theta_s – theta_r)/(1 + (alpha*h)^n)^(1-1/n),
data = sample,
control = nls.lm.control(maxiter=200,options(warn=-1)),
start = c(theta_s = 0.45, theta_r = 0.067, alpha = 0.02, n = 1.41),
trace = TRUE)
どんな助けにも感謝します、ありがとう!