I describe a command that simultaneously solves the extended estimating equations estimator for parameters in the link and variance functions along with those of the linear predictor in a generalized linear model. The method addresses difficulties in choosing the correct link and variance functions in these models. It decouples the scale of estimation for the mean model, determined by the link function, from the scale of interest for the scientifically relevant effects. It also estimates a flexible variance structure from the data, leading to efficient estimation.
BasuA., and RathouzP.J.2005. Estimating marginal and incremental effects on health outcomes using flexible link and variance function models. Biostatistics6(1): 93–109.
2.
BloughD. K., MaddenC. W., and HornbrookM.C.1999. Modeling risk using gener-lized linear models. Journal of Health Economics18: 153–171.
3.
BoxG. E. P., and CoxD.R.1964. An analysis of transformations. Journal of the Royal Statistical Society, Series B26: 211–252.
4.
CoxN. J.2000. dm79: Yet more new matrix commands. Stata Technical Bulletin56: 4–8. In Stata Technical Bulletin Reprints, vol. 10, 17–23. College Station, TX: Stata Press.
5.
GuanW., and GutierrezR.G.2002. Programmable GLM: two user-defined links. Stata Journal2(4): 378–390.
6.
HardinJ. W., and HilbeJ.M.2001. Generalized Linear Models and Extensions.College Station, TX: Stata Press.
7.
HosmerD. W.Jr., and LemeshowS.2000. Applied Logistic Regression.2nd ed. New York: Wiley.
8.
HuberP. J.1967. The behavior of maximum likelihood estimates under nonstandard conditions. In Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability, vol. 1, 221–233. Berkeley, CA: University of California Press.
9.
ManningW. G., and MullahyJ.2001. Estimating log models: to transform or not to transform?Journal of Health Economics20: 461–494.
10.
McCullaghP., and NelderJ.A.1989. Generalized Linear Models.2nd ed. London: Chapman & Hall.
11.
NelderJ. A., and PregibonD.1987. An extended quasi-likelihood function. Biometrika74: 221–232.
12.
StataCorp.2005. Stata 9 Base Reference Manual.College Station, TX: Stata Press.
13.
WedderburnR. W. M.1974. Quasilikelihood functions, generalized linear models, and the Gauss–Newton method. Biometrika61(3): 439–447.
14.
WhiteH.1980. A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. Econometrica48: 817–838.
15.
WooldridgeJ. M.1992. Some alternatives to the Box–Cox regression model. International Economic Review33: 935–955.
