draft simulation

time_frames/frames_and_delay/frames.sif

+26/5/2018 13:54:12 10 2 1
+0 5
+5 10
+10 15
+15 20
+20 25
+25 30
+30 35
+35 40
+40 50
+50 60

time_frames/frames_and_delay/radiowater.fit

+FIT1        fit_winp 0.2.0  (c) 2010 by Turku PET Centre
+Date:       2010-09-19 19:19:34
+Data file:  mean.kbq
+Data unit:  kBq/mL
+Time unit:  sec
+Nr of VOIs: 1
+Region        Plane  Start   End     dataNr WSS       parNr Type Parameters
+000001 .      .        0.010 324.493    138 7.35E+003     6 0331 5.0 1.5E+001 5.888741E+002 7.143209E-002 8.410205E+001 1.925347E-003

time_frames/frames_and_delay/sim1.sh

+#!/bin/bash
+#: Title      : frames_and_delay/sim1
+#: Date       : 2023-06-20
+#: Author     : "Vesa Oikonen" <vesa.oikonen@utu.fi>
+#: Options    : None
+
+
+printf "\n=====================================================================\n"
+printf "  Making or emptying result folder\n"
+if [ ! -d sim1 ]; then
+  mkdir sim1
+else
+  rm -f sim1/*
+fi
+
+printf "=====================================================================\n"
+printf "Making 66 s radiowater BTAC with 0.1 s sample intervals \n"
+printf "=====================================================================\n"
+
+fit2dat -c=0,66,0.1 radiowater.fit sim1/blood.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+printf "=====================================================================\n"
+printf "IDIF from LV cavity or aorta is early as compared to tissue TACs,\n"
+printf "requiring delay correction.\n"
+printf "Move BTAC forward in time 1 s at a time \n"
+printf "to simulate delay-corrected BTACs without any errors.\n"
+printf "=====================================================================\n"
+
+tactime -nogap sim1/blood.tac +1 sim1/blood01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime -nogap sim1/blood.tac +2 sim1/blood02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime -nogap sim1/blood.tac +3 sim1/blood03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime -nogap sim1/blood.tac +4 sim1/blood04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime -nogap sim1/blood.tac +5 sim1/blood05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot the original BTACs \n"
+tac2svg -legend=n -mt="" sim1/sim_frames_and_peak_blood.svg -l sim1/blood.tac sim1/blood??.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "=====================================================================\n"
+printf "Apply PET frames to the BTACs \n"
+printf "to simulate delay-corrected IDIF BTACs without any errors.\n"
+printf "=====================================================================\n"
+
+simframe sim1/blood.tac frames.sif sim1/true_idif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe sim1/blood01.tac frames.sif sim1/true_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe sim1/blood02.tac frames.sif sim1/true_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe sim1/blood03.tac frames.sif sim1/true_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe sim1/blood04.tac frames.sif sim1/true_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe sim1/blood05.tac frames.sif sim1/true_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot the true IDIF BTACs with bars and with AIF BTACs\n"
+tac4frpl sim1/true_idif00.tac sim1/bar_true_idif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/true_idif01.tac sim1/bar_true_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/true_idif02.tac sim1/bar_true_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/true_idif03.tac sim1/bar_true_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/true_idif04.tac sim1/bar_true_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/true_idif05.tac sim1/bar_true_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif00.svg -l sim1/blood.tac sim1/bar_true_idif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif01.svg -l sim1/blood01.tac sim1/bar_true_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif02.svg -l sim1/blood02.tac sim1/bar_true_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif03.svg -l sim1/blood03.tac sim1/bar_true_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif04.svg -l sim1/blood04.tac sim1/bar_true_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif05.svg -l sim1/blood05.tac sim1/bar_true_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "=====================================================================\n"
+printf "Delay-corrected IDIFs need to have the original PET frame times to be useful in pixel-by-pixel analyses.\n"
+printf "Delay correct IDIFs using tacmove which takes into account the frame start and end times \n"
+printf "=====================================================================\n"
+
+cp -f sim1/true_idif00.tac sim1/tacmove_idif01.tac
+cp -f sim1/true_idif00.tac sim1/tacmove_idif02.tac
+cp -f sim1/true_idif00.tac sim1/tacmove_idif03.tac
+cp -f sim1/true_idif00.tac sim1/tacmove_idif04.tac
+cp -f sim1/true_idif00.tac sim1/tacmove_idif05.tac
+tacmove sim1/tacmove_idif01.tac 1
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/tacmove_idif02.tac 2
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/tacmove_idif03.tac 3
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/tacmove_idif04.tac 4
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/tacmove_idif05.tac 5
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot these with bars and with AIF BTACs\n"
+tac4frpl sim1/tacmove_idif01.tac sim1/bar_tacmove_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/tacmove_idif02.tac sim1/bar_tacmove_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/tacmove_idif03.tac sim1/bar_tacmove_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/tacmove_idif04.tac sim1/bar_tacmove_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac4frpl sim1/tacmove_idif05.tac sim1/bar_tacmove_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif01.svg -l sim1/blood01.tac sim1/bar_tacmove_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif02.svg -l sim1/blood02.tac sim1/bar_tacmove_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif03.svg -l sim1/blood03.tac sim1/bar_tacmove_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif04.svg -l sim1/blood04.tac sim1/bar_tacmove_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif05.svg -l sim1/blood05.tac sim1/bar_tacmove_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+printf "=====================================================================\n"
+printf "The simple and traditional method is to assume that frame average concentration represents \n"
+printf "the concentration at frame mid time and then do dot-to-dot interpolations and integrations. \n"
+printf "=====================================================================\n"
+
+printf "\nReplace IDIF frame start and end times with frame mid times\n"
+
+tacformat -mid -f=PMOD sim1/true_idif00.tac sim1/midtime.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nCorrect that for delay \n"
+
+cp -f sim1/midtime.tac sim1/trad_idif01.tac
+cp -f sim1/true_idif00.tac sim1/trad_idif02.tac
+cp -f sim1/true_idif00.tac sim1/trad_idif03.tac
+cp -f sim1/true_idif00.tac sim1/trad_idif04.tac
+cp -f sim1/true_idif00.tac sim1/trad_idif05.tac
+tacmove sim1/trad_idif01.tac 1
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/trad_idif02.tac 2
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/trad_idif03.tac 3
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/trad_idif04.tac 4
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tacmove sim1/trad_idif05.tac 5
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nInterpolate into original frame mid times \n"
+
+interpol -f=sim1/midtime.tac sim1/trad_idif01.tac sim1/trad_mid_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=sim1/midtime.tac sim1/trad_idif02.tac sim1/trad_mid_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=sim1/midtime.tac sim1/trad_idif03.tac sim1/trad_mid_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=sim1/midtime.tac sim1/trad_idif04.tac sim1/trad_mid_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=sim1/midtime.tac sim1/trad_idif05.tac sim1/trad_mid_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot these with AIF BTACs\n"
+
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif01.svg -l sim1/blood01.tac sim1/trad_mid_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif02.svg -l sim1/blood02.tac sim1/trad_mid_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif03.svg -l sim1/blood03.tac sim1/trad_mid_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif04.svg -l sim1/blood04.tac sim1/trad_mid_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif05.svg -l sim1/blood05.tac sim1/trad_mid_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nInterpolate into original frame durations (same result as above therefore not plotted) \n"
+
+simframe --force sim1/trad_idif01.tac frames.sif sim1/trad_frame_idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe --force sim1/trad_idif02.tac frames.sif sim1/trad_frame_idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe --force sim1/trad_idif03.tac frames.sif sim1/trad_frame_idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe --force sim1/trad_idif04.tac frames.sif sim1/trad_frame_idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe --force sim1/trad_idif05.tac frames.sif sim1/trad_frame_idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+
+printf "=====================================================================\n"
+printf "Calculate AUC from zero to frame end times. \n"
+printf "Plot its error at each frame end time.\n"
+printf "=====================================================================\n"
+
+REM interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/blood.tac sim1/auc.tac
+
+interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/blood01.tac sim1/auc01.tac
+
+tacfr2x sim1/true_idif01.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc01.tac sim1/aucerror01.tac
+
+tacfr2x sim1/tacmove_idif01.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
+tacadd sim1/aucerror01.tac sim1/temp2.tac
+
+interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/trad_mid_idif01.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
+tacadd sim1/aucerror01.tac sim1/temp2.tac
+
+REM tacfr2x sim1/trad_frame_idif01.tac sim1/steps.tac
+REM interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
+REM taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
+REM tacadd sim1/aucerror01.tac sim1/temp2.tac
+
+tacren sim1/aucerror01.tac 1 True IDIF
+tacren sim1/aucerror01.tac 2 tacmove
+tacren sim1/aucerror01.tac 3 Trad
+REM tacren sim1/aucerror01.tac 4 Trad frame
+
+tac2svg -x2=65 -legend=y -mt="" -yt="s*kBq/mL" sim1/sim_frames_and_peak_aucerror01.svg -l sim1/aucerror01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+
+
+exit 0
+
+printf "=====================================================================\n"
+printf "Assuming that BTAC is derived from PET image with frame times given in SIF \n"
+printf "=====================================================================\n"
+
+simframe blood.tac frames.sif idif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe blood01.tac frames.sif idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe blood02.tac frames.sif idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe blood03.tac frames.sif idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe blood04.tac frames.sif idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+simframe blood05.tac frames.sif idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\n Plot original BTAC and its framed version with bars \n"
+tac4frpl idif00.tac blood2frpl.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -y2=460 -legend=n -mt="" sim_frames_and_delay_btac.svg -l blood2frpl.tac blood.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "=====================================================================\n"
+printf "Move framed IDIF back in time \n"
+printf "=====================================================================\n"
+
+tactime idif01.tac -1 idif01dc.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime idif02.tac -2 idif02dc.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime idif03.tac -3 idif03dc.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime idif04.tac -4 idif04dc.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tactime idif05.tac -5 idif05dc.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "=====================================================================\n"
+printf "And interpolate to original framing \n"
+printf "=====================================================================\n"
+
+tacmove 
+
+
+
+exit 0
+
+
+printf "\nPlot delay time matched original BTAC and simulated IDIF together \n"
+tac2svg -legend=n -mt="" idif00.svg -l blood.tac -b idif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif01.svg -l blood01.tac -b idif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif02.svg -l blood02.tac -b idif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif03.svg -l blood03.tac -b idif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif04.svg -l blood04.tac -b idif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif05.svg -l blood05.tac -b idif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif06.svg -l blood06.tac -b idif06.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif07.svg -l blood07.tac -b idif07.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif08.svg -l blood08.tac -b idif08.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif09.svg -l blood09.tac -b idif09.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+tac2svg -legend=n -mt="" idif10.svg -l blood10.tac -b idif10.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot all simulated IDIFs together \n"
+tac2svg -y2=460 -legend=n -mt="" sim_frames_and_peak_idif.svg idif??.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+printf "=====================================================================\n"
+printf "Compute bias in IDIFs caused by framing but time delay is accounted for \n"
+printf "=====================================================================\n"
+
+tacformat -mid -f=PMOD idif00.tac midtime.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+interpol -f=midtime.tac blood.tac bl00.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood01.tac bl01.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood02.tac bl02.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood03.tac bl03.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood04.tac bl04.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood05.tac bl05.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood06.tac bl06.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood07.tac bl07.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood08.tac bl08.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood09.tac bl09.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+interpol -f=midtime.tac blood10.tac bl10.tac 
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+taccalc idif00.tac - bl00.tac dif00.tac
+taccalc idif01.tac - bl01.tac dif01.tac
+taccalc idif02.tac - bl02.tac dif02.tac
+taccalc idif03.tac - bl03.tac dif03.tac
+taccalc idif04.tac - bl04.tac dif04.tac
+taccalc idif05.tac - bl05.tac dif05.tac
+taccalc idif06.tac - bl06.tac dif06.tac
+taccalc idif07.tac - bl07.tac dif07.tac
+taccalc idif08.tac - bl08.tac dif08.tac
+taccalc idif09.tac - bl09.tac dif09.tac
+taccalc idif10.tac - bl10.tac dif10.tac
+
+taccalc dif00.tac div bl00.tac rdif00.tac
+taccalc dif01.tac div bl01.tac rdif01.tac
+taccalc dif02.tac div bl02.tac rdif02.tac
+taccalc dif03.tac div bl03.tac rdif03.tac
+taccalc dif04.tac div bl04.tac rdif04.tac
+taccalc dif05.tac div bl05.tac rdif05.tac
+taccalc dif06.tac div bl06.tac rdif06.tac
+taccalc dif07.tac div bl07.tac rdif07.tac
+taccalc dif08.tac div bl08.tac rdif08.tac
+taccalc dif09.tac div bl09.tac rdif09.tac
+taccalc dif10.tac div bl10.tac rdif10.tac
+
+taccalc rdif00.tac x 100 pdif00.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif01.tac x 100 pdif01.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif02.tac x 100 pdif02.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif03.tac x 100 pdif03.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif04.tac x 100 pdif04.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif05.tac x 100 pdif05.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif06.tac x 100 pdif06.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif07.tac x 100 pdif07.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif08.tac x 100 pdif08.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif09.tac x 100 pdif09.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+taccalc rdif10.tac x 100 pdif10.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+printf "\nPlot percent bias curves \n"
+tac2svg -legend=n -yt="% bias" -mt="" sim_frames_and_peak_idif_bias.svg -l pdif??.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+
+
+printf "\n===================================================================\n"
+printf "   All passed! \n"
+printf "===================================================================\n"
+exit 0
+