Comprehensive Exercises on Parameter Estimation

Exercises

Exercise 1

Let the population be X∼N(μ,σ2)X \sim N(\mu, \sigma^2), with sample mean X‾\overline{X} and sample variance S2S^2. Write the unbiased estimators of μ\mu and σ2\sigma^2.

point estimationunbiased estimationnormal population

Approach: Use the sample mean and the sample variance as unbiased estimators.

Steps:

  1. Unbiased estimator of μ\mu: μ^=X‾=1n∑i=1nXi\hat{\mu} = \overline{X} = \frac{1}{n}\sum_{i=1}^n X_i
  2. Unbiased estimator of σ2\sigma^2: σ^2=S2=1n−1∑i=1n(Xi−X‾)2\hat{\sigma}^2 = S^2 = \frac{1}{n-1}\sum_{i=1}^n (X_i - \overline{X})^2

Answer: μ^=X‾\hat{\mu} = \overline{X}, σ^2=S2\hat{\sigma}^2 = S^2

Exercise 2

Estimate the parameter by the method of moments: given that the sample x1,x2,…,xnx_1, x_2, \dots, x_n comes from X∼Poisson(λ)X \sim Poisson(\lambda).

parameter estimationmethod of momentssample moments

Approach: Use the method of moments, equating the sample moments to the population moments.

Steps:

  1. First population moment: E(X)=λE(X) = \lambda
  2. First sample moment: 1n∑i=1nxi=x‾\frac{1}{n}\sum_{i=1}^n x_i = \overline{x}
  3. Set x‾=λ\overline{x} = \lambda
  4. Solving gives: λ^=x‾\hat{\lambda} = \overline{x}

Answer: λ^=x‾\hat{\lambda} = \overline{x}

Exercise 3

Estimate the parameters by maximum likelihood estimation: given that the sample x1,x2,…,xnx_1, x_2, \dots, x_n comes from X∼N(μ,σ2)X \sim N(\mu, \sigma^2).

parameter estimationmaximum likelihood estimationlikelihood function

Approach: Use maximum likelihood estimation: construct the likelihood function and differentiate.

Steps:

  1. Likelihood function: L(μ,σ2)=∏i=1n12πσe−(xi−μ)22σ2L(\mu, \sigma^2) = \prod_{i=1}^n \frac{1}{\sqrt{2\pi}\sigma} e^{-\frac{(x_i-\mu)^2}{2\sigma^2}}
  2. Log-likelihood function: ln⁡L(μ,σ2)=−n2ln⁡(2π)−n2ln⁡(σ2)−12σ2∑i=1n(xi−μ)2\ln L(\mu, \sigma^2) = -\frac{n}{2}\ln(2\pi) - \frac{n}{2}\ln(\sigma^2) - \frac{1}{2\sigma^2}\sum_{i=1}^n (x_i-\mu)^2
  3. Differentiate with respect to μ\mu: ∂ln⁡L∂μ=1σ2∑i=1n(xi−μ)=0\frac{\partial \ln L}{\partial \mu} = \frac{1}{\sigma^2}\sum_{i=1}^n (x_i-\mu) = 0
  4. Solving gives: μ^=x‾\hat{\mu} = \overline{x}
  5. Differentiate with respect to σ2\sigma^2: ∂ln⁡L∂σ2=−n2σ2+12σ4∑i=1n(xi−μ)2=0\frac{\partial \ln L}{\partial \sigma^2} = -\frac{n}{2\sigma^2} + \frac{1}{2\sigma^4}\sum_{i=1}^n (x_i-\mu)^2 = 0
  6. Solving gives: σ^2=1n∑i=1n(xi−x‾)2\hat{\sigma}^2 = \frac{1}{n}\sum_{i=1}^n (x_i-\overline{x})^2

Answer: μ^=x‾\hat{\mu} = \overline{x}, σ^2=1n∑i=1n(xi−x‾)2\hat{\sigma}^2 = \frac{1}{n}\sum_{i=1}^n (x_i-\overline{x})^2

Exercise 4

Write the 1−α1-\alpha confidence interval for the mean μ\mu of a normal population (with the variance known).

interval estimationpivotal quantityconfidence interval

Approach: Construct the confidence interval by the pivotal quantity method.

