R语言实战——线性回归分析和相关矩阵可视化

发布时间：2025-05-21 05:35

多元回归分析

多元回归分析预测法，是指通过对两个或两个以上的自变量与一个因变量的相关分析，建立预测模型进行预测的方法。当自变量与因变量之间存在线性关系时，称为多元线性回归分析。

多元线性回归的适用条件：

（1）自变量对应变量的变化具有显著影响。

（2）自变量与应变量间的线性相关必须是真实的，而非形式上的。

（3）自变量之间需有一定的互斥性。

（4）应具有完整的统计数据。

下面，我们就选取一个例子，用R实现多元线性回归模型。

训练数据：txt格式（行为样本，列为食物）。此例子中，我们选取10类食物作为特征，即自变量为10类食物。因变量为产生的蛋白质含量（Protein）。

1检查两两因变量之间的关系

x=read.table(test.txt,header=T,sep=\t,row.names=1)

x=as.matrix(x)

z=x[,1:10]

cor-cor(z,use=pairwise,method=pearson)

cor

corrplot(cor)

相关性图的展示使用的是corrplot包。具体步骤参见R语言绘图-相关矩阵的可视化。

从图中可以看出，各个自变量之间的线性程度不高，可以直接作为回归参数输入。因为如果自变量之间的线性关系超过自变量与因变量之间的线性关系，则回归模型稳定性受到破坏。

1多元回归模型建立

接下来，我们用10个特征作为自变量，蛋白质含量作为因变量。建立多元线性回归模型。

y=x[,11]

lm.result-lm(y~z)

summary(lm.result)

***表示极为显著，**表示高度显著，*表示显著。以上是模型检验的结果。因为我们的例子数据不全，所以大家可以用自己数据试一试，看看结果好不好。

1残差分析

par(mfrow=c(2,2))

plot(lm.results)

以上就是使用R实现多元线性回归的过程。