我在R中创建一个情节,使用:
plot(IQ, isAtheist)
abline(lm(isAtheist~IQ))
IQ为numeric
,isAtheist
为布尔值,值为TRUE
或FALSE
。
我试着写:
cor(IQ, isAtheist)
但它给了我一个错误:
Error in cor(IQ, isAtheist) : 'x' must be numeric
如何确定这两个变量之间的相关性?
答案 0 :(得分:5)
在这种情况下,我真的不知道你想如何解释相关性,但你可以试试cor(IQ, as.numeric(isAtheist))
。
在这种情况下,TRUE将为1且为0。
答案 1 :(得分:2)
这是我认为您可能想要的(显示叠加在箱线图上的平均IQ值的差异):
plot(IQ~isAtheist)
lines(x=c(1,2), y=predict( lm(IQ~isAtheist),
newdata=list(isAtheist=c("NO","YES") ) ) ,
col="red", type="b")
plot.formula默认情况下的X位置为as.numeric(factor(isAtheist))
,即1和2而不是0和1,这是您在使用abline
时所假设的。推断超出这些值是没有意义的,所以我选择绘制为有界线段。我将添加一个有效的例子和输出。
set.seed(123)
isAtheist=factor(c("NO","YES")[1+rep( c(0,1), 50 )])
plot(IQ~isAtheist)
lines(x=c(1,2), y=predict( lm(IQ~isAtheist),
newdata=data.frame(isAtheist=c("NO","YES") ) ) ,
col="red", type="b")