Я рассматриваю следующую линейную модель: .Y= Xβ+ ϵ
Вектор невязок оценивается как
ε^= у- Хβ^= ( Я- Х( X'Икс)- 1Икс') у= Q у= Q ( Xβ+ ϵ ) = Q ϵ
где .Q = я- Х( X'Икс)- 1Икс'
Заметим, что tr ( Q ) = n - p (след инвариантен относительно циклической перестановки) и что Q'= Q = Q2 . Следовательно, собственные значения Q равны 0 и 1 (некоторые детали приведены ниже). Следовательно, существует унитарная матрица В такая, что ( матрицы диагонализируются унитарными матрицами тогда и только тогда, когда они нормальны ).
В'Q V= Δ = диаг ( 1 , … , 1n - p раз, 0 , … , 0р раз)
Теперь пусть .К= V'ε^
Поскольку , мы имеем и, следовательно, . таким образомK~N(ε^∼ N( 0 , σ2Q )К п - р + 1 = ... = К п = 0К∼ N( 0 , σ2Δ )Кn - p + 1= … = КN= 0
∥ К∥2σ2= ∥ К⋆∥2σ2∼ χ2н - р
с .К⋆= ( К1, … , Кн - р)'
Кроме того, поскольку является унитарной матрицей, мы также имеемВ
∥ ϵ^∥2= ∥ К∥2= ∥ К⋆∥2
таким образом
RSSσ2∼ χ2н - р
Наконец, обратите внимание, что этот результат подразумевает, что
Е( RSSн - р) = σ2
Поскольку , минимальный многочлен от делит многочлен . Итак, собственные значения находятся среди и . Поскольку также является суммой собственных значений, умноженных на их кратность, мы обязательно имеем, что является собственным значением с кратностью а ноль является собственным значением с кратностью .Q z 2 - z Q 0 1 tr ( Q ) = n - p 1 n - p pQ2- Q = 0QZ2- зQ01tr ( Q ) = n - p1н - рп