Internal function: generate probability matrix for ordinal outcomes given thresholds, linear predictor and link function
estimate_p_mat(theta, z, link)numeric vector of thresholds
numeric vector of linear predictor
The link function used to model ordinal outcomes. Default is
"logit" for proportional odds model. Other options are
"cloglog" (proportional hazards model) and "probit".