nlsLM应用程序需要初始参数

时间:2013-09-23 07:27:44

标签: r regression

我是新来的,我需要你的帮助,我正在努力通过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)

感谢任何帮助,谢谢!

0 个答案:

没有答案