/*
Figure 1.2 and Figure 1.3, pp. 4-5

GCNG    Personal consumption expenditure on gasoline and oil
        (billions of $, SAAR)
GCNGQ   Real personal consumption expenditure on gasoline and oil
        (billions of 87$, SAAR)
GYDQ    Real disposable personal income (billions of 87$, SAAR)
GDC     Implicit deflator for personal consumption expenditure
        (1987 = 100, SA)
The derived series are:
POP16 = The sum of Citibase PM16 and PF16, converted from a monthly
        to a quarterly frequency
      = Civilian nonsinstitutional population, 16 years and over
        (thousands, NSA)
Y     = log(GCNGQ/POP16)
      = Log of per capita real expenditure on gasoline and oil
X2    = log {(GCNG*10^4)/(GCNGQ*GDC)}
      = Log of the real price of gasoline and oil (10 to the power of 4
        is an arbitrary scaling factor)
X3    = log (GYDQ/POP16)
      = Log of per capita real disposable income
Logarithms are to base e.  These series are also analyzed in
Chapter 4 and, with an additional series, in Chapters 8 and 9.
*/

open data auto1.asc
calendar(q) 1959:1
data(format=prn,org=columns) 1959:1 1992:1 gcng gcngq gdc gydq pop16 x2 x3 y
table
*
* Give y, x2 and x3 more descriptive labels
*
label y x2 x3
# "GAS" "PRICE" "INCOME"
graph(key=upleft,Header="Gasoline Consumption and Price",overlay=line,vlabel="Gas",ovlabel="Price") 2
# y
# x2
graph(key=upleft,Header="Gasoline Consumption and Income",overlay=line,vlabel="Gas",ovlabel="Income") 2
# y
# x3
*
spgraph(vfields=2,hfields=2,footer="Figure 1.3 Scatter plots of price and gasoline consumption")
scatter(vlabel=%l(x2),hlabel=%l(y),subheader="1959:1-1992:1",style=line) 1
# y x2 1959:1 1992:1
scatter(vlabel=%l(x2),hlabel=%l(y),subheader="1959:1-1973:3",style=line) 1
# y x2 1959:1 1973:3
scatter(vlabel=%l(x2),hlabel=%l(y),subheader="1973:4-1981:4",style=line) 1
# y x2 1973:4 1981:4
scatter(vlabel=%l(x2),hlabel=%l(y),subheader="1982:1-1992:1",style=line) 1
# y x2 1982:1 1992:1
spgraph(done)