!  Power calculation for IBD statuses sibships size 4

#define nsib 4
#define nsibm1 3
#define nvar 1
#define nvarnsib 4

G1: Model parameters
Calc NGroups=37
Begin Matrices;
 U Unit nsib nsib       ! for model for C 
 I Iden nsib nsib       ! for model for E 
 W Lower nvar nvar      ! QTL path
 X Lower nvar nvar      ! A path
 Y Lower nvar nvar      ! C path
 Z Lower nvar nvar      ! E path
 H Stan nsib nsib       ! for A of sib pairs,
End Matrices;

Value .5 H 2 1 to H nsib nsibm1
 Matrix X 0.547722557505166
 Matrix Z 0.774596669241483
 Matrix W 0.316227766016838

Begin Algebra;
 A= X*X';
 C= Y*Y';
 E= Z*Z';
 Q= W*W';
End Algebra:
End


IBD 222222 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1 1 1 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=222222.cov
End 

IBD 222111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1  1 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=222111.cov
End 

IBD 222000 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1  1 0 0 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=222000.cov
End 


IBD 211221 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 1 1 .5
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=211221.cov
End 

IBD 211112 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 .5 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=211112.cov
End 

IBD 211110 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 .5 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=211110.cov
End 

IBD 211001 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 0 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=211001.cov
End 

IBD 200220 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 1 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=200220.cov
End 

IBD 200111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=200111.cov
End 

IBD 200002 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 0 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=200002.cov
End 

IBD 121212 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 1 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=121212.cov
End 

IBD 121121 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 .5 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=121121.cov
End 

IBD 121101 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 .5 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=121101.cov
End 

IBD 121010 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 0 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=121010.cov
End  

IBD 112211 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 1 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=112211.cov
End 

IBD 112122 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 .5 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=112122.cov
End 

IBD 112100 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 .5 0 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=112100.cov
End 

IBD 112011 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 0 .5 .5
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=112011.cov
End 

IBD 110211 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 1 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=110211.cov
End 

IBD 110120 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 .5 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=110120.cov
End 

IBD 110102 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 .5 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=110102.cov
End 

IBD 110011 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 0 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=110011.cov
End  

IBD 101210 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 1 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=101210.cov
End 

IBD 101121 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 .5 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=101121.cov
End 0

IBD 101101 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 .5 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=101101.cov
End  

IBD 101012 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 0 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=101012.cov
End 

IBD 020202 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 1 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=020202.cov
End 

IBD 020111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=020111.cov
End 

IBD 020020 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 0 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=020020.cov
End 

IBD 011201 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 1 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=011201.cov
End 

IBD 011112 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 .5 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=011112.cov
End 

IBD 011110 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 .5 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=011110.cov
End 

IBD 011021 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 0 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=011021.cov
End  

IBD 002200  pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 1 0 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=002200.cov
End 

IBD 002111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=002111.cov
End 

IBD 002022 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=dummy.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 0 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option mx%e=002022.cov
End 


!  Power calculation for IBD statuses sibships size 4
!  Fit false model

#define nsib 4
#define nsibm1 3
#define nvar 1
#define nvarnsib 4

G1: Model parameters
Data Calc NGroups=37
Begin Matrices;
 U Unit nsib nsib       ! for model for C 
 I Iden nsib nsib       ! for model for E 
 W Lower nvar nvar Free       ! QTL path
 X Lower nvar nvar Free       ! A path
 Y Lower nvar nvar            ! C path
 Z Lower nvar nvar Free       ! E path
 H Stan nsib nsib       ! for A of sib pairs,
End Matrices;

Value .5 H 1 2 to H nsib nsibm1

Begin Algebra;
 A= X*X';
 C= Y*Y';
 E= Z*Z';
 Q= W*W';
End Algebra:
End


IBD 222222 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=222222.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1  1 1 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 222111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=222111.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1  1 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 222000 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=222000.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 1  1 0 0 0
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 


IBD 211221 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=211221.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 1 1 .5
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 211112 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=211112.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 .5 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 211110 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=211110.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 .5 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 211001 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=211001.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 .5  .5 0 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 200220 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=200220.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 1 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 200111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=200111.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 200002 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=200002.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 1 0  0 0 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 121212 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=121212.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 1 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 121121 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=121121.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 .5 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 121101 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=121101.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 .5 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 121010 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=121010.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 1  .5 0 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End  

IBD 112211 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=112211.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 1 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 112122 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=112122.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 .5 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 112100 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=112100.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 .5 0 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 112011 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=112011.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  1 0 .5 .5
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 110211 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=110211.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 1 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 110120 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=110120.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 .5 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 110102 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=110102.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 .5 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 110011 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=110011.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 .5  0 0 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End  

IBD 101210 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=101210.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 1 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 101121 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=101121.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 .5 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 0

IBD 101101 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=101101.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 .5 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End  

IBD 101012 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=101012.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S .5 0  .5 0 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 020202 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=020202.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 1 0 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 020111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=020111.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 020020 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=020020.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 1  0 0 1 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 011201 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=011201.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 1 0 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 011112 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=011112.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 .5 .5 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 011110 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=011110.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 .5 .5 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 011021 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=011021.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 .5  .5 0 1 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End  

IBD 002200  pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=002200.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 1 0 0 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 002111 pairs
Data NInput_vars=nvarnsib NObservations=80
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=002111.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 .5 .5 .5 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
End 

IBD 002022 pairs
Data NInput_vars=nvarnsib NObservations=40
Labels sib1 sib2 sib3 sib4
CMatrix Symmetric File=002022.cov
Matrices= Group 1
 S Stan nsib nsib       ! for pihat        
End Matrices;
 Matrix S 0 0  1 0 1 1 
Cov H@A + U@C + I@E + S@Q /
Options RSidual
Option Multiple
End 

 Drop W 1 1 1
 Option power=2.0166445E-4 1
End