安装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 ...
最新文章
- vc中常用文件操作（二） Ini文件操作
- 5年没有工资收入，他如何支撑世界上最大的免费编程社区？
- Python 内建函数 max/min的高级用法
- boost::mp11::mp_partial_sum相关用法的测试程序
- SAP UI5 component.js createContent
- php四种标量,php的标量数据类型和复合数据类型,php四种标量
- datetime报错 sql脚本_SQLSERVER 中datetime 和 smalldatetime类型分析说明
- python贪吃蛇源码_Python：游戏：贪吃蛇（附源码）
- 深度学习难？学完它，拿高薪绝对不是问题
- Eclipse快捷键与使用技巧总结
- PHP人民币金额数字转中文大写的函数
- 基于Sendmail、Dovecot和Ssl的电子邮件服务器搭建——从不安全通信到安全通信
- 全球最值得模仿的500个网站（扫描版pdf）
- OSChina 周二乱弹 ——流川枫与苍井空
- 学习笔记(01):程序员的数学：微积分-常用导数（一）：最常用到的技巧
- Adobe Photoshop 2021 22.4.2 绿色精简版
- Final Cut Pro X Guru: Advanced Trimming Final Cut Pro X Guru：高级修剪 Lynda课程中文字幕
- php 数组课件,php学习 数组课件第1/2页
- ＜script＞标签通过src属性调用js文件
- 凝眸笑靥，又是一年寂寥春