Steps:

  1. Pivotal quantity: Z=X‾−μσ/n∼N(0,1)Z = \frac{\overline{X} - \mu}{\sigma/\sqrt{n}} \sim N(0,1)
  2. P(−zα/2≤X‾−μσ/n≤zα/2)=1−αP(-z_{\alpha/2} \leq \frac{\overline{X} - \mu}{\sigma/\sqrt{n}} \leq z_{\alpha/2}) = 1 - \alpha
  3. Solving gives: μ∈[X‾−zα/2σn,X‾+zα/2σn]\mu \in [\overline{X} - z_{\alpha/2}\frac{\sigma}{\sqrt{n}}, \overline{X} + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}]

Answer: [X‾−zα/2σn,X‾+zα/2σn][\overline{X} - z_{\alpha/2}\frac{\sigma}{\sqrt{n}}, \overline{X} + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}]

Exercise 5

Consider the estimators θ^1,θ^2\hat{\theta}_1, \hat{\theta}_2. If Var(θ^1)<Var(θ^2)Var(\hat{\theta}_1)<Var(\hat{\theta}_2), which one is more efficient?

point estimationefficiencycomparison of estimators

Approach: Compare the variances of the two estimators.

Steps:

  1. If Var(θ^1)<Var(θ^2)Var(\hat{\theta}_1) < Var(\hat{\theta}_2)
  2. and both estimators are unbiased
  3. then θ^1\hat{\theta}_1 is more efficient than θ^2\hat{\theta}_2

Answer: θ^1\hat{\theta}_1 is more efficient

Exercise 6

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼Poisson(λ)X \sim Poisson(\lambda). Prove that the sample mean X‾\overline{X} is an unbiased estimator of λ\lambda.

mathematical expectationlinearity of expectationcomputation

Approach: Use the linearity of expectation.

Steps:

  1. E(X‾)=E(1n∑i=1nXi)=1n∑i=1nE(Xi)E(\overline{X}) = E(\frac{1}{n}\sum_{i=1}^n X_i) = \frac{1}{n}\sum_{i=1}^n E(X_i)
  2. Since Xi∼Poisson(λ)X_i \sim Poisson(\lambda), we have E(Xi)=λE(X_i) = \lambda
  3. E(X‾)=1n⋅nλ=λE(\overline{X}) = \frac{1}{n} \cdot n\lambda = \lambda
  4. Therefore X‾\overline{X} is an unbiased estimator of λ\lambda

Answer: the proof is complete

Exercise 7

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼N(μ,σ2)X \sim N(\mu, \sigma^2). Compare 1n∑i=1n(Xi−X‾)2\frac{1}{n}\sum_{i=1}^n (X_i - \overline{X})^2 and 1n−1∑i=1n(Xi−X‾)2\frac{1}{n-1}\sum_{i=1}^n (X_i - \overline{X})^2, and determine which one is the unbiased estimator of σ2\sigma^2.

point estimationunbiased estimationproperties of estimators

Approach: Compute the expectations of the two estimators.

Steps:

  1. Let S2=1n−1∑i=1n(Xi−X‾)2S^2 = \frac{1}{n-1}\sum_{i=1}^n (X_i - \overline{X})^2
  2. It is known that E(S2)=σ2E(S^2) = \sigma^2, that is, 1n−1∑i=1n(Xi−X‾)2\frac{1}{n-1}\sum_{i=1}^n (X_i - \overline{X})^2 is unbiased
  3. Let σ^2=1n∑i=1n(Xi−X‾)2\hat{\sigma}^2 = \frac{1}{n}\sum_{i=1}^n (X_i - \overline{X})^2
  4. E(σ^2)=E(n−1nS2)=n−1nσ2≠σ2E(\hat{\sigma}^2) = E(\frac{n-1}{n}S^2) = \frac{n-1}{n}\sigma^2 \neq \sigma^2
  5. Therefore 1n∑i=1n(Xi−X‾)2\frac{1}{n}\sum_{i=1}^n (X_i - \overline{X})^2 is not an unbiased estimator

Answer: 1n−1∑i=1n(Xi−X‾)2\frac{1}{n-1}\sum_{i=1}^n (X_i - \overline{X})^2 is the unbiased estimator

Exercise 8

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼U(0,θ)X \sim U(0, \theta). Estimate the parameter θ\theta by the method of moments.

parameter estimationmethod of momentssample moments

Approach: Use the method of moments, equating the sample moments to the population moments.

