首页 > 学术百科

用SAS实现结构方程模型分析

SAS Library  Solving Examples from EQS Manual Using Proc Calis 
Example 1: A Simple Regression Model on page 15
 
EQS code: manul1.eqs
SAS code
/TITLE
877uu  STABILITY OF POWERLESSNESS
  (EXAMPLE IN EQS MANUAL P.15)
/SPECIFICATIONS
  CAS=932; VAR=6; ME=LS;
/EQUATIONS
  V4 = 1*V2 - 1*V5 + E4;
/VARIANCES
  V2 = 9*; V5 = 9*; E4 = 2*;
/COVARIANCES
  V5,V2 = -4*;
/MATRIX
  11.834
    6.947  9.364
    6.819  5.091 12.532
    4.783  5.028  7.495  9.986
  -3.839 -3.889 -3.841 -3.625 9.610
  -2.189 -1.883 -2.175 -1.878 3.552 4.503
/END
data power(TYPE=COV);
  _type_ = 'cov';
  input _name_ $ v1-v6;
  datalines;                               
  v1  11.834    .        .        .      .        .
  v2    6.947    9.364    .        .      .        .
  v3    6.819    5.091  12.532    .      .        .
  v4    4.783    5.028    7.495    9.986    .        .
  v5  -3.839  -3.889  -3.841  -3.625  9.610    .
  v6  -2.189  -1.883  -2.175  -1.878  3.552  4.503
;
run;
proc calis cov data = power method = ls nobs = 932 ;
      Lineqs 
        V4 = b1 v2 + b2 v5 + E4;
      Std                     
        v2 = ev2(9),
        v5 = ev5(9),
    e4 = ee4(2);
      Cov             
        v2 v5 = Theta (-4);   
run;
Example 2: A Two-Equation Path Model on page 22
 EQS code: manul2.eqs
SAS code
/TITLE
  PATH ANALYSIS MODEL
  (EXAMPLE IN EQS MANUAL P.22)
/SPECIFICATIONS
  CAS=932; VAR=6; ME=ML;
/LABEL
  V1=ANOMIE67; V2=POWRLS67;
  V3=ANOMIE71; V4=POWRLS71;
/EQUATIONS
  V3 = 1*V1 + 1*V2 + E3;
  V4 = 1*V1 + 1*V2 + E4;
/VARIANCES
  V1 TO V2 = 10*;
  E3 TO E4 =  2*;
/COVARIANCES
  V2,V1= 7*;
/MATRIX
  11.834
    6.947  9.364
    6.819  5.091 12.532
    4.783  5.028  7.495  9.986
  -3.839 -3.889 -3.841 -3.625 9.610
  -2.189 -1.883 -2.175 -1.878 3.552 4.503
/diagram
title='Diagram from Manul2';
orientation=landscape;
border=yes;
tlocation=bottom;
layout=  v1 b v3 e3&
        v2 b v4 e4;
/end
/*********************************************
This example uses the same covariance matrix
from the example above.
The labels can be created in a data step
which we omitted here.
********************************************/
proc calis cov data=power method = ml nobs = 932;
      Lineqs
        v3 = a3 v1 + b3 v2 + e3,
        v4 = a4 v1 + b4 v2 + e4;
      Std 
          v1 - v2 = ev1 - ev2,
          e3 - e4 = ee3 - ee4;
      Cov                   
        v1 v2 = Theta (.7);
run;
Example 3: A Factor Analysis Model on page 29
博达大桥广告公司 
  EQS code: manul3.eqs
SAS code
/TITLE
  FACTOR ANALYSIS MODEL
  (EXAMPLE IN EQS MANUAL P.29)
/SPECIFICATIONS
  CAS=932; VAR=6; ME=GLS;
/LABEL
  V1=ANOMIA67; V2=POWRLS67;
  V3=ANOMIA71; V4=POWRLS71;
/EQUATIONS
  V1 = 2*F1 + E1;
  V2 = 2*F1 + E2;
  V3 = 2*F2 + E3;
  V4 = 2*F2 + E4;
/VARIANCES
  F1 TO F2 = 1.0;
  E1 TO E4 =  3*;
/COVARIANCES
  F2,F1= .3*;
  e3,e1 = *;
  e4,e2 = *;
/MATRIX
  11.834
    6.947  9.364
    6.819  5.091 12.532
    4.783  5.028  7.495  9.986
  -3.839 -3.889 -3.841 -3.625 9.610
  -2.189 -1.883 -2.175 -1.878 3.552 4.503
/constraints
  (v1,f1) = (v3,f2);
  (v2,f1) = (v4,f2);
  (e1,e1) = (e3,e3);
  (e2,e2) = (e4,e4);
  (e3,e1) = (e4,e2);
/DIAGRAM
title='Diagram from Manul3';
tlocation=bottom;
orientation=landscape;
border=yes;
postscript='MANUL3.PS';
layout=E1  b  E2  E3  b    E4 &
      V1  b  V2  V3  b    V4 &
      b    F1  b  b  F2;
