我是新来的,我需要你的帮助,我正在努力通过van genuchten模型拟合保水曲线,我正在尝试使用nls
的修改版本:nlsLM
方法来自R中的minpack.lm
包,但我有一个关于如何使用'start'选项的问题,我的意思是,如何确定此模型的初始参数?,现在我设置随机起始值,并且运行良好(没有警告或错误代码)但这是样本的相对数据集,所以,我的困境是,如何计算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)
感谢任何帮助,谢谢!