Hello.

Using the polychor function
> polychor(data[c(s1,s2)] )
for polychoric correlations of two ordinal variables in R takes a long time
for N=7000 (20 minutes+) and significantly slows down my computer.

Now, I have a pretty old computer, but it takes about 20 seconds for MPLUS
to print out the complete polychoric correlation matrix for all 16
variables, while I am running the R function only for 2 of the 16 variables.

Why is that? Can that process be speeded up? What makes a polychoric
correlation so much more computationally intensive than a pearson
correlation?

Thank you

	[[alternative HTML version deleted]]

