11
$\begingroup$

I am reading Ruppert's Statistics and Data Analysis for Financial Engineering, which contains the following theorem:

Let $Y_1$, $...$ , $Y_n$ be an i.i.d. sample with a CDF $F$. Suppose that $F\;$ has a density $f\;$ that is continuous and positive at $F^{-1}(q)$, 0 < $q$ < 1. Then for large $n\;$, the $q^{th}$ sample quantile is approximately normally distributed with mean equal to the population quantile $F^{-1}(q)$ and variance equal to:

$$ \frac{q(1-q)}{n \left[f\;\left(F^{-1}(q)\right)\right]^2} $$

However, consider a Normal density with mean $\mu$ and variance $\sigma^2$. The distribution of the mean of $n\;$ samples will have mean $\mu$ and variance $\sigma^2/n$. If I am understanding it correctly, for q = 0.5 the above formula gives the variance as $$ \frac{\pi\sigma^2}{2n} $$

What am I missing here? Any insights are appreciated.

$\endgroup$
4
  • 5
    $\begingroup$ Although the population mean may equal the population median, the sample mean and sample median are not always the same and do not have the same sampling distributions. $\endgroup$
    – whuber
    Commented Aug 27, 2011 at 14:53
  • $\begingroup$ Actually it's just the opposite: var(median) > var(mean). $\endgroup$
    – Ringold
    Commented Aug 27, 2011 at 19:31
  • $\begingroup$ you're right; I deleted my previous comment $\endgroup$ Commented Aug 27, 2011 at 20:31
  • $\begingroup$ If n is somewhat smaller like 10 runs, then do these theorems apply ? $\endgroup$
    – SriK
    Commented Jun 25, 2018 at 18:03

2 Answers 2

4
$\begingroup$

As @whuber said in a comment, the formula you cite for the variance of a sample quantile is correct, as is that you give for the sample median in the normal case. The sample mean and sample median have different distributions, and in fact the median has higher variance in the normal case.

In R, you could estimate the SD of the sample median for the normal case as follows:

n <- 100; n.rep <- 1000
x <- matrix(rnorm(n*n.rep), ncol=n)
sd(apply(x, 1, median)) # 0.126
$\endgroup$
1
$\begingroup$

For $n=1$ the variance is of course 1. For $n=3$ we can show $1-\sqrt{3}/\pi$. For larger $n$ your estimate $\pi/(2n)$ seems very good from my experiments: variance

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.