matrixPval {GIGSEA} | R Documentation |
matrixPval calculates the p values for the correlation coefficients based on t-statistics
matrixPval(r, df)
r |
a vector or a matrix of Pearson correlation coefficients taking values in [-1,+1] |
df |
the degree of freedom |
a vector or matrix of p values taking values in [0,1]
r <- cor(USArrests) df <- nrow(USArrests) - 2 pval1 <- matrixPval(r,df) pval2 <- matrix(ncol=ncol(USArrests),nrow=ncol(USArrests),data=0) for(i in 1:ncol(USArrests)) { for(j in 1:ncol(USArrests)) { pval2[i,j] <- cor.test(USArrests[,i],USArrests[,j])$p.val } } head(pval1) head(pval2)