Scientific Question Determine whether the breastfeeding of Nepalese children varies with child age and/or sex of child.
|
|
- Dwight Hamilton
- 6 years ago
- Views:
Transcription
1 Longtudnal Logstc Regresson: Breastfeedng of Nepalese Chldren PART II GEE models (margnal, populaton average) covered last lab Random Intercept models (subject specfc) Transton models Scentfc Queston Determne whether the breastfeedng of Nepalese chldren vares wth chld age and/or sex of chld. Data: Nepal Data (nepal.dta) as modfed n Lab 10 Outcome: Y j =I(breastfeedng j ) for ndvdual at vst number j We use vst number as our tme. We prepare our dataset lke we dd n Lab 10.. use "nepal.dta", clear ************************************** ** Dataset prep work from lab 10 **** ************************************** ** drop extra varables **. drop age age3 age4 t ** drop other varables we're not usng n ths analyss. drop wt ht arm day month year ded alve mage lt ** gen vst number varable to use as our tme varable **. sort d age. by d: gen vst=_n. tab vst vst Freq. Percent Cum Total 1, xtset d vst panel varable: d (strongly balanced) tme varable: vst, 1 to 5 delta: 1 unt Drop observatons wth mssng values on our outcome. drop f bf==. (53 observatons deleted). xtdes 1
2 d: 1,,..., 00 n = 199 vst: 1,,..., 5 T = 5 Delta(vst) = 1 unt Span(vst) = 5 perods (d*vst unquely dentfes each observaton) Dstrbuton of T_: mn 5% 5% 50% 75% 95% max Freq. Percent Cum. Pattern XXXXX We create a bnary ndcator of breastfeedng for our outcome varable. *** combne groups to gen 0-1 ndcator of breast feedng ***. gen bfbn=1 f bf==1 bf== (564 mssng values generated). replace bfbn=0 f bfbn==. (564 real changes made) For our covarates of nterest n the analyss, we create a centered age varable, an ndcator for male gender and an nteracton between male gender and centered age. ** create a centered age varable. gen agec = age * generate a bnary ndcator for male gender *. gen male=(sex==1). drop sex * defne centered age and male gender nteracton term *. gen agemale=agec*male Subject-specfc models accountng for correlaton (Random Intercept) We run a random ntercept model for breastfeedng status (y j ) where we account for correlaton of the repeated observatons on chldren by ncludng a random ntercept for each chld () and control for chld s age (centered). logtp U ( y = 1 U ) j ~ N( 0, σ ) u = β 0 = ( β 0 + β 1 agec +U ) + β j 1 +U agec j
3 (At least) two ways to ft a logstc random ntercept model n Stata 1. xtlogt. gllamm (wll be used a lot n the multlevel modelng class) xtlogt We ll ft the xtlogt model frst snce xtlogt s a more smple command Random-effects logstc regresson Number of obs = 947 Group varable (): d Number of groups = 199 Random effects u_ ~ Gaussan Obs per group: mn = 1 avg = 4.8 max = 5 Wald ch(1) = Log lkelhood = Prob > ch = bfbn Coef. Std. Err. z P> z [95% Conf. Interval] agec _cons /lnsgu sgma_u rho Lkelhood-rato test of rho=0: chbar(01) = Prob >= chbar = Report results on the OR scale: bfbn OR Std. Err. z P> z [95% Conf. Interval] agec /lnsgu sgma_u rho Coeffcent estmates have subject-specfc nterpretaton! Interpretaton of OR for agec: For a gven chld, the odds of breastfeedng decreases by 8% for a one month ncrease n that chld s age. If ths result was from a GEE (populaton average model), the nterpretaton would be: On average, for the Nepalese chldren n our study, the odds of breastfeedng decreases by 40% for each one month ncrease n age. 3
4 Intra-class correlaton n our model Two equvalent nterpretatons: 1. The proporton of the total varance n breastfeedng status that s due to dfferences between chldren (.e., the varance n the random ntercept for chld) after controllng for age.. The correlaton of the repeated measurements of breastfeedng status on the same chld after controllng for age. In our model, rho = 0.63 rh o σ u σ u + σ ε The varance of the random ntercept s (sgma_u)^ = 5.61 Dgresson on usng xtlogt check for adequate fttng of the model xtlogt uses quadrature methods to approxmate the lkelhood functon snce there s no closed form soluton. Ths means xtlogt and other methods that use quadrature can run slowly and we need to check that the specfed quadrature has adequately approxmated the lkelhood.. quadchk, nooutput Refttng model ntponts() = 8 Refttng model ntponts() = 16 Quadrature check Ftted Comparson Comparson quadrature quadrature quadrature 1 ponts 8 ponts 16 ponts Log lkelhood Dfference Relatve dfference bfbn: agec Dfference Relatve dfference bfbn: _cons Dfference Relatve dfference lnsgu: _cons Dfference Relatve dfference Relatve dfference > 10^- (1%) (arbtrary STATA rule of thumb) so reft usng > # of ntegraton ponts (default was 1). 4
5 . xtlogt bfbn agec, (d) or ntponts(18) nolog Random-effects logstc regresson Number of obs = 947 Group varable (): d Number of groups = 199 Random effects u_ ~ Gaussan Obs per group: mn = 1 avg = 4.8 max = 5 Wald ch(1) = Log lkelhood = Prob > ch = bfbn OR Std. Err. z P> z [95% Conf. Interval] agec /lnsgu sgma_u rho Lkelhood-rato test of rho=0: chbar(01) = Prob >= chbar = quadchk, nooutput Refttng model ntponts() = 14 Refttng model ntponts() = Quadrature check Ftted Comparson Comparson quadrature quadrature quadrature 18 ponts 14 ponts ponts Log lkelhood Dfference Relatve dfference bfbn: agec Dfference Relatve dfference bfbn: _cons Dfference Relatve dfference lnsgu: _cons Dfference Relatve dfference We can ncrease ntponts agan for even more assurance that the lkelhood s approprately approxmated.. xtlogt bfbn agec, (d) or ntponts(100) nolog Random-effects logstc regresson Number of obs = 947 Group varable (): d Number of groups = 199 Random effects u_ ~ Gaussan Obs per group: mn = 1 avg = 4.8 5
6 max = 5 Wald ch(1) = 57.0 Log lkelhood = Prob > ch = bfbn OR Std. Err. z P> z [95% Conf. Interval] agec /lnsgu sgma_u rho Lkelhood-rato test of rho=0: chbar(01) = Prob >= chbar = Lkelhood-rato test of rho=0: chbar(01) = Prob >= chbar = quadchk, nooutput Refttng model ntponts() = 67 Refttng model ntponts() = 133 Quadrature check Ftted Comparson Comparson quadrature quadrature quadrature 100 ponts 67 ponts 133 ponts Log lkelhood.63e-06.45e-07 Dfference -1.53e e-09 Relatve dfference bfbn: agec e e-07 Dfference 1.33e e-06 Relatve dfference bfbn: _cons e e-06 Dfference 1.847e e-06 Relatve dfference lnsgu: _cons.863e e-06 Dfference 9.159e e-07 Relatve dfference Note the dfference! Now the OR for agec s 0.6 and the estmated heterogenety s 0.9. GLLAMM generalzed lnear latent and mxed models (gllamm can do a lot more than just logstc models wth a random ntercept!) Install gllamm (you have to be connected to the nternet) ssc descrbe gllamm or update gllamm (replace prevous verson) ssc nstall gllamm, replace 6
7 Ft the same random ntercept model where we control for age. logtp U ( y = 1 U ) j ~ N( 0, σ ) u = β = ( β. gllamm bfbn agec, (d) l(logt) f(bnom) number of level 1 unts = 947 number of level unts = β 0 1 agec +U ) + β j 1 +U agec bfbn Coef. Std. Err. z P> z [95% Conf. Interval] agec _cons Varances and covarances of random effects ***level (d) var(1): ( ) ** get OR nstead of log(or) **. gllamm, eform bfbn exp(b) Std. Err. z P> z [95% Conf. Interval] agec Varances and covarances of random effects ***level (d) var(1): ( ) gllamm reports the varance of the random ntercept, σ u =15.38 We ll cover gllamm n much more detal 4 th term. Now, we are gong to compare the GEE models (lab 10) to the random ntercept model. Frst, we obtan the predcted log odds from the random ntercept model (gllamm). * get the ftted log odds. predct fttedlo, xb (xb wll be stored n fttedlo) j Next, we assgn values to the random ntercept for each chld usng Emprcal Bayes estmates (we ll cover ths more n 4 th term).. gllapred ebri, u (means and standard devatons wll be stored n ebrim1 ebris1) Produce a varable contanng the ftted probablty of breastfeedng for an average chld (a chld wth a value of the random ntercept = 0).. gen ftted_avgp = exp(fttedlo)/( 1 + exp(fttedlo) ) 7
8 Produce a varable contanng the ftted probablty of breastfeedng for each chld.. gen ftted_ndp = exp(fttedlo + ebrim1)/( 1 + exp(fttedlo + ebrim1) ) Now we reft GEE models wth three correlaton structures (from lab 10) and store the ftted values.. quetly xtgee bfbn agec, f(bn) lnk(logt) corr(ar1) robust. predct ftted_geear1, mu. label var ftted_geear1 "gee_ar1". ****Ft GEE wth unform corr model ****. quetly xtgee bfbn agec, f(bn) lnk(logt) corr(exc) robust. predct ftted_geeunf, mu. label var ftted_geeunf "gee_unf". ****Ft GEE wth ndependence corr model ****. quetly xtgee bfbn agec, f(bn) lnk(logt) corr(nd) robust. predct ftted_geend, mu. label var ftted_geend "gee_nd" Smooth the observed breastfeedng data versus agec. ksm bfbn agec, lowess bw(.4) ylab(0(.)1) lwdth(10) gen(bfbnsm) Plot the GEE and RI model fts. twoway (scatter bfbn agec, jtter(4)) (lne bfbnsm ftted_geear1 ftted_geeunf ftted_geend ftted_avgp agec, sort pstyle (p1)), ylab(0(.)1) agec bfbn gee_ar1 gee_nd lowess: bfbn gee_unf ftted_avgp 8
9 The ft from random ntercept model s qute dfferent from those from margnal model. Ths s to be expected because they are estmatng dfferent thngs and have dfferent nterpretatons. The OR from the random ntercept model s larger n absolute value (the logstc curve s steeper), as s guaranteed by theory. (Check the book!) twoway (lne ftted_ndp agec, c(l) pstyle(p15)) (lne ftted_avgp agec, sort clwdth(thck)) (lne ftted_geear1 agec, sort clwdth(thck)) agec ftted_ndp gee_ar1 ftted_avgp Transton model Basc dea: nclude the outcome varable at a prevous tme pont as a fxed effect covarate. We condton the response at tme j on the response at tme j-1 or j-, etc Generate a varable that represents breastfeedng status durng the prevous month. We often call ths lag 1 breastfeedng status.. sort d vst. by d: gen bfbn_lag1 = bfbn[_n-1] (199 mssng values generated) A very smple transton model 9
10 We ll assume condtonal ndependence once we condton on the outcome at the prevous month (frst order Markov chan) ( y = 1 U ) logt P j = β + β agecj + β y( j 1) 0 1. logt bfbn agec bfbn_lag1 Logstc regresson Number of obs = 748 LR ch() = Prob > ch = Log lkelhood = Pseudo R = bfbn Coef. Std. Err. z P> z [95% Conf. Interval] agec bfbn_lag _cons logstc bfbn agec bfbn_lag1 Logstc regresson Number of obs = 748 LR ch() = Prob > ch = Log lkelhood = Pseudo R = bfbn Odds Rato Std. Err. z P> z [95% Conf. Interval] agec bfbn_lag Interpretaton of OR on agec (populaton average) Comparng chldren of the same breastfeedng status at the prevous month, a one month ncrease n age s assocated wth a 10% decrease n the odds of breastfeedng. Interpretaton of OR on bfbn_lag1 (populaton average) Comparng chldren of the same age, those chldren who were breastfeedng at the prevous month have an odds rato of breastfeedng that s 101 tmes greater than the odds rato of breastfeedng for those chldren who were not breastfeedng at the prevous month. Huge!! Note that the 95% CI for the lag 1 bfbn varable does not nclude 1, so we should keep ths varable n our model. Less smple transton model We ll allow for a correlaton structure on the responses even after we condton on the outcome at the prevous month 10
11 . xtgee bfbn agec bfbn_lag1, nolog f(bn) l(logt) corr(unst) robust GEE populaton-averaged model Number of obs = 748 Group and tme vars: d vst Number of groups = 194 Lnk: logt Obs per group: mn = 1 Famly: bnomal avg = 3.9 Correlaton: unstructured max = 4 Wald ch() = Scale parameter: 1 Prob > ch = (Std. Err. adjusted for clusterng on d) Sem-robust bfbn Coef. Std. Err. z P> z [95% Conf. Interval] agec bfbn_lag _cons xtcorr Estmated wthn-d correlaton matrx R: c1 c c3 c4 r r r r The workng correlaton matrx estmates correlatons that are all relatvely close to zero, We ll keep the unstructured correlaton wth robust as our fnal model.. xtgee bfbn agec male agemale bfbn_lag1, nolog f(bn) l(logt) corr(unst) robust eform GEE populaton-averaged model Number of obs = 748 Group and tme vars: d vst Number of groups = 194 Lnk: logt Obs per group: mn = 1 Famly: bnomal avg = 3.9 Correlaton: unstructured max = 4 Wald ch(4) = Scale parameter: 1 Prob > ch = (Std. Err. adjusted for clusterng on d) Sem-robust bfbn Odds Rato Std. Err. z P> z [95% Conf. Interval] agec male agemale bfbn_lag xtcorr Estmated wthn-d correlaton matrx R: c1 c c3 c4 r
12 r r r Test for a gender effect. test male agemale ( 1) male = 0 ( ) agemale = 0 ch( ) = 1.58 Prob > ch = We stll do not see an effect of gender on the breastfeedng status of the Nepalese chldren. Fnal note on the transtons observed n breastfeedng status:. xttrans bfbn bfbn bfbn 0 1 Total Total Maybe our model stll doesn t adequately represent the data. There are a class of models called mover-stayer models developed n the 1950 s that attempt to represent data of ths sort (we don t cover these models n ths class). You can fnd tons of nformaton on these models wth a Google search. 1
BIO Lab 2: TWO-LEVEL NORMAL MODELS with school children popularity data
Lab : TWO-LEVEL NORMAL MODELS wth school chldren popularty data Purpose: Introduce basc two-level models for normally dstrbuted responses usng STATA. In partcular, we dscuss Random ntercept models wthout
More informationLab 4: Two-level Random Intercept Model
BIO 656 Lab4 009 Lab 4: Two-level Random Intercept Model Data: Peak expratory flow rate (pefr) measured twce, usng two dfferent nstruments, for 17 subjects. (from Chapter 1 of Multlevel and Longtudnal
More informationLecture 6: Introduction to Linear Regression
Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Recall: man dea of lnear regresson Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 8 Lnear regresson can be used to study an
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 008 Recall: man dea of lnear regresson Lnear regresson can be used to study
More informationsince [1-( 0+ 1x1i+ 2x2 i)] [ 0+ 1x1i+ assumed to be a reasonable approximation
Econ 388 R. Butler 204 revsons Lecture 4 Dummy Dependent Varables I. Lnear Probablty Model: the Regresson model wth a dummy varables as the dependent varable assumpton, mplcaton regular multple regresson
More informationLecture 3 Stat102, Spring 2007
Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture
More information1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands
Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of
More informationβ0 + β1xi. You are interested in estimating the unknown parameters β
Ordnary Least Squares (OLS): Smple Lnear Regresson (SLR) Analytcs The SLR Setup Sample Statstcs Ordnary Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals) wth OLS
More informationChapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where y + = β + β e for =,..., y and are observable varables e s a random error How can an estmaton rule be constructed for the
More informatione i is a random error
Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationSTATISTICS QUESTIONS. Step by Step Solutions.
STATISTICS QUESTIONS Step by Step Solutons www.mathcracker.com 9//016 Problem 1: A researcher s nterested n the effects of famly sze on delnquency for a group of offenders and examnes famles wth one to
More informationx i1 =1 for all i (the constant ).
Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by
More informationProperties of Least Squares
Week 3 3.1 Smple Lnear Regresson Model 3. Propertes of Least Squares Estmators Y Y β 1 + β X + u weekly famly expendtures X weekly famly ncome For a gven level of x, the expected level of food expendtures
More informationβ0 + β1xi. You are interested in estimating the unknown parameters β
Revsed: v3 Ordnar Least Squares (OLS): Smple Lnear Regresson (SLR) Analtcs The SLR Setup Sample Statstcs Ordnar Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals)
More informationDepartment of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6
Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationProblem 3.1: Error autotocorrelation and heteroskedasticity Standard variance components model:
ECON 510: Panel data econometrcs Semnar 3: October., 007 Problem 3.1: Error autotocorrelaton and heteroskedastcy Standard varance components model: (0.1) y = k+ x β + + u, ε = + u, IID(0, ), u Rewrng the
More informationReminder: Nested models. Lecture 9: Interactions, Quadratic terms and Splines. Effect Modification. Model 1
Lecture 9: Interactons, Quadratc terms and Splnes An Manchakul amancha@jhsph.edu 3 Aprl 7 Remnder: Nested models Parent model contans one set of varables Extended model adds one or more new varables to
More informationStatistics for Business and Economics
Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear
More information/ n ) are compared. The logic is: if the two
STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence
More informationIntroduction to Regression
Introducton to Regresson Dr Tom Ilvento Department of Food and Resource Economcs Overvew The last part of the course wll focus on Regresson Analyss Ths s one of the more powerful statstcal technques Provdes
More informationApplications of GEE Methodology Using the SAS System
Applcatons of GEE Methodology Usng the SAS System Gordon Johnston Maura Stokes SAS Insttute Inc, Cary, NC Abstract The analyss of correlated data arsng from repeated measurements when the measurements
More informationInterpreting Slope Coefficients in Multiple Linear Regression Models: An Example
CONOMICS 5* -- Introducton to NOT CON 5* -- Introducton to NOT : Multple Lnear Regresson Models Interpretng Slope Coeffcents n Multple Lnear Regresson Models: An xample Consder the followng smple lnear
More information9. Binary Dependent Variables
9. Bnar Dependent Varables 9. Homogeneous models Log, prob models Inference Tax preparers 9.2 Random effects models 9.3 Fxed effects models 9.4 Margnal models and GEE Appendx 9A - Lkelhood calculatons
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear
More informationPredictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore
Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.
More informationβ0 + β1xi and want to estimate the unknown
SLR Models Estmaton Those OLS Estmates Estmators (e ante) v. estmates (e post) The Smple Lnear Regresson (SLR) Condtons -4 An Asde: The Populaton Regresson Functon B and B are Lnear Estmators (condtonal
More informationSTAT 405 BIOSTATISTICS (Fall 2016) Handout 15 Introduction to Logistic Regression
STAT 45 BIOSTATISTICS (Fall 26) Handout 5 Introducton to Logstc Regresson Ths handout covers materal found n Secton 3.7 of your text. You may also want to revew regresson technques n Chapter. In ths handout,
More informationSee Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)
Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes
More informationBasically, if you have a dummy dependent variable you will be estimating a probability.
ECON 497: Lecture Notes 13 Page 1 of 1 Metropoltan State Unversty ECON 497: Research and Forecastng Lecture Notes 13 Dummy Dependent Varable Technques Studenmund Chapter 13 Bascally, f you have a dummy
More informationY = β 0 + β 1 X 1 + β 2 X β k X k + ε
Chapter 3 Secton 3.1 Model Assumptons: Multple Regresson Model Predcton Equaton Std. Devaton of Error Correlaton Matrx Smple Lnear Regresson: 1.) Lnearty.) Constant Varance 3.) Independent Errors 4.) Normalty
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationECON 351* -- Note 23: Tests for Coefficient Differences: Examples Introduction. Sample data: A random sample of 534 paid employees.
Model and Data ECON 35* -- NOTE 3 Tests for Coeffcent Dfferences: Examples. Introducton Sample data: A random sample of 534 pad employees. Varable defntons: W hourly wage rate of employee ; lnw the natural
More informationMultilevel Logistic Regression for Polytomous Data and Rankings
Outlne Multlevel Logstc Regresson for Polytomous Data and Rankngs 1. Introducton to Applcaton: Brtsh Electon Panel 2. Logstc Models as Random Utlty Models 3. Independence from Irrelevant Alternatves (IIA)
More informationThe Ordinary Least Squares (OLS) Estimator
The Ordnary Least Squares (OLS) Estmator 1 Regresson Analyss Regresson Analyss: a statstcal technque for nvestgatng and modelng the relatonshp between varables. Applcatons: Engneerng, the physcal and chemcal
More informationDO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR. Introductory Econometrics 1 hour 30 minutes
25/6 Canddates Only January Examnatons 26 Student Number: Desk Number:...... DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR Department Module Code Module Ttle Exam Duraton
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationUniversity of California at Berkeley Fall Introductory Applied Econometrics Final examination
SID: EEP 118 / IAS 118 Elsabeth Sadoulet and Daley Kutzman Unversty of Calforna at Berkeley Fall 01 Introductory Appled Econometrcs Fnal examnaton Scores add up to 10 ponts Your name: SID: 1. (15 ponts)
More informationHere is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)
Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,
More informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationBasic Business Statistics, 10/e
Chapter 13 13-1 Basc Busness Statstcs 11 th Edton Chapter 13 Smple Lnear Regresson Basc Busness Statstcs, 11e 009 Prentce-Hall, Inc. Chap 13-1 Learnng Objectves In ths chapter, you learn: How to use regresson
More informationSIMPLE LINEAR REGRESSION
Smple Lnear Regresson and Correlaton Introducton Prevousl, our attenton has been focused on one varable whch we desgnated b x. Frequentl, t s desrable to learn somethng about the relatonshp between two
More informationMaximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models
ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models
More informationJAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger
JAB Chan Long-tal clams development ASTIN - September 2005 B.Verder A. Klnger Outlne Chan Ladder : comments A frst soluton: Munch Chan Ladder JAB Chan Chan Ladder: Comments Black lne: average pad to ncurred
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationIntroduction to Dummy Variable Regressors. 1. An Example of Dummy Variable Regressors
ECONOMICS 5* -- Introducton to Dummy Varable Regressors ECON 5* -- Introducton to NOTE Introducton to Dummy Varable Regressors. An Example of Dummy Varable Regressors A model of North Amercan car prces
More informationMultinomial logit regression
07/0/6 Multnomal logt regresson Introducton We now turn our attenton to regresson models for the analyss of categorcal dependent varables wth more than two response categores: Y car owned (many possble
More informationDepartment of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution
Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationSTAT 3008 Applied Regression Analysis
STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationHow its computed. y outcome data λ parameters hyperparameters. where P denotes the Laplace approximation. k i k k. Andrew B Lawson 2013
Andrew Lawson MUSC INLA INLA s a relatvely new tool that can be used to approxmate posteror dstrbutons n Bayesan models INLA stands for ntegrated Nested Laplace Approxmaton The approxmaton has been known
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More informationBasic R Programming: Exercises
Basc R Programmng: Exercses RProgrammng John Fox ICPSR, Summer 2009 1. Logstc Regresson: Iterated weghted least squares (IWLS) s a standard method of fttng generalzed lnear models to data. As descrbed
More informationJanuary Examinations 2015
24/5 Canddates Only January Examnatons 25 DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR STUDENT CANDIDATE NO.. Department Module Code Module Ttle Exam Duraton (n words)
More informationAddressing Alternative. Multiple Regression Spring 2012
Addressng Alternatve Explanatons: Multple Regresson 7.87 Sprng 0 Dd Clnton hurt Gore example Dd Clnton hurt Gore n the 000 electon? Treatment s not lkng Bll Clnton How would you test ths? Bvarate regresson
More informationAddressing Alternative Explanations: Multiple Regression
Addressng Alternatve Explanatons: Multple Regresson 7.87 Dd Clnton hurt Gore example Dd Clnton hurt Gore n the 000 electon? Treatment s not lkng Bll Clnton How would you test ths? Bvarate regresson of
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationChapter 8 Indicator Variables
Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n
More informationStatistical analysis using matlab. HY 439 Presented by: George Fortetsanakis
Statstcal analyss usng matlab HY 439 Presented by: George Fortetsanaks Roadmap Probablty dstrbutons Statstcal estmaton Fttng data to probablty dstrbutons Contnuous dstrbutons Contnuous random varable X
More informationTHE APPLICATION OF LINEAR MIXED-EFFECTS MODEL TO THE EFFECT OF MICROCURRENT ON DECUBITUS WOUNDS. A STUDY IN LIMBURG PROVINCE OF BELGIUM
The 3 rd Internatonal Conference on Mathematcs and Statstcs (ICoMS-3) Insttut Pertanan Bogor, Indonesa, 5-6 August 008 THE APPLICATIO OF LIEAR MIXED-EFFECTS MODEL TO THE EFFECT OF MICROCURRET O DECUBITUS
More informationChapter 5 Multilevel Models
Chapter 5 Multlevel Models 5.1 Cross-sectonal multlevel models 5.1.1 Two-level models 5.1.2 Multple level models 5.1.3 Multple level modelng n other felds 5.2 Longtudnal multlevel models 5.2.1 Two-level
More informationMarginal Effects in Probit Models: Interpretation and Testing. 1. Interpreting Probit Coefficients
ECON 5 -- NOE 15 Margnal Effects n Probt Models: Interpretaton and estng hs note ntroduces you to the two types of margnal effects n probt models: margnal ndex effects, and margnal probablty effects. It
More informationModule Contact: Dr Susan Long, ECO Copyright of the University of East Anglia Version 1
UNIVERSITY OF EAST ANGLIA School of Economcs Man Seres PG Examnaton 016-17 ECONOMETRIC METHODS ECO-7000A Tme allowed: hours Answer ALL FOUR Questons. Queston 1 carres a weght of 5%; Queston carres 0%;
More informationMaximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models
ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Mamum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models for
More informationRegression Analysis. Regression Analysis
Regresson Analyss Smple Regresson Multvarate Regresson Stepwse Regresson Replcaton and Predcton Error 1 Regresson Analyss In general, we "ft" a model by mnmzng a metrc that represents the error. n mn (y
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationLearning Objectives for Chapter 11
Chapter : Lnear Regresson and Correlaton Methods Hldebrand, Ott and Gray Basc Statstcal Ideas for Managers Second Edton Learnng Objectves for Chapter Usng the scatterplot n regresson analyss Usng the method
More informationStatistics MINITAB - Lab 2
Statstcs 20080 MINITAB - Lab 2 1. Smple Lnear Regresson In smple lnear regresson we attempt to model a lnear relatonshp between two varables wth a straght lne and make statstcal nferences concernng that
More informationOutline. Zero Conditional mean. I. Motivation. 3. Multiple Regression Analysis: Estimation. Read Wooldridge (2013), Chapter 3.
Outlne 3. Multple Regresson Analyss: Estmaton I. Motvaton II. Mechancs and Interpretaton of OLS Read Wooldrdge (013), Chapter 3. III. Expected Values of the OLS IV. Varances of the OLS V. The Gauss Markov
More informationDiagnostics in Poisson Regression. Models - Residual Analysis
Dagnostcs n Posson Regresson Models - Resdual Analyss 1 Outlne Dagnostcs n Posson Regresson Models - Resdual Analyss Example 3: Recall of Stressful Events contnued 2 Resdual Analyss Resduals represent
More informationDummy variables in multiple variable regression model
WESS Econometrcs (Handout ) Dummy varables n multple varable regresson model. Addtve dummy varables In the prevous handout we consdered the followng regresson model: y x 2x2 k xk,, 2,, n and we nterpreted
More information[The following data appear in Wooldridge Q2.3.] The table below contains the ACT score and college GPA for eight college students.
PPOL 59-3 Problem Set Exercses n Smple Regresson Due n class /8/7 In ths problem set, you are asked to compute varous statstcs by hand to gve you a better sense of the mechancs of the Pearson correlaton
More informationLecture 3: Probability Distributions
Lecture 3: Probablty Dstrbutons Random Varables Let us begn by defnng a sample space as a set of outcomes from an experment. We denote ths by S. A random varable s a functon whch maps outcomes nto the
More informationChapter 14: Logit and Probit Models for Categorical Response Variables
Chapter 4: Logt and Probt Models for Categorcal Response Varables Sect 4. Models for Dchotomous Data We wll dscuss only ths secton of Chap 4, whch s manly about Logstc Regresson, a specal case of the famly
More informationBiostatistics. Chapter 11 Simple Linear Correlation and Regression. Jing Li
Bostatstcs Chapter 11 Smple Lnear Correlaton and Regresson Jng L jng.l@sjtu.edu.cn http://cbb.sjtu.edu.cn/~jngl/courses/2018fall/b372/ Dept of Bonformatcs & Bostatstcs, SJTU Recall eat chocolate Cell 175,
More informationTesting for seasonal unit roots in heterogeneous panels
Testng for seasonal unt roots n heterogeneous panels Jesus Otero * Facultad de Economía Unversdad del Rosaro, Colomba Jeremy Smth Department of Economcs Unversty of arwck Monca Gulett Aston Busness School
More informationLinear regression. Regression Models. Chapter 11 Student Lecture Notes Regression Analysis is the
Chapter 11 Student Lecture Notes 11-1 Lnear regresson Wenl lu Dept. Health statstcs School of publc health Tanjn medcal unversty 1 Regresson Models 1. Answer What Is the Relatonshp Between the Varables?.
More informationSTAT 511 FINAL EXAM NAME Spring 2001
STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte
More informationThis column is a continuation of our previous column
Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard
More informationLaboratory 1c: Method of Least Squares
Lab 1c, Least Squares Laboratory 1c: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly
More informationLaboratory 3: Method of Least Squares
Laboratory 3: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly they are correlated wth
More information1 Binary Response Models
Bnary and Ordered Multnomal Response Models Dscrete qualtatve response models deal wth dscrete dependent varables. bnary: yes/no, partcpaton/non-partcpaton lnear probablty model LPM, probt or logt models
More informationAn R implementation of bootstrap procedures for mixed models
The R User Conference 2009 July 8-10, Agrocampus-Ouest, Rennes, France An R mplementaton of bootstrap procedures for mxed models José A. Sánchez-Espgares Unverstat Poltècnca de Catalunya Jord Ocaña Unverstat
More informationStatistical Evaluation of WATFLOOD
tatstcal Evaluaton of WATFLD By: Angela MacLean, Dept. of Cvl & Envronmental Engneerng, Unversty of Waterloo, n. ctober, 005 The statstcs program assocated wth WATFLD uses spl.csv fle that s produced wth
More information28. SIMPLE LINEAR REGRESSION III
8. SIMPLE LINEAR REGRESSION III Ftted Values and Resduals US Domestc Beers: Calores vs. % Alcohol To each observed x, there corresponds a y-value on the ftted lne, y ˆ = βˆ + βˆ x. The are called ftted
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationInterval Estimation in the Classical Normal Linear Regression Model. 1. Introduction
ECONOMICS 35* -- NOTE 7 ECON 35* -- NOTE 7 Interval Estmaton n the Classcal Normal Lnear Regresson Model Ths note outlnes the basc elements of nterval estmaton n the Classcal Normal Lnear Regresson Model
More information3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. CDS Mphil Econometrics Vijayamohan. 3-Mar-14. CDS M Phil Econometrics.
Dummy varable Models an Plla N Dummy X-varables Dummy Y-varables Dummy X-varables Dummy X-varables Dummy varable: varable assumng values 0 and to ndcate some attrbutes To classfy data nto mutually exclusve
More informationU-Pb Geochronology Practical: Background
U-Pb Geochronology Practcal: Background Basc Concepts: accuracy: measure of the dfference between an expermental measurement and the true value precson: measure of the reproducblty of the expermental result
More information17 - LINEAR REGRESSION II
Topc 7 Lnear Regresson II 7- Topc 7 - LINEAR REGRESSION II Testng and Estmaton Inferences about β Recall that we estmate Yˆ ˆ β + ˆ βx. 0 μ Y X x β0 + βx usng To estmate σ σ squared error Y X x ε s ε we
More informationResource Allocation and Decision Analysis (ECON 8010) Spring 2014 Foundations of Regression Analysis
Resource Allocaton and Decson Analss (ECON 800) Sprng 04 Foundatons of Regresson Analss Readng: Regresson Analss (ECON 800 Coursepak, Page 3) Defntons and Concepts: Regresson Analss statstcal technques
More informationLogistic Regression. CAP 5610: Machine Learning Instructor: Guo-Jun QI
Logstc Regresson CAP 561: achne Learnng Instructor: Guo-Jun QI Bayes Classfer: A Generatve model odel the posteror dstrbuton P(Y X) Estmate class-condtonal dstrbuton P(X Y) for each Y Estmate pror dstrbuton
More informationDr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur
Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models
More information