R语言验证及协方差的计算公式
作者：Leonopteryxw
这篇文章主要介绍了R语言验证及协方差的计算公式的相关资料,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
协方差的计算公式及R语言进行验证
首先附上协方差公式：
来设5个样本点：（3,9），（2,7），（4,12），（5,15），（6,17）
用R绘制出散点图，大概是这样：
要求这5个点的协方差，首先样本点为5个，n=5，X依次取3,2,4,5,6，Y依次取9,7,12,15,17。X的均值为4，带入公式可得：
不难计算出结果为6.5
现在用R语言进行验证：
已知R语言里边协方差函数为cov(x,y)
我们分别用cov()函数和上述公式来进行仿真结果，代码如下：
a <- c(3,2,4,5,6)
b <- c(9,7,12,15,17)
COV=0
EX=mean(a)
EY=mean(b)
for(j in 1:5){
COV <- COV+(a[j]-EX)*(b[j]-EY)/4
}
COV
cov(a,b)
输出结果如下：
> COV
[1] 6.5
> cov(a,b)
[1] 6.5
由此可得，计算公式得出的结果完全正确
到此这篇关于R语言验证及协方差的计算公式的文章就介绍到这了,更多相关R语言协方差计算内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！