cumulative_logit() defines a one-response ordinal model for ordered
categories. The first implemented path uses a location formula mu ~ ...
and ordered cutpoints with a fixed latent logistic scale. The location
intercept is dropped internally, as in standard cumulative-link models,
because a free location intercept and free cutpoints are not jointly
identifiable.
Details
The implemented contract is
Pr(y_i <= k) = logit^-1(theta_k - mu_i), with
mu_i = X_mu[i, ] beta_mu and
theta[1] < theta[2] < ... < theta[K - 1]. fitted() returns the expected
ordered-category score, sum_k k * Pr(y_i = k). Ordinal scale or
discrimination formulas are planned but not exposed in this first
implementation.