R语句的小问题~

3205

收藏 2014-08-11

states<-as.data.frame(state.x77[,c("Murder","Population","Illiteracy","Income","Frost")])
fit<-lm(Murder~Population+Illiteracy+Income+Frost,data=states)
install.packages("car")
library(car)
vif(fit)
sqrt(vif(fit))>2 # problem?

因为刚接触R语言，没看到过这样的语句:sqrt(vif(fit))>2 # problem?，主要是# problem? 有人能给我讲解一下么？

附件列表

QQ图片20140811095512.jpg

原图尺寸 10.31 KB

QQ图片20140811095512.jpg

原图尺寸 10.31 KB

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

crystal8832

2014-8-11 10:14:58

#problem？
是注释，R中是不会执行的，如果是别人的程序，可能是认为该命令存在问题，所以自己加了这样的一句话。
你这个我没看错的话应该是在求方差膨胀因子，就是在判断多重共线问题，后面那句是认为当VIF大于4认为有多重共线问题

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

一路嘿嘿

2014-8-11 10:17:23

sqrt(vif(fit))是一个向量，

复制代码

返回相同长度向量，如果元素大于2就返回TRUE,否则就返回FALSE

#号后面代表注释，建议你看一些R的入门书吧，这种问题最好不要问别人，最基础的问题

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

liqian0614

2014-8-11 10:25:46

crystal8832 发表于 2014-8-11 10:14
#problem？
是注释，R中是不会执行的，如果是别人的程序，可能是认为该命令存在问题，所以自己加了这样的一 ...

这样啊，我还以为#有新的含义，不只是注释

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

liqian0614

2014-8-11 10:26:25

一路嘿嘿发表于 2014-8-11 10:17
sqrt(vif(fit))是一个向量，返回相同长度向量，如果元素大于2就返回TRUE,否则就返回FALSE

#号后面代表注 ...

你说的这个我懂的，我也知道#号后面代表注释，但是他写在代码里，我以为有别的含义

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

liqian0614

2014-8-11 13:39:48

crystal8832 发表于 2014-8-11 10:14
#problem？
是注释，R中是不会执行的，如果是别人的程序，可能是认为该命令存在问题，所以自己加了这样的一 ...

在检测强影响点的时候有这样一段代码：
states<-as.data.frame(state.x77[,c("Murder","Population","Illiteracy","Income","Frost")])
fit<-lm(Murder~Population+Illiteracy+Income+Frost,data=states)
cutoff<-4/(nrow(states)-length(fit$coefficients)-2)
plot(fit,which=4,cook.levels=cutoff)
abline(h=cutoff,lty=2,col="red")

请问您能帮我解析下第五条语句么？我把cook.levels=cutoff删掉，发现结果也没有任何变化，那这语句又起怎样的作用呢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群