码迷,mamicode.com
首页 > 其他好文 > 详细

R vif() 以及 step()

时间:2020-07-01 00:07:17      阅读:74      评论:0      收藏:0      [点我收藏+]

标签:desktop   main   图片   plot   com   mic   选择   loading   car   

library(package="car")
data <- read.csv("C:\\users\\mike1\\desktop\\data\\RData\\hald水泥数据.csv", header=T, sep=",",fileEncoding = "utf-8", stringsAsFactors = FALSE)
data1 <- data[2:dim(data)[2]]
data1

cor(data1)
scatterplotMatrix(data1, main="correlation")

result <- lm(y~., data=data1)
summary(result)

vif(result)  # we can directly use the function not like python class
step(result,direction="forward")
step(result,direction="backward")
step(result, direction="both")

myStep <- step(result, direction="both")
myStep1 <- step(result,direction="backward",trace=2)


summary(myStep)
summary(myStep1)



vif(myStep)

 

共线性,显示各变量之间有强相关,vif()函数在 car包中, 而step() 函数内置。

技术图片

 

 

偏相关图

技术图片

 

相关系数图:

技术图片

 

逐步回归图

技术图片

 

 

技术图片

 

 

以上只截取了部分图,但是结果与书上的不一样。最后虽然使用的是逐步回归,但是有一个系数不显著,不知道为什么没把 x4.1这个变量去掉,只能手动选择x1, x2 作为变量。

 

R vif() 以及 step()

标签:desktop   main   图片   plot   com   mic   选择   loading   car   

原文地址:https://www.cnblogs.com/zijidefengge/p/13216852.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!