double LS::draw_z(bool y, double eta)const{ double trun_prob = plogis(0, eta); double u = y ? runif(trun_prob,1) : runif(0,trun_prob); return qlogis(u,eta); }
double OrdinalLogitImputer::impute( RNG &rng, double eta, double lower_cutpoint, double upper_cutpoint) { return eta + qlogis(runif_mt( rng, plogis(lower_cutpoint - eta), plogis(upper_cutpoint - eta))); }
double F77_SUB(invcdflogis)(double *p, double *location, double *scale, int *lower_tail, int *log_p) { return qlogis(*p, *location, *scale, *lower_tail, *log_p); }
// Convert from the probabliity scale to the logit (log odds) scale. inline double logit(double prob) { return qlogis(prob); }