安装AIC准则使用前进法后退法和逐步回归法进行变量选择的r语言代码

setwd("C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业") #设定当前的工作目录
shuju=read.table("shuju.txt",header=T)
shuju #读取数据
#采用AIC原则自动选择模型-前进法
shuju.reg1<- lm(y~.,data=shuju[,-1])
shuju.regforward2 <- step(shuju.reg1,direction="forward")#按照AIC原则自动选择模型
summary(shuju.regforward2)
#采用AIC原则自动选择模型-后退法
shuju.reg2<- lm(y~.,data=shuju[,-1])
shuju.regbackward2 <- step(shuju.reg2,direction="backward")#按照AIC原则自动选择模型
summary(shuju.regbackward2)
#采用AIC原则自动选择模型-逐步回归法
shuju.reg3<- lm(y~.,data=shuju[,-1])
shuju.regboth <- step(shuju.reg3,direction="both")#按照AIC原则自动选择模型
summary(shuju.regboth)

#计算方差扩大因子VIF
shuju.reg=lm(y~x1+x2+x3+x4+x5+x6+x7,data=shuju)
library(car)
vif(shuju.reg)#计算得方差扩大因子

#计算条件数condition index
X3<-cbind(shuju$x1,shuju$x2,shuju$x3,shuju$x4,shuju$x5,shuju$x6,shuju$x7)
CX<-cor(X3)#变换后所得的X'X其实就是相关系数阵
sqrt(kappa(CX,exact=T))#与SPSS的结果略有区别
eigen(CX)#求特征根和特征向量

#剔除一些不重要的解释变量
#先剔除x4
drop1(shuju.reg,scope="x4",test='F')
shuju.reg2 <- update(shuju.reg,~.-x4)
vif(shuju.reg2)
#再剔除x5
drop1(shuju.reg2,scope="x5",test='F')
shuju.reg3 <- update(shuju.reg2,~.-x5)
vif(shuju.reg3)

summary(lm(y~x1+x2+x3+x6+x7,data=shuju))

附录，数据如下、years y x1 x2 x3 x4 x5 x6 x7
1974 172.9 11246 681 105.9 10183 4110 11242 9
1975 352.94 10335 791 107.4 10414 3996 12693 6.5
1976 447.67 13156 607 114.4 13134 4689 16681 6
1977 404.02 6127 714 110.8 15033 6876 22131 4.75
1978 409.51 27419 911 99.4 17389 8636 31353 4.75
1979 619.71 25633 1231 91.4 21715 12339 43528 9.5
1980 1121.17 95684 2760 90.8 27075 16623 70752 10
1981 1506.94 105987 2651 86.3 31827 19937 125989 16
1982 1105.79 46230 2105 125.3 35393 24787 99468 10.5
1983 933.03 37165 3030 107.4 38823 25112 82478 10.5
1984 1008.54 48787 2810 106.6 46079 24414 54936 8.5
1985 1567.56 75808 2649 115.7 47871 22970 87135 6
1986 1960.06 123128 3031 110.1 54372 24403 129884 6.5
1987 2884.88 371406 3644 105.8 65602 30531 153044 5
1988 2556.72 198569 3690 101.6 74917 37861 215033 5.25

安装AIC准则使用前进法后退法和逐步回归法进行变量选择的r语言代码相关推荐

大数据----------------R语言下依赖库与依赖包的安装
由于博主最近在学习大数据的基础,避免不了要搭建以hadoop,hbase,hive等软件为基础的环境,这一路的bug可谓是层出不穷啊!在历经万苦后终于将前面的都安装好了,顺利了一会儿,没想到在数据可视 ...
R语言安装包，安装Github包的三种方法
以安装Github上的Achilles包为例进行安装说明. 描述 1.安装包路径:https://github.com/OHDSI/Achilles#getting-started 2.该包功能:进行 ...
【数据分析R语言系列】R和RStudio的下载和安装， R在 Ubuntu 和CentOS 系统下的安装
文章目录准备工作 R 的下载和安装 RStudio 的下载和安装配置(可选) 常见问题与方案 R 在 Linux 系统下的安装 Ubuntu 从源安装 R CentOS 从源安装 R Rtools ...
R语言运行环境安装配置
想要运行R语言代码,首先需要安装好R语言的运行环境--RGui,很多人包括我自己喜欢另一个运行环境--Rstudio.下面分别就安装RGui.Rstudio做分别说明. 1.安装RGui 去R语言官方 ...
R 语言 4.2.2安装 WGCNA
文章目录 1 WGCNA库介绍 2 安装踩坑还得是官方文档这样安装我出现的问题参考 Appendix A. 安装R B. 配置环境 C. 修改镜像 1 WGCNA库介绍 WGCNA是用于加权相 ...
Wine的安装和配置（实现中文输入法和中文的复制粘贴）(转)
Wine的安装和配置(实现中文输入法和中文的复制粘贴)(转) liuspider 最后修订于 2003年4月14日版本 0.7 本文档最新版可以在 http://xbkconfp.cosoft.or ...
［更新问题］无法在安装新的版本前，为“./boot/vmlinuz-2.6.24-19-generic”做一个符号链接备份...
2019独角兽企业重金招聘Python工程师标准>>> 关于更新中出现的问题:我在升级的时候,出现一个问题: E: /var/cache/apt/archives/linux-ima ...
故障转移集群服务器重装系统,安装故障转移群集前的准备工作
安装故障转移群集前的准备工作 08/24/2016 本文内容适用于: SQL Server(所有支持的版本) 安装 SQL Server 故障转移群集之前,必须选择运行 SQL Server 的硬件 ...
［转］安装N卡驱动前的准备和一点点小建议
安装N卡驱动前的准备和一点点小建议: 1.购买N卡时,为了安装配制简单省事,请尽量购买做工出色,用料扎实且采用公版设计的N卡,谨慎购买做工用料缩水,并且兼容性不好的N卡,驱动安装不上也不排除与主机板 ...
用U盘装机大师安装GHOST WIN8系统前准备
用U盘装机大师安装GHOST WIN8系统前准备 http://www.upandashi.com/win8/ 1.首先准备一个U盘,容量最好大于4G,然后下载U盘装机大师U盘启动盘制作工具,制作成U ...

安装AIC准则使用前进法后退法和逐步回归法进行变量选择的r语言代码

安装AIC准则使用前进法后退法和逐步回归法进行变量选择的r语言代码相关推荐

最新文章

热门文章