Hat Matrix - Blockwise Formula

Blockwise Formula

Suppose the design matrix can be decomposed by columns as . Define the Hat operator as . Similarly, define the residual operator as . Then the Hat matrix of can be decomposed as follows:


H(C) = H(A) + H(M(A) B)

There are a number of applications of such a partitioning. The classical application has a column of all ones, which allows one to analyze the effects of adding an intercept term to a regression. Another use is in the fixed effects model, where is a large sparse matrix of the dummy variables for the fixed effect terms. One can use this partition to compute the hat matrix of without explicitly forming the matrix, which might be too large to fit into computer memory.

Read more about this topic:  Hat Matrix

Famous quotes containing the word formula:

    The formula for achieving a successful relationship is simple: you should treat all disasters as if they were trivialities but never treat a triviality as if it were a disaster.
    Quentin Crisp (b. 1908)