#delimit; version 2.1; use example.dta; logit choice story1 story2 story3 story4 lvratio lcratio,tab nocon; predict predict; gen logodds = log(predict/(1-predict)); lab var logodds "PREDICTED LOGIT FUNCTION"; sort story logodds; graph choice predict logodds, symbol(oi) connect(.L) jitter(2) xline(0) yline(.5) l1("CHOICE") t1("GRAPH OF LOGIT FUNCTION VS. ACTUAL CHOICE");