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 集合论