r - 2 个向量的样本协方差

Question

我正在尝试计算这两个向量之间的样本协方差。我定义了一个带有两个输入变量的函数。不知道对不对？我的样本协方差公式也不会运行。谁能帮我用R写出来？

  xv = c(1., 5.5, 7.8, 4.2, -2.7, -5.4, 8.9)
  yv = c(0.1, 1.5, 0.8, -4.2, 2.7, -9.4, -1.9)
  sampleCov= function(x,y){ 
    cov(xv,yv) = frac{sum_{i=1}^{n}(x_i-\mu_x)(y_i-\mu_y)}{n-1}].
    return (Cov(xv,yv)
  }

Answer 1

R 中有一个基函数叫做 cov这正是您想要的，但是如果您想编写一个函数(无需这样做)，您可以尝试以下操作:

COV<- function(x,y) {
  if(length(x)!=length(y)) {stop('x must have the same length as y ')}
  x.bar <- mean(x)
  y.bar <- mean(y)
  N <- length(x)

  Cov <- (sum((x-x.bar)*(y-y.bar))) / (N-1)
  return(Cov)
}

COV(xv, yv)
[1] 8.697381

cov(xv, yv)
[1] 8.697381

如您所见 COV给出与 cov 相同的结果所以你不必为此编写函数。

此外，您的函数体没有 R 语法，而是您编写了不同的 LaTex 语法。