2. [R] Estimation of the parameter for exponential distribution (Jan5,7)

2. [R] Estimation of the parameter for exponential distribution (Jan5,7)

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
# Estimation of the parameter for exponential distribution # E(X) = 1/lam lam = .5 n = 40 x = rexp(n, lam) # 2 of many possible estimators lam.hat.1 = 1/mean(x) lam.hat.2 = 1/sqrt(var(x)) lam.hat.1 lam.hat.2 # Which is better? What is their sampling distribution? # sampling distribution # numerical aproximation to obtain this based on Monte Carlo simulation experiment M = 1000 lh.1.rep = rep(0, M) lh.2.rep = rep(0, M) lambda = lam for(i.rep in 1:M){ x.sim = rexp(n, lambda) lh.1.rep[i.rep] = 1/mean(x.sim) lh.2.rep[i.rep] =1/sqrt(var(x.sim)) } # the simulated pairs (lam.hat.1 , lam.hat.2) allows us to obtain approximations to various properties # of the sampling distribution of (lam.hat.1 , lam.hat.2) # Below we explore some of these in order to have some understanding of which simple summary measures # are useful in order to compare these two estimators. # scatter plot of lam.hat.1 versus lam.hat.2 # sometimes lam.hat.1 < lam and lam.hat.2 > lam plot( lh.1.rep, lh.2.rep) # based on this simulated values of lam.hat.1 and lam.hat.2 we can approximate many properties of
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Image of page 2
This is the end of the preview. Sign up to access the rest of the document.

Page1 / 3

2. [R] Estimation of the parameter for exponential distribution (Jan5,7)

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online