Problems

Filters
Clear Filters

2 problems found

2004 Paper 3 Q14
D: 1700.0 B: 1488.4

In this question, \(\Phi(z)\) is the cumulative distribution function of a standard normal random variable. A random variable is known to have a Normal distribution with mean \(\mu\) and standard deviation either \(\sigma_0\) or \(\sigma_1\), where \(\sigma_0 < \sigma_1\,\). The mean, \(\overline{X}\), of a random sample of \(n\) values of \(X\) is to be used to test the hypothesis \(\mathrm{H}_0: \sigma = \sigma_0\) against the alternative \(\mathrm{H}_1: \sigma = \sigma_1\,\). Explain carefully why it is appropriate to use a two sided test of the form: accept \(\mathrm{H}_0\) if \(\mu - c < \overline{X} < \mu+c\,\), otherwise accept \(\mathrm{H}_1\). Given that the probability of accepting \(\mathrm{H}_1\) when \(\mathrm{H}_0\) is true is \(\alpha\), determine \(c\) in terms of \(n\), \(\sigma_0\) and \(z_{\alpha}\), where \(z_\alpha \) is defined by \(\displaystyle\Phi(z_{\alpha}) = 1 - \tfrac{1}{2}\alpha\). The probability of accepting \(\mathrm{H}_0\) when \(\mathrm{H}_1\) is true is denoted by \(\beta\). Show that \(\beta\) is independent of \(n\). Given that \(\Phi(1.960)\approx 0.975\) and that \(\Phi(0.063) \approx 0.525\,\), determine, approximately, the minimum value of \(\displaystyle \frac{\sigma_1}{\sigma_0}\) if \(\alpha\) and \(\beta\) are both to be less than \(0.05\,\).


Solution: If \(\sigma\) is smaller we should expect our sample to have a mean closer to the true mean. Therefore we should use a two sided test which accepts \(\mathrm{H}_0\) if the mean is very close to the true mean. Suppose \(\textrm{H}_0\) is true, ie \(\sigma = \sigma_0\), then note that \(X \sim N(\mu, \frac{\sigma_0^2}{n})\) \begin{align*} && 1-\alpha &= \mathbb{P}(\mu - c < X < \mu + c) \\ &&&= \mathbb{P}(\mu - c < \frac{\sigma_0}{\sqrt{n}} Z + \mu < \mu + c) \\ &&&= \mathbb{P}(- \frac{c\sqrt{n}}{\sigma_0} < Z<\frac{\sqrt{n}c}{\sigma_0}) \\ &&&= \mathbb{P}(Z<\frac{\sqrt{n}c}{\sigma_0}) -\mathbb{P}( Z<-\frac{\sqrt{n}c}{\sigma_0}) \\ &&&= \mathbb{P}(Z<\frac{\sqrt{n}c}{\sigma_0}) -(1-\mathbb{P}( Z<\frac{\sqrt{n}c}{\sigma_0})) \\ &&&= 2\mathbb{P}(Z<\frac{\sqrt{n}c}{\sigma_0})-1 \\ \Rightarrow && \Phi(\frac{\sqrt{n}c}{\sigma_0})&=1 - \tfrac12 \alpha \\ \Rightarrow && \frac{\sqrt{n}c}{\sigma_0} &= z_{\alpha} \\ && c &= \frac{\sigma_0 z_{\alpha}}{\sqrt{n}} \end{align*} Under \(\mathrm{H}_1\), \(\sigma = \sigma_1\) so \begin{align*} && \beta &= \mathbb{P}(\mu - c < X < \mu + c) \\ &&&= \mathbb{P}(-\frac{c\sqrt{n}}{\sigma_1} < Z < \frac{\sqrt{n}c}{\sigma_1}) \\ &&&= \mathbb{P}(-\frac{\sigma_0}{\sigma_1} z_{\alpha}< Z < \frac{\sigma_0}{\sigma_1} z_{\alpha}) \\ &&&= 2\Phi(\frac{\sigma_0}{\sigma_1} z_{\alpha})-1 \end{align*} which does not depend on \(n\). Suppose both \(\alpha<0.05\) and \(\beta<0.05\), then \(z_{\alpha} > 1.96\) and \(\Phi(\frac{\sigma_0}{\sigma_1}1.96)<0.525 \Rightarrow \frac{\sigma_0}{\sigma_1}1.96 < 0.063 \Rightarrow \frac{\sigma_1}{\sigma_0} > \frac{1.96}{0.063} = 31.1 \) so the ratio of variances needs to be larger than \(31.1\).

