##################################################################
#   Reference Design:  F12, F14, F16, F18 ===> HIGH vs MEDIUM    #
#                      F13, F15, F19, F20 ===>  LOW vs MEDIUM    #
#                                                                #
#   All Pairs Design:  F46-F48 HvL                               #
#                      F49-F50 MvL                               #
#                      F71-F72 HvM                               #
##################################################################

loc=/mnt/windows/UNE_data/Diets/GPR_Files

for N in F12 F14 F16 F18 F13 F15 F19 F20 F46 F48 F49 F50 F71 F72
do
   sed "s/\"//g" $loc/$N.gpr  |                    \
   awk 'NR>30 && (($9-$12)/$14)>1 {print $1"<<>>"$5, $9, $12, $9, $10}' | \
   sed 's/\./ /g' | awk '{print $1, $(NF-3), $(NF-2), $(NF-1), $NF}' | \
   awk '{print $1, log($2-$3)/log(2), log($4)/log(2), log($5)/log(2)}' | \
   awk '$3>$4 {$5=$4/$3}; $3<=$4 {$5=$3/$4}; {print $0}' | \
   awk '$5>0.85 {print $1, $2}' | sed "s/<<>>/ /g" | \
   awk 'length($1)==1 {$1="0"$1}; {print $2, $1, $3}' | sort > $N.red

   sed "s/\"//g" $loc/$N.gpr  |                    \
   awk 'NR>30 && (($18-$21)/$23)>1 {print $1"<<>>"$5, $18, $21, $18, $19}' | \
   sed 's/\./ /g' | awk '{print $1, $(NF-3), $(NF-2), $(NF-1), $NF}' | \
   awk '{print $1, log($2-$3)/log(2), log($4)/log(2), log($5)/log(2)}' | \
   awk '$3>$4 {$5=$4/$3}; $3<=$4 {$5=$3/$4}; {print $0}' | \
   awk '$5>0.85 {print $1, $2}' | sed "s/<<>>/ /g" | \
   awk 'length($1)==1 {$1="0"$1}; {print $2, $1, $3}' | sort > $N.green
done

########################
#  Make Diets_Tot.dat  #
########################

awk '{print $1, "REFE_01_" $2 "_Grn_H", $3}' F12.green | sort >  Diets_Tot.dat
awk '{print $1, "REFE_02_" $2 "_Grn_H", $3}' F14.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_03_" $2 "_Grn_H", $3}' F16.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_04_" $2 "_Grn_H", $3}' F18.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_01_" $2 "_Red_M", $3}' F12.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_02_" $2 "_Red_M", $3}' F14.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_03_" $2 "_Red_M", $3}' F16.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_04_" $2 "_Red_M", $3}' F18.red | sort >> Diets_Tot.dat

awk '{print $1, "REFE_05_" $2 "_Grn_H", $3}' F13.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_07_" $2 "_Grn_H", $3}' F15.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_06_" $2 "_Grn_H", $3}' F19.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_08_" $2 "_Grn_H", $3}' F20.green | sort >> Diets_Tot.dat
awk '{print $1, "REFE_05_" $2 "_Red_L", $3}' F13.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_07_" $2 "_Red_L", $3}' F15.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_06_" $2 "_Red_L", $3}' F19.red | sort >> Diets_Tot.dat
awk '{print $1, "REFE_08_" $2 "_Red_L", $3}' F20.red | sort >> Diets_Tot.dat

awk '{print $1, "ALLP_09_" $2 "_Red_M", $3}' F71.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_09_" $2 "_Grn_H", $3}' F71.green | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_10_" $2 "_Red_H", $3}' F72.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_10_" $2 "_Grn_M", $3}' F72.green | sort >> Diets_Tot.dat

awk '{print $1, "ALLP_11_" $2 "_Red_M", $3}' F50.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_11_" $2 "_Grn_L", $3}' F50.green | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_12_" $2 "_Red_L", $3}' F49.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_12_" $2 "_Grn_M", $3}' F49.green | sort >> Diets_Tot.dat

awk '{print $1, "ALLP_13_" $2 "_Red_L", $3}' F46.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_13_" $2 "_Grn_H", $3}' F46.green | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_14_" $2 "_Red_H", $3}' F48.red | sort >> Diets_Tot.dat
awk '{print $1, "ALLP_14_" $2 "_Grn_L", $3}' F48.green | sort >> Diets_Tot.dat

awk '{print $1, substr($2,16,1), substr($2,12,3)}' Diets_Tot.dat | sort | \
awk 'BEGIN{ gene=""; diet=""; dye=""
            nobs=1; ndiets=1; ndyes=1 }
     { if($1==gene){ nobs++
                     if($2!=diet){ ndiets++; diet=$2 }
                     if($3!=dye){ ndyes++; dye=$3 }
                   }
      else{ if(gene!="")print gene, nobs, ndiets, ndyes/ndiets
            gene=$1; diet=$2; dye=$3
            nobs=1; ndiets=1; ndyes=1
          }
     }END{ print gene, nobs, diets, ndyes/ndiets}' > Diets_Tot.genes

awk '$2>=28 && $3==3 {print $0}' Diets_Tot.genes > Diets.genes
sort Diets_Tot.dat | join - Diets.genes | \
awk '{print $1, $2, $3}' > Diets.dat

rm F??.red F??.green