st: simulation to derive power for GLM multiple regression

[Jake <daedalus702@yahoo.com>]

Hello,

I am interested in learning details about how to conduct simulations to calculate power for a test of a single coefficient in a GLM (negative binomial) multiple regression model.  I am somewhat familiar with the relevant methods outlined by Feiveson (2002).  However, what I am curious about which Feiveson does not discuss is how best to simulate when there are k covariates, in addition to the coefficient of interest.  

I assume I would treat the covariates as fixed -- using the real (already collected) data in the simulation, in other words.  Then I imagine that I would simulate only the variable for the coefficient of interest.  I'm not certain of how to generate this random variable so that it would be (asymptotically) correlated with all of the other variables in specified ways.  I imagine I would use the cholesky method?  But many of the variables are not normally distributed -- some are dichotomous, etc.

I would appreciate any help that you might provide.



Thanks,

Jacob Felson
Assistant Professor
Dept. of Sociology
William Paterson University

Reference

Feiveson, A.H. 2002 "Power by simulation."  The Stata Journal 2: 107-124.
http://www.stata.com/support/faqs/stat/power.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/