<>
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/