MML estimator (part 3a, prepare for gradient); bug fix polychor with identical thresholds; avoid NULL arguments skewness/kurtosis in simulateData()