Next: 7.22 World-wide seismicity the Up: 7. Creating GMT Graphics Previous: 7.20 Custom plot symbols Contents Index

#
7.21 Time-series of RedHat stock price

As discussed in Section 4.4.3, the annotation of time-series is generally more complicated due to the extra degrees of freedom afforded by the dual annotation system. In this example we will display the trend of the stock price of RedHat (RHAT) from their initial public offering until late 2006. The data file is a comma-separated table and the records look like this:

    Date,Open,High,Low,Close,Volume,Adj.Close*
    12-Mar-04,17.74,18.49,17.67,18.02,4827500,18.02
    11-Mar-04,17.60,18.90,17.37,18.09,7700400,18.09

Hence, we have a single header record and various prices in USD for each day of business. We will plot the trend of the opening price as a red line superimposed on a yellow envelope representing the low-to-high fluctuation during each day. We also indicate when and at what cost Paul Wessel bought a few shares, and zoom in on the developments since 2004; in the inset we label the time-axis in Finnish in honor of Linus Thorvalds. Because the time coordinates are Y2K-challenged and the order is backwards (big units of years come after smaller units like days) we must change the default input/output formats used by GMT. Finally, we want to prefix prices with the $ symbol to indicate the currency. Here is how it all comes out:

    #!/bin/sh
    #               GMT EXAMPLE 21
    #
    # Purpose:      Plot a time-series
    # GMT progs:    gmtset, gmtconvert, minmax, psbasemap, psxy
    # Unix progs:   cut, echo
    #
    ps=example_21.ps

    # File has time stored as dd-Mon-yy so set input format to match it

    gmtset INPUT_DATE_FORMAT dd-o-yy PLOT_DATE_FORMAT o ANNOT_FONT_SIZE_PRIMARY +10p
    gmtset TIME_FORMAT_PRIMARY abbreviated CHAR_ENCODING ISOLatin1+

    # Pull out a suitable region string in yyy-mm-dd format

    minmax -fT -I50 -C -H RHAT_price.csv > RHAT.info
    w=`cut -f1 RHAT.info`
    e=`cut -f2 RHAT.info`
    s=`cut -f3 RHAT.info`
    n=`cut -f4 RHAT.info`
    R="-R$w/$e/$s/$n"

    # Lay down the basemap:

    psbasemap $R -JX9i/6i -Glightgreen -K -U"Example 21 in Cookbook" -Bs1Y/WSen \
       -Bpa3Of1o/50WSen:=\$::."RedHat (RHAT) Stock Price Trend since IPO": > $ps

    # Plot main window with open price as red line over yellow envelope of low/highs

    gmtset OUTPUT_DATE_FORMAT dd-o-yy
    gmtconvert -F0,2 -f0T -Hi RHAT_price.csv > RHAT.env
    gmtconvert -F0,3 -f0T -I -Hi RHAT_price.csv >> RHAT.env
    psxy -R -J -Gyellow -O -K RHAT.env >> $ps
    psxy -R -J RHAT_price.csv -H -Wthin,red -O -K >> $ps

    # Draw P Wessel's purchase price as line and label it.  Note we temporary switch
    # back to default yyyy-mm-dd format since that is what minmax gave us.

    echo "05-May-00 0" > RHAT.pw
    echo "05-May-00 300" >> RHAT.pw
    psxy -R -J RHAT.pw -Wthinner,- -O -K >> $ps
    echo "01-Jan-99 25" > RHAT.pw
    echo "01-Jan-07 25" >> RHAT.pw
    psxy -R -J RHAT.pw -Wthick,- -O -K >> $ps
    gmtset INPUT_DATE_FORMAT yyyy-mm-dd
    echo "$w 25 12 0 17 LB Wessel purchase price" | pstext -R -J -O -K -D2i/0.05i -N >> $ps
    gmtset INPUT_DATE_FORMAT dd-o-yy

    # Get smaller region for insert for trend since 2004

    R="-R2004T/$e/$s/40"

    # Lay down the basemap, using Finnish annotations and place the insert in the upper right:

    gmtset TIME_LANGUAGE fi
    psbasemap $R -JX6i/3i -Bpa3Of3o/10:=\$:ESw -Bs1Y/ -Glightblue -O -K -X3i -Y3i >> $ps
    gmtset TIME_LANGUAGE us

    # Again, plot close price as red line over yellow envelope of low/highs

    psxy -R -J -Gyellow -O -K RHAT.env >> $ps
    psxy -R -J RHAT_price.csv -H -Wthin/red -O -K >> $ps

    # Draw P Wessel's purchase price as dashed line

    psxy -R -J RHAT.pw -Wthick,- -O >> $ps

    # Clean up after ourselves:

    rm -f RHAT.* .gmt*

which produces the plot in Figure 7.21, suggesting Wessel has missed a few trains if he had hoped to cash in on the Internet bubble...

Figure 7.21: Time-series of RedHat stock price since IPO.

Next: 7.22 World-wide seismicity the Up: 7. Creating GMT Graphics Previous: 7.20 Custom plot symbols Contents Index