1995 Paper 2 Q14
D: 1600.0 B: 1500.0

Suppose \(X\) is a random variable with probability density \[ \mathrm{f}(x)=Ax^{2}\exp(-x^{2}/2) \] for \(-\infty < x < \infty.\) Find \(A\). You belong to a group of scientists who believe that the outcome of a certain experiment is a random variable with the probability density just given, while other scientists believe that the probability density is the same except with different mean (i.e. the probability density is \(\mathrm{f}(x-\mu)\) with \(\mu\neq0\)). In each of the following two cases decide whether the result given would shake your faith in your hypothesis, and justify your answer.

  1. A single trial produces the result 87.3.
  2. 1000 independent trials produce results having a mean value \(0.23.\)
{[}Great weight will be placed on clear statements of your reasons and none on the mere repetition of standard tests, however sophisticated, if unsupported by argument. There are several possible approaches to this question. For some of them it is useful to know that if \(Z\) is normal with mean 0 and variance 1 then \(\mathrm{E}(Z^{4})=3.\){]}


Solution: Let \(Z \sim N(0,1)\), with a pdf of \(f(x) = \frac{1}{\sqrt{2\pi}} \exp(-x^2/2)\) \begin{align*} && 1 &= \int_{-\infty}^\infty Ax^2 \exp(-x^2/2) \d x \\ &&&= A\sqrt{2\pi} \int_{-\infty}^\infty x^2 \frac{1}{\sqrt{2\pi}} \exp(-x^2/2) \d x \\ &&&= A\sqrt{2\pi} \E[Z^2] = A\sqrt{2\pi} \\ \Rightarrow && A &= \frac{1}{\sqrt{2\pi}} \end{align*}

  1. The probability of seeing a result as extreme as \(87.3\) is \begin{align*} \mathbb{P}(X > 87.3) &= \frac{1}{\sqrt{2\pi}}\int_{87.3}^{\infty} x^2 \exp(-x^2/2) \d x \\ &= \left [ -\frac{1}{\sqrt{2\pi}}x \exp(-x^2/2)\right]_{87.3}^{\infty}+\int_{87.3}^{\infty}\frac{1}{\sqrt{2\pi}} \exp(-x^2/2) \d x \\ &\approx 0 +(1- \Phi(87.3)) \\ &\approx 0 \end{align*} It is very unlikely this data point has come from our distribution rather than one with a higher mean, therefore our faith is very shaken.
  2. If there are 1000 trials of this, we would expect the sample mean to be distributed according to the CLT. Each sample has mean \(0\) and variance \(\E[X^2] = \int_{-\infty}^\infty x^4 \frac{1}{\sqrt{2\pi}} \exp(-x^2/2) \d x = \E[Z^4] = 3\), therefore the sample mean is \(N(0, 3/1000)\). Therefore the probability of being \(0.23\) away is \begin{align*} && \mathbb{P}(S > 0.23) &= \mathbb{P}\left (Z > \frac{0.23}{\sqrt{3/1000}} \right) \\ &&&= \mathbb{P}\left (Z > \frac{0.23}{\sqrt{30}/100} \right) \\ &&&\approx \mathbb{P}\left (Z > \frac{0.23}{0.055} \right) \\ &&& \approx 0 \end{align*} again our faith should be shaken