[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

From |
rgutierrez@stata.com (Roberto G. Gutierrez, StataCorp.) |

To |
statalist@hsphsun2.harvard.edu |

Subject |
st: Re: -_diparm- (and its tanh option) |

Date |
Thu, 24 Oct 2002 12:06:40 -0500 |

Stephen P. Jenkins <stephenj@essex.ac.uk> asks: > -_diparm- is a very useful utility for displaying ancillary parameters. > These parameters are often estimated in a transformed metric, but we are > interested in the estimate in the original metric. Handily, -_diparm- > includes several standard transformations as options. My questions are: > 1. In the example below, why does the value for the standard error of the > ancillary parameter "rho" that gets pumped out by -_diparm- with the tanh > option differ from the value that I calculated by hand. [details omitted] You had g(t) = (exp(2*t)-1) / (exp(2*t)+1) g'(t) = 2*sqrt(exp(2*t))/((exp(2*t)+1) Actually, g'(t) = 4*exp(2*`t')/((exp(2*`t')+1)^2) Even thought you are calculating a standard error, you don't take the square root of g'(t). The variance of the transformation is (old var.)*{g'(t)}^2. . di " std error of rho = " `set'*4*exp(2*`t')/((exp(2*`t')+1)^2) std error of rho = .22367534 which agrees with your output from -biprobit-. > 2. More generally, why does -_diparm- /not/ produce the "z" and "P > z" > values? ... and yet it produces a c.i.? [I, for one, would like to see "z" > and "P > z" reported as well.] Because, generally when you transform a parameter from the estimation using _diparm, a test of H_o:transformed parameter = 0 is not generally of any interest. For example, if you exponentiate a parameter, a test of H_o:exp(beta) = 0 is not as interesting as say a test of H_o:exp(beta)=1, which is asymptotically equivalent to a test of H_o:beta=0, and that test is already provided. In the case of -biprobit-, a test of H_o: rho = 0 is asymptotically equivalent to a test of H_o: athrho = 0, provided in the line just above -rho-. In this case, had we wanted to display a test statistic for H_o:rho = 0, we would just be displaying the "same" information twice, where by "same" I mean asymptotically the same. Which of the two tests has the better small sample properties? Who knows? We merely display the version of the test done in the original estimation metric. Programmers should note, however, that specifying the -prob- option to -_diparm- will override the default behavior of not displaying the z-test for transformed coefficients. --Bobby rgutierrez@stata.com * * For searches and help try: * http://www.stata.com/support/faqs/res/findit.html * http://www.stata.com/support/statalist/faq * http://www.ats.ucla.edu/stat/stata/

- Prev by Date:
**Re: st: -_diparm- (and its tanh option)** - Next by Date:
**RE: st: -_diparm- (and its tanh option)** - Previous by thread:
**st: Re: statalist-digest V4 #1047** - Next by thread:
**RE: st: -_diparm- (and its tanh option)** - Index(es):

© Copyright 1996–2014 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |