为什么不是命令
给出的相关系数cor(t,g)
并且由命令
给出summary(tgmodel, correlation=TRUE)
运行后相同:
t<-c(0,1.2,2.3,3,4,5.2,6.3,7,8)
g<-c(12,10,8,11,6,7,2,3,3)
tgmodel<-lm(g~t)
答案 0 :(得分:1)
他们之所以不同,是因为他们在不同的事物之间存在相关性:
cor()
显示输入变量t
和g
之间的相关性。summary(lm(...), correlation=TRUE)
显示估计参数之间的相关性,即斜率和截距。如果仔细检查summary()
的输出,您会发现它显示t
和g
之间的相关系数的平方为Multiple R-squared
:
> summary(lm(g~t))
...
Multiple R-squared: 0.8357, Adjusted R-squared: 0.8122
...
> cor(t,g)**2
[1] 0.8356938