magic.magic.post.proc(X,object,w)magic after fitting the
model with model matrix X.magic, if one was.)object contains rV (${\bf V}$, say), and
scale ($\phi$, say) which can be
used to obtain the require quantities as follows. The covariance matrix of
the parameters is ${\bf VV}^\prime \phi$. The vector of
estimated degrees of freedom for each parameter is the leading diagonal of
${\bf VV}^\prime {\bf X}^\prime {\bf W}^\prime {\bf W}{\bf X}$
where $\bf{W}$ is either the
weight matrix w or the matrix diag(w). The
hat/influence matrix is given by
${\bf WX}{\bf VV}^\prime {\bf X}^\prime {\bf W}^\prime$
.magic