/*priklad na TDM Data Rodokmen cow parity testday dim milk fat protein animal sire dam A 1 a 10 12 0.42 0.36 S 0 0 A 1 b 40 23 0.805 0.69 D1 0 0 A 1 c 70 34 1.19 1.02 D2 0 0 A 2 d 20 13 0.455 0.39 D3 0 0 A 2 e 30 32 1.12 0.96 A S D1 B 1 a 80 21 0.735 0.63 B S D2 B 2 b 40 14 0.49 0.42 C S D3 B 2 c 70 25 0.875 0.75 B 3 d 10 21 0.735 0.63 B 3 e 40 28 0.98 0.84 C 2 a 50 14 0.49 0.42 C 2 b 70 11 0.385 0.33 C 3 e 60 25 0.875 0.75 model y=Hc+Xb+Wp+Za+e y - milk c- testdays b - fixed regression p - permanent environment a - animal effect e - residuum MME H - matice spojuje skupiny HTD s u§itkovost¡ HTD je v tomto pý¡padØ dan‚ test-dnem cow/parity a1 b2 c3 b1 c2 d2 e2 d3 e3 A1 1 1 1 0 0 0 0 0 0 B1 1 0 0 0 0 0 0 0 0 A2 0 0 0 0 0 1 1 0 0 B2 0 0 0 1 1 0 0 0 0 C2 0 0 0 1 0 0 0 0 0 B3 0 0 0 0 0 0 0 1 0 C3 0 0 0 0 0 0 0 1 1 X,W,Z matice regresnich promennych zahrnuj¡c¡ch pocet dim spojenych s kravou na danem DIM a odpovidaj¡m pozorovani cow/parity intercept DIM exp(-0.05DIM) A1 1 10 0.60653066 A1 1 40 0.135335283 A1 1 70 0.030197383 A2 1 20 0.367879441 A2 1 30 0.22313016 B1 1 80 0.018315639 B2 1 40 0.135335283 B2 1 70 0.030197383 B3 1 10 0.60653066 B3 1 40 0.135335283 C2 1 50 0.082084999 C2 1 70 0.030197383 C3 1 60 0.049787068 Matice R - residu ln¡ matice mezi znaky pro ka§dì DIM a znak a poýad¡ laktace milk 1 var e 5 milk 2 var e 7.5 milk 3 var e 8.5 Matice P kovariancn¡ matice PE 20.55 -0.06 -22.64 0.2 0.02 0.5 0.5 0.4 0.6 -0.06 0.003 0.05 0.2 0.02 0.5 0.5 0.4 0.6 -22.64 0.05 59.38 0.2 0.02 0.5 0.5 0.4 0.6 0.2 0.2 0.2 40.27 -0.15 -22.64 0.5 0.4 0.5 0.02 0.02 0.02 -0.15 0.003 0.05 0.5 0.4 0.5 0.5 0.5 0.5 -22.64 0.05 59.38 0.5 0.4 0.5 0.5 0.5 0.5 0.5 0.5 0.5 47.03 -0.06 -22.64 0.4 0.4 0.4 0.4 0.4 0.4 -0.06 0.003 0.05 0.6 0.6 0.6 0.5 0.5 0.5 -22.64 0.05 59.38 Matice G kovariancn¡ matice genetickych efektu 20.55 -0.06 -22.64 0.2 0.02 0.5 0.5 0.4 0.6 -0.06 0.003 0.05 0.2 0.02 0.5 0.5 0.4 0.6 -22.64 0.05 59.38 0.2 0.02 0.5 0.5 0.4 0.6 0.2 0.2 0.2 40.27 -0.15 -22.64 0.5 0.4 0.5 0.02 0.02 0.02 -0.15 0.003 0.05 0.5 0.4 0.5 0.5 0.5 0.5 -22.64 0.05 59.38 0.5 0.4 0.5 0.5 0.5 0.5 0.5 0.5 0.5 47.03 -0.06 -22.64 0.4 0.4 0.4 0.4 0.4 0.4 -0.06 0.003 0.05 0.6 0.6 0.6 0.5 0.5 0.5 -22.64 0.05 59.38 Pý¡buzensk  matice A S D1 D2 D3 A B C S 1 0 0 0 0.5 0.5 0.5 D1 0 1 0 0 0.5 0 0 D2 0 0 1 0 0 0.5 0 D3 0 0 0 1 0 0 0.5 A 0.5 0.5 0 0 1 0.5 0.5 B 0.5 0 0.5 0 0.5 1 0.5 C 0.5 0 0 0.5 0.5 0.5 1 vektor b odhadovane efekty testdays a1 b1 c1 d1 e1 f1 a2 b2 c2 d2 e2 f2 a3 b3 c3 d3 e3 f3 FRC b0_1 b1_1 b2_1 b0_2 b1_2 b2_2 b0_3 b1_3 b2_3 ` PERRC p0_1_A p1_1_A p2_1_A p0_2_A p1_2_A p2_2_A p0_3_A p1_3_A p2_3_A p0_1_B p1_1_B p2_1_B p0_2_B p1_2_B p2_2_B p0_3_B p1_3_B p2_3_B p0_1_C p1_1_C p2_1_C p0_2_C p1_2_C p2_2_C p0_3_C p1_3_C p2_3_C ARRC a0_1_S a1_1_S a2_1_S a0_2_S a1_2_S a2_2_S a0_3_S a1_3_S a2_3_S a0_1_D1 a1_1_D1 a2_1_D1 a0_2_D1 a1_2_D1 a2_2_D1 a0_3_D1 a1_3_D1 a2_3_D1 a0_1_D2 a1_1_D2 a2_1_D2 a0_2_D2 a1_2_D2 a2_2_D2 a0_3_D2 a1_3_D2 a2_3_D2 a0_1_D3 a1_1_D3 a2_1_D3 a0_2_D3 a1_2_D3 a2_2_D3 a0_3_D3 a1_3_D3 a2_3_D3 a0_1_A a1_1_A a2_1_A a0_2_A a1_2_A a2_2_A a0_3_A a1_3_A a2_3_A a0_1_B a1_1_B a2_1_B a0_2_B a1_2_B a2_2_B a0_3_B a1_3_B a2_3_B a0_1_C a1_1_C a2_1_C a0_2_C a1_2_C a2_2_C a0_3_C a1_3_C a2_3_C */; Proc Iml; * a1 b1 c1 a2 b2 c2 d2 e2 d3 e3; H={ 1 0 0 0 0 0 0 0 0 0, 0 1 0 0 0 0 0 0 0 0, 0 0 1 0 0 0 0 0 0 0, 0 0 0 0 0 0 1 0 0 0, 0 0 0 0 0 0 0 1 0 0, 1 0 0 0 0 0 0 0 0 0, 0 0 0 0 1 0 0 0 0 0, 0 0 0 0 0 1 0 0 0 0, 0 0 0 0 0 0 0 0 1 0, 0 0 0 0 0 0 0 0 0 1, 0 0 0 1 0 0 0 0 0 0, 0 0 0 0 1 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 1}; * b0_1 b1_1 b2_1 b0_1 b1_1 b2_1 b0_2 b1_2 b3_2 b0_3 b1_3 b3_3; X={ 1 10 0.60653066 0 0 0 0 0 0 , 1 40 0.135335283 0 0 0 0 0 0, 1 70 0.030197383 0 0 0 0 0 0, 0 0 0 1 20 0.367879441 0 0 0, 0 0 0 1 30 0.22313016 0 0 0, 1 80 0.018315639 0 0 0 0 0 0, 0 0 0 1 40 0.135335283 0 0 0, 0 0 0 1 70 0.030197383 0 0 0, 0 0 0 0 0 0 1 10 0.60653066 , 0 0 0 0 0 0 1 40 0.135335283, 0 0 0 1 50 0.082084999 0 0 0, 0 0 0 1 70 0.030197383 0 0 0, 0 0 0 0 0 0 1 60 0.049787068 }; * Ap0_1 Ap1_1 Ap2_1 Ap0_2 Ap1_2 Ap2_2 Ap0_3 Ap1_3 Ap3_3 Bp0_1 Bp1_1 Bp2_1 Bp0_2 Bp1_2 Bp2_2 Bp0_3 Bp1_3 Bp3_3 Cp0_1 Cp1_1 Cp2_1 Cp0_2 Cp1_2 Cp2_2 Cp0_3 Cp1_3 Cp3_3; w= { 1 10 0.60 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 1 40 0.13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 1 70 0.03 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 1 20 0.36 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 1 30 0.22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 1 80 0.018 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 1 40 0.13 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 1 70 0.03 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 10 0.61 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 40 0.14 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 50 0.082 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 70 0.03 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 60 0.05} ; * Sp0_1 Sp1_1 Sp2_1 Sp0_2 Sp1_2 Sp2_2 Sp0_3 Sp1_3 Sp3_3 D1p0_1 D1p1_1 D1p2_1 D1p0_2 D1p1_2D1Ap2D12 Ap0_3 AD11_3 Ap3_3 D2p0_1 D2p1_1 D2p2_1 D2p0_2 D2p1_2D2Bp2_D2 BpD2_3 D2p1_3 Bp3_3 D3p0_1 D3p1_1 D3p2_1 D3p0_2 D3p1_2 Cp2_2D3Cp0D33 CD31_3 Cp3_3 Ap0_1 Ap1_1 Ap2_1 Ap0_2 Ap1_2 Ap2_2 Ap0_3 Ap1_3 Ap3_3 Bp0_1 Bp1_1 Bp2_1 Bp0_2 Bp1_2 Bp2_2 Bp0_3 Bp1_3 Bp3_3 Cp0_1 Cp1_1 Cp2_1 Cp0_2 Cp1_2 Cp2_2 Cp0_3 Cp1_3 Cp3_3; z={ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 10 0.60 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 40 0.13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 70 0.03 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 20 0.36 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 30 0.22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 80 0.018 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 40 0.13 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 70 0.03 0 0 0 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 10 0.61 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 40 0.14 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 50 0.082 0 0 0 0 0 0 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 70 0.03 0 0 0, 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 60 0.05}; A={ 1 0 0 0 0.5 0.5 0.5, 0 1 0 0 0.5 0 0 , 0 0 1 0 0 0.5 0 , 0 0 0 1 0 0 0.5, 0.5 0.5 0 0 1 0.5 0.5, 0.5 0 0.5 0 0.5 1 0.5, 0.5 0 0 0.5 0.5 0.5 1 }; P={20.55 -0.06 -22.64 0.2 0.02 0.5 0.5 0.4 0.6 , -0.06 0.003 0.05 0.2 0.02 0.5 0.5 0.4 0.6 , -22.64 0.05 59.38 0.2 0.02 0.5 0.5 0.4 0.6 , 0.2 0.2 0.2 40.27 -0.15 -22.64 0.5 0.4 0.5 , 0.02 0.02 0.02 -0.15 0.003 0.05 0.5 0.4 0.5 , 0.5 0.5 0.5 -22.64 0.05 59.38 0.5 0.4 0.5 , 0.5 0.5 0.5 0.5 0.5 0.5 47.03 -0.06 -22.64, 0.4 0.4 0.4 0.4 0.4 0.4 -0.06 0.003 0.05 , 0.6 0.6 0.6 0.5 0.5 0.5 -22.64 0.05 59.38 }; I3=i(3); p0=p@i3; pi=INV(p0); G={20.55 -0.06 -22.64 0.2 0.02 0.5 0.5 0.4 0.6 , -0.06 0.003 0.05 0.2 0.02 0.5 0.5 0.4 0.6 , -22.64 0.05 59.38 0.2 0.02 0.5 0.5 0.4 0.6 , 0.2 0.2 0.2 40.27 -0.15 -22.64 0.5 0.4 0.5 , 0.02 0.02 0.02 -0.15 0.003 0.05 0.5 0.4 0.5 , 0.5 0.5 0.5 -22.64 0.05 59.38 0.5 0.4 0.5 , 0.5 0.5 0.5 0.5 0.5 0.5 47.03 -0.06 -22.64, 0.4 0.4 0.4 0.4 0.4 0.4 -0.06 0.003 0.05 , 0.6 0.6 0.6 0.5 0.5 0.5 -22.64 0.05 59.38 }; G0=a@g; *PRINT g0I; g0i=INV(G0); R={ 5 0 0 0 0 0 0 0 0 0 0 0 0, 0 5 0 0 0 0 0 0 0 0 0 0 0, 0 0 5 0 0 0 0 0 0 0 0 0 0, 0 0 0 7.5 0 0 0 0 0 0 0 0 0, 0 0 0 0 7.5 0 0 0 0 0 0 0 0, 0 0 0 0 0 5 0 0 0 0 0 0 0, 0 0 0 0 0 0 7.5 0 0 0 0 0 0, 0 0 0 0 0 0 0 7.5 0 0 0 0 0, 0 0 0 0 0 0 0 0 8.5 0 0 0 0, 0 0 0 0 0 0 0 0 0 8.5 0 0 0, 0 0 0 0 0 0 0 0 0 0 7.5 0 0, 0 0 0 0 0 0 0 0 0 0 0 7.5 0, 0 0 0 0 0 0 0 0 0 0 0 0 8.5}; RI=inv(R); *print RI; *LHS; d11=H`* RI*H; d12=H`* RI*X; d13=H`* RI*W; d14=H`* RI*Z; d21=X`* RI*H; d22=X`* RI*X; d23=X`* RI*W; d24=X`* RI*Z; d31=W`* RI*H; d32=W`* RI*X; d33=(W`* RI*W)+INV(PI); d34=W`* RI*Z; d41=Z`* RI*H; d42=Z`* RI*X; d43=Z`* RI*W; d44=(z`* RI*Z)+G0I; *RHS; y={12,23,34,13,32,21,14,25,21,28,14,11,25}; y11=H`* RI*y; y22=X`* RI*y; y33=(W`* RI*y); y44=(z`* RI*y); print d11 d22 d33 d44; print y11 y22 y33 y44;