= Example of gebvtest.log = {{{ 2013-Jan-25 09:15:24 gebvtest.py: start gebvtest.py version=2013.01.24 absolute DATADIR: /fs6c/service/jette/GTEST/sample_data Processing brd=HOL pop=ABC cou=ABC datadir=../sample_data stored 426 records from file file732 stored 426 records from file bdate_ABCHOL ------------------------------------------------------------------------------- processing trt=mil h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N ------------------------------------------------------------------------------- stored 379 records for trait mil from file file300Cf_ABCHOL stored 290 records for trait mil from file file300Df_ABCHOL stored 403 records for trait mil from file file300Cr_ABCHOL stored 349 records for trait mil from file file300Gr_ABCHOL Summary statistics on candidate bulls (CB) and test bulls (TB) -------------------------------------------------------------- Trait Variable N Mean Std Min Max -------------------------------------------------------------- mil CB EBV 177 573.367 813.411 -1908.00 2962.00 mil TB EBV 164 544.957 804.675 -1908.00 2962.00 mil TB DPGM(y) 164 513.530 778.346 -1874.00 2608.00 mil TB GEBV(x1) 164 1047.909 599.771 -726.00 2624.00 mil TB EBVr(x2) 164 1244.341 548.582 -564.00 2616.00 -------------------------------------------------------------- Details of GEBVtest calculations mil i_est = (544.957 - 573.367) / 813.411 = -0.035 mil p=0.987 x=-2.213 i=0.035 k=0.079 R2b=0.642 E(b1)=0.970 mil b1=1.022 se=0.063 E(b1)=0.970 t=0.83 R2_1=62.3 R2_2=28.8 mil passes t-test=Y bio-test=Y R2-test=Y ------------------------------------------------------------------------------- processing trt=fat h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N ------------------------------------------------------------------------------- stored 379 records for trait fat from file file300Cf_ABCHOL stored 290 records for trait fat from file file300Df_ABCHOL stored 403 records for trait fat from file file300Cr_ABCHOL stored 349 records for trait fat from file file300Gr_ABCHOL Summary statistics on candidate bulls (CB) and test bulls (TB) -------------------------------------------------------------- Trait Variable N Mean Std Min Max -------------------------------------------------------------- fat CB EBV 176 24.284 29.839 -65.00 95.00 fat TB EBV 164 22.939 28.980 -65.00 95.00 fat TB DPGM(y) 164 21.799 29.369 -68.00 98.00 fat TB GEBV(x1) 164 37.720 20.138 -19.00 86.00 fat TB EBVr(x2) 164 44.335 20.490 -4.00 88.00 -------------------------------------------------------------- Details of GEBVtest calculations fat i_est = (22.939 - 24.284) / 29.839 = -0.045 fat p=0.982 x=-2.097 i=0.045 k=0.097 R2b=0.543 E(b1)=0.953 fat b1=1.047 se=0.079 E(b1)=0.953 t=1.18 R2_1=51.8 R2_2=24.2 fat passes t-test=Y bio-test=Y R2-test=Y ------------------------------------------------------------------------------- processing trt=pro h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N ------------------------------------------------------------------------------- stored 379 records for trait pro from file file300Cf_ABCHOL stored 290 records for trait pro from file file300Df_ABCHOL stored 403 records for trait pro from file file300Cr_ABCHOL stored 349 records for trait pro from file file300Gr_ABCHOL Summary statistics on candidate bulls (CB) and test bulls (TB) -------------------------------------------------------------- Trait Variable N Mean Std Min Max -------------------------------------------------------------- pro CB EBV 176 20.926 22.004 -50.00 89.00 pro TB EBV 164 20.335 21.066 -50.00 89.00 pro TB DPGM(y) 164 19.110 20.544 -51.00 81.00 pro TB GEBV(x1) 164 35.177 15.715 -12.00 72.00 pro TB EBVr(x2) 164 40.768 15.368 -8.00 76.00 -------------------------------------------------------------- Details of GEBVtest calculations pro i_est = (20.335 - 20.926) / 22.004 = -0.027 pro p=0.990 x=-2.327 i=0.027 k=0.063 R2b=0.604 E(b1)=0.974 pro b1=1.000 se=0.066 E(b1)=0.974 t=0.40 R2_1=58.8 R2_2=25.4 pro passes t-test=Y bio-test=Y R2-test=Y ------------------------------------------------------------------------------- processing trt=scs h2=0.2500 evdate=20120801 depvar=GM byear=2004 type2x=N ------------------------------------------------------------------------------- stored 378 records for trait scs from file file300Cf_ABCHOL stored 290 records for trait scs from file file300Df_ABCHOL stored 378 records for trait scs from file file300Cr_ABCHOL stored 326 records for trait scs from file file300Gr_ABCHOL Summary statistics on candidate bulls (CB) and test bulls (TB) -------------------------------------------------------------- Trait Variable N Mean Std Min Max -------------------------------------------------------------- scs CB EBV 172 2.946 0.212 2.44 3.72 scs TB EBV 163 2.944 0.212 2.44 3.72 scs TB DPGM(y) 163 2.947 0.253 2.33 3.83 scs TB GEBV(x1) 163 2.919 0.167 2.38 3.28 scs TB EBVr(x2) 163 2.938 0.140 2.61 3.20 -------------------------------------------------------------- Details of GEBVtest calculations scs i_est = (2.944 - 2.946) / 0.212 = -0.005 scs p=0.998 x=-2.945 i=0.005 k=0.015 R2b=0.399 E(b1)=0.991 scs b1=0.950 se=0.093 E(b1)=0.991 t=-0.44 R2_1=39.6 R2_2=19.4 scs passes t-test=Y bio-test=Y R2-test=Y 2013-Jan-25 09:15:24 gebvtest.py: end }}}