在R中使用非线性函数的分段非线性回归

Question

我希望使用具有多个断点的非线性函数分段执行回归。 我已经完成了分段线性回归，但是当指定任何类型的非线性函数时，我们如何在R中设置？

具体来说，我对使用两个断点的3个函数线性，指数和指数感兴趣。 请指教

KARTHIK

Answer 1

使用nls() （非线性最小二乘）会解决你的问题吗？ 我使用了类似于此的公式，为每个“片段”添加了True / False语句：

reg = nls( y ~ (Z < 0.33) * a + (Z < 0.33) * Z * b +
        (Z >= 0.33 & Z < 0.67) * Z ^ a2 +
        (Z >= 0.67) * a3 + (Z >= 0.67) * Z * a4,
        start = list(a = 0, b = 50, a2 = 100, a3 = 150, a4 = 80),
        data = yourdata)

在上面的程式化示例中，断点位于Z = 0.33且Z = 0.67。 如果可以更具体，或者分别提供三个回归的代码，我可以使答案更具体。

Answer 2

我的建议是加载“ splines”程序包，然后在help(bs)运行示例。 您可以使用线性回归机制获得分段立方（但在结点处连续）拟合。 Harrell在他的“ rms”软件包中使用了此策略，效果极佳。 加载“ rms”并查看help(rcs) 。 该页面上的示例使用了逻辑回归的实现，但rcs()术语也在ols()和cph()中工作。