AMOS输出解读惠顿研究惠顿数据文件在各种结构方程模型中被当作经典案例,包括AMOS 和LISREL。
本文以惠顿的社会疏离感追踪研究为例详细解释AMOS的输出结果。
AMOS同样能处理与时间有关的自相关回归。
惠顿研究涉及三个潜变量,每个潜变量由两个观测变量确定。
67疏离感由67无力感(在1967年无力感量表上的得分)和67无价值感(在1967年无价值感量表上的得分)确定。
71疏离感的处理方式相同,使用1971年对应的两个量表的得分。
第三个潜变量,SES(社会经济地位)是由教育(上学年数)和SEI(邓肯的社会经济指数)确定。
解读步骤1.导入数据。
AMOS在文件ex06-a.amw中提供惠顿数据文件。
使用File/Open,选择这个文件。
在图形模式中,文件显示如下。
虽然这里是预定义模式,图形模式允许你给变量添加椭圆,方形,箭头等元素建立新模型2.模型识别。
潜变量的方差和与它关联的回归系数取决于变量的测量单位,但刚开始谁知道呢。
比如说要估计误差的回归系数同时也估计误差的方差,就好像说“我买了10块钱的黄瓜,然后你就推测有几根黄瓜,每根黄瓜多少钱”,这是不可能实现的,因为没有足够的信息。
如何告诉你“我买了10块钱的黄瓜,有5根”,你便可以推出每根黄瓜2块钱。
对潜变量,必须给它们指定一个数值,要么是与潜变量有关的回归系数,要么是它的方差。
对误差项的处理也是一样。
一旦做完这些处理,其它系数在模型中就可以被估计。
在这里我们把与误差项关联的路径设为1,再从潜变量指向观测变量的路径中选一条把它设为1。
这样就给每个潜变量设置了测量尺度,如果没有这个测量尺度,模型是不确定的。
有了这些约束,模型就可以识别了。
注释:设置的数值可以是1,也可以是其它数,这些数对回归系数没有影响,但对误差有影响,在标准化的情况下,误差项的路径系数平方等于它的测量方差。
3.解释模型。
模型设置完毕后,在图形模式中点击工具栏中计算估计按钮运行分析。
点击浏览文本按钮。
输出如下。
蓝色字体用于注解,不是AMOS输出的一部分。
TitleExample6,Model A:Exploratory analysis Stability of alienation, mediated by ses.Correlations,standard deviations and means from Wheaton et al.(1977).以上是标题,全是英文,自己翻译去吧,没有什么价值,一堆垃圾。
Notes for Group(Group number1)The model is recursive.Sample size=932各组注释:Group number1是模型内定的模型名称,因为你还没有给模型取名。
它告诉你模型为递归模型,样本量为932。
Variable Summary(Group number1)Your model contains the following variables(Group number1)Observed,endogenous variablesanomia67powles67anomia71powles71educatioSEIUnobserved,endogenous variables71_alienation67_alienationUnobserved,exogenous variableseps1eps2eps3eps4sesdelta1zeta1zeta2delta2变量汇总:对模型中的变量作一些概括,内生观测变量:67无力感,67无价值感,71无力感,71无价值感,教育和SEI。
内生非观测变量:67疏离感,71疏离感。
外生非观测变量:各种误差和社会经济地位。
注释:观测变量与非观测变量的区别:一个用方形表示,一个用椭圆表示。
内生和外生的区别:箭头指向自己的就是内生,发送箭头的就是外生。
注意区分测量模式和结构模式。
Variable counts(Group number1)Number of variables in your model:17Number of observed variables:6Number of unobserved variables:11变量计数:数数模型中的变量,变量总数为17,其中观测变量有6个,非观测变量有11个;外生变量有9个,内生变量有8个。
Parameter summary (Group number 1)模型的参数概括:固定系数11个,就是模型识别中固定的11个1。
还有6个自由的系数,9个方差对应着前面外生非观测变量。
Computation of degrees of freedom (Default model)(内定模型)的自由度计算:21"样本矩"是6个观测变量的6个样本方差加上15个协方差构成(也就是6中取2的组合数)。
15个参数是模型的6个回归系数和9个被估计的方差。
样本矩与估计参数的差为6个自由度。
Number of exogenous variables:9Number of endogenous variables:8Weights CovariancesVariancesMeans InterceptsTotal Fixed 11000011Labeled 000000Unlabeled 6090015Total17090026Number of distinct sample moments:21Number of distinct parameters to be estimated:15Degrees of freedom (21-15):6(内定模型)迭代过程:极大似然估计是一个迭代过程。
这里给出迭代历史。
这个输出是可选的,你不必直接使用它。
基本上没有什么用。
Result(Default model)Minimum was achievedChi-square=71.544Degrees of freedom=6Probability level=.000卡方拟合指数:这是所有软件都使用的最普通的拟和检验。
AMOS和LISREL把它称为卡方统计量,其它软件称为卡方拟和优度和卡方拟和劣度。
卡方拟合指数检验选定的模型协方差矩阵与观察数据协方差矩阵相匹配的假设。
原假设是模型协方差阵等于样本协方差阵。
如果模型拟合的好,卡方值应该不显著。
在这种情况下,数据拟和不好的模型被拒绝。
卡方检验的问题是样本越大,越可能拒绝模型,越可能犯第一类错误。
卡方拟和指数对违反多变量正态假设也是非常敏感。
这由卡方拟和指数的计算公式可以看出:卡方统计量=(N-1)x FN是样本量,F是模型协方差阵和样本协方差阵的最小适配函数。
这个函数比较复杂,也不知道是哪个天才搞出来的,它的计算公式中包含行列式,矩阵的迹,还要取对数,再经过一些加减运算把多维数据压缩为一个数值。
从卡方统计量的计算中可以看出,如果适配函数减少的速度没有样本量增加的速度快,即使模型协方差阵与样本协方差阵拟和的很好,但样本量的增加也会导致拒绝原假设。
这种拒绝正确建议的行为就是犯了第一类错误。
如果不服从正态分布,卡方统计量会更多地拒绝真实模型。
不过好在ML 估计比较稳健,所以即使违背了正态分布的假定,模型也能对付着用。
Maximum Likelihood EstimatesSEM使用最大似然法估计模型,而不是通常的最小二乘法。
OLS寻找数据点到回归线距离的最小平方和。
MLE寻找最大的对数似然,它反映从自变量观测值预测因变量观测值的可能性有多大。
Regression Weights:(Group number1-Default model)Estimate S.E. C.R.P Label67_alienation<---ses-.614.056-10.912***par_671_alienation<---67_alienation.705.05313.200***par_471_alienation<---ses-.174.054-3.213.001par_5 powles71<---71_alienation.849.04220.427***par_1anomia71<---71_alienation 1.000<---67_alienation.888.04320.577***par_2 anomia67<---67_alienation 1.000<---ses 1.000<---ses 5.331.43112.370***par_3回归系数是模型中带箭头的路径系数。
为了识别模型,部分系数在模型识别中已固定为1(例如,潜变量67疏离感到观测变量67无力感的路径)。
也给出路径系数的标准误。
"C.R."是临界比,它是回归系数的估计值除以它的标准误(-0.614/0.056=-10.912)。
临界比与原假设有关,在这个案例中对67疏离感和社会经济地位的原假设是回归系数为0。
如果我们处理近似标准正态分布的随机变量,在0.05的显著性水平上,临界比估计的绝对值大于 1.96称之为显著。
这样67疏离感和社会经济地位的回归系数-10.912的绝对值大于 1.96,可以说这个回归系数在0.05显著性水平上显著地不等于0。
P值给出检验原假设总体中参数是0的近似双尾概值。
它表示67疏离感和社会经济地位的回归系数显著地不等于0,p=0.001。
P值的计算假定参数估计是正态分布,它只是对大样本正确。
Variances:(Group number1-Default model)Estimate S.E. C.R.P Labelses 6.656.64110.379***par_7zeta1 5.301.48310.967***par_8zeta2 3.737.3889.623***par_9eps1 4.010.35811.186***par_10eps2 3.187.28411.242***par_11eps3 3.696.3919.443***par_12eps4 3.622.30411.915***par_13Estimate S.E. C.R.P Labeldelta1 2.944.501 5.882***par_14delta2260.63018.25614.277***par_15方差的估计,标准误和临界比和P值的解释同上。
用表格看数据总是让人眼花缭乱,还是看图示舒服些,这是上面表格数字的图形显示。
Modification Indices(Group number1-Default model)Covariances:(Group number1-Default model)M.I.Par Changeeps2<-->delta1 5.905-.424eps2<-->eps426.545.825eps2<-->eps332.071-.988eps1<-->delta1 4.609.421eps1<-->eps435.367-1.069M.I.Par Changeeps1<-->eps340.911 1.253Variances:(Group number1-Default model)M.I.Par ChangeRegression Weights:(Group number1-Default model)M.I.Par Changepowles71<---powles67 5.457.057powles71<---anomia679.006-.065anomia71<---powles67 6.775-.069anomia71<---anomia6710.352.076powles67<---powles71 5.612.054powles67<---anomia717.278-.054anomia67<---powles717.706-.070anomia67<---anomia719.065.068修正指数(MI)。