多元线性回归——《应用回归分析R语言版》

#第3章—多元线性回归data3.1<-read.table('clipboard',header=T)#导入数据
lm3.1<-lm(y~x1+x2+x3+x4+x5+x6+x7+x8+x9,data=data3.1) #建立回归方程
summary(lm3.1)#计算协方差阵与相关阵
X1<-as.matrix(data3.1[,2:10])   #将data3.1中自变量部分的数据提取出来并转换成矩阵形式进行存储
X<-cbind(1,X1)               #将元素全为1的列向量和X1合并生成矩阵X
XX<-t(X)%*%X                #计算矩阵X'Xsigma<-389.4                 #残差的标准差为  的估计值
covBeta<-sigma^2*solve(XX)   #根据公式（3.31）计算协方差矩阵，函数solve()计算矩阵的逆
covBetar<-matrix(nrow = 10, ncol = 10) #建立10行10列的矩阵，元素为空for(i in 1:10){for(j in 1:10)
r[i,j]<-covBeta[i,j]/(sqrt(covBeta[i,i])*sqrt(covBeta[j,j]))}
#根据公式（3.34）计算相关阵中每个元素的值
r
#r<-cor(data3.1)#方差分析
#install.packages("car")
library(car)
Anova(lm3.1,type="III")       #输出方差分析表，typeII用于回归分析以及不含交互作用的方差分析；typeIII用于含交互作用的方差分析；lm3.1_drop4<-lm(y~x1+x2+x3+x5+x6+x7+x8+x9,data=data3.1) #lm3.1_drop4<-update(lm3.1,.~.-x4)
summary(lm3.1_drop4)lm3.1_drop49<-lm(y~x1+x2+x3+x5+x6+x7+x8,data=data3.1)
summary(lm3.1_drop49)data3.2<-read.table('clipboard',header=T)#导入数据
lm3.2<-lm(y~x1+x2,data=data3.2)
summary(lm3.2)
r<-cor(data3.2)                       #计算相关系数阵
r#install.packages("corpcor")
library(corpcor)
pcor3.2<-cor2pcor(r) #由相关系数阵计算偏相关系数阵
pcor3.2e<-resid(lm3.2,digits=5)   #计算残差值
ZRE<-rstandard(lm3.2)    #计算学生化残差
ZRE1<-e/475.8           #计算标准化残差data3.3<-read.table('clipboard',header=T)#导入数据
cor3.3<-cor(data3.3[,-1]) #用除去第一列年份数据后剩余的样本数据计算相关阵
lm3.3<-lm(y~x1+x2+x3+x4+x5,data=data3.3)
summary(lm3.3)

多元线性回归——《应用回归分析R语言版》相关推荐

多元线性回归分析c语言,多元线性回归公式推导及R语言实现
多元线性回归多元线性回归模型实际中有很多问题是一个因变量与多个自变量成线性相关,我们可以用一个多元线性回归方程来表示. 为了方便计算,我们将上式写成矩阵形式: Y = XW 假设自变量维度为N W ...
多元线性回归公式推导及R语言实现
多元线性回归多元线性回归模型实际中有很多问题是一个因变量与多个自变量成线性相关,我们可以用一个多元线性回归方程来表示. 为了方便计算,我们将上式写成矩阵形式: Y = XW 假设自变量维度为N W ...
自变量选择与逐步回归——《应用回归分析R语言版》
data5.2<-read.csv("F:/R/应用回归分析/data/li3-1.csv",header = T) install.packages("leaps ...
【视频】什么是梯度下降？用线性回归解释和R语言估计GARCH实例
全文链接:http://tecdat.cn/?p=23606 梯度下降是一种优化算法,能够为各种问题找到最佳解决方案(点击文末"阅读原文"获取完整代码数据). 梯度下降是什么? 梯 ...
tushare 金融数据获取（R语言版）
在上次 tushare正确爬取指数数据文章后,看到后台有人留言说是希望能分享一个R语言版,那么好,你们要的R语言版今天来了. 首先,R语言只支持tushare pro,不支持tushare,因此在使 ...
带南海九段线分位数地图可视化（R语言版）
今天带来一篇承诺虾神的可视化博客.内容是使用R语言进行带南海九段线分位数地图可视化.虾神的原博文地址如下(Python版). Python实现带南海九段线分位数地图完整可视化版本(附代码及数据) 19 ...
回归模型选择(R语言版)
回归模型选择(R语言版) 选择准测: With Cp , AIC and BIC, smaller values are better, but for adjusted R2 , larger va ...
论文中常见的拟合散点验证图（R语言版）
论文中常见的拟合散点验证图(R语言版) 如上图所示,是论文中常见的validation图,python也能实现相似的图绘. 今天先介绍R语言版,python改期再介绍吧这张图需要依次实现下列功能: ...
r读取shape文件可视化_R语言读取空间数据以及ArcGIS中OLS工具回归结果可视化R语言版...
前面已经介绍过R语言读取excel的方法了,当然读取数据来说,个人还是推荐csv或txt存储(针对小数据量).大数据量的数据的话建议还是用数据库,此外也可以考虑data.table包读取,这个包也是个 ...

多元线性回归——《应用回归分析R语言版》

多元线性回归——《应用回归分析R语言版》相关推荐

最新文章

热门文章