2025-May-15 16:01:49 gebvtest.py: start

Version: 2025.05.15
Options:  -v -m hol abc ../sample_data --outdir ../ITBC_results

absolute DATADIR: /interbull/GEBV_test/web/gebvtest20250515/sample_data
Processing brd=HOL pop=ABC cou=ABC datadir=../sample_data
stored   1014 records from file file736_ABCHOL
opening output file /interbull/GEBV_test/web/gebvtest20250515/ITBC_results/file735_ABCHOL ...
opening output file /interbull/GEBV_test/web/gebvtest20250515/ITBC_results/file735adj_ABCHOL ...
reading file300Cr_ABCHOL ...
stored    991 records for trait mil from file file300Cr_ABCHOL
reading file300Gr_ABCHOL ...
stored    785 records for trait mil from file file300Gr_ABCHOL
reading file300Cf_ABCHOL ...
stored    967 records for trait mil from file file300Cf_ABCHOL
reading file300Gf_ABCHOL ...
stored    742 records for trait mil from file file300Gf_ABCHOL

-------------------------------------------------------------------------------
processing trt=mil h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N
-------------------------------------------------------------------------------
derived    401 records for trait mil for a DGEBV file300 equivalent

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         171   577.825   811.031 -1908.00  2962.00
mil   TB  VAL(y)     171  1046.220   904.221 -1436.40  4147.10 Base=Current
mil   TB GEBV(x1)    171  1068.807   599.782  -726.00  2624.00 Base=Old
mil   TB EBVr(x2)    171  1261.292   545.417  -564.00  2616.00 Base=Old
mil   TB GEBV(x1)    171   729.321   624.825 -1140.42  2349.45 Base=Current
mil   TB EBVr(x2)    171   929.844   568.190  -971.66  2341.11 Base=Current
mil   BB  EBVf(y)    570   -97.361   780.743 -2580.00  1884.00
mil   BB  EBVr(x)    570   310.270   745.217 -1984.00  2318.00
mil   BB  EBV(wt)    570    48.554    42.090     0.00    97.00
--------------------------------------------------------------
mil : DGEBV : Base 0.25: Full = -384.11 + 1.042 * Red  b1 0.903 ~> 0.867  b2 1.058 ~> 1.016

Details of GEBVtest calculations
mil i_est = (577.825 - 573.367) / 813.411 = 0.005

After Base adjustments:
mil p=0.998 x=-2.929 i=0.005 k=0.016 R2b=0.369 E(b1)=0.990 b1=0.867
mil b1-E(b1) = -0.123  t = -1.29  bootstrap P(b=Eb) = 19.8  b1_test = Pass
mil R2_1 = 36.5  R2_2 = 41.2   bootstrap P(R2_gain) = 19.4  R2_test = NS-
mil   stat_test=Y   pract_test=Y   bio_test=N   R2_test=-   overall = PASS

2025-May-15 16:01:49 gebvtest.py: finished mil

reading file300Cr_ABCHOL ...
stored    991 records for trait fat from file file300Cr_ABCHOL
reading file300Gr_ABCHOL ...
stored    785 records for trait fat from file file300Gr_ABCHOL
reading file300Cf_ABCHOL ...
stored    967 records for trait fat from file file300Cf_ABCHOL
reading file300Gf_ABCHOL ...
stored    742 records for trait fat from file file300Gf_ABCHOL

-------------------------------------------------------------------------------
processing trt=fat h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N
-------------------------------------------------------------------------------
derived    401 records for trait fat for a DGEBV file300 equivalent

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         171    24.409    29.816   -65.00    95.00
fat   TB  VAL(y)     171    36.121    32.424   -35.31   156.08 Base=Current
fat   TB GEBV(x1)    171    38.661    20.519   -19.00    86.00 Base=Old
fat   TB EBVr(x2)    171    45.135    20.639    -4.00    92.00 Base=Old
fat   TB GEBV(x1)    171    30.057    21.734   -31.02    80.20 Base=Current
fat   TB EBVr(x2)    171    36.914    21.860   -15.13    86.55 Base=Current
fat   BB  EBVf(y)    570   -12.807    32.849  -109.00    65.00
fat   BB  EBVr(x)    570    -0.105    30.713   -83.00    81.00
fat   BB  EBV(wt)    570    48.554    42.090     0.00    97.00
--------------------------------------------------------------
fat : DGEBV : Base 0.25: Full = -10.89 + 1.059 * Red  b1 0.809 ~> 0.764  b2 0.771 ~> 0.728

Details of GEBVtest calculations
fat i_est = (24.409 - 24.284) / 29.839 = 0.004

After Base adjustments:
fat p=0.999 x=-3.018 i=0.004 k=0.013 R2b=0.270 E(b1)=0.991 b1=0.764
fat b1-E(b1) = -0.227  t = -2.30  bootstrap P(b=Eb) =  1.6  b1_test = Fail
fat R2_1 = 26.7  R2_2 = 24.6   bootstrap P(R2_gain) = 64.8  R2_test = NS+
fat   stat_test=N   pract_test=N   bio_test=N   R2_test=-   overall = FAIL

