maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
14

	
15

	
16

	
17

	
18

	
19

	
20

	
21

	
22

	
23

	
24

	
25

	
26

	
27

	
28

	
29

	
            #!/usr/bin/env bash
ps=distributions.ps
# LL: Histogram in counts with 3 normal distribution curves based on L2, L1, and LMS statistsics
gmt pshistogram -Baf -BWSne+glightblue @v3206_06.txt -R-6000/0/0/1300 -JX3i/2.0i -Gorange -L1p -Z0 -W200 -P -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -Y0.75i  > $ps
echo "0 1300 Counts" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# LR: Same, but in percent
gmt pshistogram -Bxaf -Byaf+u" %" -BWSne+glightblue @v3206_06.txt -R-6000/0/0/20 -JX3i/2.0i -Gorange -L1p -Z1 -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X3.75i >> $ps
echo "0 20 Percent" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# ML: Same, but in log 1+counts
gmt pshistogram -Baf -BWSne+glightblue @v3206_06.txt -R-6000/0/0/8 -JX3i/2.0i -Gorange -L1p -Z2 -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X-3.75i -Y2.5i >> $ps
echo "0 8 log1p(counts)" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# MR: Same, but in log 1+percent
gmt pshistogram -Baf -BWSne+glightblue @v3206_06.txt -R-6000/0/0/3.5 -JX3i/2.0i -Gorange -L1p -Z3 -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X3.75i >> $ps
echo "0 3.5 log1p(percent)" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# TL: Same, but in log10 1+counts
gmt pshistogram -Baf -BWSne+glightblue @v3206_06.txt -R-6000/0/0/3.5 -JX3i/2.0i -Gorange -L1p -Z4 -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X-3.75i -Y2.5i >> $ps
echo "0 3.5 log101p(counts)" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# TR: Same, but in log10 1+percent
gmt pshistogram -Baf -BWSne+glightblue @v3206_06.txt -R-6000/0/0/1.5 -JX3i/2.0i -Gorange -L1p -Z5 -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X3.75i >> $ps
echo "0 1.5 log101p(percent)" | gmt pstext -R -J -O -K -Dj0.1i -F+jTR+f12p >> $ps
# L: Cumulative histogram
gmt pshistogram -Bxaf -Byafg100+u" %" -BWSne+glightblue @v3206_06.txt -R-6000/0/0/105 -JX3i/2.0i -Glightred -L1p -Z1 -Q -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X-3.75i -Y2.5i >> $ps
# R: Cumulative histogram staircase
gmt pshistogram -Bxaf -Byafg100+u" %" -BWSne+glightblue @v3206_06.txt -R-6000/0/0/105 -JX3i/2.0i -S -L1p -Z1 -Q -W200 -O -K -N0+p1p,red -N1+p1p,green -N2+p1p,blue -X3.75i >> $ps
gmt pslegend -Dg0/0+w1.2i+jBR+o0.1i/0.1i -R -J -O -F+p2p+gwhite << EOF >> $ps
S 0.3i - 0.45i - 1p,red 0.7i L@-2@-
S 0.3i - 0.45i - 1p,green 0.7i L@-1@-
S 0.3i - 0.45i - 1p,blue 0.7i LMS
EOF