DATA: FILE is C:\Users\donnel59\Desktop\parcelsdec.dat;
!Model 4 in Donnellan et al. JRP
VARIABLE: NAMES ARE se89a se89b se89c se90a se90b se90c se91a se91b se91c
se92a se92b se92c se94a se94b se94c
se95a se95b se95c se97a se97b se97c se99a se99b se99c
se05a se05b se05c se07a se07b se07c;
USEVAR = se89a-se07c; missing is se89a-se07c (9999);
analysis: estimator = ml;
Model = NOCOVARIANCES;
model:
SE1 by se89a (e)
se89b (f)
se89c (g);
SE2 by se90a (e)
se90b (f)
se90c (g);
SE3 by se91a (e)
se91b (f)
se91c (g);
SE4 by se92a (e)
se92b (f)
se92c (g);
SE5 by se94a (e)
se94b (f)
se94c (g);
SE6 by se95a (e)
se95b (f)
se95c (g);
SE7 by se97a (e)
se97b (f)
se97c (g);
SE8 by se99a (e)
se99b (f)
se99c (g);
SE9 by se05a (e)
se05b (f)
se05c (g);
SE10 by se07a (e)
se07b (f)
se07c (g);
se90a with se89a; se90b with se89b; se90c with se89c;
se91a with se89a; se91b with se89b; se91c with se89c;
se92a with se89a; se92b with se89b; se92c with se89c;
se94a with se89a; se94b with se89b; se94c with se89c;
se95a with se89a; se95b with se89b; se95c with se89c;
se97a with se89a; se97b with se89b; se97c with se89c;
se99a with se89a; se99b with se89b; se99c with se89c;
se05a with se89a; se05b with se89b; se05c with se89c;
se07a with se89a; se07b with se89b; se07c with se89c;
se91a with se90a; se91b with se90b; se91c with se90c;
se92a with se90a; se92b with se90b; se92c with se90c;
se94a with se90a; se94b with se90b; se94c with se90c;
se95a with se90a; se95b with se90b; se95c with se90c;
se97a with se90a; se97b with se90b; se97c with se90c;
se99a with se90a; se99b with se90b; se99c with se90c;
se05a with se90a; se05b with se90b; se05c with se90c;
se07a with se90a; se07b with se90b; se07c with se90c;
se92a with se91a; se92b with se91b; se92c with se91c;
se94a with se91a; se94b with se91b; se94c with se91c;
se95a with se91a; se95b with se91b; se95c with se91c;
se97a with se91a; se97b with se91b; se97c with se91c;
se99a with se91a; se99b with se91b; se99c with se91c;
se05a with se91a; se05b with se91b; se05c with se91c;
se07a with se91a; se07b with se91b; se07c with se91c;
se94a with se92a; se94b with se92b; se94c with se92c;
se95a with se92a; se95b with se92b; se95c with se92c;
se97a with se92a; se97b with se92b; se97c with se92c;
se99a with se92a; se99b with se92b; se99c with se92c;
se05a with se92a; se05b with se92b; se05c with se92c;
se07a with se92a; se07b with se92b; se07c with se92c;
se95a with se94a; se95b with se94b; se95c with se94c;
se97a with se94a; se97b with se94b; se97c with se94c;
se99a with se94a; se99b with se94b; se99c with se94c;
se05a with se94a; se05b with se94b; se05c with se94c;
se07a with se94a; se07b with se94b; se07c with se94c;
se97a with se95a; se97b with se95b; se97c with se95c;
se99a with se95a; se99b with se95b; se99c with se95c;
se05a with se95a; se05b with se95b; se05c with se95c;
se07a with se95a; se07b with se95b; se07c with se95c;
se99a with se97a; se99b with se97b; se99c with se97c;
se05a with se97a; se05b with se97b; se05c with se97c;
se07a with se97a; se07b with se97b; se07c with se97c;
se05a with se99a; se05b with se99b; se05c with se99c;
se07a with se99a; se07b with se99b; se07c with se99c;
se07a with se05a; se07b with se05b; se07c with se05c;
Trait by se1-se10@1;
ar1 by se1@1;
ar2 by se2@1;
ar3 by se3@1;
ar4 by se4@1;
ar6 by se5@1;
ar7 by se6@1;
ar9 by se7@1;
ar11 by se8@1;
ar17 by se9@1;
ar19 by se10@1;
s1 by se1@1;
s2 by se2@1;
s3 by se3@1;
s4 by se4@1;
s6 by se5@1;
s7 by se6@1;
s9 by se7@1;
s11 by se8@1;
s17 by se9@1;
s19 by se10@1;
ar2 on ar1 *.5 (b1);
ar3 on ar2 *.5 (b2);
ar4 on ar3 *.5 (b3);
ar6 on ar4 *.5 (b4);
ar7 on ar6 *.5 (b5);
ar9 on ar7 *.5 (b6);
ar11 on ar9 *.5 (b7);
ar17 on ar11 *.5 (b8);
ar19 on ar17 *.5 (b9);
ar1 (v);
ar2 (v1);
ar3 (v2);
ar4 (v3);
ar6 (v4);
ar7 (v5);
ar9 (v6);
ar11 (v7);
ar17 (v8);
ar19 (v9);
s1-s19 (b);
se1-se10@0;
model constraint:
v1 = v - v*b1*b1;
v2 = v - v*b2*b2;
v3 = v - v*b3*b3;
v4 = v - v*b4*b4;
v5 = v - v*b5*b5;
v6 = v - v*b6*b6;
v7 = v - v*b7*b7;
v8 = v - v*b8*b8;
v9 = v - v*b9*b9;
new (k * 9.8 , bx * .8, b13t14,
b14t15,
b15t16,
b16t17,
b17t18,
b18t19,
b19t20,
b20t21,
b21t22,
b22t23,
b23t24,
b24t22,
b25t26,
b26t27,
b27t28,
b28t29,
b29t30,
b30t31,
b31t32,
b13t32);
b1=bx**(((Log(14.09-k)-log(13.18-k))/(Log(22.0-k)-log(21.0-k))));
b2=bx**(((Log(15.09-k)-log(14.09-k))/(Log(22.0-k)-log(21.0-k))));
b3=bx**(((Log(16.07-k)-log(15.09-k))/(Log(22.0-k)-log(21.0-k))));
b4=bx**(((Log(18.14-k)-log(16.07-k))/(Log(22.0-k)-log(21.0-k))));
b5=bx**(((Log(19.18-k)-log(18.14-k))/(Log(22.0-k)-log(21.0-k))));
b6=bx**(((Log(21.15-k)-log(19.18-k))/(Log(22.0-k)-log(21.0-k))));
b7=bx**(((Log(23.79-k)-log(21.15-k))/(Log(22.0-k)-log(21.0-k))));
b8=bx**(((Log(29.61-k)-log(23.79-k))/(Log(22.0-k)-log(21.0-k))));
b9=bx**(((Log(31.70-k)-log(29.61-k))/(Log(22.0-k)-log(21.0-k))));
b13t14 = bx**(((Log(14-k)-log(13-k))/(Log(22.0-k)-log(21.0-k))));
b14t15 = bx**(((Log(15-k)-log(14-k))/(Log(22.0-k)-log(21.0-k))));
b15t16 = bx**(((Log(16-k)-log(15-k))/(Log(22.0-k)-log(21.0-k))));
b16t17 = bx**(((Log(17-k)-log(16-k))/(Log(22.0-k)-log(21.0-k))));
b17t18 = bx**(((Log(18-k)-log(17-k))/(Log(22.0-k)-log(21.0-k))));
b18t19 = bx**(((Log(19-k)-log(18-k))/(Log(22.0-k)-log(21.0-k))));
b19t20 = bx**(((Log(20-k)-log(19-k))/(Log(22.0-k)-log(21.0-k))));
b20t21 = bx**(((Log(21-k)-log(20-k))/(Log(22.0-k)-log(21.0-k))));
b21t22 = bx**(((Log(22-k)-log(21-k))/(Log(22.0-k)-log(21.0-k))));
b22t23 = bx**(((Log(23-k)-log(22-k))/(Log(22.0-k)-log(21.0-k))));
b23t24 = bx**(((Log(24-k)-log(23-k))/(Log(22.0-k)-log(21.0-k))));
b24t22 = bx**(((Log(25-k)-log(24-k))/(Log(22.0-k)-log(21.0-k))));
b25t26 = bx**(((Log(26-k)-log(25-k))/(Log(22.0-k)-log(21.0-k))));
b26t27 = bx**(((Log(27-k)-log(26-k))/(Log(22.0-k)-log(21.0-k))));
b27t28 = bx**(((Log(28-k)-log(27-k))/(Log(22.0-k)-log(21.0-k))));
b28t29 = bx**(((Log(29-k)-log(28-k))/(Log(22.0-k)-log(21.0-k))));
b29t30 = bx**(((Log(30-k)-log(29-k))/(Log(22.0-k)-log(21.0-k))));
b30t31 = bx**(((Log(31-k)-log(30-k))/(Log(22.0-k)-log(21.0-k))));
b31t32 = bx**(((Log(32-k)-log(31-k))/(Log(22.0-k)-log(21.0-k))));
b13t32 = bx**(((Log(32-k)-log(13-k))/(Log(22.0-k)-log(21.0-k))));
output: tech4; stand;