代码之家  ›  专栏  ›  技术社区  ›  Ian

如何找到R中每n行的协方差

  •  0
  • Ian  · 技术社区  · 8 年前

    我有两大组数据,每个都有2000多个数据,并试图找到每5行的协方差。

    x=c(1,2,3,4,5)
    y=c(6,7,8,9,10)
    df=data.frame(x,y)
    group=rep(1:length(df),each=2,length=length(df))
    

    cov(x[1:2,],y[1:2,])
    

    cov(x[3:4,],y[3:4,])
    
    1 回复  |  直到 8 年前
        1
  •  0
  •   AK88    8 年前
    library(zoo)    
    x = c(1,2,3,4,5)
    y = c(6,7,8,9,10)
    rows = 2
    out = rollapply(data.frame(x,y), rows, function(x) cov(x[,1],x[,2]),
                        by.column=FALSE)
    out