I just used the "margeff" module to compute the average partial effect (APE), and found the following problem:
Suppose a probit mode is: 
Y* = a1*x1 + a2*x2 + e;
y = 1 if Y* > 0; and = 0 if otherwise.
I would think that a1/a2 must be equal to APE(x1)/APE(x2) because 
APE(x1) = 1/n *sum(a1*normpdf(a1*x1+a2*x2))
APE(x2) = 1/n *sum(a2*normpdf(a1*x1+a2*x2))
But the APE estimates produced by margeff does not satisfy this property.  Anyone can advise?
Many thanks.
Wei Jiang
      
*
*   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/