xp (X = x) = P (X = 1) = θ. Hence, the method of moments estimator of θ is

Similar documents
Solutions to selected problems from the midterm exam Math 222 Winter 2015

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

Department of Mathematical and Statistical Sciences University of Alberta

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

F D D D D F. smoothed value of the data including Y t the most recent data.

Statistical Theory MT 2009 Problems 1: Solution sketches

Moment Generating Function

Statistical Theory MT 2008 Problems 1: Solution sketches

OLS bias for econometric models with errors-in-variables. The Lucas-critique Supplementary note to Lecture 17

14.02 Principles of Macroeconomics Fall 2005

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

Modified Ratio and Product Estimators for Estimating Population Mean in Two-Phase Sampling

DETERMINATION OF PARTICULAR SOLUTIONS OF NONHOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS BY DISCRETE DECONVOLUTION

Math 2414 Homework Set 7 Solutions 10 Points

Institute of Actuaries of India

Calculus BC 2015 Scoring Guidelines

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

SUMMATION OF INFINITE SERIES REVISITED

Summary. Recap. Last Lecture. Let W n = W n (X 1,, X n ) = W n (X) be a sequence of estimators for

Lecture 7: Properties of Random Samples

Pure Math 30: Explained!

Review Questions, Chapters 8, 9. f(y) = 0, elsewhere. F (y) = f Y(1) = n ( e y/θ) n 1 1 θ e y/θ = n θ e yn

Direction: This test is worth 250 points. You are required to complete this test within 50 minutes.

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

Math 6710, Fall 2016 Final Exam Solutions

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013

Solutions to Problems 3, Level 4

Theorem. Let H be a class of functions from a measurable space T to R. Assume that for every ɛ > 0 there exists a finite set of brackets

LIMITS OF FUNCTIONS (I)

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

STK4080/9080 Survival and event history analysis

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

Stat410 Probability and Statistics II (F16)

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 17, 2013

N! AND THE GAMMA FUNCTION

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Big O Notation for Time Complexity of Algorithms

MATH 472 / SPRING 2013 ASSIGNMENT 2: DUE FEBRUARY 4 FINALIZED

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

Sampling Example. ( ) δ ( f 1) (1/2)cos(12πt), T 0 = 1

ME 3210 Mechatronics II Laboratory Lab 6: Second-Order Dynamic Response

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Research Design - - Topic 2 Inferential Statistics: The t-test 2010 R.C. Gardner, Ph.D. Independent t-test

Affine term structure models

COS 522: Complexity Theory : Boaz Barak Handout 10: Parallel Repetition Lemma

ECE-314 Fall 2012 Review Questions

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017

Lecture 11 and 12: Basic estimation theory

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

12 Getting Started With Fourier Analysis

12 th Mathematics Objective Test Solutions

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

Review Answers for E&CE 700T02

K3 p K2 p Kp 0 p 2 p 3 p

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

1 Notes on Little s Law (l = λw)

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

The Moment Approximation of the First Passage Time For The Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

Extremal graph theory II: K t and K t,t

Problem Set 4 Due Oct, 12

1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.

10.3 Autocorrelation Function of Ergodic RP 10.4 Power Spectral Density of Ergodic RP 10.5 Normal RP (Gaussian RP)

Additional Tables of Simulation Results

Supplementary Information for Thermal Noises in an Aqueous Quadrupole Micro- and Nano-Trap

Exercise 3 Stochastic Models of Manufacturing Systems 4T400, 6 May

Lecture 9: Polynomial Approximations

Delta Method on Bootstrapping of Autoregressive Process. Abstract

Mathematical Issues. Written by Okito Yamashita, Last revised 2009/07/27

Section 5: Chain Rule

Lecture 12: September 27

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

Section 8 Convolution and Deconvolution

Inference of the Second Order Autoregressive. Model with Unit Roots

Asymptotics. Hypothesis Testing UMP. Asymptotic Tests and p-values

APPM 5720 Solutions to Problem Set Five

CSE 202: Design and Analysis of Algorithms Lecture 16

MCR3U FINAL EXAM REVIEW (JANUARY 2015)

Answers to QUIZ

Homework 2 Solutions

The Eigen Function of Linear Systems

Unbiased Estimation. February 7-12, 2008

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

S n. = n. Sum of first n terms of an A. P is

Advanced Integration Techniques: Integration by Parts We may differentiate the product of two functions by using the product rule:

L-functions and Class Numbers

Last Lecture. Biostatistics Statistical Inference Lecture 16 Evaluation of Bayes Estimator. Recap - Example. Recap - Bayes Estimator

On triangular billiards

Modal Analysis of a Tight String

F.Y. Diploma : Sem. II [AE/CH/FG/ME/PT/PG] Applied Mathematics

Lecture 6 Ecient estimators. Rao-Cramer bound.

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Stability. Outline Stability Sab Stability of Digital Systems. Stability for Continuous-time Systems. system is its stability:

