R语言随手记-线性回归模型诊断

gvlma()

gvlma函数可以用来检测线性拟合模型的假设是否成立，并对峰度、偏度进行验证。

install.packages("gvlma")
library(gvlma)
gvlma(fit)

结果如下，原假设是假设成立，所以当p>0.05时，才符合线性模型

• Global Stat [全局统计]：代表X个预测变量和Y个之间是否大致呈线性关系。拒绝空值(p <.05)表示X与Y中的一个或多个之间存在非线性关系。
• Skewness [偏度]：分布是正偏还是负偏，是否需要进行转换才能满足正态性假设。拒绝null(p <.05)表示应该转换数据。
• Kurtosis [峰度]：分布是峰态(高度峰值还是非常浅的峰值)，是否需要进行转换才能满足正态性假设。拒绝null(p <.05)表示应该转换数据。
• Link Function [链接功能]：因变量是真正连续的还是分类的。拒绝null(p <.05)表示应使用广义线性模型的另一种形式(例如，逻辑回归或二项式回归)。
• Heteroscedasticity [异方差性]：检验模型残差的方差在X范围内(同质性假设)是否恒定。拒绝空值(p <.05)表示残差是异方差的，因此在X范围内是非恒定的，模型在预测X比例尺的某些范围时更好/更差。