1#!/usr/bin/env bash 2# Test Reweighted L-2 measure using Reduced Major Axis area misfit 3# Make sure the Reweighted path (-Nw) yields the same regression as 4# regression just the good points after removing the outliers. 5# Test data from Dietmar Muller. 6ps=rls_on_rma.ps 7 8echo "0 red 1 red" > code.cpt 9echo "1 blue 2 blue" >> code.cpt 10gmt regress noisy_data.txt -Nr -Er > rma.txt 11gmt regress noisy_data.txt -Nr -Er -S > raw.txt 12gmt regress noisy_data.txt -Nw -Er > model.txt 13gmt psxy rma.txt -R-200/-30/200/800 -JX6i/9i -P -Baf -W2p,green -i0,2 -K -Xc > $ps 14gmt psxy model.txt -R -J -O -K -W4p -i0,2 >> $ps 15gmt regress raw.txt -i0,1 -Er | gmt psxy -R -J -O -K -W1p,orange -i0,2 >> $ps 16gmt psxy model.txt -i0,1,6 -Ccode.cpt -R -J -O -Sc0.25c -K >> $ps 17gmt pslegend -R -J -O -DjBR+w3.5i+jBR+o0.1i -F+p1p << EOF >> $ps 18S 0.2i c 0.25c blue - 0.5i Good point according to @;green;RMA via LMS@;; 19S 0.2i c 0.25c red - 0.5i Outlier according to @;green;RMA via LMS@;; 20D 0.1i 0.25p 21S 0.2i - 0.25i - 2p,green 0.5i RMA via LMS 22S 0.2i - 0.25i - 4p 0.5i Reweighted RMA 23S 0.2i - 0.25i - 1p,orange 0.5i RMA via LS on good points only 24EOF 25