範 例 三： 選拔指數的制定
 
 
程式檔名稱 : linEAR3.htm

參考文獻   : 種豬選拔指數與本省豬種改良方向。池雙慶 。1980。
             中國畜牧協會會誌　9:55-69。
 
 
程式檔指令


PROC IML;
                                           隻日增重、飼料效率及背脂厚度的遺傳相關
rG={   1.00    -0.70       0.25,   
      -0.70     1.00       0.30,
       0.25     0.30       1.00};
                                           隻日增重、飼料效率及背脂厚度的表型相關
rP={   1.00    -0.50       0.25,
      -0.50     1.00       0.15,
       0.25     0.15       1.00};


sP={0.09, 0.26, 0.46};               外表型標準偏差(ADG,  FE,      BF)?
H2={0.30, 0.35, 0.50};               遺傳率( 隻日增重、飼料效率 及 背脂厚度)
sG=SQRT(sP#H2#sP);

V0=Hdir(sP, I(3))* rP * Hdir(sP, I(3));   表型變方與共變方矩陣
G0=Hdir(sG, I(3))* rG * Hdir(sG, I(3));   遺傳變方與共變方矩陣

                                                三個性狀的經濟加權值
W= { 573,                               
    -696,
    -475};

                    PRINT V0;PRINT G0;PRINT w;
b = INV(V0)*G0*w;   PRINT b;                    解出指數的 b 值

V_I = b`*V0*b;                   PRINT V_I;     指數的變方
V_H = w`*G0*w;                   PRINT V_H;     遺傳值的變方
R_IH  = SQRT (V_I/V_H);          PRINT R_IH;    指數準確度

DLT_A = G0*b*(1 / SQRT(V_I) );   PRINT DLT_A;   遺傳改進量
                                                      (需再乘以指數的選拔強度)   
b=b*(1/SQRT(V_I/(25*25)));                      調整指數的 SD=25, 
PRINT b;                                        新的 b 值(性狀間比重不變)

RUN;



輸出 3-1  表型、遺傳的變方與共變方矩陣(V0, G0) 與經濟加權值(w)

       V0
          0.0081   -0.0117   0.0104
         -0.0117    0.0676   0.0179
          0.0104    0.0179   0.2116
        
        
       G0
          0.0024    -0.0053   0.0040
         -0.0053     0.0237   0.0150
          0.0040     0.0150   0.1058
        
        
            W
          573
         -696
         -475
        
        
輸出 3-2  原始指數係數(b)與變方(V_I)及準確度(R_IH)
        
             b
        367.60301
        -258.2219
        -272.1032
        
        
              V_I
        23940.635
        
        
              V_H
        48106.128
                
                
             R_IH
        0.7054523
        

        
輸出 3-3  遺傳改進量(DLT_A, 需再乘以選拔強度 i) 及 調整指數的 SD=25 
           時新的 b 值(性狀間比重不變)
        
          DLT_A
         0.007582 
        -0.078492      *i
        -0.201585


           B          
          59.39518
         -41.72201    
         -43.96487