《第11章 协方差结构模型课件.ppt》由会员分享，可在线阅读，更多相关《第11章 协方差结构模型课件.ppt（62页珍藏版）》请在三一办公上搜索。
1、第十章 协方差结构模型,10.1 验证性因素分析10.2 协方差结构模型10.3 LISREL8.20的使用10.4 AMOS4.0的使用,10.1 验证性因素分析1、探索性因素分析和验证性因素分析的区别（1）两者的理论假设不同：探索性因素分析的假设：所有的公共因素都相关（或都无关）；所有的公共因素都直接影响所有的观测变量；特殊因素之间相互独立；所有观测变量只受一个特殊因素的影响；公共因素和特殊因素相互独立。,验证性因素分析的基本假设：公共因素之间可以相关也可以无关；观测变量可以只受某一个或几个公共因素的影响而不必受所有公共因素的影响；特殊因素之间可以有相关。还可以出现不存在误差因素的观测变量
2、；公共因素和特殊因素相互独立。,（2）前提条件不同 当因素的结构、因素之间的关系、因素的数量未知时，应使用探索性因素分析；验证性因素分析是对已有的理论模型与数据拟合程度的一种验证。进行验证性因素分析时，必须明确指明：公共因素的个数、观测变量的个数、观测变量和公共因素之间的关系、观测变量和特殊因素之间的关系以及特殊因素之间的关系。,x1,x2,x3,x4,x5,1,2,探索型元素分析,验证性因素分析,2、验证性因素分析的数学模型验证性因素分析是在对研究问题有所了解的基础之上进行的，这种了解可以是建立在理论、实验研究或是两者结合的基础上。模型假设为：（1）在总体中，模型所有变量（观测变量、潜变量、
3、误差）都设定其平均值为零；（2）公共因素与误差项之间相互独立；（3）各独立因素之间相互独立（这一条件有时得到放宽）；（4）观测变量数大于公共因素数。,模型的估计就是求解上面协方差方程中的各个参数的估计值，以便使模型更好地重新产生观测变量的协方差矩阵。,3、验证性因素分析的步骤（1）模型定义 根据理论假设，定义观测变量和潜变量之间的关系，潜变量之间的关系以及特殊因素（误差项）之间的关系。,观测变量用矩形表示，潜变量用椭圆形表示；单箭头是从潜变量指向观测变量，表明潜变量引起了观测变量的变化；用双箭头表示相关，可以指定潜变量之间、误差项之间的相关。,（2）模型识别 识别问题就是协方差结构方程有唯一解
4、的问题。,1）恰好识别：指方程式的个数等于要估计的参数的个数，因此每个参数能求得唯一解。2）超识别：指方程式的个数多于参数估计所需要方程的个数。3）不可识别：模型中方程式的个数少于待估参数的个数，无法确定模型参数。,模型可识别的必要条件：模型中待估参数的个数要小于或等于q(q+1)/2,其中q为观测变量的个数。模型可识别的充分条件：如果潜变量之间的协方差矩阵，并且因素负荷矩阵的k列中至少有k-1列是规定的元素，则模型可识别。如果不是对角矩阵，但对角线上的元素是相同的，若因素负荷矩阵中的每一列中至少有k-1个值为规定的元素，则模型是可识别的。（k为公共因素数目）。三指标原则：a.每个潜变量有三个
5、或更多测量变量；b.因素负荷矩阵每一行有且只有一个非零值；c.残差的协方差矩阵为对角阵。如果满足上述三个条件，模型可识别。,4）两指标法则：法则1：a.因素负荷矩阵的每一行有且仅有一个非零值;b.至少有两个潜变量相关；c.每个潜变量至少有两个非零的测量变量；d.残差的协方差矩阵为对角阵。同时满足上述四个条件，则模型可识别。法则2：a.因素负荷矩阵的每一行有且仅有一个非零值;b.对于每个潜在变量，至少有一个潜变量与之相关；c.每个潜变量至少有两个非零的测量变量。同时满足上述三个条件，模型可识别。,（3）验证性因素分析的参数估计,（4）验证性因素分析模型的评价 得到了参数的估计值后，需要对模型与数
6、据间是否拟合进行评价。拟合指数主要反映了与S拟合的程度。常用的模型拟合指数有：绝对拟合指数有：2统计量，p0.05拟合优度指数GFI 0.90调整的拟合优度指数AGFI 0.90近似均方根误差RMSEA 0.90标准拟合指数NFI 0.90Tucker-Lewis 指数NNFI 0.90递增拟合指数 IFI 0.90,(5)模型修正 如果模型不能很好地拟合数据，就需要对模型进行修正或重新设定。研究者需要决定如何删除、增加和修改参数，通过修正模型可以增进模型的拟合程度。另外，对于模型的选取应该遵循省俭原则。即当两个模型同样吻合数据时，应当取两个模型中比较简单的一个。,4、验证性因素分析在测量上的
7、应用（1）构想效度和项目信度 通过数据与理论假设模型之间的吻合程度来表示一个测验构想效度的高低。信度定义为公共因素与观测变量相关的平方，它表示在一个观测变量的总方差中，能够由公共因素所解释的比例。如x1=111+1,则x1的信度为：REL(x1)=211(2)多质多法（multitrait-multimethod.MMMT）,(3)高阶因素分析,5、验证性因素分析的其他问题（1）样本容量被试数与自由参数的比例至少达到10:1（2）变量数量有的研究者主张每个潜在因子应该用三个或四个指标表示。,10.2 协方差结构模型,10.2.1 模型的产生和发展1、几个概念协方差结构模型(Covariance
8、 Structure Models,CSM);结构方程模型(Structural Equation Modeling,SEM);协方差结构分析(the analysis of covariance structure);线性结构分析(the linear structural relations models);矩结构模型(the moments structure models);结构化线性模型中的潜变量方程系统(Latent variable equation system linear model)LISREL模型,2、SEM与一般线性模型的关系（1）SEM实际是一般线性模型(Gener
9、al Linear Models,GLM)的扩展。一般线性模型包括：路径分析、典型相关、因素分析、判别分析、多元方差分析以及多元回归分析。它们都是协方差结构方程模型的特例，许多模型均可用SEM来处理和评价；（2）SEM是路径分析和因素分析的有机结合，它比传统的回归分析有更突出的优点。传统的回归分析的方法不考虑自变量之间的测量误差，会导致高估简单模型变量的真正变异量和低估相关系数。测量误差越大，用传统的回归方法导致的错误越大。,3、结构方程的优点（1）可处理多个因变量；（2）允许自变量和因变量含有测量误差；（3）容许潜在变量由多个外源指标变量组成，并可同时估计指标变量的信度和效度；（4）可采用比
10、传统方法更有弹性的测量模型，如某一观测变量或项目在SEM内可以同时从属于两个潜在变量。（5）可以考虑潜在变量之间的关系，并估计整个模型是否与数据相吻合。,10.2.2 数学模型的设定,模型中的变量：X外源观测变量；Y内源观测变量；外源潜在变量；内源潜在变量；外源观测变量的独立因子；内源观测变量的独立因子。外源变量因果关系中不受其他变量影响的变量；内源变量因果关系中受其他变量影响的变量。,协方差结构模型由两部分组成：（1）测量模型(Measurement Model)即验证性因素分析(Confirmatory Factor Analysis)（2）结构模型(Structural Equation
11、 Model)即潜变量的因果关系模型。测量模型主要用于表示观测变量与潜变量之间的关系；而结构模型主要用来表示潜变量之间的关系。,1、测量模型测量模型又称验证性因素分析模型，它是在对研究问题有所了解的基础上进行的，这种了解可以是建立在理论、实验研究或是两者相结合的基础上。模型假设为：（1）在总体中，模型中所有的变量（观测变量、潜变量、误差）都设定其平均值为零；（2）公共因子与误差项之间相互独立；（3）各独立因子之间相互独立。模型的数学表达式为：,观测变量的协方差矩阵与估计参数之间的关系：,2、结构方程模型结构方程模型用来描述潜变量之间的关系。该模型的假设为：（1）在总体中，模型所有的潜变量都是平
12、均数为零；（2）方程中的外源变量与误差之间的相关为零；（3）模型中潜变量关系不存在多余的方程。数学表达式为：+协方差方程为:,3、完整协方差结构模型模型假定：（1）所有变量来自均数为零的总体；（2）没有多余的方程存在；（3）方程中的误差项与外源变量的误差之间不相关。数学表达式：Xx+Y=y+=+协方差结构方程:,X1为高中平均成绩；X2为高考成绩；X3为高中教师评价；X4为学生性别。Y1为是否在重点院校；Y2为高校第一年平均成绩；Y3为教师对其高校第一年表现的评价。ksi1为高中成绩；ksi2为性别；eta1为是否重点院校；eta2为高校第一年成绩。,10.2.3 协方差结构模型的识别一、自由
13、度识别(必要条件)根据所定义模型的自由度与最大可能自由度之间的关系来简单判断模型是否有可能识别，这是一个必要条件。1、最大可能自由度的概念 设模型中有p个外源观测变量，q个内源观测变量，那么模型最大可能的自由度为(p+q)(p+q+1)/2.2、自由度识别方法 在协方差方程=()中包含(p+q)(p+q+1)/2个方程。如果待估参数的个数t大于方程的个数，则模型不可能识别。所以,t(p+q)(p+q+1)/2。模型的自由度df=(p+q)(p+q+1)/2t0。df 0,模型有可能识别；df0,模型一定不能识别。,P=5,q=6,(p+q)(p+q+1)/2=1112/2=66,需要估计的参数
14、:由内源变量指向测量变量的路径系数。本例为11个；测量变量的残差。11个；外源变量指向内源变量的路径系数。4个；潜变量之间的协方差和残差。家庭地位与学生能力的协方差1个，同伴关系和学业成就的2个残差，两个残差之间的协方差1个，共4个。外源潜变量的方差，2个。为了使模型可以识别，对于每个潜变量，一定要对其设立一条路径系数，通常为1，这些系数当然不需要估计，本例为4个。模型的自由度为：df=66-11-11-4-4-2+4=38,二、两阶段判别法（充分条件）两阶段判别法将完整的协方差结构模型分为验证性因素分析（测量部分）和潜变量因果关系模型（结构部分）两个部分，然后分别判断两部分模型是否可识别的方
15、法。如果测量部分和结构部分都可识别，那么，可以断定整个协方差结构模型可以识别。用两阶段判别法分析上例，,本图是上例模型的测量部分，由于符合两指标法则，所以可识别,上图为模型的结构部分，该模型为递归模型，所以可识别。协方差结构模型的两部分都可识别，所以，整个模型可以识别。,10.2.4 模型的估计 所谓模型估计就是求解模型中各个参数的估计值，使模型再生协方差矩阵尽可能好的拟合观测变量的协方差矩阵。有几种方法：1、极大似然法；2、非加权的最小二乘估计；3、广义最小二乘估计；10.2.5 模型的评价与修正拟合指数：1、绝对拟合指数（1）2=(n-1)F,一般希望得到不显著的2值。拟合优度的卡方检验受
16、样本容量的影响，常常不能很好地判断模型的拟合度。在样本容量很小时，几乎接受所有的拟合劣势的模型，而在大样本容量时，几乎拒绝所有的拟合较优的模型。,（2）、拟合优度指数(Goodness-of-fit index,GFI)和调整的拟合优度指数(Adjusted goodness-of-fit index,AGFI),（3）、近似误差均方根(Root mean square error of approximation,RMSEA),2、相对拟合指数（1）本特勒与波内特正态化拟合指数（Bentler-Bonett normed fit index,NFI）,（2）Bollen的IFI指数,（3）非
17、正态拟合指数（Nonormal fit index,NNFI）,（4）本特勒比较拟合指数(Comparative fit index,CFI),（5）信息标准指数(Information criteria indexes)Akaike Information Criterion,AIC AIC=c+2tConsistent Akaike Information Criterion,CAIC CAIC=c+1+ln(n)tExpected cross-validation index,ECVI ECVI=c/(n-1)+2t/(n-1)这些指标越小说明模型越简单并且拟合越好。,模型的修正 结构方
18、程是用于估计潜变量之间的关系，并验证某一模型是否与研究者所得的数据吻合。在实际中，研究者常根据一些统计分析结果如残差、模型修正指数，进而放宽(free)、固定(fix)或改动模型，以便使模型更拟合数据，整个过程就是模型的修正。几个概念：简约原则(principle of parsimony);等统模式(equivalent models);嵌套模型(nested model)修正指数：指一个固定参数或限定参数恢复为自由时卡方值可能减少最小的量。,10.3 LISREL8.20的使用,1、LISREL模型中个变量的符号及有关的LISREL命令,2、8个待估协方差矩阵在LISREL程序中的形式,说
19、明：ZE零矩阵；ID单位阵；IZ单位和零矩阵；ZI零和单位阵，DI对角阵；SY对称矩阵。,3、LISREL模型中的三类参数：固定参数规定了数值的参数（通常取值为0或1）；约束参数等于其他未知参数的的参数，多一个约束参数就减少了一个待估参数；自由参数其值未加任何限制的参数，是要估计的未知参数。,解决模型识别的常用方法：固定方差将潜变量标准化，令其方差为零；固定负荷将各潜变量的一个观测变量的负荷设定为1。,4、路径图绘制路径图的规则：（1）观测变量x和y要画在方框里；（2）潜变量和要画在圆圈里;（3）误差项、不画在圆圈里也不画在方框里；（4）单向箭头表示一个变量对另一个变量的直接影响，要画成直线；
20、（5）双向箭头表示相关关系而不是因果关系，要画成曲线；潜变量是作为原因变量，所以，没有箭头指向它们，是结果变量，所有指向的箭头都来自、，或结构方程的误差。路径图中的各个系数都有两个下标，第一个下标代表的是箭头终点变量的标号，第二个下标表示的是箭头起点变量的标号。,5、用PRELIS准备数据文件例1TITLE:covariance matrix prepared for confirmatory factor analysisDA NI=6 NO=73LAVisperc cubs lozenges paragraph sentence wordmeanRA FI=d:lisrel8data.da
21、tCO ALLOU MA=CM CM=d:lisrel8data.cm6、创建LISREL命令文件Lisrel命令文件包括三部分：数据输入：由DA定义模型构建：由MO定义结果输出：由OU定义,例2LISREL MODEL EXAMPLE 1=(FILENAME:STUDENT.LS8)DA NI=7 NO=1000 MA=CM CM FI=D:STUDENTSTUDENT.CMLA X1 X2 X3 X4 Y2 Y3 Y1SE Y1 Y2 Y3 X1 X2 X3 X4MO NY=3 NE=2 NX=4 NK=2 LY=FI LX=FI BE=SD,FI FR LY 3 2 FR LX 2 1
22、LX 3 1FR BE 2 1VA 1.0 LY 1 1 LY 2 2VA 1.0 LX 1 1 LX 4 2FI TE 1 1FI TD 4 4LK HIGH_SCH SUC GENDERLE KEYUNIV COLL ACHIPATH DIAGRAMOU ME=ML MI AD=OFF ND=4 IT=100 NS,（1）数据输入DA命令行功能：规定Lisrel分析的数据的类型。语法：DA NG=a NI=b NO=c MA=(矩阵类型)XMd RP=e关键词：NG组数 NI观测变量数 NO个案数 MA矩阵类型 XM规定数据的缺失值LA命令行功能：规定观测变量的标签语法：LA FI文件名
23、FO RE（变量格式说明）*,RA命令行功能：规定观测变量的标签语法：RA FI文件名 FO RE（变量格式说明）*例3 DA NI=5 RA FI=data.dat fo(5F3.1)例4 DA NI=5 NO=20 RA(5F3.1)35 42 68 93 101 37 51 68 93 101 32 51 72 34 91,CM,KM,MM,OM,PM,RM,TM命令行功能：读取汇总矩阵（summary matrix）,如协方差，相关系数矩阵。语法；如CM FI文件名 SY FO RE FU（文件格式说明）*如:CM FU(3F5.2)1.13-0.87 1.08-0.87 2.17 1
24、.83 1.08 1.83 3.25,CM SY1.13-0.87 2.17 1.08 1.83 3.25,SE命令行功能：将观测变量重新排序，Y变量在前，X变量在后。（2）模型设定MO命令行功能：模型设定语法：MO NY=e NX=f NE=g NK=h LY=mf,mm LX=mf,mm BE=mf,mm GA=mf,mm TD=mf,mm.,LE和LK指令功能:为和规定标签.语法:,LE FI=文件名 FO RELK变量格式说明*数据(标签),FI和FR命令行功能:将矩阵元素固定或释放为自由语法:FI 矩阵元素 FR,例如:FI LX(2,1)或FI LX 2 1使LX矩阵第二行第一列的
25、元素固定为0.,EQ命令行功能:规定矩阵两个或两个以上的元素相等。语法：EQ 矩阵元素列表例如：EQ LX 2 2 LY 3 3 GA 2 1将这三个矩阵相应的元素相等。CO命令行功能：规定一个参数等于其它参数的函数语法：CO 参数A其它参数的表达式例如：CO BE（1，2）BE(2，1)IR命令行功能：规定一个参数取值的区间语法：IR 参数Aa-1 1,PA命令行功能:将LISREL 的参数矩阵规定为由若干个0和1组成的参数矩阵的形式。语法：PA FI文件名 FO RE 矩阵名 变量格式说明*数据（整数1和0的模式）例如：,VA和ST命令行功能：VA的功能是对固定参数赋值，ST的功能是对自由
26、参数赋以迭代初始值。语法：VA 数值 参数矩阵的元素 ST ALL例如对以下两矩阵规定初始值,MA命令行功能：在ULS、GLS、ML、WLS和DWLS估计方法时作自由参数或固定参数的拟合函数图。NF命令行功能：规定不需要输出修正指数的自由和约束参数PD命令行功能：输出路径图语法：PATH DIAGRAM PD（3）结果输出OU命令行（1）功能：规定参数估计的过程,语法：OU ME=RC=c SL=100a NS RO AM SO,IVTSULGLMLWLDW,关键词：ME 规定估计方法，默认为极大似然法MLRC 规定岭常数，为使矩阵变为正定，需要在矩阵对角线上加一个常数。SL 模型修正指数的置
27、信水平NS 不设待估参数的初始值AM 自动修正。程序将根据修正指数的大小逐个将修正指数显著的固定参数恢复为自由估计。,OU命令行（2）功能：选择LISREL的打印结果语法：OU SE TV PC PT RS EF MR MI XM XI FS SS SC ALL TO WPND=D选项：SE 标准误TV t值PC 估计参数的相关矩阵PT 打印技术性结果RS 残差与标准残差，Q图和拟合矩阵EF 总效应和间接效应MR 多方面结果MI 修正指数XM 不计算修正指数FS 因子得分回归SS 标准化解 TO 每行打印80个字符SC 完全标准化解 WP 每行打印132个字符ALL 打印所有结果 ND 输出结
28、果大小数位数，默认值为2,OU 命令行（3）功能：程序将输出结果存入指定的文件中。语法：OU 矩阵1文件名1 矩阵2文件名2.关键词：矩阵i 替换要保存的文件名,可能的名称为LY,LX,BE,GA,PH PS,TE,TD TH AL,KA,TX,TY或MA 重新选择或变量定序后的文件名SI 拟合的方差、协方差或相关矩阵RM 潜变量在显变量上的因子负荷矩阵，使用时要有PS指令。EC 参数估计的协方差矩阵GF 所有的拟合度量PV 估计的自由参数的向量SV 相应的标准误差的向量TV 相应的t值的向量,OU命令行（4）功能：选择LISREL的打印结果语法：OU TMf IT=n AD=m EP=off关键词:TM 显示运行程序所用的CPU时间(time)IT 迭代次数（iteration）,默认次数为自由参数的三倍AD 检查参数解的可接受性EP 收敛标准，缺省值是EPS=0.000001。,LISREL的应用实例,
