```Date: Wed, 28 Jul 1999 14:01:50 -0400 Reply-To: "Hamani, Elmaache" Sender: "SAS(r) Discussion" From: "Hamani, Elmaache" Subject: Can someone help? Content-Type: text/plain; charset="iso-8859-1" Greetings all, Somebody who used model ACE in genetics could help me to write the genetic equation: Yij= m + Cj +R_1j* MZj + R_2j*DZj + R_2ij(DZ) + e_ij With two constraints : 2*Var(R_1j)=Var(R_2j) and Var(R_1j)=Var(R_2ij(DZ) ) ? THANK YOU . I' v tried to use PROC MIXED and I solved the problem partially with LDATA and LIN(q): A part of the: Family twin sexe Y MZ DZ 1 1 0 13 0 1 1 2 1 16 0 1 2 1 0 17 1 0 2 2 0 17 1 0 3 1 0 16 1 0 3 2 0 17 1 0 4 1 0 14 0 1 4 2 0 17 0 1 5 1 0 12 1 0 5 2 0 14 1 0 6 2 0 13 0 1 6 1 1 13 0 1 7 1 1 22 0 1 7 2 1 13 0 1 ; Yij= m + Cj +R_1j* MZj + R_2j*DZj + R_2ij(DZ) + e_ij With two constraints : 1- at the 2 level (elements of matrix G): 2*Var(R_1j)=Var(R_2j) 2- relationship betwin the firet level (matrix R) and the second level (matrix G): Var(R_1j)=Var(R_2ij(DZ) ). Here R_2ij(DZ) is an element of R as like e_ij. I solved problem 1, but not the 2. Here my program and part of the output: data con; input parm row col1-col3; datalines; 1 1 1 0 0 2 2 0 1 0 2 3 0 0 2 ; run; PROC MIXED data=un noclprint covtest ; class famille ; model nmis_agr = / noint; random int dz mz /type=lin(2) ldata=con sub=famille G ; repeated / r ; run; ouput: R Matrix for FAMILLE 1001 Row COL1 1 0.53656972 G Matrix Effect FAMILY Row COL1 COL2 COL3 INTERCEPT 1001 1 0.26700391 DZ 1001 2 0.09326706 MZ 1001 0.18653412 Covariance Parameter Estimates (REML) Cov Parm Subject Estimate Std Error Z Pr > |Z| LIN(1) FAMILLE 0.26700391 0.15804857 1.69 0.0911 LIN(2) FAMILLE 0.09326706 0.11151294 0.84 0.4029 DIAG 0.53656972 0.04241956 12.65 0.0001 It misses here DIAG (DZ) =? for the seconde relationship:Var(R_1j)=Var(R_2ij(DZ) ). Thanks a lot. best regards. ```

Back to: Top of message | Previous page | Main SAS-L page