我正在尝试使用8个预测变量来拟合GAM模型。 8中的一个是a * exp(b * X)形式的指数衰减函数,其中b <0。其他预测因子是线性的。
我可以使用nls找到a,b:
out <- nls(Y~a*exp(b*X1),data=dat1,start=list(a=-1.5,b=1e-4))
summary(out)
现在我想要拟合多元回归模型并找到符合以下一般形式的最佳a,b:
out <- gam(Y~nls(a*exp(b*X1)) +X2+X3+X4+X5+X6+X7+X8, data=dat1)
有没有办法在R中实现这一目标? Ilik