Math 70 Homework 8 Michael Downs. and, given n iid observations, has log-likelihood function: k i n (2) = 1 n
|
|
- Aubrey Stafford
- 6 years ago
- Views:
Transcription
1 1. (a) The poisson distribution has density function f(k; λ) = λk e λ k! and, given n iid observations, has log-likelihood function: l(λ; k 1,..., k n ) = ln(λ) k i nλ ln(k i!) (1) and score equation: dl dλ = 1 λ k i n (2) which gives the maximum likelihood estimator ˆλ = 1 n n k i = k. The fisher information can be determined from taking the variance of the score equation for only one observation k: I(λ) = var( k var(k) 1) = = 1 (3) λ λ 2 λ so the inverse Fisher information is I(λ) 1 = λ. Thus the Wald test for the hypothesis H 0 : λ = λ 0 given n iid observations is: n(ˆλ λ0 ) N (0, 1) (4) λ0 For the Likelihood Ratio test, the restricted maximum likelihood value is: l(λ 0 ; k) = ln(λ 0 ) k i nλ 0 ln(k i!) (5) and the maximum likelihood value is: l(ˆλ; k) = ln(ˆλ) k i nˆλ ln(k i!) (6) ( ) n λ Their difference is ln 0 ˆλ k i + n(ˆλ λ 0 ) giving the likelihood ratio test for the hypothesis H 0 : λ = λ 0 : ( ( ) ) λ0 2 ln k i + n(ˆλ λ 0 ) χ ˆλ 2 (1) (7) (b) Code: 1 hw8 = f u n c t i o n ( job =1.2, lambda = 10, nexp = 1000, alpha =. 0 5 ) { 3 i f ( job ==1.2) { 5 n = c ( 2, 5, 1 0, 2 0 ) ln = l e n g t h ( n ) 7 alphaexp=matrix ( ncol =2,nrow=ln ) f o r ( i in 1 : ln ) 9 { 1
2 pvwald=pvlr=rep (NA, nexp) 11 f o r ( j in 1 : nexp) { 13 obs = r p o i s ( n [ i ], lambda = lambda ) mle = mean( obs ) 15 waldstat = s q r t ( n [ i ] ) ( mle lambda ) / s q r t ( lambda ) pvwald[ j ] = 2 (1 pnorm ( abs ( waldstat ) ) ) 17 LRstat = 2 ( l o g ( lambda/mle ) sum( obs ) + n [ i ] ( mle lambda ) ) 19 pvlr [ j ] = 1 pchisq ( LRstat, df =1) 21 alphaexp [ i,1]= mean(pvwald<alpha ) alphaexp [ i,2]= mean(pvlr<alpha, na. rm=t) 23 matplot (n, alphaexp, ylim=c ( 0,. 2 ), type= b, xlab= Sample s i z e, n, ylab = S i g n i f i c a n c e l e v e l ) 25 segments ( 1, alpha, 1 0 0, alpha, lwd=3) t i t l e ( Wald vs LR f o r o b t a i n i n g e m p i r i c a l s i g n i f i c a n c e l e v e l ) 27 legend ( 1 0,. 2, c ( 1=Wald Z t e s t, 2=LR chi square 1 t e s t ), l t y =1:2, c o l =1:2) r eturn ( cbind (n, alphaexp ) ) 29 > hw8(1.2, 10, ,.05) n [1,] [2,] [3,] [4,] Plot: Wald vs LR for obtaining empirical significance level Significance level Sample size, n 1=Wald Z test 2=LR chi square 1 te (c) Code: 2
3 obs = c ( 9, 1 4, 1 7, 1 2, 1 0 ) 2 mle = mean( obs ) n u l l = n = l e n g t h ( obs ) 6 waldstat = s q r t ( n ) ( mle n u l l ) / s q r t ( n u l l ) pwald = 2 (1 pnorm ( abs ( waldstat ) ) ) 8 LRstat = 2 ( l o g ( n u l l /mle ) sum( obs ) + n ( mle n u l l ) ) 10 p l r = 1 pchisq ( LRstat, df =1) 12 p r i n t ( p value f o r wald t e s t : ) p r i n t ( pwald ) 14 p r i n t ( p value f o r l r t e s t : ) p r i n t ( p l r ) 16 p r i n t ( in both c a s e s we cannot r e j e c t the n u l l ) > hw8(1.3) [1] "p-value for wald test: " [1] [1] "p-value for lr test: " [1] [1] "in both cases we cannot reject the null" 2. (a) Using the derivation for the Wald test from the lecture notes, n(ˆp p0 ) p0 (1 p 0 ) N (0, 1) where ˆp is the MLE m where m is the number of successes in n trials we can derive n (using a method similar to the one from the lecture notes) the two power functions under null H 0 : p = p 0 at the α significance level: ( ) ( ) n(p Power 1 (p) = Φ Z 1 α p0 ) n(p + 1 Φ Z 2 1 α 2 p0 (1 p 0 ) p0 ) (8) p0 (1 p 0 ) Power 2 (p) = Φ ( Z 1 α 2 n(p p0 ) p(1 p) ) + 1 Φ ( ) n(p Z 1 α p0 ) 2 p(1 p) (9) These gives the approximate probability of rejecting the null when the alternative hypothesis is correct. (b) Code: 3
4 hw8 = f u n c t i o n ( job =1.2, lambda = 10, nexp = 1000, alpha =. 0 5, p0 =. 2 5, n = 10) 2 { i f ( job ==2.2) 4 { z = qnorm (1 alpha / 2) 6 power1 = f u n c t i o n ( pr ) 8 { d = s q r t ( n ) ( pr p0 ) / s q r t ( p0 (1 p0 ) ) 10 r eturn ( pnorm( z d ) + 1 pnorm ( z d ) ) 12 power2 = f u n c t i o n ( pr ) 14 { d = s q r t ( n ) ( pr p0 ) / s q r t ( pr (1 pr ) ) r eturn ( pnorm( z d ) + 1 pnorm ( z d ) ) 20 x = seq ( from =.001, to =. 9 9, l ength = 100) y1 = power1 ( x ) 22 y2 = power2 ( x ) p l o t ( x, y1, type = l, c o l = 2, xlab = a l t e r n a t i v e p, ylab = p r o b a b i l i t y o f r e j e c t i n g n u l l ) 24 l i n e s ( x, y2, type = l, c o l = 3) 26 emp. pow1=rep ( 0, l e n g t h ( x ) ) emp. pow2=rep ( 0, l e n g t h ( x ) ) 28 counter = 1 f o r ( i in x ) 30 { obs = rbinom (nexp, n, i ) 32 obs = obs /n s = s q r t ( n ) ( obs p0 ) / s q r t ( p0 (1 p0 ) ) 34 emp. pow1 [ counter ] = mean( abs ( s )>z ) s2 = s q r t ( n ) ( obs p0 ) / s q r t ( i (1 i ) ) 36 emp. pow2 [ counter ] = mean( abs ( s2 )>z ) 38 counter = counter l i n e s ( x, emp. pow1, type = l, c o l = 4) l i n e s ( x, emp. pow2, type = l, c o l = 5) 42 segments ( 100, alpha, 1 0 0, alpha, l t y =2, c o l =2) segments ( p0, 1,p0, 1, c o l =3) 44 t i t l e ( Four power f u n c t i o n s ) 46 legend (. 6,. 5, c ( I at n u l l, I at a l t e r n a t i v e, Empirical n u l l, Empirical a l t ), l t y = 1, c o l =2:5) 48 p r i n t ( sum o f d i f f e r e n c e o f squares f o r e m p i r i c a l 1 and power1 : ) p r i n t (sum ( ( emp. pow1 y1 ) ˆ2) ) 50 p r i n t ( sum o f d i f f e r e n c e o f squares f o r e m p i r i c a l 2 and power2 : ) 4
5 52 p r i n t (sum ( ( emp. pow2 y2 ) ˆ2) ) p r i n t ( power f u n c t i o n s at p=.35 ) 54 p r i n t ( power1 (. 3 5 ) ) p r i n t ( power2 (. 3 5 ) ) 56 p r i n t ( e m p i r i c a l value 1 at p=.35 ) obs = rbinom (nexp, n,. 3 5 ) 58 obs = obs /n s = s q r t ( n ) ( obs p0 ) / s q r t ( p0 (1 p0 ) ) 60 emp. pow = mean( abs ( s )>z ) p r i n t (emp. pow) 62 p r i n t ( e m p i r i c a l value 2 at p=.35 ) 64 s = s q r t ( n ) ( obs p0 ) / s q r t (. 3 5 (1.35) ) emp. pow = mean( abs ( s )>z ) 66 p r i n t (emp. pow) > hw8(2.2, nexp = , n =10) [1] "sum of difference of squares for empirical1 and power1: " [1] [1] "sum of difference of squares for empirical2 and power2: " [1] [1] "power functions at p=.35" [1] [1] [1] "empirical value 1 at p=.35" [1] [1] "empirical value 2 at p=.35" [1]
6 Four power functions probability of rejecting null I at null I at alternative Empirical null Empirical alt alternative p The version where the Fisher information is evaluated at the alternative hypothesis appears to be better because it expresses the fact that p can only takes values from 0 to 1. (c) The derivation for the formula is similar to the derivation from the lecture notes. Code: p =. 4 2 p0 =. 2 5 d e l t a = p p0 4 z = qnorm ( / 2) z2 = qnorm (. 1 ) 6 n1 = p0 (1 p0 ) / d e l t a ˆ2 ( z z2 ) ˆ2 8 n2 = p (1 p ) / d e l t a ˆ2 ( z z2 ) ˆ2 10 p r i n t ( number o f mice r e q u i r e d using n u l l f i s h e r power f u n c t i o n ) p r i n t ( c e i l i n g ( n1 ) ) 12 p r i n t ( number o f mice r e q u i r e d using a l t f i s h e r power f u n c t i o n ) p r i n t ( c e i l i n g ( n2 ) ) [1] "number of mice required using null fisher power function" [1] 88 [1] "number of mice required using alt fisher power function" 6
7 [1] 113 7
Econ 583 Homework 7 Suggested Solutions: Wald, LM and LR based on GMM and MLE
Econ 583 Homework 7 Suggested Solutions: Wald, LM and LR based on GMM and MLE Eric Zivot Winter 013 1 Wald, LR and LM statistics based on generalized method of moments estimation Let 1 be an iid sample
More informationLoglikelihood and Confidence Intervals
Stat 504, Lecture 2 1 Loglikelihood and Confidence Intervals The loglikelihood function is defined to be the natural logarithm of the likelihood function, l(θ ; x) = log L(θ ; x). For a variety of reasons,
More informationLikelihood Construction, Inference for Parametric Survival Distributions
Week 1 Likelihood Construction, Inference for Parametric Survival Distributions In this section we obtain the likelihood function for noninformatively rightcensored survival data and indicate how to make
More informationStatistics for laboratory scientists II
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike License. Your use of this material constitutes acceptance of that license and the conditions of use of materials on this
More informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
More informationInference for Single Proportions and Means T.Scofield
Inference for Single Proportions and Means TScofield Confidence Intervals for Single Proportions and Means A CI gives upper and lower bounds between which we hope to capture the (fixed) population parameter
More informationSection 4.6 Simple Linear Regression
Section 4.6 Simple Linear Regression Objectives ˆ Basic philosophy of SLR and the regression assumptions ˆ Point & interval estimation of the model parameters, and how to make predictions ˆ Point and interval
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationProbability theory and inference statistics! Dr. Paola Grosso! SNE research group!! (preferred!)!!
Probability theory and inference statistics Dr. Paola Grosso SNE research group p.grosso@uva.nl paola.grosso@os3.nl (preferred) Roadmap Lecture 1: Monday Sep. 22nd Collecting data Presenting data Descriptive
More informationMaximum Likelihood Large Sample Theory
Maximum Likelihood Large Sample Theory MIT 18.443 Dr. Kempthorne Spring 2015 1 Outline 1 Large Sample Theory of Maximum Likelihood Estimates 2 Asymptotic Results: Overview Asymptotic Framework Data Model
More informationLikelihoods. P (Y = y) = f(y). For example, suppose Y has a geometric distribution on 1, 2,... with parameter p. Then the pmf is
Likelihoods The distribution of a random variable Y with a discrete sample space (e.g. a finite sample space or the integers) can be characterized by its probability mass function (pmf): P (Y = y) = f(y).
More informationTUTORIAL 8 SOLUTIONS #
TUTORIAL 8 SOLUTIONS #9.11.21 Suppose that a single observation X is taken from a uniform density on [0,θ], and consider testing H 0 : θ = 1 versus H 1 : θ =2. (a) Find a test that has significance level
More informationEconomics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1,
Economics 520 Lecture Note 9: Hypothesis Testing via the Neyman-Pearson Lemma CB 8., 8.3.-8.3.3 Uniformly Most Powerful Tests and the Neyman-Pearson Lemma Let s return to the hypothesis testing problem
More informationInstitute of Actuaries of India
Institute of Actuaries of India Subject CT3 Probability & Mathematical Statistics May 2011 Examinations INDICATIVE SOLUTION Introduction The indicative solution has been written by the Examiners with the
More informationSome General Types of Tests
Some General Types of Tests We may not be able to find a UMP or UMPU test in a given situation. In that case, we may use test of some general class of tests that often have good asymptotic properties.
More informationModeling Overdispersion
James H. Steiger Department of Psychology and Human Development Vanderbilt University Regression Modeling, 2009 1 Introduction 2 Introduction In this lecture we discuss the problem of overdispersion in
More informationChapter 4: Constrained estimators and tests in the multiple linear regression model (Part III)
Chapter 4: Constrained estimators and tests in the multiple linear regression model (Part III) Florian Pelgrin HEC September-December 2010 Florian Pelgrin (HEC) Constrained estimators September-December
More information" M A #M B. Standard deviation of the population (Greek lowercase letter sigma) σ 2
Notation and Equations for Final Exam Symbol Definition X The variable we measure in a scientific study n The size of the sample N The size of the population M The mean of the sample µ The mean of the
More informationMultiple Sample Categorical Data
Multiple Sample Categorical Data paired and unpaired data, goodness-of-fit testing, testing for independence University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html
More informationMIT Spring 2015
Assessing Goodness Of Fit MIT 8.443 Dr. Kempthorne Spring 205 Outline 2 Poisson Distribution Counts of events that occur at constant rate Counts in disjoint intervals/regions are independent If intervals/regions
More informationExercises and Answers to Chapter 1
Exercises and Answers to Chapter The continuous type of random variable X has the following density function: a x, if < x < a, f (x), otherwise. Answer the following questions. () Find a. () Obtain mean
More informationOutline. Unit 3: Inferential Statistics for Continuous Data. Outline. Inferential statistics for continuous data. Inferential statistics Preliminaries
Unit 3: Inferential Statistics for Continuous Data Statistics for Linguists with R A SIGIL Course Designed by Marco Baroni 1 and Stefan Evert 1 Center for Mind/Brain Sciences (CIMeC) University of Trento,
More informationStat 135, Fall 2006 A. Adhikari HOMEWORK 6 SOLUTIONS
Stat 135, Fall 2006 A. Adhikari HOMEWORK 6 SOLUTIONS 1a. Under the null hypothesis X has the binomial (100,.5) distribution with E(X) = 50 and SE(X) = 5. So P ( X 50 > 10) is (approximately) two tails
More informationsimple if it completely specifies the density of x
3. Hypothesis Testing Pure significance tests Data x = (x 1,..., x n ) from f(x, θ) Hypothesis H 0 : restricts f(x, θ) Are the data consistent with H 0? H 0 is called the null hypothesis simple if it completely
More informationGOV 2001/ 1002/ E-2001 Section 3 Theories of Inference
GOV 2001/ 1002/ E-2001 Section 3 Theories of Inference Solé Prillaman Harvard University February 11, 2015 1 / 48 LOGISTICS Reading Assignment- Unifying Political Methodology chs 2 and 4. Problem Set 3-
More informationEstimation MLE-Pandemic data MLE-Financial crisis data Evaluating estimators. Estimation. September 24, STAT 151 Class 6 Slide 1
Estimation September 24, 2018 STAT 151 Class 6 Slide 1 Pandemic data Treatment outcome, X, from n = 100 patients in a pandemic: 1 = recovered and 0 = not recovered 1 1 1 0 0 0 1 1 1 0 0 1 0 1 0 0 1 1 1
More informationSTAT 536: Genetic Statistics
STAT 536: Genetic Statistics Tests for Hardy Weinberg Equilibrium Karin S. Dorman Department of Statistics Iowa State University September 7, 2006 Statistical Hypothesis Testing Identify a hypothesis,
More informationExercises Chapter 4 Statistical Hypothesis Testing
Exercises Chapter 4 Statistical Hypothesis Testing Advanced Econometrics - HEC Lausanne Christophe Hurlin University of Orléans December 5, 013 Christophe Hurlin (University of Orléans) Advanced Econometrics
More informationThe Chi-Square and F Distributions
Department of Psychology and Human Development Vanderbilt University Introductory Distribution Theory 1 Introduction 2 Some Basic Properties Basic Chi-Square Distribution Calculations in R Convergence
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationAnswer Key for STAT 200B HW No. 7
Answer Key for STAT 200B HW No. 7 May 5, 2007 Problem 2.2 p. 649 Assuming binomial 2-sample model ˆπ =.75, ˆπ 2 =.6. a ˆτ = ˆπ 2 ˆπ =.5. From Ex. 2.5a on page 644: ˆπ ˆπ + ˆπ 2 ˆπ 2.75.25.6.4 = + =.087;
More informationBIOSTATS Intermediate Biostatistics Spring 2017 Exam 2 (Units 3, 4 & 5) Practice Problems SOLUTIONS
BIOSTATS 640 - Intermediate Biostatistics Spring 2017 Exam 2 (Units 3, 4 & 5) Practice Problems SOLUTIONS Practice Question 1 Both the Binomial and Poisson distributions have been used to model the quantal
More informationLet us first identify some classes of hypotheses. simple versus simple. H 0 : θ = θ 0 versus H 1 : θ = θ 1. (1) one-sided
Let us first identify some classes of hypotheses. simple versus simple H 0 : θ = θ 0 versus H 1 : θ = θ 1. (1) one-sided H 0 : θ θ 0 versus H 1 : θ > θ 0. (2) two-sided; null on extremes H 0 : θ θ 1 or
More informationTopic 15: Simple Hypotheses
Topic 15: November 10, 2009 In the simplest set-up for a statistical hypothesis, we consider two values θ 0, θ 1 in the parameter space. We write the test as H 0 : θ = θ 0 versus H 1 : θ = θ 1. H 0 is
More informationMaximum-Likelihood Estimation: Basic Ideas
Sociology 740 John Fox Lecture Notes Maximum-Likelihood Estimation: Basic Ideas Copyright 2014 by John Fox Maximum-Likelihood Estimation: Basic Ideas 1 I The method of maximum likelihood provides estimators
More informationParametric Modelling of Over-dispersed Count Data. Part III / MMath (Applied Statistics) 1
Parametric Modelling of Over-dispersed Count Data Part III / MMath (Applied Statistics) 1 Introduction Poisson regression is the de facto approach for handling count data What happens then when Poisson
More informationHorse Kick Example: Chi-Squared Test for Goodness of Fit with Unknown Parameters
Math 3080 1. Treibergs Horse Kick Example: Chi-Squared Test for Goodness of Fit with Unknown Parameters Name: Example April 4, 014 This R c program explores a goodness of fit test where the parameter is
More informationPart 1.) We know that the probability of any specific x only given p ij = p i p j is just multinomial(n, p) where p k1 k 2
Problem.) I will break this into two parts: () Proving w (m) = p( x (m) X i = x i, X j = x j, p ij = p i p j ). In other words, the probability of a specific table in T x given the row and column counts
More informationStatistical Computing with R MATH , Set 6 (Monte Carlo Methods in Statistical Inference)
Statistical Computing with R MATH 6382 1, Set 6 (Monte Carlo Methods in Statistical Inference) Tamer Oraby UTRGV tamer.oraby@utrgv.edu 1 Based on textbook. Last updated November 14, 2016 Tamer Oraby (University
More informationHYPOTHESIS TESTING: THE CHI-SQUARE STATISTIC
1 HYPOTHESIS TESTING: THE CHI-SQUARE STATISTIC 7 steps of Hypothesis Testing 1. State the hypotheses 2. Identify level of significant 3. Identify the critical values 4. Calculate test statistics 5. Compare
More informationMunitions Example: Negative Binomial to Predict Number of Accidents
Math 37 1. Treibergs Munitions Eample: Negative Binomial to Predict Number of Accidents Name: Eample July 17, 211 This problem comes from Bulmer, Principles of Statistics, Dover, 1979, which is a republication
More informationWald s theorem and the Asimov data set
Wald s theorem and the Asimov data set Eilam Gross & Ofer Vitells ATLAS statistics forum, Dec. 009 1 Outline We have previously guessed that the median significance of many toy MC experiments could be
More informationTesting and Model Selection
Testing and Model Selection This is another digression on general statistics: see PE App C.8.4. The EViews output for least squares, probit and logit includes some statistics relevant to testing hypotheses
More informationPump failure data. Pump Failures Time
Outline 1. Poisson distribution 2. Tests of hypothesis for a single Poisson mean 3. Comparing multiple Poisson means 4. Likelihood equivalence with exponential model Pump failure data Pump 1 2 3 4 5 Failures
More informationScore Test and Likelihood Ratio Test for Zero- Inflated Binomial Distribution and Geometric Distribution
Western Kentucky University TopSCHOLAR Masters Theses & Specialist Projects Graduate School Spring 2018 Score Test and Likelihood Ratio Test for Zero- Inflated Binomial Distribution and Geometric Distribution
More informationCh 7 : One-Sample Test of Hypothesis
Summer 2017 UAkron Dept. of Stats [3470 : 461/561] Applied Statistics Ch 7 : One-Sample Test of Hypothesis Contents 1 Preliminaries 3 1.1 Test of Hypothesis...............................................................
More information2.6.3 Generalized likelihood ratio tests
26 HYPOTHESIS TESTING 113 263 Generalized likelihood ratio tests When a UMP test does not exist, we usually use a generalized likelihood ratio test to verify H 0 : θ Θ against H 1 : θ Θ\Θ It can be used
More informationThe Multinomial Model
The Multinomial Model STA 312: Fall 2012 Contents 1 Multinomial Coefficients 1 2 Multinomial Distribution 2 3 Estimation 4 4 Hypothesis tests 8 5 Power 17 1 Multinomial Coefficients Multinomial coefficient
More informationA3. Statistical Inference Hypothesis Testing for General Population Parameters
Appendix / A3. Statistical Inference / General Parameters- A3. Statistical Inference Hypothesis Testing for General Population Parameters POPULATION H 0 : θ = θ 0 θ is a generic parameter of interest (e.g.,
More informationLecture 32: Asymptotic confidence sets and likelihoods
Lecture 32: Asymptotic confidence sets and likelihoods Asymptotic criterion In some problems, especially in nonparametric problems, it is difficult to find a reasonable confidence set with a given confidence
More informationCSE 312 Final Review: Section AA
CSE 312 TAs December 8, 2011 General Information General Information Comprehensive Midterm General Information Comprehensive Midterm Heavily weighted toward material after the midterm Pre-Midterm Material
More informationComparing two independent samples
In many applications it is necessary to compare two competing methods (for example, to compare treatment effects of a standard drug and an experimental drug). To compare two methods from statistical point
More informationName: SOLUTIONS Exam 01 (Midterm Part 2 take home, open everything)
Name: SOLUTIONS Exam 01 (Midterm Part 2 take home, open everything) To help you budget your time, questions are marked with *s. One * indicates a straightforward question testing foundational knowledge.
More informationSession 3 The proportional odds model and the Mann-Whitney test
Session 3 The proportional odds model and the Mann-Whitney test 3.1 A unified approach to inference 3.2 Analysis via dichotomisation 3.3 Proportional odds 3.4 Relationship with the Mann-Whitney test Session
More informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationV. Introduction to Likelihood
GEOS 33000/EVOL 33000 17 January 2005 Modified July 25, 2007 Page 1 V. Introduction to Likelihood 1 Likelihood basics 1.1 Under assumed model, let X denote data and Θ the parameter(s) of the model. 1.2
More informationUnit 9 ONE Sample Inference SOLUTIONS
BIOSTATS 540 Fall 017 Introductory Biostatistics Solutions Unit 9 ONE Sample Inference SOLUTIONS 1. This exercise gives you practice calculating a confidence interval for the mean of a Normal distribution
More informationML Testing (Likelihood Ratio Testing) for non-gaussian models
ML Testing (Likelihood Ratio Testing) for non-gaussian models Surya Tokdar ML test in a slightly different form Model X f (x θ), θ Θ. Hypothesist H 0 : θ Θ 0 Good set: B c (x) = {θ : l x (θ) max θ Θ l
More informationFinal Exam. 1. (6 points) True/False. Please read the statements carefully, as no partial credit will be given.
1. (6 points) True/False. Please read the statements carefully, as no partial credit will be given. (a) If X and Y are independent, Corr(X, Y ) = 0. (b) (c) (d) (e) A consistent estimator must be asymptotically
More information1/24/2008. Review of Statistical Inference. C.1 A Sample of Data. C.2 An Econometric Model. C.4 Estimating the Population Variance and Other Moments
/4/008 Review of Statistical Inference Prepared by Vera Tabakova, East Carolina University C. A Sample of Data C. An Econometric Model C.3 Estimating the Mean of a Population C.4 Estimating the Population
More informationMAS361. MAS361 1 Turn Over SCHOOL OF MATHEMATICS AND STATISTICS. Medical Statistics
t r r t r r t r t s s MAS361 SCHOOL OF MATHEMATICS AND STATISTICS Medical Statistics Autumn Semester 2015 16 2 hours t s 2 r t t 1 t t r t t r s t rs t2 r t s q st s r t r t r 2 t st s rs q st s rr2 q
More informationSummary of Chapter 7 (Sections ) and Chapter 8 (Section 8.1)
Summary of Chapter 7 (Sections 7.2-7.5) and Chapter 8 (Section 8.1) Chapter 7. Tests of Statistical Hypotheses 7.2. Tests about One Mean (1) Test about One Mean Case 1: σ is known. Assume that X N(µ, σ
More informationBinary response data
Binary response data A Bernoulli trial is a random variable that has two points in its sample space. The two points may be denoted success/failure, heads/tails, yes/no, 0/1, etc. The probability distribution
More informationV. Introduction to Likelihood
GEOS 36501/EVOL 33001 20 January 2012 Page 1 of 24 V. Introduction to Likelihood 1 Likelihood basics 1.1 Under assumed model, let X denote data and Θ the parameter(s) of the model. 1.2 Probability: P (X
More informationST417 Introduction to Bayesian Modelling. Conjugate Modelling (Poisson-Gamma)
ST417 Introduction to Bayesian Modelling Conjugate Modelling (Poisson-Gamma) Slides based on the source provided courtesy of Prof. D.Draper, UCSC ST417 Intro to Bayesian Modelling 1 Integer-Valued Outcomes
More informationMath 181B Homework 1 Solution
Math 181B Homework 1 Solution 1. Write down the likelihood: L(λ = n λ X i e λ X i! (a One-sided test: H 0 : λ = 1 vs H 1 : λ = 0.1 The likelihood ratio: where LR = L(1 L(0.1 = 1 X i e n 1 = λ n X i e nλ
More informationPubHlth Intermediate Biostatistics Spring 2015 Exam 2 (Units 3, 4 & 5) Study Guide
PubHlth 640 - Intermediate Biostatistics Spring 2015 Exam 2 (Units 3, 4 & 5) Study Guide Unit 3 (Discrete Distributions) Take care to know how to do the following! Learning Objective See: 1. Write down
More informationNormal (Gaussian) distribution The normal distribution is often relevant because of the Central Limit Theorem (CLT):
Lecture Three Normal theory null distributions Normal (Gaussian) distribution The normal distribution is often relevant because of the Central Limit Theorem (CLT): A random variable which is a sum of many
More informationChapter 11. Hypothesis Testing (II)
Chapter 11. Hypothesis Testing (II) 11.1 Likelihood Ratio Tests one of the most popular ways of constructing tests when both null and alternative hypotheses are composite (i.e. not a single point). Let
More informationIntroduction to Survey Analysis!
Introduction to Survey Analysis! Professor Ron Fricker! Naval Postgraduate School! Monterey, California! Reading Assignment:! 2/22/13 None! 1 Goals for this Lecture! Introduction to analysis for surveys!
More informationSuggested solutions to written exam Jan 17, 2012
LINKÖPINGS UNIVERSITET Institutionen för datavetenskap Statistik, ANd 73A36 THEORY OF STATISTICS, 6 CDTS Master s program in Statistics and Data Mining Fall semester Written exam Suggested solutions to
More informationConfidence Intervals, Testing and ANOVA Summary
Confidence Intervals, Testing and ANOVA Summary 1 One Sample Tests 1.1 One Sample z test: Mean (σ known) Let X 1,, X n a r.s. from N(µ, σ) or n > 30. Let The test statistic is H 0 : µ = µ 0. z = x µ 0
More informationQuantitative Introduction ro Risk and Uncertainty in Business Module 5: Hypothesis Testing
Quantitative Introduction ro Risk and Uncertainty in Business Module 5: Hypothesis Testing M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu October
More informationProbability and Statistics Notes
Probability and Statistics Notes Chapter Seven Jesse Crawford Department of Mathematics Tarleton State University Spring 2011 (Tarleton State University) Chapter Seven Notes Spring 2011 1 / 42 Outline
More informationGoodness of Fit Goodness of fit - 2 classes
Goodness of Fit Goodness of fit - 2 classes A B 78 22 Do these data correspond reasonably to the proportions 3:1? We previously discussed options for testing p A = 0.75! Exact p-value Exact confidence
More informationz and t tests for the mean of a normal distribution Confidence intervals for the mean Binomial tests
z and t tests for the mean of a normal distribution Confidence intervals for the mean Binomial tests Chapters 3.5.1 3.5.2, 3.3.2 Prof. Tesler Math 283 Fall 2018 Prof. Tesler z and t tests for mean Math
More informationBIOS 625 Fall 2015 Homework Set 3 Solutions
BIOS 65 Fall 015 Homework Set 3 Solutions 1. Agresti.0 Table.1 is from an early study on the death penalty in Florida. Analyze these data and show that Simpson s Paradox occurs. Death Penalty Victim's
More informationSimple logistic regression
Simple logistic regression Biometry 755 Spring 2009 Simple logistic regression p. 1/47 Model assumptions 1. The observed data are independent realizations of a binary response variable Y that follows a
More informationHypothesis testing: theory and methods
Statistical Methods Warsaw School of Economics November 3, 2017 Statistical hypothesis is the name of any conjecture about unknown parameters of a population distribution. The hypothesis should be verifiable
More informationInference for Proportions, Variance and Standard Deviation
Inference for Proportions, Variance and Standard Deviation Sections 7.10 & 7.6 Cathy Poliak, Ph.D. cathy@math.uh.edu Office Fleming 11c Department of Mathematics University of Houston Lecture 12 Cathy
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 15: Examples of hypothesis tests (v5) Ramesh Johari ramesh.johari@stanford.edu 1 / 32 The recipe 2 / 32 The hypothesis testing recipe In this lecture we repeatedly apply the
More informationTopic 3: Sampling Distributions, Confidence Intervals & Hypothesis Testing. Road Map Sampling Distributions, Confidence Intervals & Hypothesis Testing
Topic 3: Sampling Distributions, Confidence Intervals & Hypothesis Testing ECO22Y5Y: Quantitative Methods in Economics Dr. Nick Zammit University of Toronto Department of Economics Room KN3272 n.zammit
More informationLogistic Regression Analysis
Logistic Regression Analysis Predicting whether an event will or will not occur, as well as identifying the variables useful in making the prediction, is important in most academic disciplines as well
More informationLecture 10: Generalized likelihood ratio test
Stat 200: Introduction to Statistical Inference Autumn 2018/19 Lecture 10: Generalized likelihood ratio test Lecturer: Art B. Owen October 25 Disclaimer: These notes have not been subjected to the usual
More informationEcon 325: Introduction to Empirical Economics
Econ 325: Introduction to Empirical Economics Lecture 6 Sampling and Sampling Distributions Ch. 6-1 Populations and Samples A Population is the set of all items or individuals of interest Examples: All
More informationTheory of Inference: Homework 4
Theory of Inference: Homework 4 1. Here is a slightly technical question about hypothesis tests. Suppose that Y 1,..., Y n iid Poisson(λ) λ > 0. The two competing hypotheses are H 0 : λ = λ 0 versus H
More informationPsychology 282 Lecture #4 Outline Inferences in SLR
Psychology 282 Lecture #4 Outline Inferences in SLR Assumptions To this point we have not had to make any distributional assumptions. Principle of least squares requires no assumptions. Can use correlations
More informationThe outline for Unit 3
The outline for Unit 3 Unit 1. Introduction: The regression model. Unit 2. Estimation principles. Unit 3: Hypothesis testing principles. 3.1 Wald test. 3.2 Lagrange Multiplier. 3.3 Likelihood Ratio Test.
More informationUNIVERSITY OF TORONTO SCARBOROUGH Department of Computer and Mathematical Sciences FINAL EXAMINATION, APRIL 2013
UNIVERSITY OF TORONTO SCARBOROUGH Department of Computer and Mathematical Sciences FINAL EXAMINATION, APRIL 2013 STAB57H3 Introduction to Statistics Duration: 3 hours Last Name: First Name: Student number:
More informationPoisson Regression. Ryan Godwin. ECON University of Manitoba
Poisson Regression Ryan Godwin ECON 7010 - University of Manitoba Abstract. These lecture notes introduce Maximum Likelihood Estimation (MLE) of a Poisson regression model. 1 Motivating the Poisson Regression
More informationLecture 15. Hypothesis testing in the linear model
14. Lecture 15. Hypothesis testing in the linear model Lecture 15. Hypothesis testing in the linear model 1 (1 1) Preliminary lemma 15. Hypothesis testing in the linear model 15.1. Preliminary lemma Lemma
More informationDrawing Inferences from Statistics Based on Multiyear Asset Returns
Drawing Inferences from Statistics Based on Multiyear Asset Returns Matthew Richardson ames H. Stock FE 1989 1 Motivation Fama and French (1988, Poterba and Summer (1988 document significant negative correlations
More informationCollated responses from R-help on confidence intervals for risk ratios
Collated responses from R-help on confidence intervals for risk ratios Michael E Dewey November, 2006 Introduction This document arose out of a problem assessing a confidence interval for the risk ratio
More informationLogistic Regression. Interpretation of linear regression. Other types of outcomes. 0-1 response variable: Wound infection. Usual linear regression
Logistic Regression Usual linear regression (repetition) y i = b 0 + b 1 x 1i + b 2 x 2i + e i, e i N(0,σ 2 ) or: y i N(b 0 + b 1 x 1i + b 2 x 2i,σ 2 ) Example (DGA, p. 336): E(PEmax) = 47.355 + 1.024
More information1/15. Over or under dispersion Problem
1/15 Over or under dispersion Problem 2/15 Example 1: dogs and owners data set In the dogs and owners example, we had some concerns about the dependence among the measurements from each individual. Let
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Yan Lu Jan, 2018, week 3 1 / 67 Hypothesis tests Likelihood ratio tests Wald tests Score tests 2 / 67 Generalized Likelihood ratio tests Let Y = (Y 1,
More informationTentamentsskrivning: Statistisk slutledning 1. Tentamentsskrivning i Statistisk slutledning MVE155/MSG200, 7.5 hp.
Tentamentsskrivning: Statistisk slutledning 1 Tentamentsskrivning i Statistisk slutledning MVE155/MSG200, 7.5 hp. Tid: tisdagen den 17 mars, 2015 kl 14.00-18.00 Examinator och jour: Serik Sagitov, tel.
More informationH 2 : otherwise. that is simply the proportion of the sample points below level x. For any fixed point x the law of large numbers gives that
Lecture 28 28.1 Kolmogorov-Smirnov test. Suppose that we have an i.i.d. sample X 1,..., X n with some unknown distribution and we would like to test the hypothesis that is equal to a particular distribution
More informationF79SM STATISTICAL METHODS
F79SM STATISTICAL METHODS SUMMARY NOTES 9 Hypothesis testing 9.1 Introduction As before we have a random sample x of size n of a population r.v. X with pdf/pf f(x;θ). The distribution we assign to X is
More informationExperimental Design and Statistical Methods. Workshop LOGISTIC REGRESSION. Jesús Piedrafita Arilla.
Experimental Design and Statistical Methods Workshop LOGISTIC REGRESSION Jesús Piedrafita Arilla jesus.piedrafita@uab.cat Departament de Ciència Animal i dels Aliments Items Logistic regression model Logit
More information