How to get polychoric corr. with significance

Hi all
I am new to this forum.

I have been desperately looking for a way to compute a polychoric correlation matrix, with significance in R. If that is very hard then polychoric correlation between two variables with significance would be sufficient.

What I have tried so far:
poly <- polychor(var1,var2)
poly <- polychor(DatM) #where DatM is a DF converted to matrix

hetcor(Dat2) #I am however uncertain hetcor is something I would want if I am looking for polychoric correlation.

polychoric(Dat$for2a,smooth=TRUE,global=TRUE,polycor=FALSE, ML = FALSE, std.err=TRUE)
None of these give significance. I have read that the Permutation test [could give me significance][1] for the corr=0 hypothesis for each correlation pair. The packages coin and lmPerm do allow one to compute the Permutation test. However, I am unsure how.

thanks in advance, Helgi