Date:         Wed, 28 Jul 1999 14:01:50 -0400
Reply-To:     "Hamani, Elmaache" <elmaachh@ERE.UMONTREAL.CA>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         "Hamani, Elmaache" <elmaachh@ERE.UMONTREAL.CA>
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.