#
# Computation of deregressed proof for sires.
#
# usage: awk -v h2=0.25 -f pvr_drp.awk sol_and_acc > drp.txt
#
# You can change the relitability as -v h2=value.
# The default heritability is 0.25.
#
BEGIN{
    # default h2=0.25 equiv. kd=14
    if(h2<=0){ h2=0.25 }
    kd=(4-2*h2)/h2
    print "h2=",h2,"; kd=",kd > "/dev/stderr"
}
NR>1{
    DE_EBV = kd*$5/(1 - $5)
    DE_PA  = kd*$10/(1 - $10)
    DE_R   = DE_EBV - DE_PA
    R = DE_R/(DE_R + kd)
    if(R>0){
        DRP = $6 + ($4-$6)/R
    } else {
        R = 0.0
        DRP = 0.0
    }
    print $0,DRP,R
}