通用线性模型使我们可以将ANOVA模型编写为回归模型。假设我们有两组,每个组都有两个观测值,即向量中有四个观测值。然后,原始的,overparametrized模型是È (Ý )= X ⋆ β ⋆,其中X ⋆是预测值的矩阵,即,虚设编码表示变量:
(μ 1 μ 1 μ 2 μ 2) = (1 1 0 1 1 0 1 0 1 1 0 1ÿË(y)= X⋆β⋆X⋆
⎛⎝⎜⎜⎜μ1个μ1个μ2μ2⎞⎠⎟⎟⎟= ⎛⎝⎜⎜⎜1个1个1个1个1个1个00001个1个⎞⎠⎟⎟⎟⎛⎝⎜β⋆0β⋆1个β⋆2⎞⎠⎟
的参数不能被识别为因为X ⋆具有秩2( (X ⋆)' X ⋆不可逆)。要改变这种状况,我们引入约束β ⋆ 1 = 0(治疗对比),这为我们提供了新的模式Ë (Ÿ )= X β:
(((X⋆)′X⋆)−1(X⋆)′E(y)X⋆(X⋆)′X⋆β⋆1=0E(y)=Xβ
⎛⎝⎜⎜⎜μ1μ1μ2μ2⎞⎠⎟⎟⎟=⎛⎝⎜⎜⎜11110011⎞⎠⎟⎟⎟(β0β2)
μ1=β0β0μ2=β0+β2β2μ2−μ1
tψ=∑cjβjψ0c=(0,1)′β2=0μ2−μ1=0ψ^=∑cjβ^jβ^=(X′X)−1X′yψ
t=ψ^−ψ0σ^c′(X′X)−1c−−−−−−−−−√
σ^2=∥e∥2/(n−Rank(X))∥e∥2Rank(X)=2(X′X)−1X′=(.5−.5.5−.50.50.5)β^0=0.5y1+0.5y2=M1β^2=−0.5y1−0.5y2+0.5y3+0.5y4=M2−M1c′(X′X)−1c
t=M2−M1−0σ^=M2−M1∥e∥2/(n−2)−−−−−−−−−−√
ttn−Rank(X)n−2t(M2−M1)2/1∥e∥2/(n−2)=SSb/dfbSSw/dfw=FFbwFn−Rank(X)
βj1≤jψ