<>
Nick pointed out invlogit(), which does not appear in the function  
list of density functions, but of math functions. Quite so, but it  
doesn't save much typing:
local cum1  1/(1+exp(-((_b[_cons]+ _b[trunk]*mut + _b[mpg]*mum +  
_b[d]*1) - mu)/s))
scalar onea = invlogit((_b[_cons]+ _b[trunk]*mut + _b[mpg]*mum +  
_b[d]*1 - mu)/s)
will indeed generate the same values, but you still have to compute mu  
and s and include them in the function argument.
In my prior example I had the line
gen d = (runiform()>0.7)
which did not appear in the posting.
Kit Baum   |   Boston College Economics & DIW Berlin   |   http://ideas.repec.org/e/pba1.html
                              An Introduction to Stata Programming   
|   http://www.stata-press.com/books/isp.html
   An Introduction to Modern Econometrics Using Stata  |   http://www.stata-press.com/books/imeus.html
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/