Question

In: Statistics and Probability

The exponential distribution with rate λ has mean μ = 1/λ. Thus the method of moments...

The exponential distribution with rate λ has mean μ = 1/λ. Thus the method of moments estimator of λ is 1/X. Use the following steps to verify that X is unbiased, but 1/X is biased.

a) Generate 10000 samples of size n = 5 from the standard exponential distribution (i.e. λ = 1) using rexp(50000) and arranging the 50000 random numbers in a matrix with 5 rows.

b) Use the apply() function to compute the 10000 sample means and store them in the object means. The 10000 estimators of λ can be stored in the object lambdas by lambdas = 1/means

c) Compute the sample mean of the object means, and sample mean of the object lambdas. What can you say about the bias of X and of 1/X?

d) Repeat with a sample of size n = 10, using rexp(100000), and report your estimate of the bias of 1/X. Has the bias decreased?

Solutions

Expert Solution

 rm(list=ls(all=TRUE)) > #a) Given information > n=5; > lambda=1 ; > #Random Number Generation > samples=rexp(50000,lambda) > x=matrix(samples,nrow=5,ncol=10000) > dim(x) [1] 5 10000 > ##b) > means=apply(x,2,mean) > lambdas=1/means > ##c) > smeans=mean(means) > cat("Sample mean of object means=",smeans) Sample mean of object means= 1.010265 > slambdas=mean(lambdas) > cat("Sample mean of object lambdas=",slambdas) Sample mean of object lambdas= 1.234539 > cat("Bias of X is=",abs(lambda-smeans)) Bias of X is= 0.01026538 > cat("Bias of 1/X is=",abs(lambda-slambdas)) Bias of 1/X is= 0.2345389 > cat("Bias of X which is insignificant less than bias of 1/X, so that X is good estimator than 1/X, And hence X is Undbiased Estimator and 1/X is biased") Bias of X which is insignificant less than bias of 1/X, so that X is good estimator than 1/X, And hence X is Undbiased Estimator and 1/X is biased > ##d ) Procedure for n=10 > n=10; > lambda=1 > #Random Number Generation > samples=rexp(100000,lambda) > x=matrix(samples,nrow=10,ncol=10000) > dim(x) [1] 10 10000 > means=apply(x,2,mean) > lambdas=1/means > smeans=mean(means) > cat("Sample mean of object means=",smeans) Sample mean of object means= 0.9975897 > slambdas=mean(lambdas) > cat("Sample mean of object lambdas=",slambdas) Sample mean of object lambdas= 1.114355 > cat("Bias of X is=",abs(lambda-smeans)) Bias of X is= 0.002410267 > cat("Bias of 1/X is=",abs(lambda-slambdas)) Bias of 1/X is= 0.1143547 > cat("When n=5 the bias of 1/X is greater than the bias of 1/X when n=10,So we conclude that as sample size(n) increases bias of estimator decreases.") When n=5 the bias of 1/X is greater than the bias of 1/X when n=10,So we conclude that as sample size(n) increases bias of estimator decreases.

Related Solutions

The exponential distribution with rateλhas meanμ= 1/λ. Thus the method of moments estimator of λ is...
The exponential distribution with rateλhas meanμ= 1/λ. Thus the method of moments estimator of λ is 1/X. Use the following steps to verify that X is unbiased, but 1/X is biased. USING R: a) Generate 10000 samples of size n= 5 from the standard exponential distribution (i.e.λ= 1) using rexp(50000) and arranging the 50000 random numbers in a matrix with 5 rows. b) Use the apply() function to compute the 10000 sample means and store them in the object means....
If Y is a random variable with exponential distribution with mean 1/λ, show that E (Y...
If Y is a random variable with exponential distribution with mean 1/λ, show that E (Y ^ k) = k! / (λ^k), k = 1,2, ...
For a random sample of sizenfrom an Exponential distribution with rate parameter λ (so that the...
For a random sample of sizenfrom an Exponential distribution with rate parameter λ (so that the density is fY(y) =λe−λy), derive the maximum likelihood estimator, the methods of moments estimator, and the Bayes estimator (that is, the posterior mean) using a prior proportional to λe−λ, for λ >0. (Hint: the posterior distribution will be a Gamma.)
Consider a two-server queue with Exponential arrival rate λ. Suppose servers 1 and 2 have exponential...
Consider a two-server queue with Exponential arrival rate λ. Suppose servers 1 and 2 have exponential rates μ1 and μ2, with μ1 > μ2. If server 1 becomes idle, then the customer being served by server 2 switches to server 1. a) Identify a condition on λ,μ1,μ2 for this system to be stable, i.e., the queue does not grow indefinitely long. b) Under that condition, and the long-run proportion of time that server 2 is busy.
Consider a two-server queue with Exponential arrival rate λ. Suppose servers 1 and 2 have exponential...
Consider a two-server queue with Exponential arrival rate λ. Suppose servers 1 and 2 have exponential rates µ1 and µ2, with µ1 > µ2. If server 1 becomes idle, then the customer being served by server 2 switches to server 1. a) Identify a condition on λ, µ1, µ2 for this system to be stable, i. e., the queue does not grow infinitely long. b) Under that condition, find the long-run proportion of time that server 2 is busy.
The number of imperfections in an object has a Poisson distribution with a mean λ =...
The number of imperfections in an object has a Poisson distribution with a mean λ = 8.5. If the number of imperfections is 4 or less, the object is called "top quality." If the number of imperfections is between 5 and 8 inclusive, the object is called "good quality." If the number of imperfections is between 9 and 12 inclusive, the object is called "normal quality." If the number of imperfections is 13 or more, the object is called "bad...
X and Y are independent Exponential random variables with mean=4, λ = 1/2. 1) Find the...
X and Y are independent Exponential random variables with mean=4, λ = 1/2. 1) Find the joint CDF of the random variables X, Y and  Find the probability that 4X > Y . 2) Find the expected value of X^3 + X*Y .
Exponential distribution with λ = 0.5 Apply Chi-Square test to see whether the your hypothesized distribution...
Exponential distribution with λ = 0.5 Apply Chi-Square test to see whether the your hypothesized distribution fits your data. using excel
Suppose the time to failure follows an exponential distribution with parameter λ. Based on type II...
Suppose the time to failure follows an exponential distribution with parameter λ. Based on type II censoring determine the likelihood and the maximum likelihood estimator for λ. Assume that n = 3 and the experiment ends after 2 failures.
a) Find the analytic MLE formula for exponential distribution exp(λ). Show that MLE is the same...
a) Find the analytic MLE formula for exponential distribution exp(λ). Show that MLE is the same as MoM estimator here. (b) A random sample of size 6 from the exp(λ) distribution results in observations: 1.636, 0.374, 0.534, 3.015, 0.932, 0.179. Find the MLE on this data set in two ways: by numerical optimization of the likelihood and by the analytic formula. For (b): please give both values from the analytic MLE formula and numerical MLE solution on this data set....
ADVERTISEMENT
ADVERTISEMENT
ADVERTISEMENT