FAQ/covary - CBU statistics Wiki
location: FAQ / covary

# Testing the value of a mean adjusted for a covariate in a single group

Suppose we have an outcome, y, related to a variable, x, whose mean we wish to compare to zero, adjusted for x. For a linear regression with coefficient, B and centering x by subtracting its mean we have

Predicted Y = intercept + B(x - mean of x)

When x is equal to its mean the predicted value of Y is the value of the intercept so we can just read off the t-value for the intercept or constant term in the regression coefficients output to test whether the adjusted y mean equals zero.

This is analogous to a (multi-group) ancova where the covariate group mean is replaced by the overall pooled covariate group mean. The group means are adjusted depending upon how far each covariate group mean is from the overall covariate mean. The adjusted response mean is simply the predicted response evaluated at the overall mean of the covariate.

In the case of one sample the overall mean and the group mean are identical so we obtain the adjusted mean at the covariate mean. The y mean, as measured by the intercept, is independent of the covariate when it has its mean subtracted from all values (see below) and, therefore, equals the (unadjusted) y mean.

What does change when adding a covariate is the standard error of the mean because we have removed variance in the outcome attributed to the covariate so the covariate adjusted standard error of the mean will be less than the unadjusted one.

More particularly the standard error of both means = $$\mbox{MSE}/sqrt{\mbox{N}}}$$ but the MSE (Mean Square Error) of unexplained response variance with the covariate present will be less. The t-test of the intercept will have N-2 degrees of freedom because two terms (intercept and slope) are estimated.

When the variable, x, is centered by subtracting its overall mean

 Matrix X 1, ..., 1 $$x_text{1}-\bar{x}$$, ..., $$x_text{N}- \bar{x}$$

$$X^text{T}X$$ =

 N 0 0 $$\sum_text{i}(x_text{i}-\bar{x})^text{2}$$

$$(XTX)-1$$ =

 1/N 0 0 $$\frac{1}{\sum_text{i}(x_text{i}-\bar{x})^text{2}}$$

$$X^text{T}Y$$ =

 $$\sum_text{i}y_text{i}$$ $$(\sum_text{i}(x_text{i}-\bar{x})y_text{i})$$

Since the diagonals in the inverse are both zero it follows since the least squares estimates are $$(XTX)-1XTY$$ that the intercept equals the response mean regardless of the presence of the centred covariate.

None: FAQ/covary (last edited 2013-05-15 15:50:50 by PeterWatson)