The model still conditioned on the information from the data with mu=mu[y_argsort], combining with ordered transformation (the latent mu need to satisfy the ordering condition after indexing) is how it fit to the data.
1 Like