! Phenotypic Simplex Model Neale & Cardon (1992) Example #define NVAR 6 G1: Simplex structure Data NInput_vars=6 NObservations=66 NGroups=1 Labels weight1 weight2 weight3 weight4 weight5 weight6 CMatrix 51.470 53.867 58.150 55.399 59.709 63.658 58.573 63.124 67.071 72.820 56.597 60.644 65.104 70.421 71.634 52.222 55.930 60.485 66.239 67.360 66.855 Begin Matrices; A Diag NVAR NVAR Free ! innovations B Lower NVAR NVAR Free ! beta transmissions U Diag NVAR NVAR Free ! measurement error variances I Iden NVAR NVAR ! identity matrix End Matrices; Specify B 0 101 0 0 102 0 0 0 103 0 0 0 0 104 0 0 0 0 0 105 0 EQ u 1 1 u 2 2 u 3 3 u 4 4 u 5 5 u 6 6 !Measurement error equated across occasions Begin Algebra; C = (I-B)~*(A*A')*(I-B)~'+U*U'; !Phenotypic covariance matrix D = (I-B)~*(A*A')*(I-B)~'; !Covariance due to simplex process E = U*U'; !Measurement error variance F = \d2v(C); !Phenotypic variance at each time point G = \d2v(D); !Variance due to simplex process at each time point H = \d2v(E); !Error variance J = \d2v(A*A'); !Innovation variance at each time point End Algebra; Labels Row F total_variance Labels Col F var1 var2 var3 var4 var5 var6 Labels Row G latent_variance Labels Col G var(eta1) var(eta2) var(eta3) var(eta4) var(eta5) var(eta6) Labels Row H measurement_error_variance Labels Col H error1 error2 error3 error4 error5 error6 Labels Row J innovation_variance Labels Col J innov1 innov2 innov3 innov4 innov5 innov6 Covariances C / Options rs Options Multiple Start 7 All End