UA STAT687 线性模型理论I 线性模型概述
UA STAT687 線性模型理論I 線性模型概述
- 線性回歸
- One-way ANOVA
- Two-way ANOVA
- Nested Design
- Cross Design
- ANCOVA
線性模型是統計中一類模型的總稱,包括線性回歸模型、ANOVA模型、線性時間序列模型等,因此線性模型理論研究的是這一類模型共同的性質。稱解釋變量XXX與被解釋變量yyy之間滿足
y=Xb+ey=Xb+ey=Xb+e
的模型為線性模型,其中bbb與X,yX,yX,y無關,eee是隨機誤差。
第I部分介紹一些常見的線性模型,這些模型在UA MATH571A與UA MATH571B這兩個系列中都分別介紹過,只是這里再綜合一下,讓大家對線性模型具體指哪一類模型有一個總覽性的認識,并且掌握用矩陣來描述分量形式表示的模型。
線性回歸
XXX是一個N×kN \times kN×k的矩陣,表示解釋變量,每nnn行表示第nnn組樣本,每kkk列表示第kkk個解釋變量的所有樣本;yyy是NNN維列向量,表示被解釋變量的所有樣本,線性回歸模型可以表示為:
y=Xβ+ey = X \beta + ey=Xβ+e
其中eee表示隨機誤差,通常假設Ee=0,Cov(e)=σ2INEe=0,Cov(e) = \sigma^2I_NEe=0,Cov(e)=σ2IN?(Gauss-Markov假設);另外,線性回歸中的線性具體指的是β\betaβ關于yyy是線性的。線性回歸模型要寫成矩陣的形式是非常直觀的。
One-way ANOVA
考慮一般情況,不平衡的RCD:
yij=μ+αi+eij,i=1,?,a;j=1,?,niy_{ij}=\mu+\alpha_i+e_{ij},i=1,\cdots,a; j=1,\cdots,n_iyij?=μ+αi?+eij?,i=1,?,a;j=1,?,ni?
這個模型可以改寫成
需要注意的是yyy與eee是一個二階的array,改寫成矩陣的時候寫成一個列向量,排序方式是先iii后jjj。αi\alpha_iαi?表示treatment factor第iii個factor level的treatment effect,它與yi.y_{i.}yi.?以及1ni\textbf{1}_{n_i}1ni??對應。1\textbf{1}1填充的位置表示對應位置的response由哪些effect構成。
Two-way ANOVA
Nested Design
假設有兩個factor X,YX,YX,Y,YYY嵌套在XXX中:XXX的取值為F,PF,PF,P,YYY的取值為A,B,D,E,MA,B,D,E,MA,B,D,E,M,其中D,MD,MD,M嵌套在FFF中,A,B,EA,B,EA,B,E嵌套在PPP中,則模型可以寫為
yXYk=μ+αX+βY(X)+eXYky_{XYk}=\mu+\alpha_X+\beta_{Y(X)}+e_{XYk}yXYk?=μ+αX?+βY(X)?+eXYk?
改寫成矩陣為:
Cross Design
考慮模型
yijk=μ+αi+βj+eijk,i=1,?,a,j=1,?,by_{ijk}=\mu+\alpha_i+\beta_j+e_{ijk},i=1,\cdots,a,j=1,\cdots,byijk?=μ+αi?+βj?+eijk?,i=1,?,a,j=1,?,b
改寫為矩陣:
ANCOVA
在RCD中,如果試驗單位存在重要的covariate時,可以用ANCOVA,把covariate也加入到模型中,與One-way ANOVA,模型要修正為
yij=μ+αi+βxij+eij,i=1,?,a;j=1,?,niy_{ij}=\mu+\alpha_i+\beta x_{ij}+e_{ij},i=1,\cdots,a; j=1,\cdots,n_iyij?=μ+αi?+βxij?+eij?,i=1,?,a;j=1,?,ni?
改寫成矩陣:
總結
以上是生活随笔為你收集整理的UA STAT687 线性模型理论I 线性模型概述的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: UA MATH523A 实分析1 集合论
- 下一篇: UA MATH523A 实分析1 集合论