用R中的nlsList进行非线性回归
以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)
}
我想通过使用plot(mycars2 [[i]])进行散点图,然后使用tanh回归。 不幸的是我总是有这样的错误:
dim(x) must have positive length
你可以帮帮我吗? 我可以直接使用吗?
nlme() or other functions to do regressions and xyplot results?
非常感谢你!
链接地址: http://www.djcxy.com/p/30949.html