来自图书《MPlus中介调节模型》
模型方程如下:
Y = b0 + b1M + b2MV + b3MQ + b4MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程,得到:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b3(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
模型方程如下:
Y = b0 + b1M + b2MV + b3MQ + b4MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程,得到:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b3(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
展开括号,得到:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + a0b2V + a1b2XV + a2b2WV + a3b2ZV + a4b2XWV + a5b2XZV + a6b2WZV + a7b2XWZV + a0b3Q + a1b3XQ + a2b3WQ + a3b3ZQ + a4b3XWQ + a5b3XZQ + a6b3WZQ + a7b3XWZQ + a0b4VQ + a1b4XVQ + a2b4WVQ + a3b4ZVQ + a4b4XWVQ + a5b4XZVQ + a6b4WZVQ + a7b4XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
模型方程如下:
Y = b0 + b1M + b2MV + b3MQ + b4MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程,得到:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b3(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
展开括号,得到:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + a0b2V + a1b2XV + a2b2WV + a3b2ZV + a4b2XWV + a5b2XZV + a6b2WZV + a7b2XWZV + a0b3Q + a1b3XQ + a2b3WQ + a3b3ZQ + a4b3XWQ + a5b3XZQ + a6b3WZQ + a7b3XWZQ + a0b4VQ + a1b4XVQ + a2b4WVQ + a3b4ZVQ + a4b4XWVQ + a5b4XZVQ + a6b4WZVQ + a7b4XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
将包含X的项合并,整理成`Y = a + bX`的形式:
Y = (b0 + a0b1 + a2b1W + a3b1Z + a6b1WZ + a0b2V + a2b2WV + a3b2ZV + a6b2WZV + a0b3Q + a2b3WQ + a3b3ZQ + a6b3WZQ + a0b4VQ + a2b4WVQ + a3b4ZVQ + a6b4WZVQ + c2'W + c3'Z + c6'WZ + c8'V + c9'Q + c12'VQ) + (a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b2V + a4b2WV + a5b2ZV + a7b2WZV + a1b3Q + a4b3WQ + a5b3ZQ + a7b3WZQ + a1b4VQ + a4b4WVQ + a5b4ZVQ + a7b4WZVQ + c1' + c4'W + c5'Z + c7'WZ + c10'V + c11'Q + c13'VQ)X
模型方程如下:
Y = b0 + b1M + b2MV + b3MQ + b4MVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
M = a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ
将第二个方程代入第一个方程,得到:
Y = b0 + b1(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ) + b2(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)V + b3(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)Q + b4(a0 + a1X + a2W + a3Z + a4XW + a5XZ + a6WZ + a7XWZ)VQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
展开括号,得到:
Y = b0 + a0b1 + a1b1X + a2b1W + a3b1Z + a4b1XW + a5b1XZ + a6b1WZ + a7b1XWZ + a0b2V + a1b2XV + a2b2WV + a3b2ZV + a4b2XWV + a5b2XZV + a6b2WZV + a7b2XWZV + a0b3Q + a1b3XQ + a2b3WQ + a3b3ZQ + a4b3XWQ + a5b3XZQ + a6b3WZQ + a7b3XWZQ + a0b4VQ + a1b4XVQ + a2b4WVQ + a3b4ZVQ + a4b4XWVQ + a5b4XZVQ + a6b4WZVQ + a7b4XWZVQ + c1'X + c2'W + c3'Z + c4'XW + c5'XZ + c6'WZ + c7'XWZ + c8'V + c9'Q + c10'XV + c11'XQ + c12'VQ + c13'XVQ
将包含X的项合并,整理成`Y = a + bX`的形式:
Y = (b0 + a0b1 + a2b1W + a3b1Z + a6b1WZ + a0b2V + a2b2WV + a3b2ZV + a6b2WZV + a0b3Q + a2b3WQ + a3b3ZQ + a6b3WZQ + a0b4VQ + a2b4WVQ + a3b4ZVQ + a6b4WZVQ + c2'W + c3'Z + c6'WZ + c8'V + c9'Q + c12'VQ) + (a1b1 + a4b1W + a5b1Z + a7b1WZ + a1b2V + a4b2WV + a5b2ZV + a7b2WZV + a1b3Q + a4b3WQ + a5b3ZQ + a7b3WZQ + a1b4VQ + a4b4WVQ + a5b4ZVQ + a7b4WZVQ + c1' + c4'W + c5'Z + c7'WZ + c10'V + c11'Q + c13'VQ)X
因此,在给定W, Z, V, Q的条件下:
X对Y的间接效应:
(a1 + a4W + a5Z + a7WZ)(b1 + b2V + b3Q + b4VQ)
X对Y的直接效应:
c1' + c4'W + c5'Z + c7'WZ + c10'V + c11'Q + c13'VQ
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z + a1*b2*LOW_V + a4*b2*LOW_W*LOW_V + a5*b2*LOW_Z*LOW_V + a7*b2*LOW_W*LOW_Z*LOW_V + a1*b3*LOW_Q + a4*b3*LOW_W*LOW_Q + a5*b3*LOW_Z*LOW_Q + a7*b3*LOW_W*LOW_Z*LOW_Q + a1*b4*LOW_V*LOW_Q + a4*b4*LOW_W*LOW_V*LOW_Q + a5*b4*LOW_Z*LOW_V*LOW_Q + a7*b4*LOW_W*LOW_Z*LOW_V*LOW_Q;
...(其余的ILLLL到IHHHH, DLLLL到DHHHH, TLLLL到THHHH的定义)...
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z + a1*b2*LOW_V + a4*b2*LOW_W*LOW_V + a5*b2*LOW_Z*LOW_V + a7*b2*LOW_W*LOW_Z*LOW_V + a1*b3*LOW_Q + a4*b3*LOW_W*LOW_Q + a5*b3*LOW_Z*LOW_Q + a7*b3*LOW_W*LOW_Z*LOW_Q + a1*b4*LOW_V*LOW_Q + a4*b4*LOW_W*LOW_V*LOW_Q + a5*b4*LOW_Z*LOW_V*LOW_Q + a7*b4*LOW_W*LOW_Z*LOW_V*LOW_Q;
...(其余的ILLLL到IHHHH, DLLLL到DHHHH, TLLLL到THHHH的定义)...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z + a1*b2*LOW_V + a4*b2*LOW_W*LOW_V + a5*b2*LOW_Z*LOW_V + a7*b2*LOW_W*LOW_Z*LOW_V + a1*b3*LOW_Q + a4*b3*LOW_W*LOW_Q + a5*b3*LOW_Z*LOW_Q + a7*b3*LOW_W*LOW_Z*LOW_Q + a1*b4*LOW_V*LOW_Q + a4*b4*LOW_W*LOW_V*LOW_Q + a5*b4*LOW_Z*LOW_V*LOW_Q + a7*b4*LOW_W*LOW_Z*LOW_V*LOW_Q;
...(其余的ILLLL到IHHHH, DLLLL到DHHHH, TLLLL到THHHH的定义)...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
...(PLLLL到PHHHH的定义)...
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z + a1*b2*LOW_V + a4*b2*LOW_W*LOW_V + a5*b2*LOW_Z*LOW_V + a7*b2*LOW_W*LOW_Z*LOW_V + a1*b3*LOW_Q + a4*b3*LOW_W*LOW_Q + a5*b3*LOW_Z*LOW_Q + a7*b3*LOW_W*LOW_Z*LOW_Q + a1*b4*LOW_V*LOW_Q + a4*b4*LOW_W*LOW_V*LOW_Q + a5*b4*LOW_Z*LOW_V*LOW_Q + a7*b4*LOW_W*LOW_Z*LOW_V*LOW_Q;
...(其余的ILLLL到IHHHH, DLLLL到DHHHH, TLLLL到THHHH的定义)...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
...(PLLLL到PHHHH的定义)...
PLOT:
TYPE = plot2;
USEVARIABLES = X M W Z V Q Y XW XZ WZ XV XQ VQ MV MQ XWZ XVQ MVQ;
DEFINE:
MQ = M*Q;
MV = M*V;
XW = X*W;
XZ = X*Z;
XQ = X*Q;
XV = X*V;
WZ = W*Z;
VQ = V*Q;
MVQ = M*V*Q;
XWZ = X*W*Z;
XVQ = X*V*Q;
ANALYSIS:
TYPE = GENERAL;
ESTIMATOR = ML;
BOOTSTRAP = 10000;
MODEL:
[Y] (b0);
Y ON M (b1);
Y ON MV (b2);
Y ON MQ (b3);
Y ON MVQ (b4);
Y ON X (cdash1);
Y ON W (cdash2);
Y ON Z (cdash3);
Y ON XW (cdash4);
Y ON XZ (cdash5);
Y ON WZ (cdash6);
Y ON XWZ (cdash7);
Y ON V (cdash8);
Y ON Q (cdash9);
Y ON XV (cdash10);
Y ON XQ (cdash11);
Y ON VQ (cdash12);
Y ON XVQ (cdash13);
[M] (a0);
M ON X (a1);
M ON W (a2);
M ON Z (a3);
M ON XW (a4);
M ON XZ (a5);
M ON WZ (a6);
M ON XWZ (a7);
MODEL CONSTRAINT:
NEW(LOW_W MED_W HIGH_W LOW_Z MED_Z HIGH_Z LOW_V MED_V HIGH_V LOW_Q MED_Q HIGH_Q
ILLLL IMLLL IHLLL ILMLL IMMLL IHMLL ILHLL IMHLL IHHLL
ILLML IMLML IHLML ILMML IMMML IHMML ILHML IMHML IHHML
ILLHL IMLHL IHLHL ILMHL IMMHL IHMHL ILHHL IMHHL IHHHL
ILLLM IMLLM IHLLM ILMLM IMMLM IHMLM ILHLM IMHLM IHHLM
ILLMM IMLMM IHLMM ILMMM IMMMM IHMMM ILHMM IMHMM IHHMM
ILLHM IMLHM IHLHM ILMHM IMMHM IHMHM ILHHM IMHHM IHHHM
ILLLH IMLLH IHLLH ILMLH IMMLH IHMLH ILHLH IMHLH IHHLH
ILLMH IMLMH IHLMH ILMMH IMMMH IHMMH ILHMH IMHMH IHHMH
ILLHH IMLHH IHLHH ILMHH IMMHH IHMHH ILHHH IMHHH IHHHH
DLLLL DMLLL DHLLL DLMLL DMMLL DHMLL DLHLL DMHLL DHHLL
DLLML DMLML DHLML DLMML DMMML DHMML DLHML DMHML DHHML
DLLHL DMLHL DHLHL DLMHL DMMHL DHMHL DLHHL DMHHL DHHHL
DLLLM DMLLM DHLLM DLMLM DMMLM DHMLM DLHLM DMHLM DHHLM
DLLMM DMLMM DHLMM DLMMM DMMMM DHMMM DLHMM DMHMM DHHMM
DLLHM DMLHM DHLHM DLMHM DMMHM DHMHM DLHHM DMHHM DHHHM
DLLLH DMLLH DHLLH DLMLH DMMLH DHMLH DLHLH DMHLH DHHLH
DLLMH DMLMH DHLMH DLMMH DMMMH DHMMH DLHMH DMHMH DHHMH
DLLHH DMLHH DHLHH DLMHH DMMHH DHMHH DLHHH DMHHH DHHHH
TLLLL TMLLL THLLL TLMLL TMMLL THMLL TLHLL TMHLL THHLL
TLLML TMLML THLML TLMML TMMML THMML TLHML TMHML THHML
TLLHL TMLHL THLHL TLMHL TMMHL THMHL TLHHL TMHHL THHHL
TLLLM TMLLM THLLM TLMLM TMMLM THMLM TLHLM TMHLM THHLM
TLLMM TMLMM THLMM TLMMM TMMMM THMMM TLHMM TMHMM THHMM
TLLHM TMLHM THLHM TLMHM TMMHM THMHM TLHHM TMHHM THHHM
TLLLH TMLLH THLLH TLMLH TMMLH THMLH TLHLH TMHLH THHLH
TLLMH TMLMH THLMH TLMMH TMMMH THMMH TLHMH TMHMH THHMH
TLLHH TMLHH THLHH TLMHH TMMHH THMHH TLHHH TMHHH THHHH);
LOW_W = #LOWW;
MED_W = #MEDW;
HIGH_W = #HIGHW;
LOW_Z = #LOWZ;
MED_Z = #MEDZ;
HIGH_Z = #HIGHZ;
LOW_V = #LOWV;
MED_V = #MEDV;
HIGH_V = #HIGHV;
LOW_Q = #LOWQ;
MED_Q = #MEDQ;
HIGH_Q = #HIGHQ;
ILLLL = a1*b1 + a4*b1*LOW_W + a5*b1*LOW_Z + a7*b1*LOW_W*LOW_Z + a1*b2*LOW_V + a4*b2*LOW_W*LOW_V + a5*b2*LOW_Z*LOW_V + a7*b2*LOW_W*LOW_Z*LOW_V + a1*b3*LOW_Q + a4*b3*LOW_W*LOW_Q + a5*b3*LOW_Z*LOW_Q + a7*b3*LOW_W*LOW_Z*LOW_Q + a1*b4*LOW_V*LOW_Q + a4*b4*LOW_W*LOW_V*LOW_Q + a5*b4*LOW_Z*LOW_V*LOW_Q + a7*b4*LOW_W*LOW_Z*LOW_V*LOW_Q;
...(其余的ILLLL到IHHHH, DLLLL到DHHHH, TLLLL到THHHH的定义)...
PLOT(PLLLL PMLLL PHLLL PLMLL PMMLL PHMLL PLHLL PMHLL PHHLL
PLLML PMLML PHLML PLMML PMMML PHMML PLHML PMHML PHHML
PLLHL PMLHL PHLHL PLMHL PMMHL PHMHL PLHHL PMHHL PHHHL
PLLLM PMLLM PHLLM PLMLM PMMLM PHMLM PLHLM PMHLM PHHLM
PLLMM PMLMM PHLMM PLMMM PMMMM PHMMM PLHMM PMHMM PHHMM
PLLHM PMLHM PHLHM PLMHM PMMHM PHMHM PLHHM PMHHM PHHHM
PLLLH PMLLH PHLLH PLMLH PMMLH PHMLH PLHLH PMHLH PHHLH
PLLMH PMLMH PHLMH PLMMH PMMMH PHMMH PLHMH PMHMH PHHMH
PLLHH PMLHH PHLHH PLMHH PMMHH PHMHH PLHHH PMHHH PHHHH);
LOOP(XVAL,1,5,0.1);
...(PLLLL到PHHHH的定义)...
PLOT:
TYPE = plot2;
OUTPUT:
STAND CINT(bcbootstrap);