Steps:

  1. First population moment: E(X)=θ2E(X) = \frac{\theta}{2}
  2. First sample moment: X‾=1n∑i=1nXi\overline{X} = \frac{1}{n}\sum_{i=1}^n X_i
  3. Set X‾=θ2\overline{X} = \frac{\theta}{2}
  4. Solving gives: θ^=2X‾\hat{\theta} = 2\overline{X}

Answer: θ^=2X‾\hat{\theta} = 2\overline{X}

Exercise 9

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼Exp(λ)X \sim Exp(\lambda). Estimate the parameter λ\lambda by maximum likelihood estimation.

parameter estimationmaximum likelihood estimationlikelihood function

Approach: Use maximum likelihood estimation: construct the likelihood function and differentiate.

Steps:

  1. Likelihood function: L(λ)=∏i=1nλe−λXi=λne−λ∑i=1nXiL(\lambda) = \prod_{i=1}^n \lambda e^{-\lambda X_i} = \lambda^n e^{-\lambda \sum_{i=1}^n X_i}
  2. Log-likelihood function: ln⁡L(λ)=nln⁡λ−λ∑i=1nXi\ln L(\lambda) = n \ln \lambda - \lambda \sum_{i=1}^n X_i
  3. Differentiate: ∂ln⁡L∂λ=nλ−∑i=1nXi=0\frac{\partial \ln L}{\partial \lambda} = \frac{n}{\lambda} - \sum_{i=1}^n X_i = 0
  4. Solving gives: λ^=n∑i=1nXi=1X‾\hat{\lambda} = \frac{n}{\sum_{i=1}^n X_i} = \frac{1}{\overline{X}}

Answer: λ^=1X‾\hat{\lambda} = \frac{1}{\overline{X}}

Exercise 10

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼N(μ,σ2)X \sim N(\mu, \sigma^2) with σ2\sigma^2 unknown. Find the 95%95\% confidence interval for μ\mu.

interval estimationt distributionconfidence interval

Approach: Construct the confidence interval using the t distribution.

Steps:

  1. Pivotal quantity: T=X‾−μS/n∼t(n−1)T = \frac{\overline{X} - \mu}{S/\sqrt{n}} \sim t(n-1)
  2. P(−t0.025(n−1)≤X‾−μS/n≤t0.025(n−1))=0.95P(-t_{0.025}(n-1) \leq \frac{\overline{X} - \mu}{S/\sqrt{n}} \leq t_{0.025}(n-1)) = 0.95
  3. Solving gives: μ∈[X‾−t0.025(n−1)Sn,X‾+t0.025(n−1)Sn]\mu \in [\overline{X} - t_{0.025}(n-1)\frac{S}{\sqrt{n}}, \overline{X} + t_{0.025}(n-1)\frac{S}{\sqrt{n}}]

Answer: [X‾−t0.025(n−1)Sn,X‾+t0.025(n−1)Sn][\overline{X} - t_{0.025}(n-1)\frac{S}{\sqrt{n}}, \overline{X} + t_{0.025}(n-1)\frac{S}{\sqrt{n}}]

Exercise 11

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼N(μ,σ2)X \sim N(\mu, \sigma^2). Find the 1−α1-\alpha confidence interval for σ2\sigma^2.

interval estimationchi-square distributionconfidence interval

Approach: Construct the confidence interval using the chi-square distribution.

Steps:

  1. Pivotal quantity: χ2=(n−1)S2σ2∼χ2(n−1)\chi^2 = \frac{(n-1)S^2}{\sigma^2} \sim \chi^2(n-1)
  2. P(χ1−α/22(n−1)≤(n−1)S2σ2≤χα/22(n−1))=1−αP(\chi^2_{1-\alpha/2}(n-1) \leq \frac{(n-1)S^2}{\sigma^2} \leq \chi^2_{\alpha/2}(n-1)) = 1 - \alpha
  3. Solving gives: σ2∈[(n−1)S2χα/22(n−1),(n−1)S2χ1−α/22(n−1)]\sigma^2 \in [\frac{(n-1)S^2}{\chi^2_{\alpha/2}(n-1)}, \frac{(n-1)S^2}{\chi^2_{1-\alpha/2}(n-1)}]

Answer: [(n−1)S2χα/22(n−1),(n−1)S2χ1−α/22(n−1)][\frac{(n-1)S^2}{\chi^2_{\alpha/2}(n-1)}, \frac{(n-1)S^2}{\chi^2_{1-\alpha/2}(n-1)}]

Exercise 12

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼N(μ,σ2)X \sim N(\mu, \sigma^2) with σ2\sigma^2 known. Find the 90%90\% confidence interval for μ\mu.

