/*..................inH1bmdo.................*/ /* v indexu procenta slozek mleka */ /* ze souhrnneho genotypu vypusteno maso, vaha kravy ponechana */ /* dosazenim EH 0 */ /* dojitelnost a zdravi vypusteno z genotypu i kontroly uzitkovosti */ /* pouze hlavni index, */ /*................vypocet indexu pro H byky...... */ /* varianta 1 */ /*.....................................27.2.2003.VUZV..................*/ dm output 'clear'; dm log 'clear';title "Varianta 1 bez masa, dojitelnosti a odolnosti"; proc iml; /*.....reset print;....*/ /*...................vlastnosti v genotypu - na ktere se slechti ....*/ /*......vlastnosti v genotypu.*/ /* 1 mlecna plazma kg 2 tuk kg 3 bilkoviny kg ................................................... 4 dojitelnost kg/min 5 zanety vemene % 6 metabolické poruchy % ................................................... 7 mezidobi primy efekt dnu 8 mezidobi maternalni efekt dnu 9 prebihani jalovic primy efekt dnu 10 prebihani jalovic maternalni efekt dnu ................................................... 11 prubeh porodu primy efekt 0,01 tridy 12 prubeh porodu maternalni efekt 0,01 tridy 13 mrtve rozena telata primy efekt % 14 mrtve rozena telata maternalni efekt % ................................................... 15 vek pri prvem oteleni dnu 16 dlouhovekost laktaci(mezidobí) ................................................... 17 vaha krav kg 18 netto prirustek g/den 19 vyteznost % 20 zmasilost EUROP 0,01 tridy 21 protucneni 0,01 tridy 22 spotreba zivin MJ NE / kg prirustu ...................................................*/ /*.....EH - Kc na jednotku..................*/ eh = { -1.13 , 13.12 , 93.62 , 0 , 0 , 0 , -39 , -23.13 , -19 , -17.86 , -66 , -33 , -200 , -100 , -7 , 1000 , -10 , 0 , 0 , 0 , 0 , 0 }; /*.......vsechny EH ................*/ /*................vsechny EH .........*/ ehc = { -1.13 , 13.12 , 93.62 , 700 , -19 , -30 , -39 , -23.13 , -19 , -17.86 , -66 , -33 , -200 , -100 , -7 , 1000 , -10 , 9.76 , 160.1 , -5.19 , -1.08 , -163.4 }; /*.............. genetická smerodatná odchylka ........*/ Sg = {501 , 21.06 , 15.61 , 0.20 , 7.5 , 4 , 7 , 7 , 5.5 , 5.5 , 0.022 , 0.013 , 2.5 , 2.0 , 30 , 0.65 , 17.50 , 40 , 0.2 , 0.5 , 0.3 , 1.5 }; /*....... geneticke korelace vlastnosti v genotypu v % ....*/ Gkor = {100 68 89 40 30 15 -6 15 -6 6 10 -5 10 -5 10 -10 10 10 -10 -8 -5 -10 , 0 100 80 40 15 15 -6 15 -10 6 10 -5 10 -5 10 -10 10 10 -2 -12 -5 -10 , 0 0 100 40 31 15 -6 15 -10 6 10 -5 10 -5 10 -10 20 10 -9 -8 -5 -10 , 0 0 0 100 25 15 -6 5 -6 5 10 -5 10 -5 10 -10 10 10 -10 -8 -5 -10 , 0 0 0 0 100 27 5 15 5 20 -1 10 -1 10 10 -20 10 -10 -5 -5 -5 -5 , 0 0 0 0 0 100 4 15 4 20 -1 10 -1 10 10 -20 10 -10 -5 -5 -5 -5 , 0 0 0 0 0 0 100 -5 32 -5 5 -5 5 -5 -15 10 0 -10 10 -2 0 0 , 0 0 0 0 0 0 0 100 -5 10 -20 20 -5 19 10 -20 10 9 -10 -7 20 0 , 0 0 0 0 0 0 0 0 100 -5 5 -5 5 -5 -15 10 0 -10 10 -2 0 0 , 0 0 0 0 0 0 0 0 0 100 -20 20 -5 19 10 -20 10 9 -10 -7 20 0 , 0 0 0 0 0 0 0 0 0 0 100 -34 50 -30 -10 -20 10 9 2 -6 5 0 , 0 0 0 0 0 0 0 0 0 0 0 100 -10 40 15 -22 -10 -10 0 -10 20 0 , 0 0 0 0 0 0 0 0 0 0 0 0 100 -30 -10 -20 10 9 2 -6 5 0 , 0 0 0 0 0 0 0 0 0 0 0 0 0 100 15 -22 -10 -10 0 -10 20 0 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 -10 15 15 0 0 0 0 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 -10 -5 3 10 10 0 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 40 10 -10 -10 -40 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 30 -30 -20 -80 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 -60 -10 -20 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 20 20 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 10 , 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 100 }; /*.....geneticka kovariancni matice......................*/ d=diag(gkor); g=gkor+gkor`-d; /* print g; */ gkov=0.01*sg#g#sg`; /* print gkov; */ /*...........vlastnosti v indexu - plemenne hodnoty....*/ /* 1 PH MLE 2 PH t% 3 PH b% ................................................... 4 vl. plodnost PH 5 pl. dcer PH ................................................... 6 dojitelnost PH .................................................... 7 hranatost 8 velikost 9 sir. hrudniku 10 hloubka tela 11 sklon zade 12 sirka zade 13 post. zad. koncetiny 14 paznehty 15 pred upnuti vemene 16 vyska zad. upnuti 17 zavesny vaz 18 hloubka vemene 19 rozmistneni pred. struku 20 delka struku ..................................................... 21 mlecny charakter 22 kapacita 23 zad 24 koncetiny 25 vemeno ................................................... 26 vyska v krizi 27 obvod hrudi 28 vyska v kohoutku ...................................................*/ /*............smerodatná odchylka plemennych hodnot ........*/ Sph = { 416.516 , 0.21 , 0.10 , 2.519 , 1.522 , 3.161 , 1.126 , 1.183 , 1.601 , 1.493 , 1.462 , 1.519 , 1.730 , 1.723 , 1.445 , 1.502 , 1.460 , 1.637 , 1.622 , 1.490 , 0.777 , 1.285 , 1.325 , 1.533 , 1.166 , 1.189 , 1.370 , 1.235 }; /*.....rozptyl plemennych hodnot.....*/ s2ph=sph#sph; /*............geneticka smerodatná odchylka vlastnosti v indexu ........*/ SgI = {451.775 , 0.234787 , 0.111803 , 2.655259 , 2.090628 , 5.77117 , 1.45366 , 1.527246 , 2.066882 , 1.927455 , 1.887434 , 1.961021 , 2.23342 , 2.224383 , 1.865487 , 1.939074 , 1.884852 , 2.113358 , 2.093993 , 1.923582 , 1.003103 , 1.658928 , 1.710568 , 1.979094 , 1.5053 , 1.534992 , 1.768662 , 1.594378 }; /*..rozptyl geneticky vlastnosti v indexu s plemennymi hodnotami....*/ s2gi=sgi#sgi; /*.......spolehlivost dilcich plemennych hodnot.......*/ ms2gi=i(28)#s2gi; ims2gi=inv(ms2gi); r2ph=s2ph#vecdiag(ims2gi); r2ph=round(r2ph,0.01); /*....print sph sgi r2ph;......*/ /*....... korelace plemennych hodnot vlastnosti v indexu v %....*/ PHkor = { 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -39 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -45 60 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -6 2 -1 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -13 -6 -7 -15 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , 21 -3 -11 8 3 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , 15 -2 -11 10 -14 19 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , 1 8 -9 -14 4 7 33 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -1 -5 2 2 -7 6 -20 19 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -6 2 -6 -1 -5 7 -1 17 75 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -4 3 -8 -17 12 -4 11 38 -36 -1 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -10 8 -1 8 1 2 6 23 41 34 -8 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , 5 0 -4 1 -19 14 8 -28 -16 -7 -8 5 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , 7 -8 -4 -1 2 3 3 31 10 -2 6 -12 -50 100 0 0 0 0 0 0 0 0 0 0 0 0 0 0 , -10 2 -1 8 -15 9 -1 13 30 23 -5 11 -8 7 100 0 0 0 0 0 0 0 0 0 0 0 0 0 , 15 -14 -16 6 0 4 43 13 1 7 -2 1 -11 1 10 100 0 0 0 0 0 0 0 0 0 0 0 0 , 4 -3 -3 1 0 1 24 20 8 3 1 6 -20 18 9 20 100 0 0 0 0 0 0 0 0 0 0 0 , -10 0 9 1 -6 -5 8 26 -14 -30 -3 0 -6 6 41 10 21 100 0 0 0 0 0 0 0 0 0 0 , -3 -1 11 4 -2 11 8 6 14 9 -9 9 -1 -2 19 8 6 17 100 0 0 0 0 0 0 0 0 0 , 4 -1 -11 -7 -21 -15 10 2 -2 5 -1 -13 -13 12 10 4 17 2 -25 100 0 0 0 0 0 0 0 0 , 15 -10 -20 -3 -17 19 88 36 -14 2 7 -2 9 1 5 48 23 9 13 6 100 0 0 0 0 0 0 0 , 6 5 -6 -2 -8 7 9 45 77 87 3 41 -19 10 27 12 10 -20 9 5 18 100 0 0 0 0 0 0 , -7 4 0 12 3 6 4 17 46 38 -33 78 -10 1 18 20 19 2 14 12 10 47 100 0 0 0 0 0 , 10 -7 -20 5 6 0 17 32 8 5 0 -6 -71 66 10 26 26 7 -3 7 19 20 21 100 0 0 0 0 , 15 -9 -12 3 -5 8 41 27 12 12 -5 2 -28 13 42 76 37 31 17 12 50 25 28 43 100 0 0 0 , 1 8 -7 -14 4 6 33 99 19 17 38 23 -30 32 14 12 22 29 6 2 35 44 17 33 27 100 0 0 , -4 1 -4 -8 -5 2 -10 39 81 75 8 40 -16 10 21 -1 6 -11 10 -1 -5 84 38 7 7 39 100 0 , 2 8 -7 -12 4 7 34 95 21 20 35 25 -21 24 13 13 21 23 4 3 37 49 19 25 24 95 44 100 }; /*......kovariancni matice plemennych hodnot ...............*/ d=diag(phkor); ph=phkor+phkor`-d; /* print phkor; */ phkov=0.01*sph#ph#sph`; /* print phkov;*/ /* .geneticke korelace vlastnosti v indexu k vlastnostem v genotypu .*/ IGkor= { 99 71 90 15 45 30 -6 15 -10 15 10 -5 10 -5 17 -15 10 10 -10 -8 -5 -10 , -39 9 -25 -6 -26 10 7 9 -3 9 -2 -5 -2 -5 10 5 -10 20 15 0 15 -10 , -55 -21 -28 -8 -24 10 4 7 4 7 -2 -5 -2 -5 10 5 -10 14 20 0 5 -10 , 8 10 7 8 10 4 -100 5 -80 5 -5 5 -5 5 15 -10 0 10 -10 2 0 0 , 7 20 20 7 -20 6 6 -100 5 -10 20 -20 5 -19 -10 39 -10 -9 10 7 -20 0 , 15 10 15 100 25 20 -6 5 -6 5 10 -5 10 -5 18 -10 10 10 -10 -8 -5 -10 , 10 10 10 15 15 15 0 -14 0 -14 0 10 0 10 10 10 -10 -10 -10 -10 -10 0 , 5 5 5 7 0 -5 -14 4 -14 4 15 0 15 0 -5 -5 50 40 5 5 5 -10 , 11 10 11 6 0 -5 -2 7 -2 7 15 0 15 0 -5 -5 81 40 5 5 5 -10 , 11 10 11 6 0 -5 -2 7 -2 7 15 0 15 0 -5 -5 70 40 5 5 5 -10 , 8 -9 5 -4 8 7 -17 12 -17 12 0 -10 0 -10 -10 0 0 0 0 0 0 0 , -12 -5 -7 2 0 0 -8 -1 -8 -1 -10 10 -10 10 0 -10 40 20 10 5 5 -10 , -5 -5 -5 -5 -5 -5 1 -19 1 -19 -10 10 -10 10 0 0 0 0 0 0 0 0 , 6 -10 6 3 5 5 -1 2 -1 2 0 0 0 0 0 17 0 0 0 0 0 0 , -10 -10 -15 9 -20 -10 8 -15 8 -15 0 0 0 0 0 20 0 0 0 0 0 0 , 10 10 10 4 3 0 6 0 6 0 0 0 0 0 0 15 0 0 0 0 0 0 , 10 10 10 0 0 0 0 0 0 0 0 0 0 0 0 15 0 0 0 0 0 0 , -10 -10 -10 -5 -20 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 , 6 1 6 11 13 0 4 -2 4 -2 0 0 0 0 0 10 0 0 0 0 0 0 , 6 1 6 -15 -12 0 -5 -10 -5 -10 0 0 0 0 0 5 0 0 0 0 0 0 , 15 15 15 15 10 5 -3 -17 -3 -17 0 0 0 0 0 20 0 0 0 0 0 0 , 6 10 10 7 2 5 -2 -8 -2 -8 10 10 10 10 10 10 80 40 20 20 20 -20 , -7 8 1 6 -2 0 -12 -3 -12 -3 5 20 5 20 10 10 30 15 5 5 5 -10 , 15 6 16 -19 -11 -20 5 -13 5 -13 0 -10 0 -10 0 15 0 0 0 0 0 0 , 15 15 15 10 20 11 3 12 3 12 0 0 0 0 0 20 0 0 0 0 0 0 , 5 5 5 7 0 -5 -14 4 -14 4 15 0 15 0 -5 -5 50 40 5 5 5 -10 , 6 10 10 7 2 5 -2 -8 -2 -8 10 10 10 10 10 10 80 40 20 20 20 20 , 5 5 5 7 0 -5 -14 4 -14 4 15 0 15 0 -0 -5 50 40 5 5 5 -10 }; /*........kovariance plemennych hodnot k vlastnostem v genotypu........................*/ igkov=0.01*r2ph#sgi#igkor#sg`; /* print igkov; */ /*........ vypocet indexu ............................................................*/ /*......korekce matic - pronasobeno koeficientem...*/ igkov=igkov*1; phkov=phkov*1; /*......vyber bloku matic do indexu, totez jako vyber vlastnosti....*/ /*.................................................................................*/ /*varianta 1*/ phkov1=phkov[1:5,1:5]; igkov1=igkov[1:5,1:22]; phkov2=phkov[1:5,7:20]; phkov12=phkov1||phkov2; phkov3=phkov[7:20,1:5]; igkov2=igkov[7:20,1:22]; phkov4=phkov[7:20,7:20]; phkov34=phkov3||phkov4; phkov=phkov12//phkov34; igkov=igkov1//igkov2; /*...........vahy v indexu ....b.......*/ inphk=inv(phkov); b=inphk*igkov*eh; /*print b r2ph ; */ /*...........rozptyl souhrnneho genotypu..........*/ s2sg=eh`*gkov*eh; /*...podil vlastnosti v genotypu podle promenlivosti v %.....*/ pvg=0*eh; do i= 1 to 22 by 1; eei=eh; eei[i]=0; pvg[i]= s2sg-eei`*gkov*eei; end; svg=sum(pvg); pvg=100*pvg/svg; pvg=round(pvg,1); /*..........rozptyl indexu..............*/ s2ind1=b`*igkov*eh; s2ind2=b`*phkov*b; /* ....... vyznam vlastnosti v indexu.....*/ zz=0*b; do i= 1 to 19 by 1; zz[i]=100-100*((s2ind1-(b[i]*b[i]/inphk[i,i]))/s2ind1)##0.5; end; szz=sum(zz); zzp=100*zz/szz; zzp=round(zzp,1); /*.........smerodatna odchylka indexu....*/ sind=s2ind1##0.5; /*..........spolehlivost indexu.........*/ r2ind=s2ind1/s2sg; print s2sg svg s2ind1 /*s2ind2*/ sind r2ind ; /*........selekcni efekt u vlastnosti pri selekcnim rozdilu 1 smerodatna odchylka ......*/ dg=b`*igkov*(1/sind); dgt=dg`; /*........selekcni efekt v penezich......*/ pen=ehc#dgt; /*........efekt podlu skupin vlastnosti...*/ mlek = pen[1]+pen[2]+pen[3]; mas = pen[17]+pen[18]+pen[19]+pen[20]+pen[21]+pen[22]; druhot = pen[4]+pen[5]+pen[6]+pen[7]+pen[8]+pen[9]+pen[10] +pen[11]+pen[12]+pen[13]+pen[14]+pen[15]+pen[16]; /*.........efekt celkem .................*/ efcelk=mlek+mas+druhot; /*........podil efektu podle vlastnosti...*/ podpen=pen*(100/efcelk); podpen=round(podpen,1); print b /*r2ph*/ zzp dgt pen podpen pvg ; podmlek=mlek*100/efcelk; podmlek=round(podmlek,1); podmas=mas*100/efcelk; podmas=round(podmas,1); poddruh=druhot*100/efcelk; poddruh=round(poddruh,1); print mlek mas druhot efcelk podmlek podmas poddruh; /*...............konec.základniho indexu...............................*/ quit;