ON THE AUTOREGRESSIVE FRACTIONAL UNIT INTEGRATED MOVING AVERAGE (ARFUIMA) PROCESS

Fermat Numbers in Multinomial Coefficients

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

A Note on Random k-sat for Moderately Growing k

Transcription:

Exercise 7 / page 356 Noe ha X i are ii from Beroulli(θ where 0 θ a Meho of momes: Sice here is oly oe parameer o be esimae we ee oly oe equaio where we equae he rs sample mome wih he rs populaio mome, m µ, a we solve i erms of θ m X a µ E(X, where E(X xp (X x P (X θ x0 Hece, he meho of momes esimaor of θ is ˆθ MM X Noe hf more ha half ou of he observaios are equal o, he X > a ˆθ MM ges a value which is ou of he rage of parameer θ Recall ha 0 θ MLE of θ: We have o eermie he likelihoo fucio rs: L(θ x P (X i x i θ θ x i ( θ x i θ x i ( θ x i a for θ (0, ] he log-likelihoo is of he form: log L(θ x x i log θ + ( x i log( θ (Noe ha θ [0, ] bu he log-likelihoo is o ee for θ 0 So, we may ee o suy he fucio L(θ x a his specic poi separaely The rs erivaive of log L(θ x wih respec o θ is log L(θ x θ θ x i θ ( θ x i θ( x i θ( θ ( x θ θ( θ For θ log L(θ x 0 we ge θ x The seco erivaive of log L(θ x wih respec o θ is ( log L(θ x [ ] x i θ θ θ θ θ θ ( θ < 0 Therefore, θ x is global mamum poi (a we o o ee o suy L(θ x a θ 0 Noe ha θ [0, ], ie we ee o he global mamum of log L(θ x wihi his rage of θ values Therefore, we coclue ha: ˆθ MLE { X if X if X > or ˆθ MLE mi { X, }

(Observe ha for θ < x, log L(θ x > 0, ie log L(θ x a subsequely L(θ x are icreasig, θ a for θ > x, log L(θ x < 0, ie log L(θ x a subsequely L(θ x are ecreasig b θ MSE(ˆθ MM [ Bias(ˆθ MM ] + V ar(ˆθmm, where Bias(ˆθ MM E(ˆθ MM θ a ˆθ MM X Moreover, E(ˆθ MM E( X E(X θ Bias(ˆθ MM 0, a sice X i are ii from Beroulli(θ Thus, V ar( X V ar(x θ ( θ MSE(ˆθ MM θ ( θ For MSE(ˆθ MLE i is preferable o work wih he eiio of MSE, ie: ( { } ( { } Xi MSE(ˆθ MLE E (ˆθMLE θ E mi X, θ E mi, θ Sice X i Beroulli(θ, 0 θ, he X i Bi(, θ wih 0 θ Le T { } X i The mi Xi, is a fucio of he raom variable T where T Bi(, θ Thus: ( { } T ( { } ( MSE(ˆθ MLE E mi, θ mi, θ θ ( θ 0 Noe hf is eve he mi {, } { for 0,, / for (/ +,, a if is o he mi {, } Le [/] { / ( / { if eve if o, he for 0,, ( / for ( + /,, { } { mi, for 0,, [/] for [/] +,, Thus, he expecaio above ca be wrie as MSE(ˆθ MLE [/] 0 ( ( θ θ ( θ + [/]+ ( θ ( θ ( θ

c We will base our aswer o he quesio which of he wo esimaors is prefere o he MSE crierio To compare he MSE of he wo esimaors i woul help if we express MSE(ˆθ MM i a similar fashio as we i for MSE(ˆθ MLE Thus, MSE(ˆθ MM E ( X θ E ( T θ [/] 0 ( ( θ 0 θ ( θ + ( ( θ [/]+ θ ( θ ( ( θ θ ( θ Sice > for [/] +,,, he ( θ < ( θ for all [/] +,, a [/]+ ( θ ( θ ( θ < [/]+ ( ( θ θ ( θ for all θ (0, ] Hece, MSE(ˆθ MLE < MSE(ˆθ MM for all θ (0, ] If θ 0 he MSE(ˆθ MLE MSE(ˆθ MM Thus, ˆθ MLE is preferre base o he MSE crierio Oe coul hik as follows as well: Sice ˆθ MM X a ˆθ MLE mi { X, } he ˆθMLE ˆθ MM Therefore, (ˆθMLE θ (ˆθMM θ a E (ˆθMLE θ E (ˆθMM θ (see Theorem 5, page 57, ie MSE(ˆθMLE MSE(ˆθ MM Thus, we woul prefer ˆθ MLE agai o be o he safe sie bu i woul o be clear i which cases equaliy hols where boh esimaors perform equally well base o he MSE crierio Exercise 738 / page 356 To aswer he quesio of he exercise we use Corollary 735 a I orer o use he corollary we have o check rs wheher f(x θ saises he coiio of Cramér-Rao iequaliy, ie wheher he ierchage of iegraio a iereiaio is allowe Noice ha f(x θ θx θ is a expoeial family isribuio (i ca be wrie i he form f(x θ h(xc(θ exp ( k w k(θ k (x Thus, he requireme of Cramér-Rao iequaliy hols So, ow we ca apply he corollary The likelihoo fucio is of he form: L(θ x f(x θ f(x i θ θx θ i ( θ θ x i a The log L(θ x log θ + (θ log x i ( θ log L(θ x θ + log x i log x i θ, 3

ie log L(θ x ca be facorize i he form log L(θ x a(θ [W (X τ(θ] Hece, accorig o Corollary 735, log X i θ θ is a bes ubiase esimaor of θ b f(x θ log θ θ θx is a expoeial family isribuio a herefore he requireme of Cramér- Rao iequaliy hols Thus, we ca apply Corollary 735 The likelihoo fucio is of he form: a The L(θ x f(x θ f(x i θ ( log θ θ θx i log θ θ x i θ ( log L(θ x log (log θ log (θ + x i log θ log L(θ x θ θ log θ θ + θ x i θ ( ( θ x θ log θ ie log L(θ x ca be facorize i he form log L(θ x a(θ [W (X τ(θ] Hece, accorig o Corollary 735, X is a bes ubiase esimaor of θ θ θ θ log θ, Exercise 740/ p 36 Oe soluio coul be usig Cramér-Rao iequaliy & Corollary 730 & Lemma 73: Firs of all, we check ha X is a ubiase esimaor of p Sice X i are ii from Beroulli(p, i,,, we ge E( X E(X p We also compue he variace of X i orer o compare i wih he Cramér-Rao bou, V ar( X V ar(x p( p Sice Beroulli(p is a expoeial family isribuio, he ierchage of iegraio a iffereiaio is allowe a V ar( X <, so he requiremes of Cramér-Rao iequaliy hol Thus, ( p E( X V ar( X [ ( ], E log f(x p pθ [ ( where E( X ] [ ( ] p a E log f(x p E log f(x p sice X p p p p i are ii (see [ ( ] ( Corollary 730 Furhermore, E log f(x p E log f(x p because Beroulli(p p p is a expoeial family isribuio a Lemma 73 hols To compue he eomiaor: f(x p p x ( p x log f(x p x log p + ( x log( p p log f(x p x ( x p p p log f(x p x ( x p ( p 4

( ( E log f(x p E X ( X E(X + E(X p p ( p p ( p p + p p( p Thus, he Cramér-Rao lower bou is ( E( X p [ ( ] E log f(x p pθ ( E log f(x p p p( p p( p which is acually equal o V ar( X Hece, X is a bes ubiase esimaor of p Base o Theorem 739 we coclue h is he bes ubiase esimaor of p Aleraively, we coul use Corollary 735: X i are ii from Beroulli(p which saises he coiio of he Cramér-Rao Theorem sice Beroulli(p is a expoeial family isribuio The likelihoo fucio is of he form: L(p x f(x p f(x i p p x i ( p x i p x i ( p x i a he log-likelihoo is log L(p x ( log p + ( x i log( p a he rs erivaive of he log-likelihoo wih respec o θ is log L(p x p p x i p p p( p ( x p p( p Sice log L(p x ca be facorize i he form log L(p x a(p [W (x τ(p], accorig o p p Corollary 735, X is a ubiase esimaor of p a i aais he Cramér-Rao Lower bou Hece, X is a bes ubiase esimaor of p a base o Theorem 739 i is he bes ubiase esimaor of p Exercise 74/ p 363 X,, X is a raom sample where E(X i µ a V ar(x i σ a ( E X i E(X i µ µ If, he E ( X i µ, ie a ubiase esimaor of µ b Firsly, we eermie he variace of a esimaor of ype X i where : ( V ar X i a i V ar(x i a i σ σ a i 5

To he esimaor wih he miimum variace we ca see he above as a fucio of 's, i,, ha we wa o miimize uer he coiio Equivalely, we ca miimize he fucio: ( f(a, a,, a a i + (Observe ha he laer is a fucio of 's, i,,, a a is expresse as a base o he resricio ha The parial erivaives of f(a, a,, a are f(a, a,, a a a ( a f(a, a,, a a ( ( f(a, a,, a a a ie f(a, a,, a ( for i,, Seig all parial erivaives equal o 0 we ge a a a ie, for i,, Sice he righ ha sie i all equaios is he same, i hols ha se a a a k Replacig i ay of he above equaios he 's wih k we ge: ( k k k ( k k k The seco parial erivaives are f(a, a,, a 4 > 0, i,, 6

a j f(a, a,, a > 0, i, j,, a i j Thus, f(a, a,, a reaches is global miimum a he poi (a, a,, a (,,, Sice a, a his poi a akes he value: a Hece, i he case of X i variables wih he same mea a variace, amog all liear ubiase esimaors of ype X i for µ, he oe wih all 's equal o /, ie X, has he miimum variace Is variace is V ar( X σ V ar(x 7