interval estimationZ distributionconfidence interval

Approach: Construct the confidence interval using the normal distribution.

Steps:

  1. Pivotal quantity: Z=X‾−μσ/n∼N(0,1)Z = \frac{\overline{X} - \mu}{\sigma/\sqrt{n}} \sim N(0,1)
  2. P(−z0.05≤X‾−μσ/n≤z0.05)=0.90P(-z_{0.05} \leq \frac{\overline{X} - \mu}{\sigma/\sqrt{n}} \leq z_{0.05}) = 0.90
  3. Solving gives: μ∈[X‾−z0.05σn,X‾+z0.05σn]\mu \in [\overline{X} - z_{0.05}\frac{\sigma}{\sqrt{n}}, \overline{X} + z_{0.05}\frac{\sigma}{\sqrt{n}}]

Answer: [X‾−z0.05σn,X‾+z0.05σn][\overline{X} - z_{0.05}\frac{\sigma}{\sqrt{n}}, \overline{X} + z_{0.05}\frac{\sigma}{\sqrt{n}}]

Exercise 13

Compare the advantages and disadvantages of the method of moments and maximum likelihood estimation.

parameter estimationmethod of momentsmaximum likelihood estimation

Approach: Compare them in terms of computational complexity, estimation performance, applicability, and so on.

Steps:

  1. Computational complexity:

    • Method of moments: simple to compute, only the sample moments need to be computed
    • Maximum likelihood estimation: relatively more complex to compute, requiring differentiation and solving equations
  2. Estimation performance:

    • Method of moments: may not yield the optimal estimator
    • Maximum likelihood estimation: yields the optimal estimator under regularity conditions
  3. Applicability:

    • Method of moments: applicable to all kinds of distributions, without needing to know the specific form
    • Maximum likelihood estimation: requires knowing the specific form of the distribution
  4. Asymptotic properties:

    • Method of moments: has consistency, but does not necessarily have asymptotic efficiency
    • Maximum likelihood estimation: has consistency, asymptotic normality, and asymptotic efficiency

Answer: the method of moments is simple to compute but its performance may not be as good as maximum likelihood estimation, whereas maximum likelihood estimation is more complex to compute but performs better

Exercise 14

Let the variances of two populations be σ12,σ22\sigma_1^2, \sigma_2^2. To test H0:σ12=σ22H_0: \sigma_1^2=\sigma_2^2, which test should be used?

hypothesis testingchoice of methodtest for equality of variances

Approach: Choose the appropriate test method according to the object being tested.

Steps:

  1. Test the ratio of the variances of two normal populations
  2. Use the F test
  3. Test statistic: F=S12S22∼F(n1−1,n2−1)F = \frac{S_1^2}{S_2^2} \sim F(n_1-1, n_2-1)

Answer: the F test

Exercise 15

Let X1,X2,…,XnX_1, X_2, \dots, X_n be a sample from the population X∼B(1,p)X \sim B(1,p). Estimate the parameter pp by maximum likelihood estimation.

parameter estimationmaximum likelihood estimationlikelihood function

Approach: Use maximum likelihood estimation: construct the likelihood function and differentiate.

Steps:

  1. Likelihood function: L(p)=∏i=1npXi(1−p)1−Xi=p∑i=1nXi(1−p)n−∑i=1nXiL(p) = \prod_{i=1}^n p^{X_i}(1-p)^{1-X_i} = p^{\sum_{i=1}^n X_i}(1-p)^{n-\sum_{i=1}^n X_i}
  2. Log-likelihood function: ln⁡L(p)=∑i=1nXiln⁡p+(n−∑i=1nXi)ln⁡(1−p)\ln L(p) = \sum_{i=1}^n X_i \ln p + (n-\sum_{i=1}^n X_i) \ln(1-p)
  3. Differentiate: ∂ln⁡L∂p=∑i=1nXip−n−∑i=1nXi1−p=0\frac{\partial \ln L}{\partial p} = \frac{\sum_{i=1}^n X_i}{p} - \frac{n-\sum_{i=1}^n X_i}{1-p} = 0
  4. Solving gives: p^=1n∑i=1nXi=X‾\hat{p} = \frac{1}{n}\sum_{i=1}^n X_i = \overline{X}

Answer: p^=X‾\hat{p} = \overline{X}

Found an issue with this content?

Tell us what could be improved.