2009-06-27

如何在SEM中处理dummy变量和显示R平方?

月出 @ 2009-06-09:

祝老师,您好!您的伯克真的很好,blogspot我一直登不上。我初学LISREL,有几个问题想请教您:

1、虚拟变量在SEM中一样用0、1处理吗?我在原先运行得出的模型里加入了性别变量,但就显示

W_A_R_N_I_N_G: PHI is not positive definite
W_A_R_N_I_N_G: PSI is not positive definite
W_A_R_N_I_N_G: The solution was found non-admissible after  50 iterations. The following solution is preliminary and is provided only for the purpose of tracing the source of the problem. Setting AD> 50 or AD=OFF may solve the problem. 这是怎么回事?

2、子模型原先得出的变量之间的关系为什么进入母模型后显著性发生变化?路径系数的正负也反了?正常吗?怎么解释?比方说总风险和各类风险之间本来是这么个关系,加入信任置于总风险前,总风险和各类风险之间的关系就变化了。可能有什么问题存在吗?

3、我怎么在LISREL中显示可决系数R平方阿?

庄主 @ 2009-06-27:

1、SEM中可以包括虚拟变量(dummy variables)。一般是先用LISREL中的PRELIS工具将含虚拟变量的相关系数矩阵转化成tetrachoric correlation矩阵(参见http://ourworld.compuserve.com/homepages/jsuebersax/tetra.htm),然后再进行分析。

你得到的出错信息由很多种原因造成,可能是你模型的设置有问题(如性别要factor loading = 1、factory error = 0)、也许是性别与某些其它观测变量之间的高度相关、等等,很难一下子知道是什么原因。如果用了上述转化矩阵,仍有类似问题,那么更有可能是模型设置的问题。

2、子模型的关系进入母模型之后显著性出现变化是常见的、但是如果系数方向也反了,说明母模型中含有一个disorter的第三者变量(具体参见我前些日子的有关帖子)。你也可以将先后两个模型的系数用表格的形式贴出来给大家看看。

3、LISREL将R平方叫做 Sqaured Multiple Correlations (SMC)。如果你是用SIMPLIS,需要加一句“LISREL OUTPUT EF“ 来显示。需要尤其注意的是,LISRE同时给出计算Structural Form 和 Reduced Form两种结果,前者中的SMC是错的,要看后者。(具体参见Joreskog的文章 What is the interpretation of R2?)

2 comments:

Anonymous said...

LISREL中的解释方差,是不是可以通过内生变量的residual的方差计算出来,即把1减去psi。

庄主 said...

如果你仔细读一下我引用的Joreskog的文章,就可以知道,如果因变量的公式右边(即path图的左边)不含有endorgeneous factor的话,structural form = reduced form,而因变量的R平方确实就是1-PSI。但是,如果因变量的预测变量中含有endorgeneous factor的话(这往往更常见3),那么structural form与reduced form不相等,后者的R平方小于前者,也不再等于1-PSI了。