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