Independence of Sample mean and Sample range of Normal Distribution 
Let $X_1,\dots,X_n$ be i.i.d. random variables with $X_1 \sim N(\mu,\sigma^2)$. Let $\bar X =\sum_{i=1}^n X_i/n$ and $R = X_{(n)}-X_{(1)}$, where $X_{(i)}$ is the $i$ the order statistic. Show that $\bar X$ and $R$ are independently distributed.

I know that sample mean and sample variance of normal distribution are independent. But this result state that sample mean and sample range of normal distribution are also independent. I know that $\bar X \sim N(\mu,\sigma^2/n)$. 
 A: This exercise is an immediate application of Basu's theorem.
To begin with, first note that $\bar{X}$ and $R$ are independent if and only if $\sigma^{-1}\bar{X}$ and $\sigma^{-1}R$ are independent, so we may assume $X_1, \ldots, X_n \text{ i.i.d.} \sim N(\mu, 1)$.
It is well-known that $\bar{X}$ is a sufficient and complete statistic for $\mu$, so according to Basu's theorem, to show that $\bar{X}$ and $R$ are independent, it remains to show $R$ is ancillary, i.e., $R$'s distribution is independent of $\mu$. This is easily seen by noticing
$$R = X_{(n)} - X_{(1)}= (X_{(n)} - \mu) - (X_{(1)} - \mu).$$
Clearly, the distribution of $(X_{(1)} - \mu, \ldots, X_{(n)} - \mu)$ is identical to the distribution of $(Z_{(1)}, \ldots, Z_{(n)})$, where $Z_1, \ldots, Z_n \text{ i.i.d.} \sim N(0, 1)$, which is distribution-constant. Consequently, the distribution of $R$ does not depend on $\mu$. This completes the proof.
A: A well-known property of the normal distribution is that the joint distribution of the differences $X_i-\overline X$, $i=1,2,\ldots,n-1$ is independent of the distribution of sample mean $\overline X$. This is also apparent from the independence of $\overline X=\frac{1}{n}\sum\limits_{i=1}^n X_i$ with the sample variance $\frac{1}{n-1}\sum\limits_{i=1}^n (X_i-\overline X)^2$, giving us $\operatorname{Cov}(X_i-\overline X,\overline X)=0$ for each $i$. The joint normality of $\overline X$ and $X_i-\overline X$ then implies their independence. Since the sample range $R=\max(X_i-\overline X)-\min(X_i-\overline X)$, $i=1,2,\ldots,n$ is a measurable function of $X_i-\overline X$ for each $i$, we can justify that $\overline X$ is independent of $R$.

An alternative proof of this result in line with @whuber's answer appears in this paper by J. Daly.
The argument provided in this article roughly goes along these lines:

We assume without loss of generality that $\mu=0$ and $\sigma^2=1$.
The joint characteristic function of $\overline X$ and the $\frac{n(n-1)}{2}$ differences $X_j-X_k$, $j<k$, is then
$$\varphi(t,t_{jk})=\frac{1}{(2\pi)^{n/2}}\int_{\mathbb{R^n}}\exp\left[-\frac{1}{2}\sum_{j=1}^nx_j^2+i\frac{t}{n}\sum_{j=1}^nx_j+i\sum_{1\le j<k\le n}t_{jk}(x_j-x_k)\right]\mathrm{d}\mathbf x$$
Completing the square in the exponent and further simplification leads to
$$\varphi(t,t_{jk})=\exp\left[-\frac{1}{2}\sum_{j=1}^n\left(\frac{t}{n}+\sum_{k=1}^n\left(t_{jk}-t_{kj}\right)\right)^2\right]$$
, which factors into the marginal characteristic functions
$$\varphi(t)\varphi(t_{jk})=\exp\left(-\frac{t^2}{2n}\right).\exp\left[-\frac{1}{2}\sum_{j=1}^n\left(\sum_{k=1}^n\left(t_{jk}-t_{kj}\right)\right)^2\right]$$
Thus the differences $X_j-X_k$ are jointly independent of $\overline X$.
Since the sample range $R=\max|X_j-X_k|$ is a measurable function of these differences, it follows that $\overline X$ and $R$ are independently distributed.

A: This is evidently a self-study question, so I do not intend to deprive you of the satisfaction of developing your own answer.  Moreover, I'm sure there are many possible solutions.  But for guidance, consider these observations:


*

*When a random variable $X$ is independent of other random variables $Y_1, \ldots, Y_m$, then $X$ is independent of any function of them, $f(Y_1, \ldots, Y_m)$.  (See Functions of Independent Random Variables for more about this.)

*Because the $X_i$ are jointly Normal, $X_1 + \cdots + X_n$ is independent of all the differences $Y_{ij} = X_i - X_j$, since their covariances are zero.
Because the range can be expressed as
$$X_{(n)} - X_{(1)} = \max_{i,j}(|X_i - X_j|) = \max_{i,j}(|Y_{ij}|)$$
you can exploit (1) and (2) to finish the proof.

For more intuition, a quick simulation might be of some help.  The following shows the marginal and joint distribution of the mean and range in the case $n=3$, using $10,000$ independent datasets.  The joint distribution clearly is not bivariate Normal, so the temptation to prove independence by means of a zero correlation--although a good idea--is bound to fail.  However, a close analysis of these results ought to suggest that the conditional distribution of the range does not vary with the mean.  (The appearance of some variation at the right and left is due to the paucity of outcomes with such extreme means.)

Here is the R code that produced these figures. It is easily modified to vary $n$, vary the simulation size, and to analyze the simulation results more extensively.
n <- 3; n.sim <- 1e4
sim <- apply(matrix(rnorm(n * n.sim), n), 2, function(y) c(mean(y), diff(range(y))))
par(mfrow=c(1,3))
hist(sim[1,], xlab="Mean", main="Histogram of Means")
hist(sim[2,], xlab="Range", main="Histogram of Ranges")
plot(sim[1,], sim[2,], pch=16, col="#00000020", xlab="Mean", ylab="Range")

