使用二进制comparsionřchisq.test（）数据帧上(R chisq.test() on

我想要做的尺寸（50x752）的数据帧一chisq.test。我想对所有可能的PAIRE比较中的所有列的p值（由多个测试调整）。最后我想回去矩阵（50×50），以产生调整CHISQ p值的热图。这是我做的时刻，但这个远远beeing理想。

第一步：做pairewise比较

function(data,p.adjust.method="holm")
{
cor.mat <- cor(data)
x<-ncol(data)#nb of column in matrix here 50
y<-nrow(data)#nb of column in matrix here 758
index<-t(combn(x, 2)) #create the matrix position of output for all possible combination
nindex <- nrow(index)
pvals <- numeric(nindex)

for (i in 1:nindex)
{
pvals[i]<-chisq.test(data[, index[i, 1]], data[, index[i,2]])$p.value   
}
pvals<-p.adjust(pvals,method = p.adjust.method)
out <- as.data.frame(cbind(index, pvals))
}

步骤2：该输出表被变换成一个矩阵使用

   dcast(df,V2~V1,fill=1) # thanx to Roland for this function!

但是，这都不尽如人意，因为我不反映在最后的矩阵p值，我不得不操纵第一函数的输出来获得充满0对角线（列与自己进行比较时）。对你的帮助表示感谢！

像这样？

#some data
set.seed(42)
df <- data.frame(a=rbinom(1000,5,0.3),
                 b=rbinom(1000,5,0.001),
                 c=rbinom(1000,5,0.1),
                 d=rbinom(1000,5,0.9))

#function to calculate the adj. p-value
fun <- function(x,y) {
  p.adjust(chisq.test(df[,x],df[,y])$p.value,method="holm",n=choose(ncol(df),2))
}

p.adj <- outer(names(df),names(df),FUN=Vectorize(fun)) #use outer to get a matrix
diag(p.adj) <- 1  #you should find out why chisq.test returns zero at the diag
rownames(p.adj) <- names(df)
colnames(p.adj) <- names(df)

p.adj
#  a         b c         d
#a 1 1.0000000 1 1.0000000
#b 1 1.0000000 1 0.6152165
#c 1 1.0000000 1 1.0000000
#d 1 0.6152165 1 1.0000000