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

From |
<Colin.Vance@dlr.de> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
st: gllapred vs gllasim take two |

Date |
Fri, 24 Jun 2005 20:29:07 +0200 |

Greetings, This question follows up on a previous inquiry about the difference between gllapred and gllasim but with a few additional details. I'm estimating a probit model on panel data and would like to use the results to simulate outcomes over different values of one of the explanatory variables, holding the other variables fixed at their mean. My question is whether gllapred or gllsim is a better tool for this (and my hunch is the latter). Specifically, I'd like to do the following. 1. Estimate a panel probit model using: gllamm y age other1 other2, i(persid) family(binom) link(probit) nip(20) 2. Holding other1 and other2 fixed at their mean values, generate predicted values of y over different integer values of age, say, ranging between 18 and 65. To implement step 2, I've created an artificial data set in which other1 and other2 are fixed at their mean values and age varies over the range of interest. So I first estimate the model on the real data. Then I open the artificial data and type: gllasim pred1, linpred fsample The above seems like the right approach given my objective, but I've noted a few interesting things: 1. the approach doesn't work using gllapred unless the dependent variable is included in the artificial data (which I guess has to do with the fact that gllapred includes empirical Bayes). 2. gllapred produces the same answer when used repeatedly, whereas gllasim always produces a different answer, suggesting that the latter may be sampling from a distribution of the parameter estimates I guess point 2 would actually be another good reason to use gllsim in my case. I could set up a simple code to implement the command, say 1000 times, and, after taking the mean, I'd have something akin to a Monte Carlo simulation of a predicted value. One problem is that there is very little documentation of gllasim in the manual or elsewhere, so it's hard to know if my idea is even in the ballpark. Any insights offered would be greatly appreciated. Many thanks, Colin Colin Vance, Ph.D. German Aerospace Center Rutherfordstrasse 2 12209 Berlin Germany

`<<winmail.dat>>`

**References**:**st: ivreg2 revision available***From:*Kit Baum <baum@bc.edu>

**Re: st: ivreg2 revision available***From:*Roger Harbord <roger.harbord@bristol.ac.uk>

**Re: st: ivreg2 revision available***From:*Kit Baum <baum@bc.edu>

- Prev by Date:
**st: Tests for Biprobit** - Next by Date:
**st: reshaping data** - Previous by thread:
**Re: st: ivreg2 revision available** - Next by thread:
**st: How to save "within-group inequality index"** - Index(es):

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