以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?
非常感谢!