Residuals from AR1xAR1 one model

From: kwstat <asremlforum_at_VSNI.CO.UK>
Date: Thu, 25 Sep 2008 19:43:53 +0100

I usually think of residuals as being the difference between the observed values and the modeled/fitted values. Residual plots can then be used as diagnostics for the fit of the model.

When I fit a model with an AR1xAR1 spatial structure in the residual covariance matrix, the above philosophy is not quite right. For example, if I look at a heat map/level plot of the residuals, I usually see a smoothish, undulating trend. This is the structure imposed by the correlation of the residuals, and is not structure in the data that that has not been adequately modeled.

There is a resid() function in ASREML-R, but these aren't really deviations from the model in the usual sense, are they?

What is the right way to think about these quantities?

Kevin Wright

-------------------- m2f --------------------

Sent using Mail2Forum (

Read this topic online here:

-------------------- m2f --------------------
Received on Sat Sep 25 2008 - 19:43:53 EST

This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)