Differences between revisions 1 and 4 (spanning 3 versions)
⇤ ← Revision 1 as of 2013-01-29 17:42:04 →
Size: 36
Editor: João
Comment:
← Revision 4 as of 2013-01-29 17:44:50 → ⇥
Size: 5921
Editor: João
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
Describe public/gebvtest_log here. gebvtest_log_ABCHOL
{{{
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
}}}

gebvtest_log_ABCHOL 

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 2024-08-13 10:39:05 by Valentina)