#################################################################
# Prepare Data file for VCE Input
#
#	Input for VCE			Location in Diets.dat
#
#  	1. Y (Intensity)		$3
#	2. Fixed (Array|Block|Dye)	$2
#	3. Gene				$1
#	4. Gene*Array|Block		$1*substr($2,1,10)
#	5. Gene*Dye        		$1*substr($2,12,3)
#	6. Gene*Diet       		$1*substr($2,16,1)
#
#  NB: All the effects need to be codified and formated
#################################################################

echo Intens Fixed Gene GxAB GxDye GxDiet > vce.dat

awk '{print $3, $2, $1, \
            $1"__"substr($2,1,10), \
            $1"__"substr($2,12,3), \
            $1"__"substr($2,16,1)}' Diets.dat | sort +1 | \
awk 'BEGIN{ fixed=""; n=0 }
            {if($2==fixed){ print $1, n, $3, $4, $5, $6 }
             else{ n++
                   print $1, n, $3, $4, $5, $6
                   fixed=$2 } }' | sort +2 | \
awk 'BEGIN{ gene=""; n=0 }
            {if($3==gene){ print $1, $2, n, $4, $5, $6 }
             else{ n++
                   print $1, $2, n, $4, $5, $6
                   gene=$3 } }' | sort +3 | \
awk 'BEGIN{ gxab=""; n=0 }
            {if($4==gxab){ print $1, $2, $3, n, $5, $6 }
             else{ n++
                   print $1, $2, $3, n, $5, $6
                   gxab=$4 } }' | sort +4 | \
awk 'BEGIN{ gxdye=""; n=0 }
            {if($5==gxdye){ print $1, $2, $3, $4, n, $6 }
             else{ n++
                   print $1, $2, $3, $4, n, $6
                   gxdye=$5 } }' | sort +5 | \
awk 'BEGIN{ gxdiet=""; n=0 }
            {if($6==gxdiet){ print $1, $2, $3, $4, $5, n }
             else{ n++
                   print $1, $2, $3, $4, $5, n
                   gxdiet=$6 } }' | \
awk '{printf"%12.5E%8d%8d%8d%8d%8d\n", $1, $2, $3, $4, $5, $6}' >> vce.dat