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

安装AIC准则使用前进法后退法和逐步回归法进行变量选择的r语言代码setwd(“C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业”) #设定当前的工作目录shuju=read.table(“shuju.txt”,header=T)shuju #读取数据#采用AIC原则自动选择模型-前进法shuju.reg1shuju.regforward2summary(shuju.regforward2)#采用A

大家好,又见面了,我是你们的朋友全栈君。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

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/142381.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • Request header field Content-Type is not allowed by Access-Control-Allow-Headers跨域

    Request header field Content-Type is not allowed by Access-Control-Allow-Headers跨域跨域错误提示:XMLHttpRequestcannotloadhttp://xxx.com.RequestheaderfieldContent-TypeisnotallowedbyAccess-Control-Allow-Headers.解决方案:例如php服务端程序设置头:header(‘Access-Control-Allow-Origin:*…

  • ip route add命令详解_ip route什么意思

    ip route add命令详解_ip route什么意思最近使用haproxy+keepalived负载均衡,发现虚拟IP有问题,配置文件里配的VIP是···.4.143,但是执行ipa发现还有一个虚拟IP···.4.134;应该是以前配置文件keepalived.conf里设置过VIP···.4.134以前的VIP一直没有消除,手动清除VIP:ipaddrdel···.4.134/32deveth0顺便总结下常用的网络管理命…

  • Windows Server 2008 防火墙开放 Oracle 的1521端口

    Windows Server 2008 防火墙开放 Oracle 的1521端口在防火墙的入站规则中,新建端口规则。过程如下例图片所示:同理可以开放EM用的1158端口。执行完后用下面命令测试telnetSERVER_IP1521参考资料[1]WindowsServer2008防火墙如何配置(5).http://www.bitscn.com/netpro/firewall/200711/118934_

  • 不能复制文字的网页文字复制怎么办_html循环粘贴

    不能复制文字的网页文字复制怎么办_html循环粘贴网页无法复制文字怎么办?当我们在电脑上需要复制某个网页上的文字时,发现我们不能选择复制粘贴文字,那这种情况该怎么解决呢,网页无法复制文字怎么办,怎么解决网页无法复制粘贴文字情况,下面就和小编一起来看看吧!1.可以使用谷歌浏览器扩展程序AllowCopy解决问题,打开谷歌浏览器的网上应用店,搜索【AllowCopy】;2.然后找到SimpleAllowCopy,点击【添加至Chrome】将其…

  • 浅谈顶级域名和WWW二级域名收录区别[通俗易懂]

    浅谈顶级域名和WWW二级域名收录区别[通俗易懂]浅谈顶级域名和WWW二级域名收录区别      谈到这个话题,大家首先要认识两点,什么是顶级域名,什么是二级域名。顶级域名:简单讲不带前缀的域名称为顶级域名,例如:g2hh.com,ghh.com二级域名:位于顶级域名的下一级域名称为二级域名,例如:www.g2hh.com,由于大家经常浏览带有www的网页,误以为带有www的域名为顶级,其实不然,像www.g2hh.com的域名

  • tomcat官网如何下载低版本的tomcat

    tomcat官网如何下载低版本的tomcat

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号