some success

notes

@@ -10,7 +10,67 @@ Ed
 
 
 
+4bin status
+─────────────
+    bin bounds:
+        0.008
+        0.02011893
+        0.05059644
+        0.12724332
+        0.32
+    
+    all converged but 17
+    1158A is suspect
+    4368A is suspect
+    4392A is suspect
 
 
-4368A is suspect
+7bin status
-4392A is suspect
+─────────────
+    bin bounds:
+        0.008
+        0.01444256
+        0.02607345
+        0.04707094
+        0.08497812
+        0.15341274
+        0.27695915
+        0.5
+
+    converged:
+        2600A
+        3471A -- suspect
+        4368A
+        4775A
+        6175A
+        6439A -- suspect
+        7647A
+        9157A
+
+    stalled:
+        1746A
+        1928A
+        5404A
+        5468A
+
+    broken:
+        1158A
+        1479A
+        2246A
+        3465A
+        4392A
+        8560A
+
+
+
+Hi Ed,
+
+
+
+I could not find any combination that allowed all models to converge. In all the combinations I've tried over the last 2 weeks, I never found 5 or 6 bins to be an advantage over 7 bins. I think it's reasonable for us to expect we can find 7-bin binnings that will work for some subsets with size about 1/3 of all lightcurves. We can tailor the bins to what we think are more important. That said, I was at least able to find one 4-bin binning that gave convergence for 17 of the 18 lightcurves. Some results are suspect because clag didn't report output for the fitting attempts for some parameters at those wavelengths. I could not get 1764A to converge under any binnings.
+
+We'll want to adjust the window spacing for the graphs, but for now I wanted to be sure the low values showed up. Again, these all have fully-developed errors for the lag measurements as reported by Abdu's program. The 7-bin overlayed result is just a best attempt to find a comprehensive setting. The boundaries and a suspect and convergence report are given below this message.
+
+I had to squeeze the frequency domain to get this to work well. I listed the boundaries below. Pushing up to 0.4 or down to 0.005 breaks many of these. I chose the outside boundaries based on two things: 0.008 is the nyquist frequency for something that occurs every ~60 days, right? So we wouldn't expect there to be much consistently observable lag on that time scale, since the smallest lightcurves aren't much longer than 100 days. At the high end, the lightcurves have a sampling rate of about 1 per day, which suggests the frequency scale should include up to 0.5. As I stated, that was a bit high, but bringing it down to 0.32 at that stage in my troubleshooting left me with the (considerable) result we see now.
+
+I will include a tar.gz archive with tables that are used to create this graph. The tables beginning with "cackett_" are tables in the format you requested corresponding to all of the results shown in the plots.

scripts/plot.sh

@@ -59,19 +59,24 @@ case $1 in
                                     sed "s|%FILEA%|$tabfile|"|
                                     sed "s|%LABEL%|$echo_band|")
             gnuplot_input="${gnuplot_input_edit}"
-        done
+            comp_tabfile=$(echo $tabfile|sed "s|${2}|${3}|")
-        for tabfile in data/tables/${3}/lag_*.tab;
-        do
-            ref_band_extracted=$(basename $tabfile|
-                                    sed 's|lag_\([0-9]\{4\}A\)_[0-9]\{4\}A.tab|\1|')
-            echo_band=$(basename $tabfile|
-                        sed 's|lag_[0-9]\{4\}A_\([0-9]\{4\}A\).tab|\1|')
-            if [[ "$echo_band" == "$ref_band" ]] ; then continue; fi
             gnuplot_input_edit=$(echo "$gnuplot_input"|
-                                    sed "s|%FILEB%|$tabfile|"|
+                                    sed "s|%FILEB%|$comp_tabfile|")
-                                    sed "s|%LABEL%|$echo_band|")
             gnuplot_input="${gnuplot_input_edit}"
+
         done
+#        for tabfile in data/tables/${3}/lag_*.tab;
+#        do
+#            ref_band_extracted=$(basename $tabfile|
+#                                    sed 's|lag_\([0-9]\{4\}A\)_[0-9]\{4\}A.tab|\1|')
+#            echo_band=$(basename $tabfile|
+#                        sed 's|lag_[0-9]\{4\}A_\([0-9]\{4\}A\).tab|\1|')
+#            if [[ "$echo_band" == "$ref_band" ]] ; then continue; fi
+#            gnuplot_input_edit=$(echo "$gnuplot_input"|
+#                                    sed "s|%FILEB%|$tabfile|"|
+#                                    sed "s|%LABEL%|$echo_band|")
+#            gnuplot_input="${gnuplot_input_edit}"
+#        done
         echo "$gnuplot_input"|perl -pe 's|␤|\n|g' > ${gnuplot_file}
         gnuplot $gnuplot_file
         rm $gnuplot_file

scripts/templates/lag_atlas_overlay.gp

@@ -5,7 +5,7 @@ set termopt enhanced
 set macros
 
     # Placement of the a,b,c,d labels in the graphs
-    POS = "at graph 0.61,0.85 font 'Times,12'"
+    POS = "at graph 0.72,0.88 font 'Times,9'"
     
     # x- and ytics for each row resp. column
     NOXNUMS = "unset xlabel;\
@@ -34,10 +34,11 @@ set macros
 unset key
 set logscale x
 set xtics auto font 'Times,9' offset 0,.5
-set ytics (-1,0,1,2,3,4,5) font 'Times,9'
+set ytics (-4,-2,0,2,4) font 'Times,9'
-set ytics add ('' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1)
+set ytics add ('' -5.5 1,'' -4.5 1,'' -3.5 1,'' -2.5 1,'' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1)
+set ytics add ('' -5 2,'' -3 2,'' -1 2,'' 1 2,'' 3 2,'' 5 2)
 set xrange [0.005:0.620];
-set yrange [-2:6]
+set yrange [-5:5]
 
 # Draw line at origin
 set arrow from 0.005,0 to 0.620,0 nohead lt 3 lc rgb 'black'