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)
上一篇 2022年5月24日 上午7:40
下一篇 2022年5月24日 上午7:40


相关推荐

  • 初学SSTI

    初学SSTISSTI 的基础学习

    2026年3月16日
    1
  • 无名汉化组官网_什么是并查集

    无名汉化组官网_什么是并查集永无乡包含 n 座岛,编号从 1 到 n ,每座岛都有自己的独一无二的重要度,按照重要度可以将这 n 座岛排名,名次用 1 到 n 来表示。某些岛之间由巨大的桥连接,通过桥可以从一个岛到达另一个岛。如果从岛 a 出发经过若干座(含 0 座)桥可以到达岛 b ,则称岛 a 和岛 b 是连通的。现在有两种操作:B x y 表示在岛 x 与岛 y 之间修建一座新桥。Q x k 表示询问当前与岛 x 连通的所有岛中第 k 重要的是哪座岛,即所有与岛 x 连通的岛中重要度排名第 k 小的岛是哪座,请你输出那

    2022年8月10日
    5
  • Java分布式锁

    Java分布式锁Java分布式锁我的理解应该叫集群锁或者跨实例锁锁的作用是在多线程情况下,控制线程同步访问变量,执行代码块、方法,例如synchronized,在单个jvm进程中,这样是奏效的。但是在分布式环境中,单个服务往往都是要部署多台实例的,在有多个jvm进程的集群里,synchronized就达不到我们的要求了。synchronized只能控制当前jvm进程中的线程,对于其它jvm进程中的线程,它无能为力。也就是说有可能一个jvm中的线程是同步执行的,在此过程中,或许会有集群里其它jvm的线程执行到

    2022年6月10日
    48
  • Android在屏幕中模拟时钟

    Android在屏幕中模拟时钟AnalogClock 组件用来在 Android 中显示模拟时钟 可以在 XML 布局文件中通过标记添加 格式如下 amp lt AnalogClock 属性列表 amp gt 布局文件代码如下 amp lt xmlversion 1 0 encoding utf 8 amp

    2026年3月17日
    1
  • 网页中部分非主流代码[通俗易懂]

    网页中部分非主流代码[通俗易懂]一:颜色代码0000000000200000400000600000800000a00000c00000ff…

    2022年7月2日
    50
  • fastjson的json字符串转List

    fastjson的json字符串转Listfastjson的json字符串转Listpom依赖<dependency><groupId>com.alibaba</groupId><artifactId>fastjson</artifactId><version>1.2.62</version></dependency>json字符串转List:JSONObject.paresArray(str,class)Str

    2022年5月15日
    46

发表回复

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

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