不同大小的向量之间的相关矩阵

时间:2016-03-17 15:07:57

标签: r cross-correlation pearson-correlation

我正在尝试为不平衡数据实现相关矩阵,这有可能吗?我的代码显示我一直在努力做到。有人可以帮我解决这个问题吗?

我的代码

变量

var1<-rnorm(200,5,0.25)##  200 values
var2<-1:500 ## 500 values
var3<-rnorm(100,5,0.25)## 100 values
var4<-500:1## 500 values

相关矩阵函数

pn <- function(X){crossprod(!is.na(X))}

cor.prob <- function(X){
     pair.SampSize <- pn(X)
     above1 <- row(pair.SampSize) < col(pair.SampSize)
     pair.df <- pair.SampSize[above1] - 2

     R <- cor(X, use="pair")
     above2 <- row(R) < col(R)
     r2 <- R[above2]^2
     Fstat <- (r2 * pair.df)/(1 - r2)
     R[above2] <- 1 - pf(Fstat, 1, pair.df)
     R
     }


correla <- round(cor.prob(var1,var2,var3,var4),4)

不起作用

由于

0 个答案:

没有答案