我想拟合一个线性模型(lm),其中残差方差明显取决于解释变量。
我知道这样做的方法是使用带有 Gamma 族的 glm 对方差进行建模,然后将其倒数放入 lm 函数中的权重中(例如: http: //nitro.biosci.arizona.edu/r/chapter31 .pdf )
我在想:
- 这是唯一的技术吗?
- 还有哪些相关的方法?
- 与这种类型的建模相关的 R 包/功能是什么?(除 glm、lm 外)
我想拟合一个线性模型(lm),其中残差方差明显取决于解释变量。
我知道这样做的方法是使用带有 Gamma 族的 glm 对方差进行建模,然后将其倒数放入 lm 函数中的权重中(例如: http: //nitro.biosci.arizona.edu/r/chapter31 .pdf )
我在想:
R包中的gls
函数nlme
可以同时估计回归和与方差的关系。请参阅weights
帮助页面上的参数和第二个示例。
使用gamlss包,您可以将响应的误差分布建模为解释变量的线性、非线性或平滑函数。这似乎是一种非常强大的方法(我学到了很多关于模型选择过程中可能出现的所有可能性),并且在上面链接中引用的几本出版物(包括书籍)中很好地解释了一切。