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
public/gebvtest_log (last edited 2013-01-29 17:46:53 by João)
