1

私はここに来たばかりで、あなたの助けが必要です。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)

どんな助けにも感謝します、ありがとう!

4

0 に答える 0