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