AIC = Log likelihood = BIC =
|
|
- Lucinda Woods
- 5 years ago
- Views:
Transcription
1 - log: /mnt/ide1/home/sschulh1/apc/apc_examplelog log type: text opened on: 21 Jul 2006, 18:08:20 *replicate table 5 and cols 7-9 of table 3 in Yang, Fu and Land (2004) *Stata can maximize GLM objective functions using two different numerical *optimization methods: Newton-Raphson (NR) and iterative reweighted least *squares (IRLS) NR is the default in Stata However, NR presents a *problem in replicating Yang, Fu and Land: The paper set the scale *parameter equal to the deviance divided by the residual degrees of *freedom, but Stata allows this scale parameter only with IRLS and not *with NR So, we show two sets of results below: * 1 NR optimization with scale parameter=pearson chi-squared/residual df * 2 IRLS optimization with scale parameter=deviance/residual df *Version 1 is basically the default in Stata Version 2 matches what was *done in the paper The results are numerically identical to the number of *decimal places shown in the paper use apc_example_datadta *first, using Newton-Raphson optimization (the default in Stata) and scale(x2 > ) * (scale parameter = Pearson chi-squared / residual degrees of freedom) #delim ; delimiter now ; apc_ie death_f if age<=90, > age(age) period(year) cohort(cohort) family(poisson) link(log) > exposure(exp_f) scale(x2); Iteration 0: log likelihood = Iteration 1: log likelihood = Iteration 2: log likelihood = Iteration 3: log likelihood = Iteration 4: log likelihood = Iteration 5: log likelihood = Intrinsic estimator of APC effects No of obs = 152 Optimization : ML Residual df = 102 Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = AIC = Log likelihood = BIC = OIM age_
2 age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ period_ period_ period_ period_ period_ period_ period_ period_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ _cons apc_cglim death_f if age<=90,
3 > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(x2) constraint("a5=a10"); Iteration 0: log likelihood = Iteration 1: log likelihood = Iteration 2: log likelihood = Iteration 3: log likelihood = Iteration 4: log likelihood = Iteration 5: log likelihood = Generalized linear models No of obs = 152 Optimization : ML Residual df = 102 Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = AIC = Log likelihood = BIC = OIM _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _P_ _P_ _P_ _P_ _P_ _P_ _P_ _C_ _C_ _C_ _C_
4 _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of Pearson X2-based dispersion) _A_5=_A_10 drop _A* _P* _C*; apc_cglim death_f if age<=90, > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(x2) constraint("p1965=p1960"); Iteration 0: log likelihood = Iteration 1: log likelihood = Iteration 2: log likelihood = Iteration 3: log likelihood = Iteration 4: log likelihood = Iteration 5: log likelihood = Generalized linear models No of obs = 152 Optimization : ML Residual df = 102 Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = AIC = Log likelihood = BIC = OIM
5 _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _P_ _P_ _P_ _P_ _P_ _P_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of Pearson X2-based dispersion) _P_1965=_P_1960
6 drop _A* _P* _C*; apc_cglim death_f if age<=90, > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(x2) constraint("c1995=c1990"); Iteration 0: log likelihood = Iteration 1: log likelihood = Iteration 2: log likelihood = Iteration 3: log likelihood = Iteration 4: log likelihood = Iteration 5: log likelihood = Generalized linear models No of obs = 152 Optimization : ML Residual df = 102 Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = AIC = Log likelihood = BIC = OIM _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _P_ _P_ _P_ _P_ _P_ _P_ _P_
7 _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of Pearson X2-based dispersion) _C_1995=_C_1990 drop _A* _P* _C*; #delim cr delimiter now cr *next, using IRLS optimization (the default in S-Plus) and scale(dev) * (scale parameter = deviance / residual degrees of freedom) #delim ; delimiter now ; apc_ie death_f if age<=90, > age(age) period(year) cohort(cohort) family(poisson) link(log) > exposure(exp_f) scale(dev) irls; Iteration 1: deviance = Iteration 2: deviance = Iteration 3: deviance = Iteration 4: deviance = Iteration 5: deviance = Iteration 6: deviance = Iteration 7: deviance = Intrinsic estimator of APC effects No of obs = 152 Optimization : MQL Fisher scoring Residual df = 102 (IRLS EIM) Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson =
8 AIC = Deviance = BIC = EIM age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ age_ period_ period_ period_ period_ period_ period_ period_ period_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ cohort_
9 cohort_ cohort_ cohort_ cohort_ cohort_ cohort_ _cons apc_cglim death_f if age<=90, > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(dev) irls constraint("a5=a10"); Iteration 1: deviance = Iteration 2: deviance = Iteration 3: deviance = Iteration 4: deviance = Iteration 5: deviance = Iteration 6: deviance = Iteration 7: deviance = Generalized linear models No of obs = 152 Optimization : MQL Fisher scoring Residual df = 102 (IRLS EIM) Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = BIC = EIM _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_
10 _P_ _P_ _P_ _P_ _P_ _P_ _P_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of deviance-based dispersion) _A_5=_A_10 drop _A* _P* _C*; apc_cglim death_f if age<=90, > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(dev) irls constraint("p1965=p1960"); Iteration 1: deviance = Iteration 2: deviance = Iteration 3: deviance = Iteration 4: deviance = Iteration 5: deviance = Iteration 6: deviance = Iteration 7: deviance = Generalized linear models No of obs = 152 Optimization : MQL Fisher scoring Residual df = 102
11 (IRLS EIM) Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = BIC = EIM _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _P_ _P_ _P_ _P_ _P_ _P_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_
12 _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of deviance-based dispersion) _P_1965=_P_1960 drop _A* _P* _C*; apc_cglim death_f if age<=90, > age(age) period(year) cohort(cohort) > agepfx("_a") periodpfx("_p") cohortpfx("_c") > family(poisson) link(log) > exposure(exp_f) scale(dev) irls constraint("c1995=c1990"); Iteration 1: deviance = Iteration 2: deviance = Iteration 3: deviance = Iteration 4: deviance = Iteration 5: deviance = Iteration 6: deviance = Iteration 7: deviance = Generalized linear models No of obs = 152 Optimization : MQL Fisher scoring Residual df = 102 (IRLS EIM) Scale parameter = 1 Deviance = (1/df) Deviance = Pearson = (1/df) Pearson = BIC = EIM _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_ _A_
13 _A_ _A_ _A_ _A_ _P_ _P_ _P_ _P_ _P_ _P_ _P_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _C_ _cons (Standard errors scaled using square root of deviance-based dispersion) _C_1995=_C_1990 drop _A* _P* _C*; #delim cr delimiter now cr log close log: /mnt/ide1/home/sschulh1/apc/apc_examplelog log type: text closed on: 21 Jul 2006, 18:08:28 -
Maximum Likelihood Estimation
Maximum Likelihood Estimation EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #6 EPSY 905: Maximum Likelihood In This Lecture The basics of maximum likelihood estimation Ø The engine that
More information############################ ### toxo.r ### ############################
############################ ### toxo.r ### ############################ toxo < read.table(file="n:\\courses\\stat8620\\fall 08\\toxo.dat",header=T) #toxo < read.table(file="c:\\documents and Settings\\dhall\\My
More informationCategorical Outcomes. Statistical Modelling in Stata: Categorical Outcomes. R by C Table: Example. Nominal Outcomes. Mark Lunt.
Categorical Outcomes Statistical Modelling in Stata: Categorical Outcomes Mark Lunt Arthritis Research UK Epidemiology Unit University of Manchester Nominal Ordinal 28/11/2017 R by C Table: Example Categorical,
More informationIntro to GLM Day 2: GLM and Maximum Likelihood
Intro to GLM Day 2: GLM and Maximum Likelihood Federico Vegetti Central European University ECPR Summer School in Methods and Techniques 1 / 32 Generalized Linear Modeling 3 steps of GLM 1. Specify the
More informationMultiple Regression and Logistic Regression II. Dajiang 525 Apr
Multiple Regression and Logistic Regression II Dajiang Liu @PHS 525 Apr-19-2016 Materials from Last Time Multiple regression model: Include multiple predictors in the model = + + + + How to interpret the
More informationStatistics 175 Applied Statistics Generalized Linear Models Jianqing Fan
Statistics 175 Applied Statistics Generalized Linear Models Jianqing Fan Example 1 (Kyhposis data): (The data set kyphosis consists of measurements on 81 children following corrective spinal surgery. Variable
More informationNegative Binomial Model for Count Data Log-linear Models for Contingency Tables - Introduction
Negative Binomial Model for Count Data Log-linear Models for Contingency Tables - Introduction Statistics 149 Spring 2006 Copyright 2006 by Mark E. Irwin Negative Binomial Family Example: Absenteeism from
More informationOrdinal Multinomial Logistic Regression. Thom M. Suhy Southern Methodist University May14th, 2013
Ordinal Multinomial Logistic Thom M. Suhy Southern Methodist University May14th, 2013 GLM Generalized Linear Model (GLM) Framework for statistical analysis (Gelman and Hill, 2007, p. 135) Linear Continuous
More informationis the bandwidth and controls the level of smoothing of the estimator, n is the sample size and
Paper PH100 Relationship between Total charges and Reimbursements in Outpatient Visits Using SAS GLIMMIX Chakib Battioui, University of Louisville, Louisville, KY ABSTRACT The purpose of this paper is
More informationproc genmod; model malform/total = alcohol / dist=bin link=identity obstats; title 'Table 2.7'; title2 'Identity Link';
BIOS 6244 Analysis of Categorical Data Assignment 5 s 1. Consider Exercise 4.4, p. 98. (i) Write the SAS code, including the DATA step, to fit the linear probability model and the logit model to the data
More informationChapter 8 Exercises 1. Data Analysis & Graphics Using R Solutions to Exercises (May 1, 2010)
Chapter 8 Exercises 1 Data Analysis & Graphics Using R Solutions to Exercises (May 1, 2010) Preliminaries > library(daag) Exercise 1 The following table shows numbers of occasions when inhibition (i.e.,
More informationProjects for Bayesian Computation with R
Projects for Bayesian Computation with R Laura Vana & Kurt Hornik Winter Semeter 2018/2019 1 S&P Rating Data On the homepage of this course you can find a time series for Standard & Poors default data
More informationModel fit assessment via marginal model plots
The Stata Journal (2010) 10, Number 2, pp. 215 225 Model fit assessment via marginal model plots Charles Lindsey Texas A & M University Department of Statistics College Station, TX lindseyc@stat.tamu.edu
More informationEvaluation of a New Variance Components Estimation Method Modi ed Henderson s Method 3 With the Application of Two Way Mixed Model
Evaluation of a New Variance Components Estimation Method Modi ed Henderson s Method 3 With the Application of Two Way Mixed Model Author: Weigang Qie; Chenfan Xu Supervisor: Lars Rönnegård June 0th, 009
More informationLogistic Regression with R: Example One
Logistic Regression with R: Example One math = read.table("http://www.utstat.toronto.edu/~brunner/appliedf12/data/mathcat.data") math[1:5,] hsgpa hsengl hscalc course passed outcome 1 78.0 80 Yes Mainstrm
More information[BINARY DEPENDENT VARIABLE ESTIMATION WITH STATA]
Tutorial #3 This example uses data in the file 16.09.2011.dta under Tutorial folder. It contains 753 observations from a sample PSID data on the labor force status of married women in the U.S in 1975.
More informationLogistic Regression. Logistic Regression Theory
Logistic Regression Dr. J. Kyle Roberts Southern Methodist University Simmons School of Education and Human Development Department of Teaching and Learning Logistic Regression The linear probability model.
More informationGeneralized Linear Models
Generalized Linear Models Scott Creel Wednesday, September 10, 2014 This exercise extends the prior material on using the lm() function to fit an OLS regression and test hypotheses about effects on a parameter.
More information> budworm$samplogit < log((budworm$y+0.5)/(budworm$m budworm$y+0.5))
budworm < read.table(file="n:\\courses\\stat8620\\fall 08\\budworm.dat",header=T) #budworm < read.table(file="c:\\documents and Settings\\dhall\\My Documents\\Dan's Work Stuff\\courses\\STAT8620\\Fall
More informationFinal Exam - section 1. Thursday, December hours, 30 minutes
Econometrics, ECON312 San Francisco State University Michael Bar Fall 2013 Final Exam - section 1 Thursday, December 19 1 hours, 30 minutes Name: Instructions 1. This is closed book, closed notes exam.
More informationAddiction - Multinomial Model
Addiction - Multinomial Model February 8, 2012 First the addiction data are loaded and attached. > library(catdata) > data(addiction) > attach(addiction) For the multinomial logit model the function multinom
More informationCreation of Synthetic Discrete Response Regression Models
Arizona State University From the SelectedWorks of Joseph M Hilbe 2010 Creation of Synthetic Discrete Response Regression Models Joseph Hilbe, Arizona State University Available at: https://works.bepress.com/joseph_hilbe/2/
More informationtm / / / / / / / / / / / / Statistics/Data Analysis User: Klick Project: Limited Dependent Variables{space -6}
PS 4 Monday August 16 01:00:42 2010 Page 1 tm / / / / / / / / / / / / Statistics/Data Analysis User: Klick Project: Limited Dependent Variables{space -6} log: C:\web\PS4log.smcl log type: smcl opened on:
More informationAnd The Winner Is? How to Pick a Better Model
And The Winner Is? How to Pick a Better Model Part 2 Goodness-of-Fit and Internal Stability Dan Tevet, FCAS, MAAA Goodness-of-Fit Trying to answer question: How well does our model fit the data? Can be
More informationDuration Models: Parametric Models
Duration Models: Parametric Models Brad 1 1 Department of Political Science University of California, Davis January 28, 2011 Parametric Models Some Motivation for Parametrics Consider the hazard rate:
More informationTransport Data Analysis and Modeling Methodologies
Transport Data Analysis and Modeling Methodologies Lab Session #14 (Discrete Data Latent Class Logit Analysis based on Example 13.1) In Example 13.1, you were given 151 observations of a travel survey
More informationStat 328, Summer 2005
Stat 328, Summer 2005 Exam #2, 6/18/05 Name (print) UnivID I have neither given nor received any unauthorized aid in completing this exam. Signed Answer each question completely showing your work where
More informationActuarial Research on the Effectiveness of Collision Avoidance Systems FCW & LDW. A translation from Hebrew to English of a research paper prepared by
Actuarial Research on the Effectiveness of Collision Avoidance Systems FCW & LDW A translation from Hebrew to English of a research paper prepared by Ron Actuarial Intelligence LTD Contact Details: Shachar
More informationContents Utility theory and insurance The individual risk model Collective risk models
Contents There are 10 11 stars in the galaxy. That used to be a huge number. But it s only a hundred billion. It s less than the national deficit! We used to call them astronomical numbers. Now we should
More informationboxcox() returns the values of α and their loglikelihoods,
Solutions to Selected Computer Lab Problems and Exercises in Chapter 11 of Statistics and Data Analysis for Financial Engineering, 2nd ed. by David Ruppert and David S. Matteson c 2016 David Ruppert and
More informationMaximum Likelihood Estimation Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised January 13, 2018
Maximum Likelihood Estimation Richard Williams, University of otre Dame, https://www3.nd.edu/~rwilliam/ Last revised January 3, 208 [This handout draws very heavily from Regression Models for Categorical
More informationYour Name (Please print) Did you agree to take the optional portion of the final exam Yes No. Directions
Your Name (Please print) Did you agree to take the optional portion of the final exam Yes No (Your online answer will be used to verify your response.) Directions There are two parts to the final exam.
More informationSTK Lecture 7 finalizing clam size modelling and starting on pricing
STK 4540 Lecture 7 finalizing clam size modelling and starting on pricing Overview Important issues Models treated Curriculum Duration (in lectures) What is driving the result of a nonlife insurance company?
More informationA Comparison of Univariate Probit and Logit. Models Using Simulation
Applied Mathematical Sciences, Vol. 12, 2018, no. 4, 185-204 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.818 A Comparison of Univariate Probit and Logit Models Using Simulation Abeer
More informationStat 401XV Exam 3 Spring 2017
Stat 40XV Exam Spring 07 I have neither given nor received unauthorized assistance on this exam. Name Signed Date Name Printed ATTENTION! Incorrect numerical answers unaccompanied by supporting reasoning
More informationCase Study: Applying Generalized Linear Models
Case Study: Applying Generalized Linear Models Dr. Kempthorne May 12, 2016 Contents 1 Generalized Linear Models of Semi-Quantal Biological Assay Data 2 1.1 Coal miners Pneumoconiosis Data.................
More informationSean Howard Econometrics Final Project Paper. An Analysis of the Determinants and Factors of Physical Education Attendance in the Fourth Quarter
Sean Howard Econometrics Final Project Paper An Analysis of the Determinants and Factors of Physical Education Attendance in the Fourth Quarter Introduction This project attempted to gain a more complete
More informationsociology SO5032 Quantitative Research Methods Brendan Halpin, Sociology, University of Limerick Spring 2018 SO5032 Quantitative Research Methods
1 SO5032 Quantitative Research Methods Brendan Halpin, Sociology, University of Limerick Spring 2018 Lecture 10: Multinomial regression baseline category extension of binary What if we have multiple possible
More informationLog-linear Modeling Under Generalized Inverse Sampling Scheme
Log-linear Modeling Under Generalized Inverse Sampling Scheme Soumi Lahiri (1) and Sunil Dhar (2) (1) Department of Mathematical Sciences New Jersey Institute of Technology University Heights, Newark,
More informationIntroduction to fractional outcome regression models using the fracreg and betareg commands
Introduction to fractional outcome regression models using the fracreg and betareg commands Miguel Dorta Staff Statistician StataCorp LP Aguascalientes, Mexico (StataCorp LP) fracreg - betareg May 18,
More informationAnalytics on pension valuations
Analytics on pension valuations Research Paper Business Analytics Author: Arno Hendriksen November 4, 2017 Abstract EY Actuaries performs pension calculations for several companies where both the the assets
More informationCreating synthetic discrete-response regression models
The Stata Journal (2010) 10, Number 1, pp. 104 124 Creating synthetic discrete-response regression models Joseph M. Hilbe Arizona State University and Jet Propulsion Laboratory, CalTech Hilbe@asu.edu Abstract.
More informationMaximum Likelihood Estimation Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised January 10, 2017
Maximum Likelihood Estimation Richard Williams, University of otre Dame, https://www3.nd.edu/~rwilliam/ Last revised January 0, 207 [This handout draws very heavily from Regression Models for Categorical
More informationProblem Set 6 ANSWERS
Economics 20 Part I. Problem Set 6 ANSWERS Prof. Patricia M. Anderson The first 5 questions are based on the following information: Suppose a researcher is interested in the effect of class attendance
More informationEcon 371 Problem Set #4 Answer Sheet. 6.2 This question asks you to use the results from column (1) in the table on page 213.
Econ 371 Problem Set #4 Answer Sheet 6.2 This question asks you to use the results from column (1) in the table on page 213. a. The first part of this question asks whether workers with college degrees
More informationBrief Sketch of Solutions: Tutorial 2. 2) graphs. 3) unit root tests
Brief Sketch of Solutions: Tutorial 2 2) graphs LJAPAN DJAPAN 5.2.12 5.0.08 4.8.04 4.6.00 4.4 -.04 4.2 -.08 4.0 01 02 03 04 05 06 07 08 09 -.12 01 02 03 04 05 06 07 08 09 LUSA DUSA 7.4.12 7.3 7.2.08 7.1.04
More informationIntroduction to Population Modeling
Introduction to Population Modeling In addition to estimating the size of a population, it is often beneficial to estimate how the population size changes over time. Ecologists often uses models to create
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Generalized Linear Models - IIIb Henrik Madsen March 18, 2012 Henrik Madsen () Chapman & Hall March 18, 2012 1 / 32 Examples Overdispersion and Offset!
More informationORDERED MULTINOMIAL LOGISTIC REGRESSION ANALYSIS. Pooja Shivraj Southern Methodist University
ORDERED MULTINOMIAL LOGISTIC REGRESSION ANALYSIS Pooja Shivraj Southern Methodist University KINDS OF REGRESSION ANALYSES Linear Regression Logistic Regression Dichotomous dependent variable (yes/no, died/
More informationDuration Models: Modeling Strategies
Bradford S., UC-Davis, Dept. of Political Science Duration Models: Modeling Strategies Brad 1 1 Department of Political Science University of California, Davis February 28, 2007 Bradford S., UC-Davis,
More informationBEcon Program, Faculty of Economics, Chulalongkorn University Page 1/7
Mid-term Exam (November 25, 2005, 0900-1200hr) Instructions: a) Textbooks, lecture notes and calculators are allowed. b) Each must work alone. Cheating will not be tolerated. c) Attempt all the tests.
More informationLecture 5a: ARCH Models
Lecture 5a: ARCH Models 1 2 Big Picture 1. We use ARMA model for the conditional mean 2. We use ARCH model for the conditional variance 3. ARMA and ARCH model can be used together to describe both conditional
More informationLogit Models for Binary Data
Chapter 3 Logit Models for Binary Data We now turn our attention to regression models for dichotomous data, including logistic regression and probit analysis These models are appropriate when the response
More informationARIMA ANALYSIS WITH INTERVENTIONS / OUTLIERS
TASK Run intervention analysis on the price of stock M: model a function of the price as ARIMA with outliers and interventions. SOLUTION The document below is an abridged version of the solution provided
More informationLecture Slides. Elementary Statistics Twelfth Edition. by Mario F. Triola. and the Triola Statistics Series. Section 7.4-1
Lecture Slides Elementary Statistics Twelfth Edition and the Triola Statistics Series by Mario F. Triola Section 7.4-1 Chapter 7 Estimates and Sample Sizes 7-1 Review and Preview 7- Estimating a Population
More informationAppendix. Table A.1 (Part A) The Author(s) 2015 G. Chakrabarti and C. Sen, Green Investing, SpringerBriefs in Finance, DOI /
Appendix Table A.1 (Part A) Dependent variable: probability of crisis (own) Method: ML binary probit (quadratic hill climbing) Included observations: 47 after adjustments Convergence achieved after 6 iterations
More informationWesVar uses repeated replication variance estimation methods exclusively and as a result does not offer the Taylor Series Linearization approach.
CHAPTER 9 ANALYSIS EXAMPLES REPLICATION WesVar 4.3 GENERAL NOTES ABOUT ANALYSIS EXAMPLES REPLICATION These examples are intended to provide guidance on how to use the commands/procedures for analysis of
More informationWesVar Analysis Example Replication C7
WesVar Analysis Example Replication C7 WesVar 5.1 is primarily a point and click application and though a text file of commands can be used in the WesVar (V5.1) batch processing environment, all examples
More informationMODEL SELECTION CRITERIA IN R:
1. R 2 statistics We may use MODEL SELECTION CRITERIA IN R R 2 = SS R SS T = 1 SS Res SS T or R 2 Adj = 1 SS Res/(n p) SS T /(n 1) = 1 ( ) n 1 (1 R 2 ). n p where p is the total number of parameters. R
More informationLabor Market Returns to Two- and Four- Year Colleges. Paper by Kane and Rouse Replicated by Andreas Kraft
Labor Market Returns to Two- and Four- Year Colleges Paper by Kane and Rouse Replicated by Andreas Kraft Theory Estimating the return to two-year colleges Economic Return to credit hours or sheepskin effects
More informationModelling the potential human capital on the labor market using logistic regression in R
Modelling the potential human capital on the labor market using logistic regression in R Ana-Maria Ciuhu (dobre.anamaria@hotmail.com) Institute of National Economy, Romanian Academy; National Institute
More informationStep 1: Load the appropriate R package. Step 2: Fit a separate mixed model for each independence claim in the basis set.
Step 1: Load the appropriate R package. You will need two libraries: nlme and lme4. Step 2: Fit a separate mixed model for each independence claim in the basis set. For instance, in Table 2 the first basis
More information6. Genetics examples: Hardy-Weinberg Equilibrium
PBCB 206 (Fall 2006) Instructor: Fei Zou email: fzou@bios.unc.edu office: 3107D McGavran-Greenberg Hall Lecture 4 Topics for Lecture 4 1. Parametric models and estimating parameters from data 2. Method
More informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationEXST7015: Multiple Regression from Snedecor & Cochran (1967) RAW DATA LISTING
Multiple (Linear) Regression Introductory example Page 1 1 options ps=256 ls=132 nocenter nodate nonumber; 3 DATA ONE; 4 TITLE1 ''; 5 INPUT X1 X2 X3 Y; 6 **** LABEL Y ='Plant available phosphorus' 7 X1='Inorganic
More informationParameter Estimation
Parameter Estimation Bret Larget Departments of Botany and of Statistics University of Wisconsin Madison April 12, 2007 Statistics 572 (Spring 2007) Parameter Estimation April 12, 2007 1 / 14 Continue
More informationSTATISTICAL METHODS FOR CATEGORICAL DATA ANALYSIS
STATISTICAL METHODS FOR CATEGORICAL DATA ANALYSIS Daniel A. Powers Department of Sociology University of Texas at Austin YuXie Department of Sociology University of Michigan ACADEMIC PRESS An Imprint of
More informationThe method of Maximum Likelihood.
Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed
More informationBradley-Terry Models. Stat 557 Heike Hofmann
Bradley-Terry Models Stat 557 Heike Hofmann Outline Definition: Bradley-Terry Fitting the model Extension: Order Effects Extension: Ordinal & Nominal Response Repeated Measures Bradley-Terry Model (1952)
More informationThe SAS System 11:03 Monday, November 11,
The SAS System 11:3 Monday, November 11, 213 1 The CONTENTS Procedure Data Set Name BIO.AUTO_PREMIUMS Observations 5 Member Type DATA Variables 3 Engine V9 Indexes Created Monday, November 11, 213 11:4:19
More informationCredit Risk Modelling
Credit Risk Modelling Tiziano Bellini Università di Bologna December 13, 2013 Tiziano Bellini (Università di Bologna) Credit Risk Modelling December 13, 2013 1 / 55 Outline Framework Credit Risk Modelling
More informationStock Market Liquidity and Firm Performance
American-Eurasian Journal of Scientific Research 12 (1): 37-46, 2017 ISSN 1818-6785 IDOSI Publications, 2017 DOI: 10.5829/idosi.aejsr.2017.37.46 Stock Market Liquidity and Firm Performance Tarika Singh
More informationİnsan TUNALI 8 November 2018 Econ 511: Econometrics I. ASSIGNMENT 7 STATA Supplement
İnsan TUNALI 8 November 2018 Econ 511: Econometrics I ASSIGNMENT 7 STATA Supplement. use "F:\COURSES\GRADS\ECON511\SHARE\wages1.dta", clear. generate =ln(wage). scatter sch Q. Do you see a relationship
More informationCredit Risk Modeling Using Excel and VBA with DVD O. Gunter Loffler Peter N. Posch. WILEY A John Wiley and Sons, Ltd., Publication
Credit Risk Modeling Using Excel and VBA with DVD O Gunter Loffler Peter N. Posch WILEY A John Wiley and Sons, Ltd., Publication Preface to the 2nd edition Preface to the 1st edition Some Hints for Troubleshooting
More informationExample 2.3: CEO Salary and Return on Equity. Salary for ROE = 0. Salary for ROE = 30. Example 2.4: Wage and Education
1 Stata Textbook Examples Introductory Econometrics: A Modern Approach by Jeffrey M. Wooldridge (1st & 2d eds.) Chapter 2 - The Simple Regression Model Example 2.3: CEO Salary and Return on Equity summ
More informationAccounting. Stock market liquidity and firm performance. 1. Introduction
Accounting 1 (2015) 29 36 Contents lists available at GrowingScience Accounting homepage: www.growingscience.com/ac/ac.html Stock market liquidity and firm performance Tarika Singh a*, Monika Gupta b and
More informationLogistic Regression Analysis
Revised July 2018 Logistic Regression Analysis This set of notes shows how to use Stata to estimate a logistic regression equation. It assumes that you have set Stata up on your computer (see the Getting
More informationGeneralized Multilevel Regression Example for a Binary Outcome
Psy 510/610 Multilevel Regression, Spring 2017 1 HLM Generalized Multilevel Regression Example for a Binary Outcome Specifications for this Bernoulli HLM2 run Problem Title: no title The data source for
More informationAttachment JEM 4 Hearing Exhibit 116 Page 1 of 11] Residential Sales 5,817,938 90,842,431 96,660,368. Normal HDD
Residential Sales Page 1 of 11] F=C*(D E)*B H=G F B C D E F G H Month HDD Coefficient Customers HDD Normal HDD Weather Impact Actual Billed WN Billed Jan-16 0.011693 1,256,043 1,183 1,102 1,193,343 18,806,509
More informationLet us assume that we are measuring the yield of a crop plant on 5 different plots at 4 different observation times.
Mixed-effects models An introduction by Christoph Scherber Up to now, we have been dealing with linear models of the form where ß0 and ß1 are parameters of fixed value. Example: Let us assume that we are
More informationECON Introductory Econometrics. Seminar 4. Stock and Watson Chapter 8
ECON4150 - Introductory Econometrics Seminar 4 Stock and Watson Chapter 8 empirical exercise E8.2: Data 2 In this exercise we use the data set CPS12.dta Each month the Bureau of Labor Statistics in the
More informationLapse Modeling for the Post-Level Period
Lapse Modeling for the Post-Level Period A Practical Application of Predictive Modeling JANUARY 2015 SPONSORED BY Committee on Finance Research PREPARED BY Richard Xu, FSA, Ph.D. Dihui Lai, Ph.D. Minyu
More informationMortality Projections Committee WORKING PAPER 91. CMI Mortality Projections Model consultation technical paper. August 2016 ISSN
ISSN 2044-3145 Mortality Projections Committee WORKING PAPER 91 CMI Mortality Projections Model consultation technical paper August 2016 NOTE: This document is being made available publicly and its use
More information1 Stat 8053, Fall 2011: GLMMs
Stat 805, Fall 0: GLMMs The data come from a 988 fertility survey in Bangladesh. Data were collected on 94 women grouped into 60 districts. The response of interest is whether or not the woman is using
More informationYou created this PDF from an application that is not licensed to print to novapdf printer (http://www.novapdf.com)
Monday October 3 10:11:57 2011 Page 1 (R) / / / / / / / / / / / / Statistics/Data Analysis Education Box and save these files in a local folder. name:
More informationGamma Distribution Fitting
Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics
More informationModule 4 Bivariate Regressions
AGRODEP Stata Training April 2013 Module 4 Bivariate Regressions Manuel Barron 1 and Pia Basurto 2 1 University of California, Berkeley, Department of Agricultural and Resource Economics 2 University of
More information11. Logistic modeling of proportions
11. Logistic modeling of proportions Retrieve the data File on main menu Open worksheet C:\talks\strirling\employ.ws = Note Postcode is neighbourhood in Glasgow Cell is element of the table for each postcode
More informationAssignment #5 Solutions: Chapter 14 Q1.
Assignment #5 Solutions: Chapter 14 Q1. a. R 2 is.037 and the adjusted R 2 is.033. The adjusted R 2 value becomes particularly important when there are many independent variables in a multiple regression
More informationA Test of the Normality Assumption in the Ordered Probit Model *
A Test of the Normality Assumption in the Ordered Probit Model * Paul A. Johnson Working Paper No. 34 March 1996 * Assistant Professor, Vassar College. I thank Jahyeong Koo, Jim Ziliak and an anonymous
More informationNon linearity issues in PD modelling. Amrita Juhi Lucas Klinkers
Non linearity issues in PD modelling Amrita Juhi Lucas Klinkers May 2017 Content Introduction Identifying non-linearity Causes of non-linearity Performance 2 Content Introduction Identifying non-linearity
More informationMaximum Likelihood Estimation
Maximum Likelihood Estimation The likelihood and log-likelihood functions are the basis for deriving estimators for parameters, given data. While the shapes of these two functions are different, they have
More informationQuantitative Techniques Term 2
Quantitative Techniques Term 2 Laboratory 7 2 March 2006 Overview The objective of this lab is to: Estimate a cost function for a panel of firms; Calculate returns to scale; Introduce the command cluster
More informationGeneralized Linear Models II: Applying GLMs in Practice Duncan Anderson MA FIA Watson Wyatt LLP W W W. W A T S O N W Y A T T.
Generalized Linear Models II: Applying GLMs in Practice Duncan Anderson MA FIA Watson Wyatt LLP W W W. W A T S O N W Y A T T. C O M Agenda Introduction / recap Model forms and model validation Aliasing
More informationHeteroskedasticity. . reg wage black exper educ married tenure
Heteroskedasticity. reg Source SS df MS Number of obs = 2,380 -------------+---------------------------------- F(2, 2377) = 72.38 Model 14.4018246 2 7.20091231 Prob > F = 0.0000 Residual 236.470024 2,377.099482551
More informationPort(A,B) is a combination of two stocks, A and B, with standard deviations A and B. A,B = correlation (A,B) = 0.
Corporate Finance, Module 6: Risk, Return, and Cost of Capital Practice Problems (The attached PDF file has better formatting.) Updated: July 19, 2007 Exercise 6.1: Minimum Variance Portfolio Port(A,B)
More informationbook 2014/5/6 15:21 page 261 #285
book 2014/5/6 15:21 page 261 #285 Chapter 10 Simulation Simulations provide a powerful way to answer questions and explore properties of statistical estimators and procedures. In this chapter, we will
More informationStudy Guide on Testing the Assumptions of Age-to-Age Factors - G. Stolyarov II 1
Study Guide on Testing the Assumptions of Age-to-Age Factors - G. Stolyarov II 1 Study Guide on Testing the Assumptions of Age-to-Age Factors for the Casualty Actuarial Society (CAS) Exam 7 and Society
More informationCREDIT RISK MODELING IN R. Logistic regression: introduction
CREDIT RISK MODELING IN R Logistic regression: introduction Final data structure > str(training_set) 'data.frame': 19394 obs. of 8 variables: $ loan_status : Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1
More informationDetermining Probability Estimates From Logistic Regression Results Vartanian: SW 541
Determining Probability Estimates From Logistic Regression Results Vartanian: SW 541 In determining logistic regression results, you will generally be given the odds ratio in the SPSS or SAS output. However,
More information