2025-May-15 16:01:50 gebvtest.py: finished fat

reading file300Cr_ABCHOL ...
stored    991 records for trait pro from file file300Cr_ABCHOL
reading file300Gr_ABCHOL ...
stored    785 records for trait pro from file file300Gr_ABCHOL
reading file300Cf_ABCHOL ...
stored    967 records for trait pro from file file300Cf_ABCHOL
reading file300Gf_ABCHOL ...
stored    742 records for trait pro from file file300Gf_ABCHOL

-------------------------------------------------------------------------------
processing trt=pro h2=0.4000 evdate=20120801 depvar=GM byear=2004 type2x=N
-------------------------------------------------------------------------------
derived    457 records for trait pro for a DGEBV file300 equivalent

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         171    21.567    21.541   -50.00    89.00
pro   TB  VAL(y)     171    36.780    23.995   -29.35   102.68 Base=Current
pro   TB GEBV(x1)    171    35.895    15.844   -12.00    72.00 Base=Old
pro   TB EBVr(x2)    171    41.368    15.377    -8.00    76.00 Base=Old
pro   TB GEBV(x1)    171    25.688    16.486   -24.15    63.26 Base=Current
pro   TB EBVr(x2)    171    31.383    15.999   -19.98    67.42 Base=Current
pro   BB  EBVf(y)    620    -4.660    23.030   -72.00    61.00
pro   BB  EBVr(x)    620     7.800    22.285   -56.00    64.00
pro   BB  EBV(wt)    620    35.936    38.719     0.00    97.00
--------------------------------------------------------------
pro : DGEBV : Base 0.25: Full = -11.66 + 1.040 * Red  b1 0.871 ~> 0.837  b2 0.825 ~> 0.793

Details of GEBVtest calculations
pro i_est = (21.567 - 20.926) / 22.004 = 0.029

After Base adjustments:
pro p=0.989 x=-2.293 i=0.029 k=0.068 R2b=0.342 E(b1)=0.954 b1=0.837
pro b1-E(b1) = -0.118  t = -1.25  bootstrap P(b=Eb) = 19.4  b1_test = Pass
pro R2_1 = 32.6  R2_2 = 27.6   bootstrap P(R2_gain) = 83.3  R2_test = NS+
pro   stat_test=Y   pract_test=Y   bio_test=N   R2_test=-   overall = PASS

2025-May-15 16:01:50 gebvtest.py: finished pro

reading file300Cr_ABCHOL ...
stored    966 records for trait scs from file file300Cr_ABCHOL
reading file300Gr_ABCHOL ...
stored    762 records for trait scs from file file300Gr_ABCHOL
reading file300Cf_ABCHOL ...
stored    966 records for trait scs from file file300Cf_ABCHOL
reading file300Gf_ABCHOL ...
stored    762 records for trait scs from file file300Gf_ABCHOL

-------------------------------------------------------------------------------
processing trt=scs h2=0.2500 evdate=20120801 depvar=GM byear=2004 type2x=N
-------------------------------------------------------------------------------
derived    530 records for trait scs for a DGEBV file300 equivalent

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         165     2.944     0.210     2.44     3.72
scs   TB  VAL(y)     165     2.739     1.983    -3.14     9.85 Base=Current
scs   TB GEBV(x1)    165     2.917     0.170     2.38     3.28 Base=Old
scs   TB EBVr(x2)    165     2.938     0.140     2.61     3.20 Base=Old
scs   TB GEBV(x1)    165     2.930     0.167     2.40     3.29 Base=Current
scs   TB EBVr(x2)    165     2.950     0.138     2.63     3.21 Base=Current
scs   BB  EBVf(y)    710     3.157     0.250     2.48     3.85
scs   BB  EBVr(x)    710     3.134     0.224     2.50     3.73
scs   BB  EBV(wt)    710     9.288    21.446     0.00    95.00
--------------------------------------------------------------
scs : DGEBV : Base 0.25: Full = 0.06 + 0.984 * Red  b1 1.147 ~> 1.166  b2 0.401 ~> 0.408

Details of GEBVtest calculations
scs i_est = (2.944 - 2.946) / 0.212 = -0.005

After Base adjustments:
scs p=0.998 x=-2.947 i=0.005 k=0.015 R2b=0.011 E(b1)=0.985 b1=1.166
scs b1-E(b1) =  0.181  t =  0.24  bootstrap P(b=Eb) = 79.8  b1_test = Pass
scs R2_1 =  1.1  R2_2 =  0.1   bootstrap P(R2_gain) = 83.5  R2_test = NS+
scs   stat_test=Y   pract_test=Y   bio_test=Y   R2_test=-   overall = PASS

2025-May-15 16:01:51 gebvtest.py: finished scs

2025-May-15 16:01:51 gebvtest.py: end
