Displaying 1 result from an estimated 1 matches for "pvaluemat".
Did you mean:
pvalue_ct
2012 Mar 28
1
resampling for correlation and testing
...dx<-sample(1:10, 10)
idx
# [1] 8 4 5 7 1 9 2 10 6 3
rho.B<-cor(B[,idx])
} ## repeat this x time (currently 500)
## in essence we then have the following :
rho.arrayB<-array(runif((10*10)*500), dim=c(10,10,500))
## Then test if rho.A[1,1] come from the distribution of rho.B[1,1]
pvalueMat[1,1]<-wilcox.test(rho.array[1,1,] , rho.A[1,1])$p.value
However, my array size would be 2300 x 2300 x 500 which R won't let me even make as an empty structure. Any suggestion are more than welcomed !!
Cheers,
Paul