%%% Title of object:	Best linear unbiased estimation in linear models
%%% Canonical Name:	BestLinearUnbiasedEstimatinInLinearModels
%%% Type:	Definition
%%% Created on:	2010-08-20 14:01:00
%%% Modified on:	2010-08-24 11:34:30
%%% Creator:	sjp
%%% Modifier:	jkimmel
%%% Author:	sjp
%%% 
%%% Classification:	msc:62J05
%%% Keywords:	Best linear unbiased, BLUE, BLUP, Gauss--Markov Theorem, Generalized inverse, Ordinary least squares, OLSE
%%% Preamble:
\documentclass[10pt]{article}
% this is the default PlanetMath preamble. as your knowledge
% of TeX increases, you will probably want to edit this, but
% it should be fine as is for beginners.

% almost certainly you want these
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{amsfonts}

% used for TeXing text within eps files
%\usepackage{psfrag}
% need this for including graphics (\includegraphics)
%\usepackage{graphicx}
% for neatly defining theorems and propositions
%\usepackage{amsthm}
% making logically defined graphics
%\usepackage{xypic}

% there are many more packages, add them here as you need them

\usepackage{amsthm}
\usepackage{mathrsfs} % CURLY C, N, M

% define commands here
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\newcommand*{\BETA}{\ensuremath{\boldsymbol{\beta}}} % bold beta
\newcommand*{\BETAH}{\ensuremath{\boldsymbol{\hat\beta}}} % bold beta hat
\newcommand*{\BETAT}{\ensuremath{\boldsymbol{\tilde\beta}}}
\newcommand{\BLUE}{\mbox{\small$\mathrm{BLUE}$}}
\newcommand{\BLUP}{\mbox{\small$\mathrm{BLUP}$}}
\newcommand*{\betat}{\tilde\beta}
\DeclareMathOperator{\C}{\ensuremath{\mathscr C}} % column space
\DeclareMathOperator{\cov}{cov} % covariance
\DeclareMathOperator{\EE}{E} % expectation
\DeclareMathOperator{\E}{E} % expectation
\newcommand*{\EPS}{\ensuremath{\boldsymbol{\varepsilon}}} % bold eps
\newcommand*{\GAMMA}{\ensuremath{\boldsymbol{\gamma}}} %
\newcommand*{\M}{\ensuremath{\mathscr M}} %
\newcommand*{\mx}[1]{\mathbf{#1}} % vectors and matrices
\DeclareMathOperator{\NS}{\ensuremath{\mathscr N}} % nullspace
\newcommand*{\OLSE}{\mbox{\small$\mathrm{OLSE}$}}
\DeclareMathOperator{\rank}{rank} % r (for rank)
\DeclareMathOperator{\tr}{trace} % r (for rank)
\newcommand*{\rz}{\ensuremath{\mathbb{R}}} % blackboard R
\newcommand*{\SIGMA}{\ensuremath{\boldsymbol{\Sigma}}} % bold cap Sigma
\DeclareMathOperator{\var}{var}

\newtheorem{theorem}{Theorem}
%%% Content:
\begin{document}
\title{Best Linear Unbiased Estimation in Linear Models}
\author{Simo Puntanen\thanks{Department of Mathematics and Statistics,
FI-33014 University of Tampere, Tampere, Finland.
Email: \texttt{simo.puntanen@uta.fi}} \\
University of Tampere, Finland
\and
George P. H. Styan\thanks{Department of Mathematics and Statistics,
McGill University, 805 ouest rue Sherbrooke
Street West, Montr\'{e}al (Qu\'{e}bec), Canada H3A 2K6.
Email: {\tt styan@math.mcgill.ca}}\\
McGill University, Montr\'{e}al, Canada
}
%\date{}
\maketitle
\noindent{\bf Keywords and Phrases:} Best linear unbiased, BLUE, BLUP, Gauss--Markov Theorem, Generalized inverse, Ordinary least squares, OLSE.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{Introduction}
In this article we consider the general linear model
(Gauss--Markov model)
\begin{equation*}
\mx y = \mx X \BETA + \EPS,
\quad \text{or in short\ }
\M = \{ \mx y, \, \mx X \BETA, \, \sigma^2 \mx V \},
\end{equation*}
where
$\mx X$ is a known $n\times p$ model matrix, the
vector $\mx y$ is an observable $n$-dimensional random vector,
$\BETA$
is a $p\times 1$ vector of unknown parameters, and
$\EPS$ is an unobservable vector of random errors
with expectation
$\EE(\EPS ) = \mx 0,$ and
covariance matrix
$\cov( \EPS) = \sigma^2 \mx V,$
where
$\sigma^2 >0$ is an unknown constant. The nonnegative
definite (possibly singular) matrix $\mx V $ is known. In our
considerations $\sigma ^2$ has no role and hence we may put
$\sigma^2=1.$


As regards the notation,
we will use the symbols
$\mx A',$
$\mx A^{-},$
$\mx A^{+},$
$\C(\mx A),$
$\C(\mx A)^{\bot},$
and
$\NS(\mx A)$
to denote, respectively,
the transpose,
a generalized inverse,
the Moore--Penrose inverse,
the column space,
the orthogonal complement of the column space,
and the null space,
of the matrix~$\mx A.$
By $(\mx A:\mx B)$ we denote the partitioned matrix with
$\mx A$ and $\mx B$ as submatrices.
By $\mx A^{\bot}$ we denote any matrix satisfying
$\C(\mx A^{\bot}) = \NS(\mx A') = \C(\mx A)^{\bot}.$
Furthermore, we will write
$\mx P_{\mx A} = \mx A\mx A^{+}
= \mx A(\mx A'\mx A)^{-}\mx A'$
to denote the orthogonal projector (with respect to the standard
inner product) onto
$\C(\mx A).$
In particular, we denote
$\mx{H} = \mx P_{\mx X}$ and $ \mx{M} = \mx I_n - \mx H.$
One choice for $\mx X^{\bot}$ is of course the projector
$\mx M.$



Let $\mx K' \BETA$ be a given vector of parametric functions specified
by $\mx K' \in \rz^{q\times p}.$
Our object is to find a (homogeneous) linear estimator $\mx A \mx y$
which would provide an unbiased and in some sense ``best'' estimator
for $\mx K' \BETA$ under the model $\M.$
However, not all parametric functions have linear unbiased
estimators; those which have are called \emph{estimable} parametric functions,
and then there exists a matrix $\mx A$ such
that
\begin{equation*}
\E(\mx{Ay}) = \mx{AX}\BETA = \mx K' \BETA
\quad \text{for all $\BETA \in \rz^p.$}
\end{equation*}
Hence
$\mx{K}' \BETA$ is estimable
if and only if
there exists a matrix $\mx A$ such that $\mx{K}' = \mx{A}\mx{X},$ i.e.,
$ \C(\mx K ) \subset \C(\mx X').$



The ordinary least squares estimator
of $\mx K' \BETA$ is defined
as
$ \OLSE(\mx K' \BETA) = \mx K' \BETAH, $
where $\BETAH$ is any solution to the normal equation
$\mx X' \mx X \BETAH = \mx X' \mx y$; hence
$\BETA = \BETAH$ minimizes $(\mx y - \mx X\BETA)' (\mx y - \mx X\BETA)$
and it can be expressed as $\BETAH = (\mx X' \mx X) ^{-}\mx X' \mx y,$
while $\mx X\BETAH = \mx H \mx y.$
Now the condition $\C(\mx K ) \subset \C(\mx X')$ guarantees that
$\mx K'\BETAH$ is unique, even though $\BETAH$ may not be unique.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{The Best Linear Unbiased Estimator (BLUE)}
%
The expectation $\mx X\BETA$ is trivially estimable
and $\mx{Gy}$ is unbiased for $\mx X\BETA$ whenever
$ \mx{G}\mx X = \mx{X}.$
An unbiased linear estimator $\mx{Gy}$
for $\mx X\BETA$ is defined to be
the {\it best} linear unbiased estimator,
$\BLUE$, for $\mx X\BETA$ under $\M$ if
\begin{equation*}
\cov( \mx{G} \mx y) \leq_\textup{L} \cov( \mx{L} \mx y) \quad
\text{for all } \mx{L} \colon
\mx{L}\mx X = \mx{X},
\end{equation*}
where ``$\leq_\text{L}$'' refers to the L\"{o}wner partial ordering.
In other words, $\mx{G} \mx y$ has the smallest covariance matrix
(in the L\"{o}wner sense) among all linear unbiased estimators.
We denote the $\BLUE$ of $\mx X\BETA$ as
$ \BLUE(\mx X\BETA) = \mx X \BETAT. $
If $\mx X$ has full column rank, then $\BETA$ is estimable
and an unbiased estimator $\mx A\mx y$ is the $\BLUE$ for $\BETA$ if
$ \mx{AVA}' \leq_\textup{L} \mx{BVB}'
$ for all $\mx{B}$ such that
$ \mx{BX} = \mx{I}_p. $
The L\"owner ordering is a very strong ordering implying for example
\begin{gather*}
\var(\betat_i) \le \var(\beta^{*}_i) \,, \quad i = 1,\dotsc,p , \\[\jot]
\tr [\cov(\BETAT)] \le \tr [\cov(\BETA^{*})] , \qquad
\det[\cov(\BETAT)] \le \det[\cov(\BETA^{*})],
\end{gather*}
for any linear unbiased estimator $\BETA^{*}$ of $\BETA$; here
$\var$ refers to the variance and
``det'' denotes
the determinant.

The following theorem gives the ``Fundamental $\BLUE$ equation'';
see, e.g.,
Rao (1967),
Zyskind (1967)
and Puntanen, Styan and Werner (2000).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{theorem}\label{thm: fundablue}%
Consider the general linear model $ \M =\{\mx y,\,\mx X\BETA,\,\mx V\}.$
Then
the estimator $\mx{Gy}$
is the $\BLUE$ for $\mx X\BETA$ if and only if $\mx G$
satisfies the equation
\begin{equation} \label{eq: 30jan09-fundablue}
\mx{G}(\mx{X} : \mx{V}\mx{X}^{\bot} ) = (\mx{X} : \mx{0}) .
\end{equation}
The corresponding
condition for $\mx{Ay}$ to be the $\BLUE$ of an estimable parametric
function $\mx{K}' \BETA$
is
$ \mx{A}(\mx{X} : \mx{V}\mx{X}^{\bot} ) = (\mx{K}' : \mx{0}). $
\end{theorem}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


It is sometimes convenient to express
\eqref{eq: 30jan09-fundablue}
in the following form, see
Rao (1971).
%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%

\begin{theorem}[Pandora's Box]
Consider the general linear model $ \M =\{\mx y,\,\mx X\BETA,\,\mx V\}.$
Then
the estimator $\mx{Gy}$
is the $\BLUE$ for $\mx X\BETA$ if and only if there exists
a matrix $\mx{L} \in \rz^{p \times n}$
so that $\mx G$ is a solution to
\begin{equation*}
\begin{pmatrix}
\mx V & \mx X \\
\mx X' & \mx 0
\end{pmatrix}
\begin{pmatrix}
\mx G' \\
\mx{L}
\end{pmatrix} =
\begin{pmatrix}
\mx 0 \\
\mx X'
\end{pmatrix}.
\end{equation*}
\end{theorem}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The equation \eqref{eq: 30jan09-fundablue} has a unique solution
for $\mx G$ if and only if $\C(\mx X : \mx V) = \rz^n.$
Notice that under $\M$ we assume that the observed value of
$\mx y$ belongs to the subspace $\C(\mx X : \mx V)$
with probability $1$; this is the consistency condition
of the linear model,
see, e.g.,
Baksalary, Rao and Markiewicz (1992).
The consistency condition means, for example, that whenever we have
some statements which involve the random vector~$\mx y$, these
statements need hold only for those values of $\mx y$ that belong
to $\C(\mx{X}:\mx{V}).$
The general solution for $\mx G$
can be expressed, for example, in
the following ways:
%%%%%%%%%%%%%%
\[
\mx G_1 = \mx{X}(\mx{X}'\mx{W}^{-}\mx{X})^{-}\mx{X}'\mx{W}^{-}
+ \mx F_{1}(\mx{I }_n - \mx W\mx W^{-} ) ,
\]
\[
\mx G_2 = \mx{H} - \mx{HVM}(\mx{MVM})^{-}\mx{M} + \mx F_{2}[\mx{I}_n -
\mx{MVM}( \mx{MVM} )^{-} ]\mx M ,
\]
where $\mx F_{1}$ and $\mx F_{2}$ are arbitrary
matrices,
$\mx {W}= \mx V + \mx X\mx U\mx X'$ and $\mx U$ is any arbitrary conformable
matrix such that $\C(\mx W) = \C(\mx X : \mx V).$
Notice that even though $\mx G$ may not be unique, the numerical value
of $\mx G\mx y$ is unique because $\mx y \in \C(\mx X : \mx V).$
If $\mx V$ is positive definite,
then
$\BLUE(\mx X\BETA) = \mx X(\mx X' \mx V^{-1} \mx X)^{-} \mx X' \mx V^{-1} \mx y.$
Clearly $\OLSE(\mx X\BETA) = \mx H\mx y$ is the $\BLUE$ under
$\{ \mx y, \, \mx X\BETA , \, \sigma^2\mx I \}.$
It is also worth noting that the matrix $\mx G$ satisfying
\eqref{eq: 30jan09-fundablue} can be interpreted as a
projector: it is a projector onto $\C(\mx X)$ along $\C(\mx V\mx X^{\bot}),$
see Rao (1974).


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{OLSE vs. BLUE}
Characterizing the equality of the Ordinary Least Squares Estimator
$(\OLSE)$ and the $\BLUE$ has received a lot
of attention in the literature,
since Anderson (1948),
but the major breakthroughs were made
by
Rao (1967) and
Zyskind (1967);
for a detailed review, see
Puntanen and Styan (1989).
For some further references from those years we may mention
Kruskal (1968),
Watson (1967),
and Zyskind and Martin (1969).


We present below six characterizations for the $\OLSE$ and
the $\BLUE$ to be equal (with probability $1$).

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{theorem}[$\OLSE$ vs. $\BLUE$]\label{cor: olseblue}%
Consider the general linear model $ \M =\{\mx y,\,\mx X\BETA,\,\mx V\}.$
Then $\OLSE(\mx{X}\BETA) = \BLUE(\mx{X}\BETA)$ if and only
if any one of the following six equivalent conditions holds.
$($Note: $\mx{V}$ may be replaced by its Moore--Penrose inverse
$\mx{V}^+$ and $\mx{H}$ and $\mx{M} = \mx I_n - \mx H$ may be
interchanged.$)$
\begin{align*}
&\textup{(1)} \; \mx{HV} = \mx{VH} , \\
&\textup{(2)} \; \mx{H}\mx{V}\mx{M} = \mx{0} \,, \\
&\textup{(3)} \; \C(\mx{V}\mx{H})\subset\C(\mx H) ,
\\
&\textup{(4)} \; \C(\mx{X}) \text{ has a basis comprising
$r= \rank(\mx X)$ orthonormal eigenvectors of $\mx V,$}
\\
&\textup{(5)}\; \mx{V} = \mx{HAH} + \mx{MBM}
\text{ for some $\mx A $ and $\mx B,$}
\\
&\textup{(6)}\; \mx{V} = \alpha\mx{I}_n +
\mx{HKH} + \mx{M}\mx L\mx M \text{ for some $\alpha \in \rz$, and $\mx K$ and
$\mx L.$}
\end{align*}
\end{theorem}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



Theorem \ref{cor: olseblue} shows at once that
under $\{ \mx y, \, \mx X\BETA, \, \mx I_n \}$ the $\OLSE$ of
$\mx X\BETA$ is trivially the $\BLUE$; this result is often called
the Gauss--Markov Theorem.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Two Linear Models}
Consider now two linear models
$ \M_{1} = \{ \mx y, \, \mx X\BETA, \, \mx V_1 \}$
and $
\; \M_{2} = \{ \mx y, \, \mx X\BETA, \, \mx V_2 \},$
which differ only in their covariance matrices.
For the proof of
the
following proposition and related discussion, see, e.g.,
Rao (1971, Th.~5.2, Th.~5.5),
and
Mitra and Moore (1973, Th.~3.3, Th.~4.1--4.2).




%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{theorem} \label{propo: 1111}
Consider the linear models
$ \M_1 = \{\mx y, \, \mx X\BETA, \, \mx V_1 \}$
and
$ \M_2 = \{ \mx y, \, \mx X\BETA, \, \mx V_2 \},$
and let the notation
$ \{\BLUE(\mx X\BETA \mid \M_1) \} \subset
\{\BLUE(\mx X\BETA \mid \M_2) \} $
mean that every representation of the $\BLUE$ for $\mx X\BETA$ under $\M_1$
remains the $\BLUE$ for $\mx X\BETA$ under $\M_2$. Then the following
statements are equivalent:
\begin{align*}
&\textup{(1)} \;
\{ \BLUE(\mx X\BETA \mid \M_1) \} \subset
\{ \BLUE(\mx X\BETA \mid \M_2) \} ,\\
&\textup{(2)} \;
\C(\mx V_2\mx X^{\bot}) \subset \C(\mx V_1 \mx X^\bot) ,\\
&\textup{(3)} \;
\mx V_2 = \mx V_1+ \mx{X} \mx N_1 \mx X' +
\mx V_1\mx M\mx N_2 \mx M \mx V_1,
\text{ for some $\mx N_1$ and $\mx N_2,$} \\
&\textup{(4)} \;
\mx V_2 =
\mx{X} \mx N_3 \mx X' + \mx V_1\mx M\mx N_4 \mx M \mx V_1,
\text{ for some $\mx N_3$ and $\mx N_4.$}
\end{align*}

\end{theorem}
%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%


Notice that obviously
\begin{align*}
\{ \BLUE(\mx X \BETA \mid \M_1) \}
=
\{ \BLUE(\mx X \BETA \mid \M_2) \}
\iff
\C(\mx V_2\mx X^{\bot}) = \C(\mx V_1 \mx X^\bot) .
\end{align*}
For the equality
between the
$\BLUE$s of $\mx X_1\BETA_1$
under two partitioned models, see
Haslett and Puntanen (2010a).



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Model with New Observations: Best Linear Unbiased Predictor (BLUP)}
Consider the model
$ \M = \{\mx y,\,\mx X\BETA,\,\mx V\},$
and
let $\mx y_f$
denote an $m\times 1$ unobservable random vector containing
\emph{new observations}. The new observations are assumed to follow
the linear model
$ \mx y_f = \mx X_f\BETA +\EPS_f ,$
where $\mx X_f$ is a known $m\times p$ model matrix associated with new
observations, $\BETA$ is the same vector of unknown parameters as
in $\M$, and $\EPS_f$ is an $m \times 1$ random
error vector associated with new observations.
Our goal is to predict the random vector $\mx y_f$
on the basis of $\mx y$.
The expectation and the covariance matrix are
\begin{equation*}
\E\begin{pmatrix}
\mx y \\
\mx y_f
\end{pmatrix} =
\begin{pmatrix}
\mx X\BETA \\
\mx X_f\BETA
\end{pmatrix} , \quad
\cov\begin{pmatrix}
\mx y \\
\mx y_f
\end{pmatrix}
=
\begin{pmatrix}
\mx V = \mx V_{11} & \mx{V}_{12} \\
\mx{V}_{21} & \mx V_{22}
\end{pmatrix} ,
\end{equation*}
which we may write as
\begin{equation*}
\M_f = \left \{
\begin{pmatrix}
\mx y \\
\mx y_f
\end{pmatrix},\,
\begin{pmatrix}
\mx X\BETA \\
\mx X _f\BETA
\end{pmatrix},\,
\begin{pmatrix}
\mx V & \mx{V}_{12} \\
\mx{V}_{21} & \mx V_{22}
\end{pmatrix} \right \}.
\end{equation*}

A linear predictor $\mx{Ay}$ is said to be unbiased for $\mx y_f$ if
$ \E(\mx{Ay}) = \E(\mx{y}_f) = \mx X_f\BETA$
for all $\BETA\in\rz^{p}.$
Then the random vector
$\mx y_f$ is said to be unbiasedly predictable.
Now an unbiased linear predictor $\mx{Ay}$ is the
best linear unbiased predictor, $\BLUP$, for $\mx y_f$
if the L\"owner ordering
\begin{equation*}
\cov(\mx{Ay}-\mx y_f) \leq_\textup{L} \cov(\mx{By}-\mx y_f)
\end{equation*}
holds for all $\mx B$ such that $\mx{By}$ is an unbiased linear
predictor for $\mx{y}_f$.

The following theorem characterizes the $\BLUP$;
see, e.g.,
Christensen (2002, p.~283),
and
Isotalo and Puntanen (2006, p.~1015).
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%


\begin{theorem}[Fundamental $\BLUP$ equation]\label{propo:BLUP-funda}
Consider the linear model
$\M_f,$ where
$\mx{X}_f\BETA$ is a given estimable parametric function.
Then the linear estimator $\mx{Ay}$ is the best
linear unbiased predictor \textup{(}$\BLUP$\textup{)} for $\mx y_f$ if
and only if $\mx{A}$ satisfies the equation
\begin{equation*}
\mx{A}(\mx{X} : \mx{V} \mx X^{\bot}) = (\mx X_f : \mx{V}_{21} \mx X^{\bot} ).
\end{equation*}
In terms of Pandora's Box (Theorem 2), $\mx{Ay}$ is the $\BLUP$
for $\mx y_f$ if and only if there exists a matrix~$\mx L$ such that
$\mx{A}$ satisfies the equation
\begin{equation*}
\begin{pmatrix}
\mx V & \mx X \\
\mx X' & \mx 0
\end{pmatrix}
\begin{pmatrix}
\mx A' \\
\mx L
\end{pmatrix} =
\begin{pmatrix}
\mx{V}_{12} \\
\mx X_{f}'
\end{pmatrix}.
\end{equation*}
\end{theorem}
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{The Mixed Model}
A mixed linear model can be presented as
\begin{equation*}
\mx y = \mx X\BETA + \mx Z \GAMMA +\EPS ,
\quad \text{or shortly } \quad
\M_{\mathrm{mix}}
= \{ \mx y,\, \mx X\BETA + \mx Z\GAMMA, \, \mx D,\,\mx R \} ,
\end{equation*}
where $\mx X \in \rz^{n \times p}$ and $\mx Z \in \rz^{n \times q}$ are
known matrices, $\BETA \in \rz^{p}$ is a vector of unknown fixed
effects, $\GAMMA$ is an unobservable vector ($q$ elements) of
\emph{random effects} with
$\cov(\GAMMA,\EPS) = \mx 0_{q \times p}$ and
\begin{equation*}
\E(\GAMMA) = \mx 0_q , \quad
\cov(\GAMMA) = \mx D_{q \times q}, \quad
\E(\EPS) = \mx 0_n \,, \quad
\cov(\EPS) = \mx R_{n\times n} .
\end{equation*}
This leads directly to:
%%%%%%%%%
%%%%%%%%%
%%%%%%%%%
%%%%%%%%%
\begin{theorem}
Consider the mixed model
$ \M_{\mathrm{mix}}
= \{ \mx y,\, \mx X\BETA + \mx Z\GAMMA, \, \mx D,\,\mx R \}.$
Then the linear estimator
$\mx B \mx y$ is the $\BLUE$ for $\mx X\BETA$ if and only if
\begin{equation*}
\mx B(\mx X : \SIGMA \mx X^{\bot}) = (\mx X : \mx{0}) ,
\end{equation*}
where $\SIGMA= \mx Z\mx D\mx Z' + \mx R$.
Moreover,
$\mx A \mx y$ is the $\BLUP$ for $\GAMMA$ if and only if
\begin{equation*}
\mx A(\mx X : \SIGMA \mx X^{\bot}) = (\mx 0 : \mx{D}\mx{Z}' \mx X^{\bot}).
\end{equation*}
In terms of
Pandora's Box (Theorem 2), $\mx A \mx y = \BLUP(\GAMMA)$ if and only if there
exists a matrix $\mx L$ such that $\mx{A}$ satisfies the equation
\begin{equation*}
\begin{pmatrix}
\SIGMA & \mx X \\
\mx X' & \mx 0
\end{pmatrix}
\begin{pmatrix}
\mx A' \\
\mx L
\end{pmatrix} =
\begin{pmatrix}
\mx Z \mx D \\
\mx 0
\end{pmatrix}.
\end{equation*}
\end{theorem}
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%

For the equality
between the
$\BLUP$s
under two mixed models, see
Haslett and Puntanen (2010b, 2010c).


\subsection{Note}
Reprinted with permission from Lovric, Miodrag (2011),
\emph{International
Encyclopedia of Statistical Science.} Heidelberg:
Springer Science+Business Media, LLC.

%%%%%%%%%
%%%%%%%%%
%%%%%%%%%
%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{References}
\begin{enumerate}



\item
Anderson, T. W. (1948).
On the theory of testing serial correlation.
\emph{Skandinavisk Aktuarietidskrift}, \textbf{31}, 88--116.
\\
\item
Baksalary, Jerzy K.; Rao, C. Radhakrishna and Markiewicz,
Augustyn (1992).
A study of the influence of the `natural restrictions' on
estimation problems in the singular {G}auss--{M}arkov model,
\textit{Journal of Statistical Planning and Inference}, \textbf{31}, 335--351.
\\
\item
Christensen, Ronald (2002).
\textit{Plane Answers to Complex Questions:
The Theory of Linear Models,} 3rd Edition. Springer, New York.
\\
\item
Haslett, Stephen J. and Puntanen, Simo (2010a).
Effect of adding regressors on the equality of the BLUEs under two linear models.
\textit{Journal of Statistical Planning and Inference}, \textbf{140}, 104--110,
\\
\item
Haslett, Stephen J. and Puntanen, Simo (2010b).
Equality of BLUEs or BLUPs under two linear models using stochastic restrictions.
\emph{Statistical Papers}, \textbf{51}, 465--475.
\\
\item
Haslett, Stephen J. and Puntanen, Simo (2010c).
On the equality of the BLUPs under two linear mixed models.
\emph{Metrika}, available online, DOI 10.1007/s00184-010-0308-6.
\\
\item
Isotalo, Jarkko and Puntanen, Simo (2006).
Linear prediction sufficiency for new observations in the general Gauss--Markov
model.
\emph{Communications in Statistics: Theory and Methods},
\textbf{35}, 1011--1023.
\\
\item
Kruskal, William (1967).
When are Gauss--Markov and least squares estimators identical?
{A} coordinate-free approach.
\emph{The Annals of Mathematical Statistics}, \textbf{39}, 70--75.
\\
\item
Mitra, Sujit Kumar and Moore, Betty Jeanne (1973).
Gauss--Markov estimation with an incorrect dispersion matrix.
\textit{Sankhy\={a}, Series~A}, \textbf{35}, 139--152.
\\
\item
Puntanen, Simo and Styan, George P. H. (1989).
The equality of the ordinary least squares estimator and the best
linear unbiased estimator [with comments by Oscar Kempthorne and by Shayle~R.
Searle and with ``Reply'' by the authors].
\textit{The American Statistician}, \textbf{43}, 153--164.
\\
\item
Puntanen, Simo; Styan, George P.~H. and Werner, Hans Joachim
(2000). Two matrix-based proofs that the linear estimator \textit{Gy} is
the best linear unbiased estimator. \textit{Journal of Statistical
Planning and Inference}, \textbf{88}, 173--179.
\\
\item
Rao, C. Radhakrishna (1967).
Least squares theory using an estimated dispersion matrix
and its application to measurement of signals. In
\textit{Proceedings of the Fifth Berkeley Symposium on Mathematical
Statistics
and Probability: Berkeley, California, 1965/1966},
vol.~1 (Eds. Lucien~M. Le Cam \& Jerzy Neyman), University of California Press,
Berkeley, pp.~355--372.
\\
\item
Rao, C. Radhakrishna (1971).
Unified theory of linear estimation.
\textit{Sankhy\=a, Series~A}, \textbf{33}, 371--394.
[Corrigenda (1972), \textbf{34}, p.~194 and p.~477.]
\\
\item
Rao, C. Radhakrishna (1974).
Projectors, generalized inverses and the BLUE's.
\textit{Journal of the Royal Statistical Society, Series~B}, \textbf{36}, 442--448.
\\
\item
Watson, Geoffrey S. (1967).
Linear least squares regression.
\emph {The Annals of Mathematical Statistics}, \textbf{38}, 1679--1699.
\\
\item
Zyskind, George (1967).
On canonical forms, non-negative covariance matrices and best
and simple least squares linear estimators in linear models.
\emph{The Annals of Mathematical Statistics}, \textbf{38}, 1092--1109.
\\
\item
Zyskind, George and Martin, Frank B. (1969).
On best linear estimation and general Gauss--Markov theorem
in linear models with arbitrary nonnegative covariance
structure.
\emph{SIAM Journal on Applied Mathematics}, \textbf{17}, 1190--1202.

\end{enumerate}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\end{document}
