occor<-WGCNA::corAndPvalue(t(x)/colSums(x),method = 'pearson') mtadj<-multtest::mt.rawp2adjp(unlist(occor$p),proc='BH') adpcor<-mtadj$adjp[order(mtadj$index),2] occor.p<-matrix(adpcor,dim(t(x)/colSums(x))[2]) ## R value occ...