AIC准则选三个变量的r语言代码

AIC准则选三个变量的r语言代码setwd(“C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业”) #设定当前的工作目录shuju=read.table(“shuju.txt”,header=T)shuju #读取数据#采用AIC原则自动选择模型-前进法library(MASS)stepAIC(lm(y~.,data=shuju[,-1])) #146.

大家好,又见面了,我是你们的朋友全栈君。

setwd(“C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业”)  #设定当前的工作目录
shuju=read.table(“shuju.txt”,header=T)
shuju   #读取数据
#采用AIC原则自动选择模型-前进法
library(MASS)
stepAIC(lm(y~.,data=shuju[,-1]))  #146.11
stepAIC(lm(y~,data=shuju[,-1]))  #181.53
stepAIC(lm(y~x3,data=shuju[,-1]))   #204.33
stepAIC(lm(y~x4,data=shuju[,-1]))   #172.52min
stepAIC(lm(y~x5,data=shuju[,-1]))  #182.18
stepAIC(lm(y~x6,data=shuju[,-1]))  #172.76
stepAIC(lm(y~x7,data=shuju[,-1])) #204.22
stepAIC(lm(y~x4+x1,data=shuju[,-1]))  #155.75min
stepAIC(lm(y~x4+x2,data=shuju[,-1]))  #174.14
stepAIC(lm(y~x4+x3,data=shuju[,-1]))  #172.9
stepAIC(lm(y~x4+x5,data=shuju[,-1]))  #173.71
stepAIC(lm(y~x4+x6,data=shuju[,-1]))  #168.67
stepAIC(lm(y~x4+x7,data=shuju[,-1]))  #174.35
stepAIC(lm(y~x4+x1+x2,data=shuju[,-1]))  #157.21
stepAIC(lm(y~x4+x1+x3,data=shuju[,-1]))  #157.5
stepAIC(lm(y~x4+x1+x5,data=shuju[,-1]))  #157.65
stepAIC(lm(y~x4+x1+x6,data=shuju[,-1]))  #148.55min
stepAIC(lm(y~x4+x1+x7,data=shuju[,-1]))  #155.63


setwd(“C:/Users/IBM/Desktop/研一课程/2.2回归分析/回归作业”)  #设定当前的工作目录
shuju=read.table(“shuju.txt”,header=T)
shuju   #读取数据
#采用AIC原则自动选择模型-后退法
library(MASS)
stepAIC(lm(y~x1+x2+x3+x4+x5+x6,data=shuju[,-1]))  #146.11min
stepAIC(lm(y~x1+x2+x3+x4+x5+x7,data=shuju[,-1]))  #159.96
stepAIC(lm(y~x1+x2+x3+x6+x5+x7,data=shuju[,-1]))  #153.58
stepAIC(lm(y~x1+x2+x6+x4+x5+x7,data=shuju[,-1]))  #149.11
stepAIC(lm(y~x1+x6+x3+x4+x5+x7,data=shuju[,-1]))  #148.08
stepAIC(lm(y~x6+x2+x3+x4+x5+x7,data=shuju[,-1]))  #165.79
stepAIC(lm(y~x1+x2+x3+x4+x5,data=shuju[,-1]))  #161.14
stepAIC(lm(y~x1+x2+x3+x4+x6,data=shuju[,-1]))  #151.88min
stepAIC(lm(y~x1+x2+x3+x5+x6,data=shuju[,-1]))  #153.83
stepAIC(lm(y~x1+x3+x4+x5+x6,data=shuju[,-1]))  #150.08
stepAIC(lm(y~x2+x3+x4+x5+x6,data=shuju[,-1]))  #163.79
stepAIC(lm(y~x1+x2+x3+x4,data=shuju[,-1]))  #159.17
stepAIC(lm(y~x1+x2+x3+x6,data=shuju[,-1]))  #152.87
stepAIC(lm(y~x1+x2+x4+x6,data=shuju[,-1]))  #150.49
stepAIC(lm(y~x1+x3+x4+x6,data=shuju[,-1]))  #150.2min
stepAIC(lm(y~x1+x3+x4,data=shuju[,-1]))  #AIC=157.5
stepAIC(lm(y~x1+x3+x6,data=shuju[,-1])) #AIC=155.12
stepAIC(lm(y~x1+x4+x6,data=shuju[,-1])) #AIC=148.55min

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

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

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/142279.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • PE结构

    PE结构PE文件是Windows操作系统下使用的可执行文件格式。它是微软在UNIX平台的COFF(通用对象文件格式)基础上制作而成。最初设计用来提高程序在不同操作系统上的移植性,但实际上这种文件格式

    2021年12月13日
    53
  • Spring Boot拦截器配置拦截登陆「建议收藏」

    Spring Boot拦截器配置拦截登陆「建议收藏」一,pom.xml的配置这里很简单,先引入spring-boot-starter-parent,parent是父模块,由父模块统一进行spring-boot版本管理,dependencies中与spring-boot启动绑定的包不需要再指定版本。<parent><groupId>org.springframework.boot<…

    2022年7月25日
    18
  • 博客作业03–栈和队列

    博客作业03–栈和队列

    2021年5月26日
    121
  • 大话IT第四期:苹果ipad2面世 竞争对手觊觎(音频)

    大话IT第四期:苹果ipad2面世 竞争对手觊觎(音频)

    2021年8月10日
    67
  • leetcode 792_leetcode5421编码

    leetcode 792_leetcode5421编码给定一个 m x n 二维字符网格 board 和一个字符串单词 word 。如果 word 存在于网格中,返回 true ;否则,返回 false 。单词必须按照字母顺序,通过相邻的单元格内的字母构成,其中“相邻”单元格是那些水平相邻或垂直相邻的单元格。同一个单元格内的字母不允许被重复使用。示例 1:输入:board = [[“A”,“B”,“C”,“E”],[“S”,“F”,“C”,“S”],[“A”,“D”,“E”,“E”]], word = “ABCCED”输出:true示例 2:输

    2022年8月9日
    5
  • dp3003打印机怎么清零_最小宽度dp

    dp3003打印机怎么清零_最小宽度dp有台奇怪的打印机有以下两个特殊要求:打印机每次只能打印由 同一个字符 组成的序列。每次可以在任意起始和结束位置打印新字符,并且会覆盖掉原来已有的字符。给你一个字符串 s ,你的任务是计算这个打印机打印它需要的最少打印次数。示例 1:输入:s = “aaabbb”输出:2解释:首先打印 “aaa” 然后打印 “bbb”。示例 2:输入:s = “aba”输出:2解释:首先打印 “aaa” 然后在第二个位置打印 “b” 覆盖掉原来的字符 ‘a’。 提示:1 <= s.le

    2022年8月9日
    5

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

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