我正在尝试为不平衡数据实现相关矩阵,这可能吗?我下面的代码显示了我一直在尝试做的事情。有人可以帮我解决这个问题吗?
我的代码
变量
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)
不工作
谢谢