0

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

我的代码

变量

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)

不工作

谢谢

4

0 回答 0