Question

In: Statistics and Probability

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. Also, please submit the R code for numerically finding the MLE.

Solutions

Expert Solution

## Following is the R code for part (b)

rm(list=ls(all=TRUE))
n <- 6;
data <- c(1.636, 0.374, 0.534, 3.015, 0.932, 0.179);
log_likelihood <- function(lambda){
lik <- n*log(lambda) - lambda*sum(data);
return(-lik);
}
parm <- 1; # initial guess
# finding mle numerically using optim function in r
numerical_mle <- optim(parm, log_likelihood, method="Brent", lower=0, upper=10)$par
analytic_mle = sum(data)/n; # finding mle analytically
cat("\n Analytic MLE is ", analytic_mle, "\n");
cat("\n Numerical MLE is ", numerical_mle, "\n");


Related Solutions

1) Xi~Bernoulli(p), MLE and MOME of p 2) Xi~Exp(λ), MLE and MOME of λ 3) Xi~Normal(μ,...
1) Xi~Bernoulli(p), MLE and MOME of p 2) Xi~Exp(λ), MLE and MOME of λ 3) Xi~Normal(μ, σ2 ), MLE and MOME of μ, σ2
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, ...
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....
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...
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.)
Find the (probability generating function) p.g.f.’s of the following distributions:[3+3=6] •P(X=x) =(exp(−λ)λ^x)/((1−exp(−λ))x!) , for x= 1,2,3,...,...
Find the (probability generating function) p.g.f.’s of the following distributions:[3+3=6] •P(X=x) =(exp(−λ)λ^x)/((1−exp(−λ))x!) , for x= 1,2,3,..., and λ >0. •P(X=x) =((pq)^x)(1−q^(N+1))^−1,for x= 0,1,..., N; where 0 < p < 1, p+q= 1.
Given the cumulative distribution of an exponential random variable find: The probability density function Show that...
Given the cumulative distribution of an exponential random variable find: The probability density function Show that it is a valid probability function The moment generating function The Expected mean The variance
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 X1,···, Xn ∼ Exp(λ) are independent. What is the distribution of X1/S where S =...
Suppose X1,···, Xn ∼ Exp(λ) are independent. What is the distribution of X1/S where S = X1+X2+···+Xn? Please show me how to do this without using the property of chi-squared dist.
Suppose X1,···, Xn ∼ Exp(λ) are independent. What is the distribution of X1/S where S =...
Suppose X1,···, Xn ∼ Exp(λ) are independent. What is the distribution of X1/S where S = X1+X2+···+Xn? Please show me how to do this without using the property of chi-squared dist.
ADVERTISEMENT
ADVERTISEMENT
ADVERTISEMENT