1#!/usr/bin/env bash 2# Testing gmtregress on the data in Draper & Smith [1998] 3# Applied Regression Analysis, 3rd Ed, Wiley. 4 5ps=draper.ps 6cat << EOF > draper.txt 7# Table 1.1 in Draper & Smith, Applied Regression Analysis 835.3 10.98 929.7 11.13 1030.8 12.51 1158.8 8.40 1261.4 9.27 1371.3 8.73 1474.4 6.36 1576.7 8.50 1670.7 7.82 1757.5 9.14 1846.4 8.24 1928.9 12.19 2028.1 11.88 2139.1 9.57 2246.8 10.94 2348.5 9.58 2459.3 10.09 2570.0 8.11 2670.0 6.83 2774.5 8.88 2872.1 7.68 2958.1 8.47 3044.6 8.86 3133.4 10.36 3228.6 11.08 33EOF 34# First plot data and basic LS fit with equation 35txt=$(gmt regress -Ey -N2 -Fxm draper.txt -T0 | awk '{printf "85 13 @!y\\303 = %.4f %.4f x\n", $17, $15}') 36gmt psbasemap -R20/85/6/13 -JX6.5i/4i -P -Xc -K -Baf -BWSne > $ps 37gmt regress -Ey -N2 -Fxym draper.txt | awk '{printf "> error\n%s %s\n%s %s\n", $1, $2, $1, $3}' | gmt psxy -R -J -O -K -W0.25p,red,- >> $ps 38gmt psxy -R -J -O -K draper.txt -Sc0.2c -Gblue >> $ps 39gmt regress -Ey -N2 -Fxm -T25/80/2+n draper.txt | gmt psxy -R -J -O -K -W2p >> $ps 40echo "$txt" | gmt pstext -R -J -O -K -F+jRT+f18p -Dj0.1i >> $ps 41# Redo plot and basic LS fit but also show 68%, 95% & 99% confidence band 42gmt psbasemap -R -J -O -K -Baf -BWSNe+t"Draper & Smith [1998] Regression" -Y4.75i >> $ps 43gmt regress -Ey -N2 -Fxmc -T25/80/1 -C99 draper.txt | gmt psxy -R -J -O -K -L+d -Glightgreen >> $ps 44gmt regress -Ey -N2 -Fxmc -T25/80/1 -C95 draper.txt | gmt psxy -R -J -O -K -L+d -Glightorange >> $ps 45gmt regress -Ey -N2 -Fxmc -T25/80/1 -C68 draper.txt | gmt psxy -R -J -O -K -L+d -Glightred -W2p >> $ps 46gmt psxy -R -J -O -K draper.txt -Sc0.2c -Gblue >> $ps 47gmt pslegend -DjTR+w1.65i+jRT+o0.1i/0.1i -R -J -O -F+p1p << EOF >> $ps 48S 0.1i s 0.15i lightgreen - 0.25i 99% Confidence 49S 0.1i s 0.15i lightorange - 0.25i 95% Confidence 50S 0.1i s 0.15i lightred - 0.25i 68% Confidence 51EOF 52