R中nlsList与组的非线性回归

时间:2016-12-19 14:01:36

标签: r nls nlme

以R中的mtcars为例:我想对每个圆柱体和每个圆群的mpg和disp之间的关系进行非线性回归 在我的情况下,我有一个比mtcars大得多的数据帧。我知道没有足够的数据。但为了简化问题,我将始终使用mtcars,但假设我们有足够的数据来进行多组非线性回归,公式为tanh类型。 我想为每个cyl和每组vs做一个tanh类型回归,所以我尝试了:

mtcars2 <-list()
reg2 <-list()
attach(mtcars)
for (i in levels(cyl))
{
    mtcars2[[i]] <-
        groupedData( mpg ~ disp | vs,
             data = subset(mtcars,cyl==i))
    reg2[[i]] <- nlsList(mpg ~ A + B * tanh(disp/C), data=mtcars2[[i]],
                   start = list(A=1, B=1, C=0.5), na.action = na.omit)
}

我想通过绘图(mycars2 [[i]])和tanh回归来做散点图。不幸的是,我总是有错误:

dim(x) must have positive length
你可以帮帮我吗?我可以直接使用

nlme() or other functions to do regressions and xyplot results?

非常感谢!

0 个答案:

没有答案