/END
/**********************************************
This example uses the same covariance matrix
as in example 1.
*********************************************/
proc calis cov data=power method = gls nobs = 932;
lineqs
  v1 = a1 f1 + e1,
  v2 = a2 f1 + e2,
  v3 = a1 f2 + e3,
  v4 = a2 f2 + e4;
std
  f1- f2 = 1.0,
  e1 e3 = ee1,
  e2 e4 = ee2;
cov
  f1 f2 = theta1,
  e1 e3 = theta2,
  e2 e4 = theta2;
run;
Example 4: A Complete Latent Variable Model on page 33
 
成都理工学院学报
 EQS code: manul4.eqs
SAS code
/TITLE
  A COMPLETE LATENT VARIABLE MODEL
  (EXAMPLE IN EQS MANUAL P.33)
/SPECIFICATIONS
  CAS=932; VAR=6; ME=ML; MAT=COV;
  data='MANUL4.DAT';
/LABEL
  V1=ANOMIA67; V2=POWRLS67;
  V3=ANOMIA71; V4=POWRLS71;
  V5=EDUCATON; V6=OCCUPATN;
/EQUATIONS
  V1 =      F1 + E1;
  V2 = .833 F1 + E2;
  V3 =      F2 + E3;
  V4 = .833 F2 + E4;
  V5 =      F3 + E5;
  V6 =  .5 *F3 + E6;
  F1 = -.5 *F3 + D1;
  F2 =  .5 *F1 - .5*F3 + D2;
/VARIANCES
  D1 TO D2 = 4*;
  F3 = 6*;
  E1 TO E6 = 3*;
/COVARIANCES
  E1,E3 = .2*;
  E2,E4 = .2*;
/CONSTRAINTS
  (E1,E1) = (E3,E3);
  (E2,E2) = (E4,E4);
  (E3,E1) = (E4,E2);
/END
/*********************************************
  This example uses the same covariance matrix
  as in example 1.
  ********************************************/
 
fifa2003
proc calis cov data=power method = ml nobs = 932;
    lineqs
  V1 =      F1 + E1,
  V2 = .833 F1 + E2,
  V3 =      F2 + E3,
  V4 = .833 F2 + E4,
  V5 =      F3 + E5,
  V6 =  a6 F3 + E6,
  F1 =  c1 F3 + D1,
  F2 =  c2 F1 + c3 F3 + D2;
std
方位角
  D1 - D2 = ed:,
  F3 = ef3,
  E1 = ee1,
  e3 = ee1,
  e2 = ee2,
  e4 = ee2,
  e5 = ee3,
  e6 = ee4;
cov
    E1 E3 = theta1,
    E2 E4 = theta1;
run;
Example 5: A Second-order Factor Analysis Model on page 38 
EQS code: manul5.eqs
SAS code
/TITLE
  A SECOND-ORDER FACTOR ANALYSIS MODEL
  (EXAMPLE IN EQS MANUAL P.38)
/SPECIFICATIONS
  CAS=932; VAR=6; ME=GLS;
/LABEL
  V1=ANOMIA67; V2=POWRLS67;
  V3=ANOMIA71; V4=POWRLS71;
  V5=EDUCATON; V6=OCCUPATN;
/EQUATIONS
  V1 =      F1 + E1;
  V2 =    2*F1 + E2;
  V3 =      F2 + E3;
  V4 =    2*F2 + E4;
  F1 =    1*F3 + D1;
  F2 =    1*F3 + D2;
/VARIANCES
  F3 = 1;
android 开发环境  F3 = 1;
  D1 TO D2 = 1*;
  E1 TO E4 = 3*;
/CONSTRAINTS
  (F1,F3) = (F2,F3);
/MATRIX
  11.834
    6.947  9.364
    6.819  5.091 12.532
    4.783  5.028  7.495  9.986
  -3.839 -3.889 -3.841 -3.625 9.610
  -2.189 -1.883 -2.175 -1.878 3.552 4.503
/END
/*********************************************
This example uses the same covariance matrix
  as in example 1.
********************************************/
proc calis cov data=power method = gls nobs=932;
lineqs
  V1 =      F1 + E1,
  V2 =  c1 F1 + E2,
  V3 =      F2 + E3,
  V4 =  c2 F2 + E4,
  F1 =  c3 F3 + D1,
  F2 =  c3 F3 + D2;
std
  F3 = 1,
  D1 - D2 = ed:,
  E1 - E4 = ee:;
run;
Example 6: A Nonstandard Model on page 104 

本文发布于:2024-09-25 14:36:28，感谢您对本站的认可！

本文链接：https://www.17tex.com/xueshu/453348.html

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系，我们将在24小时内删除。

上一篇：模型嵌岩桩试验及数值分析
下一篇：时间序列预测步骤_建立时间序列预测模型的5个简单步骤
标签：广告公司   成都   环境   大桥
留言与评论（共有 0 条评论）
   
验证码：
推荐文章
排行榜
热门标签
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议