When I use the tim.sh to generate my own TIMs, I found the tim_summed.txt output file was not for regions but for CpG sites. After debugging, I found there is a code error in this line "awk -v window="$window_size" '{print $1 "\t" $2-$window "\t" $3+$window}' $output_file"_sites" > $outputfile""$window_size". It should be no $ before window inside awk.
Hi Christa,
Thank you for developing CelFiE.
When I use the tim.sh to generate my own TIMs, I found the tim_summed.txt output file was not for regions but for CpG sites. After debugging, I found there is a code error in this line "awk -v window="$window_size" '{print $1 "\t" $2-$window "\t" $3+$window}' $output_file"_sites" > $outputfile""$window_size". It should be no $ before window inside awk.
Hope this is helpful.