DOCSLIB.ORG

Negative Binomial Regression

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text

Load more

STATISTICS COLUMN Negative Binomial Regression Shengping Yang PhD ,Gilbert Berdine MD In the hospital stay study discussed recently, it was mentioned that “in case overdispersion exists, Poisson regression model might not be appropriate.” I would like to know more about the appropriate modeling method in that case. Although Poisson regression modeling is wide- regression is commonly recommended. ly used in count data analysis, it does have a limita- tion, i.e., it assumes that the conditional distribution The Poisson distribution can be considered to of the outcome variable is Poisson, which requires be a special case of the negative binomial distribution. that the mean and variance be equal. The data are The negative binomial considers the results of a se- said to be overdispersed when the variance exceeds ries of trials that can be considered either a success the mean. Overdispersion is expected for contagious or failure. A parameter ψ is introduced to indicate the events where the first occurrence makes a second number of failures that stops the count. The negative occurrence more likely, though still random. Poisson binomial distribution is the discrete probability func- regression of overdispersed data leads to a deflated tion that there will be a given number of successes standard error and inflated test statistics. Overdisper- before ψ failures. The negative binomial distribution sion may also result when the success outcome is will converge to a Poisson distribution for large ψ. not rare. To handle such a situation, negative binomial 0 5 10 15 20 25 0 5 10 15 20 25 Figure 1. Comparison of Poisson and negative binomial distributions. Figure 1 shows that when ψ is small (e.g., ψ=5), Corresponding author: Shengping Yang a negative binomial distribution is more spread than Contact Information: [email protected]. a Poisson distribution with the same mean. However, DOI: 10.12746/swrccc2015.0310.135 when ψ is large (e.g., ψ =500), the two distributions mostly overlap. 50 The Southwest Respiratory and Critical Care Chronicles 2015;3(10) Shengping Yang et al. Negative Binomial Regression 1. The negative binomial regression model. where the conditional variance is greater than the conditional mean of the outcome variable (overdisper- Negative binomial distribution is defined as a sion). In addition, it is clear that when ψ is very large, discrete distribution of the number of successes in a µ 2 i → 0, Var ( y ; ψ, µ )~ ~ µ = E ( y ; ψ, µ ) , and the sequence of independent and identically distributed ψ i i i i i Bernoulli trials before a specified number of failures expected value or mean converges to the variance. are observed. More intuitively, it can be viewed as a Poisson distribution with parameter λ, where λ itself is 2. Application of negative binomial regression in not fixed but a random variable which follows a Gam- count data analysis. ma distribution. The Gamma distribution is a contin- uous probability function with a shape parameter ψ, Applying a negative binomial regression to a rate parameter δ, and the Gamma function Г. The model count data with overdispersion is straightfor- physical meaning of the Gamma distribution is an av- ward using available statistical software. For exam- erage or expected waiting time for a random event. ple, the SAS code (see below) for negative binomial regression modeling is very similar to that for Poisson Now, suppose that variable vi follows a Gamma dis- tribution that regression. In the hospital length of stay study, the outcome variable is LOS (length of stay), and the risk factors of interest are corticosteroid (whether or not k(v ;ψ;δ) = δ ψ v ψ-1 e -vi δ, where v > 0, δ > 0, ψ > 0, i i i a patient was treated with corticosteroids within 60 Г (ψ) minutes), race, and gender. Comparing with Poisson regression, the only change is to replace the dist (dis- 2 and E(vi) = ψ / δ , Var(vi) = ψ / δ . By setting ψ = δ, tribution) option “poisson” with “nb”, which is the ab- and transforming the one parameter Gamma distribu- breviation for “negative binomial”. tion as a function of the Poisson mean λ , we get ψ i ψ ψ-1 -λi /µ proc genmod; g(λi ;ψ; µi ) = ( ψ / µi ) λi e i . Г (ψ) class corticosteriods race gender; model LOS = corticosteriods race gender To get the marginal distribution of the outcome <other risk factors>/dist=nb; run; variable yi , we have As previously described, the class statement is used to define that corticosteriods, race, and gender- are all categorical variables. And the dist=nb option is used to indicate that the conditional distribution of LOS is assumed to be negative binomial. ψ yi = yi + ψ - 1 ψ µi , Applying a negative binomial regression using R ( ψ - 1 ) µ + ψ µ + ψ software is also straightforward. Instead of using the ( i ) ( i ) glm function previously used for applying Poisson re- which is the probability mass function of a neg- gression, the glm.nb function, which is modified ver- ative binomial distribution (NB2). Note that sion of the glm function, can be readily used. Since the glm.nb function is specifically developed for nega- E(y µ ) = µ , and i ; ψ, i i tive binomial regression modeling, we do not need to 2 include the “family=” option. Var (yi ; ψ, µi ) = µi + µi , ψ The Southwest Respiratory and Critical Care Chronicles 2015;3(10) 51 Shengping Yang et al. Negative Binomial Regression glm.nb (LOS ~ corticosteroid + race + gender + of re-admission is 0 for these patients). If this is the <other risk factors>, data=data). cause of overdispersion, then zero-inflated regression is appropriate to model these data (Lambert, 1992). 3. Assumptions of Negative binomial regression. Considering that zero-inflated count data are Negative binomial regression shares many generated by a mixture of two statistical processes, common assumptions with Poisson regression, such i.e., the first one always generates zero counts and as linearity in model parameters, independence of the second one generates both zero and nonzero individual observations, and the multiplicative effects counts. Correspondingly, a logit model can be used to of independent variables. However, comparing with determine if the outcome of an individual observation Poisson regression, negative binomial regression al- is from the always-zero or not-always-zero groups, lows the conditional variance of the outcome variable and then a Poisson or negative binomial model can to be greater than its conditional mean, which offers be used to model the counts for the not-always-zero greater flexibility in model fitting. Note that negative group. We will not discuss the details of zero-inflated binomial regression does not handle the underdis- regression here; however, the implementation of ze- persion situation, where the conditional variance is ro-inflated regression is straightforward in SAS, and smaller than the conditional mean. Fortunately, un- the example code is derdispersion is rare in practice. proc genmod; 4. Interpretation of the result of a Poisson regres- class <corticosteriods risk factors>; sion. model count = <risk factors> /dist=zip; zeromodel = <risk factors>/link=logit; A statistical test is highly recommended after run; running a Poisson regression to determine whether overdispersion exists. One such test is to apply an The dist (distribution) option is defined as “zip” ordinary least square regression without the intercept (zero-inflated Poisson), and an addition statement, term (Cameron and Trivedi, 1996). The dependent zeromodel is used to model which group (either the 2 always-zero or the not-always-zero) an individual ob- variable is defined as (yi - yi ) - yi ,where yi = exp (Xi b) servation comes from. yi is the predicted value obtained using the Poisson re- In contrast, sometimes, the data-generating gression. By using the yi as the independent variable in the ordinary least square regression, a t test can be process does not allow 0s; for example, if we are in- performed, and a small p value (p<0.05) implies that terested in modeling the number of wounds for pa- overdispersion exists. tients at a trauma center, we might find that all the pa- tients have at least one wound. Under this situation, 5. Zero-inflated and zero-truncated Poisson/nega- a zero-truncated negative binomial regression can be tive binomial regressions. used. There are many potential causes of overdis- persion; one is the excessive zeros in the data. For example, suppose that the outcome of interest is the number of hospital re-admissions, and then it is quite obvious that there might be a large number of pa- tients who do not have any re-admission (the number 52 The Southwest Respiratory and Critical Care Chronicles 2015;3(10) Shengping Yang et al. Negative Binomial Regression Author affiliation : Shengping Yang is a biostatistician in the Department of Pathology at TTHUSC. Gilbert Berdine is a pulmonary physician in the Department of Internal Medicine at TTUHSC. Submitted: 3/26/2015 Accepted: 4/7/2015 Published electronically: 4/15/2015 References 1. Cameron AC, Trivedi P K. Count data models for financial data. Handbook of Statistics 1996, 14, Statistical Methods in Finance, 363-392, Amsterdam 2. Lambert D. Zero-inflated Poisson regression, with an application to defects in manufacturing. Technometrics 1992; 34(1): 1-14. 3. North-Holland C, Greenwood M, Yule GU. An enquiry into the nature of frequency distributions representative of multiple happenings with particular reference of multiple attacks of disease or of repeated accidents. J R Statist Soc 1920; 83: 25–279 The Southwest Respiratory and Critical Care Chronicles 2015;3(10) 53.
Recommended publications
  • Lecture 22: Bivariate Normal Distribution Distribution

    Lecture 22: Bivariate Normal Distribution Distribution

    6.5 Conditional Distributions General Bivariate Normal Let Z1; Z2 ∼ N (0; 1), which we will use to build a general bivariate normal Lecture 22: Bivariate Normal Distribution distribution. 1 1 2 2 f (z1; z2) = exp − (z1 + z2 ) Statistics 104 2π 2 We want to transform these unit normal distributions to have the follow Colin Rundel arbitrary parameters: µX ; µY ; σX ; σY ; ρ April 11, 2012 X = σX Z1 + µX p 2 Y = σY [ρZ1 + 1 − ρ Z2] + µY Statistics 104 (Colin Rundel) Lecture 22 April 11, 2012 1 / 22 6.5 Conditional Distributions 6.5 Conditional Distributions General Bivariate Normal - Marginals General Bivariate Normal - Cov/Corr First, lets examine the marginal distributions of X and Y , Second, we can find Cov(X ; Y ) and ρ(X ; Y ) Cov(X ; Y ) = E [(X − E(X ))(Y − E(Y ))] X = σX Z1 + µX h p i = E (σ Z + µ − µ )(σ [ρZ + 1 − ρ2Z ] + µ − µ ) = σX N (0; 1) + µX X 1 X X Y 1 2 Y Y 2 h p 2 i = N (µX ; σX ) = E (σX Z1)(σY [ρZ1 + 1 − ρ Z2]) h 2 p 2 i = σX σY E ρZ1 + 1 − ρ Z1Z2 p 2 2 Y = σY [ρZ1 + 1 − ρ Z2] + µY = σX σY ρE[Z1 ] p 2 = σX σY ρ = σY [ρN (0; 1) + 1 − ρ N (0; 1)] + µY = σ [N (0; ρ2) + N (0; 1 − ρ2)] + µ Y Y Cov(X ; Y ) ρ(X ; Y ) = = ρ = σY N (0; 1) + µY σX σY 2 = N (µY ; σY ) Statistics 104 (Colin Rundel) Lecture 22 April 11, 2012 2 / 22 Statistics 104 (Colin Rundel) Lecture 22 April 11, 2012 3 / 22 6.5 Conditional Distributions 6.5 Conditional Distributions General Bivariate Normal - RNG Multivariate Change of Variables Consequently, if we want to generate a Bivariate Normal random variable Let X1;:::; Xn have a continuous joint distribution with pdf f defined of S.
  • An Introduction to Poisson Regression Russ Lavery, K&L Consulting Services, King of Prussia, PA, U.S.A

    An Introduction to Poisson Regression Russ Lavery, K&L Consulting Services, King of Prussia, PA, U.S.A

    NESUG 2010 Statistics and Analysis An Animated Guide: An Introduction To Poisson Regression Russ Lavery, K&L Consulting Services, King of Prussia, PA, U.S.A. ABSTRACT: This paper will be a brief introduction to Poisson regression (theory, steps to be followed, complications and interpretation) via a worked example. It is hoped that this will increase motivation towards learning this useful statistical technique. INTRODUCTION: Poisson regression is available in SAS through the GENMOD procedure (generalized modeling). It is appropriate when: 1) the process that generates the conditional Y distributions would, theoretically, be expected to be a Poisson random process and 2) when there is no evidence of overdispersion and 3) when the mean of the marginal distribution is less than ten (preferably less than five and ideally close to one). THE POISSON DISTRIBUTION: The Poison distribution is a discrete Percent of observations where the random variable X is expected distribution and is appropriate for to have the value x, given that the Poisson distribution has a mean modeling counts of observations. of λ= P(X=x, λ ) = (e - λ * λ X) / X! Counts are observed cases, like the 0.4 count of measles cases in cities. You λ can simply model counts if all data 0.35 were collected in the same measuring 0.3 unit (e.g. the same number of days or 0.3 0.5 0.8 same number of square feet). 0.25 λ 1 0.2 3 You can use the Poisson Distribution = 5 for modeling rates (rates are counts 0.15 20 per unit) if the units of collection were 8 different.
  • Overdispersion, and How to Deal with It in R and JAGS (Requires R-Packages AER, Coda, Lme4, R2jags, Dharma/Devtools) Carsten F

    Overdispersion, and How to Deal with It in R and JAGS (Requires R-Packages AER, Coda, Lme4, R2jags, Dharma/Devtools) Carsten F

    Overdispersion, and how to deal with it in R and JAGS (requires R-packages AER, coda, lme4, R2jags, DHARMa/devtools) Carsten F. Dormann 07 December, 2016 Contents 1 Introduction: what is overdispersion? 1 2 Recognising (and testing for) overdispersion 1 3 “Fixing” overdispersion 5 3.1 Quasi-families . .5 3.2 Different distribution (here: negative binomial) . .6 3.3 Observation-level random effects (OLRE) . .8 4 Overdispersion in JAGS 12 1 Introduction: what is overdispersion? Overdispersion describes the observation that variation is higher than would be expected. Some distributions do not have a parameter to fit variability of the observation. For example, the normal distribution does that through the parameter σ (i.e. the standard deviation of the model), which is constant in a typical regression. In contrast, the Poisson distribution has no such parameter, and in fact the variance increases with the mean (i.e. the variance and the mean have the same value). In this latter case, for an expected value of E(y) = 5, we also expect that the variance of observed data points is 5. But what if it is not? What if the observed variance is much higher, i.e. if the data are overdispersed? (Note that it could also be lower, underdispersed. This is less often the case, and not all approaches below allow for modelling underdispersion, but some do.) Overdispersion arises in different ways, most commonly through “clumping”. Imagine the number of seedlings in a forest plot. Depending on the distance to the source tree, there may be many (hundreds) or none. The same goes for shooting stars: either the sky is empty, or littered with shooting stars.
  • Poisson Versus Negative Binomial Regression

    Poisson Versus Negative Binomial Regression

    Handling Count Data The Negative Binomial Distribution Other Applications and Analysis in R References Poisson versus Negative Binomial Regression Randall Reese Utah State University [email protected] February 29, 2016 Randall Reese Poisson and Neg. Binom Handling Count Data The Negative Binomial Distribution Other Applications and Analysis in R References Overview 1 Handling Count Data ADEM Overdispersion 2 The Negative Binomial Distribution Foundations of Negative Binomial Distribution Basic Properties of the Negative Binomial Distribution Fitting the Negative Binomial Model 3 Other Applications and Analysis in R 4 References Randall Reese Poisson and Neg. Binom Handling Count Data The Negative Binomial Distribution ADEM Other Applications and Analysis in R Overdispersion References Count Data Randall Reese Poisson and Neg. Binom Handling Count Data The Negative Binomial Distribution ADEM Other Applications and Analysis in R Overdispersion References Count Data Data whose values come from Z≥0, the non-negative integers. Classic example is deaths in the Prussian army per year by horse kick (Bortkiewicz) Example 2 of Notes 5. (Number of successful \attempts"). Randall Reese Poisson and Neg. Binom Handling Count Data The Negative Binomial Distribution ADEM Other Applications and Analysis in R Overdispersion References Poisson Distribution Support is the non-negative integers. (Count data). Described by a single parameter λ > 0. When Y ∼ Poisson(λ), then E(Y ) = Var(Y ) = λ Randall Reese Poisson and Neg. Binom Handling Count Data The Negative Binomial Distribution ADEM Other Applications and Analysis in R Overdispersion References Acute Disseminated Encephalomyelitis Acute Disseminated Encephalomyelitis (ADEM) is a neurological, immune disorder in which widespread inflammation of the brain and spinal cord damages tissue known as white matter.
  • Section 1: Regression Review

    Section 1: Regression Review

    Section 1: Regression review Yotam Shem-Tov Fall 2015 Yotam Shem-Tov STAT 239A/ PS 236A September 3, 2015 1 / 58 Contact information Yotam Shem-Tov, PhD student in economics. E-mail: [email protected]. Office: Evans hall room 650 Office hours: To be announced - any preferences or constraints? Feel free to e-mail me. Yotam Shem-Tov STAT 239A/ PS 236A September 3, 2015 2 / 58 R resources - Prior knowledge in R is assumed In the link here is an excellent introduction to statistics in R . There is a free online course in coursera on programming in R. Here is a link to the course. An excellent book for implementing econometric models and method in R is Applied Econometrics with R. This is my favourite for starting to learn R for someone who has some background in statistic methods in the social science. The book Regression Models for Data Science in R has a nice online version. It covers the implementation of regression models using R . A great link to R resources and other cool programming and econometric tools is here. Yotam Shem-Tov STAT 239A/ PS 236A September 3, 2015 3 / 58 Outline There are two general approaches to regression 1 Regression as a model: a data generating process (DGP) 2 Regression as an algorithm (i.e., an estimator without assuming an underlying structural model). Overview of this slides: 1 The conditional expectation function - a motivation for using OLS. 2 OLS as the minimum mean squared error linear approximation (projections). 3 Regression as a structural model (i.e., assuming the conditional expectation is linear).
  • Generalized Linear Models

    Generalized Linear Models

    CHAPTER 6 Generalized linear models 6.1 Introduction Generalized linear modeling is a framework for statistical analysis that includes linear and logistic regression as special cases. Linear regression directly predicts continuous data y from a linear predictor Xβ = β0 + X1β1 + + Xkβk.Logistic regression predicts Pr(y =1)forbinarydatafromalinearpredictorwithaninverse-··· logit transformation. A generalized linear model involves: 1. A data vector y =(y1,...,yn) 2. Predictors X and coefficients β,formingalinearpredictorXβ 1 3. A link function g,yieldingavectoroftransformeddataˆy = g− (Xβ)thatare used to model the data 4. A data distribution, p(y yˆ) | 5. Possibly other parameters, such as variances, overdispersions, and cutpoints, involved in the predictors, link function, and data distribution. The options in a generalized linear model are the transformation g and the data distribution p. In linear regression,thetransformationistheidentity(thatis,g(u) u)and • the data distribution is normal, with standard deviation σ estimated from≡ data. 1 1 In logistic regression,thetransformationistheinverse-logit,g− (u)=logit− (u) • (see Figure 5.2a on page 80) and the data distribution is defined by the proba- bility for binary data: Pr(y =1)=y ˆ. This chapter discusses several other classes of generalized linear model, which we list here for convenience: The Poisson model (Section 6.2) is used for count data; that is, where each • data point yi can equal 0, 1, 2, ....Theusualtransformationg used here is the logarithmic, so that g(u)=exp(u)transformsacontinuouslinearpredictorXiβ to a positivey ˆi.ThedatadistributionisPoisson. It is usually a good idea to add a parameter to this model to capture overdis- persion,thatis,variationinthedatabeyondwhatwouldbepredictedfromthe Poisson distribution alone.
  • (Introduction to Probability at an Advanced Level) - All Lecture Notes

    (Introduction to Probability at an Advanced Level) - All Lecture Notes

    Fall 2018 Statistics 201A (Introduction to Probability at an advanced level) - All Lecture Notes Aditya Guntuboyina August 15, 2020 Contents 0.1 Sample spaces, Events, Probability.................................5 0.2 Conditional Probability and Independence.............................6 0.3 Random Variables..........................................7 1 Random Variables, Expectation and Variance8 1.1 Expectations of Random Variables.................................9 1.2 Variance................................................ 10 2 Independence of Random Variables 11 3 Common Distributions 11 3.1 Ber(p) Distribution......................................... 11 3.2 Bin(n; p) Distribution........................................ 11 3.3 Poisson Distribution......................................... 12 4 Covariance, Correlation and Regression 14 5 Correlation and Regression 16 6 Back to Common Distributions 16 6.1 Geometric Distribution........................................ 16 6.2 Negative Binomial Distribution................................... 17 7 Continuous Distributions 17 7.1 Normal or Gaussian Distribution.................................. 17 1 7.2 Uniform Distribution......................................... 18 7.3 The Exponential Density...................................... 18 7.4 The Gamma Density......................................... 18 8 Variable Transformations 19 9 Distribution Functions and the Quantile Transform 20 10 Joint Densities 22 11 Joint Densities under Transformations 23 11.1 Detour to Convolutions......................................
  • Modeling Overdispersion with the Normalized Tempered Stable Distribution

    Modeling Overdispersion with the Normalized Tempered Stable Distribution

    Modeling overdispersion with the normalized tempered stable distribution M. Kolossiatisa, J.E. Griffinb, M. F. J. Steel∗,c aDepartment of Hotel and Tourism Management, Cyprus University of Technology, P.O. Box 50329, 3603 Limassol, Cyprus. bInstitute of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. cDepartment of Statistics, University of Warwick, Coventry, CV4 7AL, U.K. Abstract A multivariate distribution which generalizes the Dirichlet distribution is intro- duced and its use for modeling overdispersion in count data is discussed. The distribution is constructed by normalizing a vector of independent tempered sta- ble random variables. General formulae for all moments and cross-moments of the distribution are derived and they are found to have similar forms to those for the Dirichlet distribution. The univariate version of the distribution can be used as a mixing distribution for the success probability of a binomial distribution to define an alternative to the well-studied beta-binomial distribution. Examples of fitting this model to simulated and real data are presented. Keywords: Distribution on the unit simplex; Mice fetal mortality; Mixed bino- mial; Normalized random measure; Overdispersion 1. Introduction In many experiments we observe data as the number of observations in a sam- ple with some property. The binomial distribution is a natural model for this type of data. However, the data are often found to be overdispersed relative to that ∗Corresponding author. Tel.: +44-2476-523369; fax: +44-2476-524532 Email addresses: [email protected] (M. Kolossiatis), [email protected] (J.E. Griffin), [email protected] (M.
  • Variance Partitioning in Multilevel Logistic Models That Exhibit Overdispersion

    Variance Partitioning in Multilevel Logistic Models That Exhibit Overdispersion

    J. R. Statist. Soc. A (2005) 168, Part 3, pp. 599–613 Variance partitioning in multilevel logistic models that exhibit overdispersion W. J. Browne, University of Nottingham, UK S. V. Subramanian, Harvard School of Public Health, Boston, USA K. Jones University of Bristol, UK and H. Goldstein Institute of Education, London, UK [Received July 2002. Final revision September 2004] Summary. A common application of multilevel models is to apportion the variance in the response according to the different levels of the data.Whereas partitioning variances is straight- forward in models with a continuous response variable with a normal error distribution at each level, the extension of this partitioning to models with binary responses or to proportions or counts is less obvious. We describe methodology due to Goldstein and co-workers for appor- tioning variance that is attributable to higher levels in multilevel binomial logistic models. This partitioning they referred to as the variance partition coefficient. We consider extending the vari- ance partition coefficient concept to data sets when the response is a proportion and where the binomial assumption may not be appropriate owing to overdispersion in the response variable. Using the literacy data from the 1991 Indian census we estimate simple and complex variance partition coefficients at multiple levels of geography in models with significant overdispersion and thereby establish the relative importance of different geographic levels that influence edu- cational disparities in India. Keywords: Contextual variation; Illiteracy; India; Multilevel modelling; Multiple spatial levels; Overdispersion; Variance partition coefficient 1. Introduction Multilevel regression models (Goldstein, 2003; Bryk and Raudenbush, 1992) are increasingly being applied in many areas of quantitative research.
  • Meta-Analysis of Count Data What Is Count Data Complicated by What Is

    Meta-Analysis of Count Data What Is Count Data Complicated by What Is

    What is count data Meta-analysis of count data • Data that can only take the non-negative integer values 0, 1, 2, 3, ……. • Examples in RCTs – Episodes of exacerbation of asthma – Falls Peter Herbison – Number of incontinence episodes Department of Preventive and • Time scale can vary from the whole study Social Medicine period to the last few (pick your time University of Otago period). Complicated by What is large? • Different exposure times .4 – Withdrawals .3 – Deaths – Different diary periods .2 • Sometimes people only fill out some days of the Probability diary as well as the periods being different .1 • Counts with a large mean can be treated 0 as normally distributed 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 n Mean 1 Mean 2 Mean 3 Mean 4 Mean 5 Analysis as rates Meta-analysis of rate ratios • Count number of events in each arm • Formula for calculation in Handbook • Calculate person years at risk in each arm • Straightforward using generic inverse • Calculate rates in each arm variance • Divide to get a rate ratio • In addition to the one reference in the • Can also be done by poisson regression Handbook: family – Guevara et al. Meta-analytic methods for – Poisson, negative binomial, zero inflated pooling rates when follow up duration varies: poisson a case study. BMC Medical Research – Assumes constant underlying risk Methodology 2004;4:17 Dichotimizing Treat as continuous • Change data so that it is those who have • May not be too bad one or more events versus those who • Handbook says beware of skewness have none • Get a 2x2 table • Good if the purpose of the treatment is to prevent events rather than reduce the number of them • Meta-analyse as usual Time to event Not so helpful • Usually time to first event but repeated • Many people look at the distribution and events can be used think that the only distribution that exists is • Analyse by survival analysis (e.g.
  • Generalized Linear Models and Point Count Data: Statistical Considerations for the Design and Analysis of Monitoring Studies

    Generalized Linear Models and Point Count Data: Statistical Considerations for the Design and Analysis of Monitoring Studies

    Generalized Linear Models and Point Count Data: Statistical Considerations for the Design and Analysis of Monitoring Studies Nathaniel E. Seavy,1,2,3 Suhel Quader,1,4 John D. Alexander,2 and C. John Ralph5 ________________________________________ Abstract The success of avian monitoring programs to effec- Key words: Generalized linear models, juniper remov- tively guide management decisions requires that stud- al, monitoring, overdispersion, point count, Poisson. ies be efficiently designed and data be properly ana- lyzed. A complicating factor is that point count surveys often generate data with non-normal distributional pro- perties. In this paper we review methods of dealing with deviations from normal assumptions, and we Introduction focus on the application of generalized linear models (GLMs). We also discuss problems associated with Measuring changes in bird abundance over time and in overdispersion (more variation than expected). In order response to habitat management is widely recognized to evaluate the statistical power of these models to as an important aspect of ecological monitoring detect differences in bird abundance, it is necessary for (Greenwood et al. 1993). The use of long-term avian biologists to identify the effect size they believe is monitoring programs (e.g., the Breeding Bird Survey) biologically significant in their system. We illustrate to identify population trends is a powerful tool for bird one solution to this challenge by discussing the design conservation (Sauer and Droege 1990, Sauer and Link of a monitoring program intended to detect changes in 2002). Short-term studies comparing bird abundance in bird abundance as a result of Western juniper (Juniper- treated and untreated areas are also important because us occidentalis) reduction projects in central Oregon.
  • Large Scale Conditional Covariance Matrix Modeling, Estimation and Testing

    Large Scale Conditional Covariance Matrix Modeling, Estimation and Testing

    Large Scale Conditional Covariance Matrix Modeling, Estimation and Testing Zhuanxin Ding1 Frank Russell Company, Tacoma, WA 98401, USA and Robert F. Engle University of California, San Diego, La Jolla, CA 92093, USA June 1996 This Revision May 2001 Abstract A new representation of the diagonal Vech model is given using the Hadamard product. Sufficient conditions on parameter matrices are provided to ensure the positive definiteness of covariance matrices from the new representation. Based on this, some new and simple models are discussed. A set of diagnostic tests for multivariate ARCH models is proposed. The tests are able to detect various model misspecifications by examing the orthogonality of the squared normalized residuals. A small Monte-Carlo study is carried out to check the small sample performance of the test. An empirical example is also given as guidance for model estimation and selection in the multivariate framework. For the specific data set considered, it is found that the simple one and two parameter models and the constant conditional correlation model perform fairly well. Keywords: conditional covariance, Multivariate ARCH, Hadamard product, M-test. 1 We would like to thank Clive Granger, Hong-Ye Gao, Doug Martin, Jeffrey Wang and Hal White for many helpful discussions. 2 1. Introduction It has now become common knowledge that the risk of the financial market, which is mainly represented by the variance of individual stock returns and the covariance between different assets or with the market, is highly forecastable. The research in this area however is far from complete and the application of different econometric models is just at its beginning.