Analysis of covariance

From Wikipedia, the free encyclopedia

Jump to: navigation, search

Analysis of covariance (ANCOVA) is a general linear model with one continuous outcome variable (quantitative) and one or more factor variables (qualitative). ANCOVA is a merger of ANOVA and regression for continuous variables. ANCOVA tests whether certain factors have an effect on the outcome variable after removing the variance for which quantitative predictors (covariates) account. The inclusion of covariates can increase statistical power because it accounts for some of the variability.

Contents

[edit] Assumptions

As any statistical procedure, ANCOVA makes certain assumptions about the data entered into the model. Only if these assumptions are met, at least approximately, will ANCOVA yield valid results. Specifically, ANCOVA, just like ANOVA, assumes that the residuals [1] are normally distributed and homoscedastic. Further, since ANCOVA is a method based in linear regression, the relationship of the dependent variable to the independent variable(s) must be linear in the parameters.

[edit] Power considerations

While the inclusion of a covariate into an ANOVA generally increases statistical power by accounting for some of the variance in the dependent variable and thus increasing the ratio of variance explained by the independent variables, adding a covariate into ANOVA also reduces the degrees of freedom (see below). Accordingly, adding a covariate which accounts for very little variance in the dependent variable might actually reduce power.

[edit] Equations

[edit] One-factor ANCOVA analysis

One factor analysis is appropriate when dealing with more than 3 populations; k populations. The single factor has k levels equal to the k populations. n samples from each population are chosen randomly from their respective population.

[edit] Calculating the sum of squared deviates for the independent variable X and the dependent variable Y

The sum of squared deviates (SS): SSTy, SSTry, and SSEy must be calculated using the following equations for the dependent variable, Y. The SS for the covariate must also be calculated, the two necessary values are SSTx and SSEx.

The total sum of squares determines the variability of all the samples. nT represents the total number of samples:

SST_y=\sum_{i=1}^n\sum_{j=1}^kY_{ij}^2-\frac{\left(\sum_{i=1}^n\sum_{j=1}^kY_{ij}\right)^2}{n_T}

The sum of squares for treatments determines the variability between populations or factors. nk represents the number of factors

SSTr_y=\sum_{i=1}^n\left(\frac{\sum_{j=1}^kY_{ij}^2}{n_k}\right)-\frac{\left(\sum_{i=1}^n\sum_{j=1}^kY_{ij}\right)^2}{n_T}

The sum of squares for error determines the variability within each population or factor. nn represents the number of samples with a given population:

SSE_y=\sum_{i=1}^n\sum_{j=1}^kY_{ij}^2-\sum_{i=1}^n\left(\frac{\sum_{j=1}^kY_{ij}^2}{n_k}\right)

The total sum of squares is equal to the sum of squares for treatments and the sum of squares for error:

SST_y=SSTr_y+SSE_y\,

[edit] Calculating the covariance of X and Y

The total sum of square covariates determines the covariance of X and Y within all the data samples:

SCT=\sum_{i=1}^n\sum_{j=1}^kX_{ij}Y_{ij}-\frac{\left(\sum_{i=1}^n\sum_{j=1}^kX_{ij}\right)\left(\sum_{i=1}^n\sum_{j=1}^kY_{ij}\right)}{n_T}
SCE=\sum_{j=1}^k\left(\sum_{i=1}^nX_{ij}Y_{ij}-\frac{\sum_{i=1}^n(X_{ij}Y_{ij})}{n_n}\right)

[edit] Adjusting SSTy

The correlation between X and Y is r_T^2.

r_T^2=\frac{SCT^2}{SST_xSST_y}
r_n^2=\frac{SCE^2}{SSE_xSSE_y}

The proportion of covariance is subtracted from the dependent, SSy values:

SST_{yadj}=SST_y-r_T^2\,
SSE_{yadj}=SSE_y-r_n^2\,
SSTryadj = SSTyadj − SSEyadj

[edit] Adjusting the means of each population k

The mean of each population is adjusted in the following manner:

M_{y_iadj}=M_{y_i}-\frac{SCE_y}{SCE_x}(M_{x_i}-M_{x_T})

[edit] Analysis using adjusted sum of squares values

Mean squares for treatments where dfTr is equal to NT − k − 1. dfTr is one less than in ANOVA to account for the covariance and dfE = k − 1:

MSTr=\frac{SSTr}{df_{Tr}}
MSE=\frac{SSE}{df_E}

The F statistic is

F_{df_E,df_\mathrm{Tr}}=\frac{\mathrm{MSTr}}{\mathrm{MSE}}.

[edit] References

  1. ^ http://www.statsoft.com/textbook/stbasic.html

[edit] External links

At Wikiversity, you can learn about: Analysis of covariance
Retrieved from "http://en.wikipedia.org/wiki/Analysis_of_covariance"
Personal tools
Languages