Time Varying Markov switching model

Discussion of models with structural breaks or endogenous switching.

Time Varying Markov switching model

Postby TomDoan » Tue Jan 22, 2008 6:46 pm

Note: An updated version of this is available at

http://www.estima.com/forum/viewtopic.php?f=30&t=513


This is (almost a) replication of parts of Filardo, "Business Cycle Phases and Their Transitional Dynamics", JBES, July 1994, vol 12, no 3, pp 299-307. I'm fairly sure that the code for estimating the time-varying Markov switching model is correct; I'm just not sure about some aspects of the data massaging.



Code: Select all
*
* The data used in the estimation are standard series:
*  IP - Total Industrial Production (1987=100) Federal Reserve Board
*  CLI - Composite Index of Eleven Leading Indicators, BEA
*  DFI - 1 month diffusion Index for CLI, BEA
*  XLI - Stock-Watson Experimental Leading Indicator, from Stock and Watson
*  SP - Standard and Poors 500 Composite Stock Index
*  FED - Federal Funds Rate, Federal Reserve Board
*  TSPREAD - the 10 year less the one year constant maturity Treasury interest rate.
*
cal(m) 1948
open data filardo.dat
data(format=prn,org=cols) 1948:1 1992:8 year month ip cli dfi xli sp fed tspread
*
set ipgr = 100.0*log(ip/ip{1})
*
* The graph in the paper appears to be standardized (probably by dividing through
* by 2.0 times the sample standard deviation).
*
graph(footer="Figure 1. Log Growth Rate of Monthly IP")
# ipgr
*
linreg ipgr * 1959:12
# constant ipgr{1 to 4}
compute stdearly=sqrt(%seesq)
linreg ipgr 1960:1 *
# constant ipgr{1 to 4}
compute stdlate=sqrt(%seesq)
*
set ipgradjust = %if(t<=1959:12,ipgr*stdlate/stdearly,ipgr)
*
set g = ipgradjust
*
* Transform the information variables
*
set cligr   = 100.0*log(cli/cli{1})
set spgr    = 100.0*log(sp/sp{1})
set xlidiff = xli-xli{1}
set ffdiff  = fed-fed{1}
set tspdiff = tspread-tspread{1}
set dfismooth = dfi+2*dfi{1}+2*dfi{2}+dfi{3}
*
* And center them to zero means
*
dofor s = cligr spgr xlidiff ffdiff tspdiff dfismooth
   diff(center) s
end dofor s
*
compute nlags=4
*****************************************************************************************
*
* This is the standard setup for a two state Hamilton model, except that a1 and a2
* are designed to be the log odds ratios for the transition probabilities.
*
compute nstates=2**(nlags+1)
dec vect mu(2) phi(nlags)
dec vect tp(5)
*
nonlin(parmset=standard) mu phi a1 a2 sigma
*
* These will be the formulas which will give the time-varying index values for the
* logistic probabilities.
*
dec frml p1frml p2frml
*
* To speed matters up, create a lookup table which gives the mapping from the
* coding for the states into the binary choice for each lag. The current state is
* 1 for 1-16 and 2 for 17-32. The first lag is 1 for 1-8 and 17-24 and 2 for 9-16
* and 25-32, etc.
*
dec rect[integer] lagstate(nstates,nlags+1)
ewise lagstate(i,j)=1+%mod((i-1)/2**(nlags+1-j),2)
*
* We also create a mapping from state to state which will pick out the transition
* probability. transit will be a number from 1 to 5. 1 is the most common value,
* which represents a zero transition probability. For instance, if we're in
* (1,1,2,2,1), the only states which can be obtained are {1,1,1,2,2} and
* {2,1,1,2,2}. The transition probability to the first of these would be p11, and
* for the second it would be 1-p11. Given our coding, it's easy to tell whether a
* state can be obtained from another, since the final four spots of the new state
* have to be the same as the lead four spots of the old state.
*
dec rect[integer] transit(nstates,nstates)
ewise transit(i,j)=fix(%if(%mod(i-1,2**nlags)==(j-1)/2,2*lagstate(i,1)+lagstate(j,1)-1,1))
dec vect pstar(nstates)
*
* The mask series is used to sum up the probability of being in state 1 at any
* point in time. The prob1 series will hold the result of this calculation.
*
dec vect mask(nstates)
ewise mask(i)=(i<=nstates/2)
*
source markov.src
*
* In order to get smoothed (as opposed to filtered) estimates of the
* probabilities, we need to save the sequence of estimates of the probabilities
* before seeing the period t data, and those after seeing the period t data.  The
* former are saved in the SERIES[VECT] pt_t1, while the latter are in pt_t.
*
dec series[vect] pt_t pt_t1 psmooth
gset pt_t    = %zeros(nstates,1)
gset pt_t1   = %zeros(nstates,1)
gset psmooth = %zeros(nstates,1)
*
*****************************************************************************************
*
* Computes the vector of probabilities (levels, not logs) of the augmented states
* at period <<time>>. G (the data being analyzed), NSTATES, LAGSTATE and the
* parameters MU, PHI and SIGMA are all global variables.
*
function HamiltonF time
type vector HamiltonF
type integer time
*
local integer i j
local real u
*
dim HamiltonF(nstates)
do i=1,nstates
   compute u=g(time)-mu(lagstate(i,1))
   do j=1,nlags
      compute u=u-phi(j)*(g(time-j)-mu(lagstate(i,j+1)))
   end do j
   compute HamiltonF(i)=%density(u/sigma)/sigma
end do i
end
*****************************************************************************************
*
* Computes the vector of probabilities for the augmented states of a Hamilton
* switching model at t given t-1. If the transition probabilities are
* time-varying, they need to be recomputed and put into the tp vector before using
* this. NSTATES, TP, TRANSIT and PSTAR are all global variables.
*
function HamiltonPt_t1
type vector HamiltonPt_t1
*
local integer i j
local real sp
*
dim HamiltonPt_t1(nstates)
*
do i=1,nstates
   compute sp=0.0
   do j=1,nstates
      compute sp=sp+tp(transit(i,j))*pstar(j)
   end do j
   compute HamiltonPt_t1(i)=sp
end do i
end
*****************************************************************************************
*
* This updates the TP array when time varying probabilities are used. It doesn't
* actually have a return value.
*
function TVPProbs time
type integer time
*
local real p11 p22
*
compute p11=%logistic(p1frml(time),1.0),p22=%logistic(p2frml(time),1.0)
compute tp=||0.0,p11,1-p22,1-p11,p22||
end
*****************************************************************************************
*
* pstar (the probability distribution of the states at the previous period) needs
* to be initialized with each function evaluation, as it gets updated at each
* period. We also copy the transition probabilities into the tp vector. The
* initial distribution is set equal to the ergodic probabilities of the states,
* given the current parameters. The solution method is shown on page 684 of
* Hamilton.
*
* The "A" matrix from 684 has 1-the transition probs in the top NxN, with a row
* of 1's below that.
*
function MarkovInit
type vector MarkovInit
*
local rect a
local integer i j
*
compute p11=%logistic(a1,1.0),p22=%logistic(a2,1.0)
compute tp=||0.0,p11,1-p22,1-p11,p22||
dim a(nstates+1,nstates)
ewise a(i,j)=%if(i>nstates,1.0,(i==j)-tp(transit(i,j)))
compute MarkovInit=%xcol(inv(tr(a)*a)*tr(a),nstates+1)
end
*****************************************************************************************
*
procedure MSSmoothTV pt_t pt_t1 psmooth
type series[vect] pt_t pt_t1 *psmooth
*
local rect pfull
local integer time
local vect pstar tp
*
dim pfull(nstates,nstates)
compute pstar=pt_t(%regstart())
gset psmooth %regstart() %regend() = pstar
*
do time=%regend()-1,%regstart(),-1
   compute p11=%logistic(p1frml(time+1),1.0),p22=%logistic(p2frml(time+1),1.0)
   compute tp=||0.0,p11,1-p22,1-p11,p22||
   ewise pfull(i,j)=tp(transit(i,j))
   gset psmooth time time = pt_t.*(tr(pfull)*(psmooth{-1}./pt_t1{-1}))
end do time
end
*
* LOGLFTP computes the log likelihood for a fixed transition probability model
* LOGLTVTP computes the log likelihood for a time-varying transition model
*
frml loglftp = f=HamiltonF(t),pt_t1=HamiltonPt_t1(),pt_t=pstar=%msupdate(f,pt_t1,fpt),log(fpt)
frml logltvtp = TVPProbs(t),f=HamiltonF(t),pt_t1=HamiltonPt_t1(),pt_t=pstar=%msupdate(f,pt_t1,fpt),log(fpt)
*
*****************************************************************************************
*
* Note that the states aren't globally identified. We would get exactly the same
* likelihood if state 1 were contraction and 2 were expansion. The choice of
* initial values for mu and the optimization method can alter whether, in the
* process of estimation, the states "trade places."
*
compute a1=log(.70/(1-.70)),a2=log(.85/(1-.85))
compute mu(1)=0.0,mu(2)=4.00/12.0
linreg g
# constant g{1 to nlags}
compute phi=%xsubvec(%beta,2,nlags+1)
compute sigma=sqrt(%seesq)
*
maximize(parmset=standard,start=(pstar=MarkovInit()),method=bfgs) loglftp %regstart() *
compute baselogl=%funcval,baseregs=%nreg
compute p11f=p11,p22f=p22
*
nonlin(parmset=tv) a11 a21
frml p1frml = a1+a11*cligr{1}
frml p2frml = a2+a21*cligr{1}
maximize(parmset=standard+tv,start=(pstar=MarkovInit()),method=bfgs) logltvtp %regstart() *
cdf(title="Likelihood Ratio Test for TVP") chisqr 2*(%funcval-baselogl) %nreg-baseregs
@MSSmoothTV pt_t pt_t1 psmooth
set p1 %regstart() %regend() = %dot(mask,psmooth)
graph(min=0.0,max=1.0,footer="Smoothed Probability of Recession")
# p1
set p11t = %logistic(p1frml(t),1.0)
set p22t = %logistic(p2frml(t),1.0)
graph(vgrid=||p11f||,min=0.0,max=1.0,footer="Transition Probability: Recession->Recession")
# p11t
graph(vgrid=||p22f||,min=0.0,max=1.0,footer="Transition Probability: Expansion->Expansion")
# p22t


Code: Select all
    YEAR     MONTH    IP       CLI       DFI      XLI      SP      FED      TSPREAD
    1948        1    23.48     37.9    0.000       NA    14.83       NA       NA
    1948        2    23.52     37.5    0.278       NA    14.10       NA       NA
    1948        3    23.26     37.8    0.667       NA    14.30       NA       NA
    1948        4    23.30     38.2    0.611       NA    15.40       NA       NA
    1948        5    23.70     37.7    0.333       NA    16.15       NA       NA
    1948        6    24.01     38.3    0.611       NA    16.82       NA       NA
    1948        7    24.01     38.1    0.200       NA    16.42       NA       NA
    1948        8    23.92     37.7    0.200       NA    15.94       NA       NA
    1948        9    23.74     37.1    0.000       NA    15.76       NA       NA
    1948       10    23.92     37.2    0.700       NA    16.19       NA       NA
    1948       11    23.61     36.9    0.150       NA    15.29       NA       NA
    1948       12    23.39     36.3    0.000       NA    15.19       NA       NA
    1949        1    23.17     35.7    0.200       NA    15.36       NA       NA
    1949        2    22.95     35.5    0.350       NA    14.77       NA       NA
    1949        3    22.51     35.3    0.250       NA    14.91       NA       NA
    1949        4    22.38     34.8    0.300       NA    14.89       NA       NA
    1949        5    22.07     34.9    0.600       NA    14.78       NA       NA
    1949        6    22.02     34.8    0.400       NA    13.97       NA       NA
    1949        7    21.98     35.5    0.900       NA    14.76       NA       NA
    1949        8    22.20     36.3    0.700       NA    15.29       NA       NA
    1949        9    22.42     37.4    0.900       NA    15.49       NA       NA
    1949       10    21.58     37.6    0.550       NA    15.89       NA       NA
    1949       11    22.16     38.3    0.800       NA    16.11       NA       NA
    1949       12    22.55     38.9    0.800       NA    16.54       NA       NA
    1950        1    22.95     39.8    1.000       NA    16.88       NA       NA
    1950        2    23.04     40.3    0.900       NA    17.21       NA       NA
    1950        3    23.79     40.7    0.650       NA    17.35       NA       NA
    1950        4    24.58     41.3    0.800       NA    17.84       NA       NA
    1950        5    25.15     42.4    0.800       NA    18.44       NA       NA
    1950        6    25.90     42.7    0.600       NA    18.74       NA       NA
    1950        7    26.74     44.4    0.800       NA    17.38       NA       NA
    1950        8    27.57     46.1    0.700       NA    18.43       NA       NA
    1950        9    27.40     45.6    0.300       NA    19.08       NA       NA
    1950       10    27.57     45.6    0.500       NA    19.87       NA       NA
    1950       11    27.53     45.4    0.150       NA    19.83       NA       NA
    1950       12    28.01     45.8    0.500       NA    19.75       NA       NA
    1951        1    28.10     46.9    0.650       NA    21.21       NA       NA
    1951        2    28.28     46.2    0.450       NA    22.00       NA       NA
    1951        3    28.41     45.9    0.400       NA    21.63       NA       NA
    1951        4    28.45     44.7    0.500       NA    21.92       NA       NA
    1951        5    28.37     44.2    0.450       NA    21.93       NA       NA
    1951        6    28.23     42.9    0.100       NA    21.55       NA       NA
    1951        7    27.79     42.0    0.200       NA    21.93       NA       NA
    1951        8    27.53     41.2    0.300       NA    22.89       NA       NA
    1951        9    27.71     41.3    0.550       NA    23.48       NA       NA
    1951       10    27.66     41.3    0.600       NA    23.36       NA       NA
    1951       11    27.88     41.0    0.400       NA    22.71       NA       NA
    1951       12    28.06     41.2    0.500       NA    23.41       NA       NA
    1952        1    28.37     41.4    0.600       NA    24.19       NA       NA
    1952        2    28.54     41.7    0.650       NA    23.75       NA       NA
    1952        3    28.63     41.8    0.600       NA    23.81       NA       NA
    1952        4    28.37     41.8    0.450       NA    23.74       NA       NA
    1952        5    28.10     41.9    0.650       NA    23.73       NA       NA
    1952        6    27.84     42.8    0.900       NA    24.38       NA       NA
    1952        7    27.40     42.8    0.600       NA    25.08       NA       NA
    1952        8    29.16     43.5    0.600       NA    25.18       NA       NA
    1952        9    30.22     44.9    0.800       NA    24.78       NA       NA
    1952       10    30.52     44.3    0.500       NA    24.26       NA       NA
    1952       11    31.14     44.3    0.600       NA    25.03       NA       NA
    1952       12    31.32     44.5    0.545       NA    26.04       NA       NA
    1953        1    31.41     44.8    0.545       NA    26.18       NA       NA
    1953        2    31.58     44.9    0.591       NA    25.86       NA       NA
    1953        3    31.85     44.8    0.591       NA    25.99       NA       NA
    1953        4    31.98     44.7    0.409       NA    24.71       NA     0.47
    1953        5    32.15     44.4    0.364       NA    24.84       NA     0.57
    1953        6    32.02     43.8    0.182       NA    23.95       NA     0.66
    1953        7    32.42     43.7    0.273       NA    24.29       NA     0.55
    1953        8    32.24     42.7    0.091       NA    24.39       NA     0.67
    1953        9    31.58     41.8    0.182       NA    23.27       NA     0.67
    1953       10    31.32     41.7    0.364       NA    23.97       NA     0.87
    1953       11    30.57     41.3    0.364       NA    24.50       NA     1.01
    1953       12    29.82     41.4    0.545       NA    24.83       NA     0.93
    1954        1    29.60     41.4    0.545       NA    25.46       NA     1.07
    1954        2    29.69     41.8    0.727       NA    26.02       NA     1.33
    1954        3    29.51     42.0    0.727       NA    26.57       NA     1.24
    1954        4    29.34     42.4    0.909       NA    27.63       NA     1.33
    1954        5    29.51     43.0    0.909       NA    28.73       NA     1.52
    1954        6    29.60     43.6    0.909       NA    28.96       NA     1.56
    1954        7    29.64     44.1    0.864       NA    30.13     0.80     1.46
    1954        8    29.60     44.3    0.636       NA    30.73     1.22     1.48
    1954        9    29.64     44.9    0.727       NA    31.45     1.06     1.35
    1954       10    30.00     45.8    1.000       NA    32.18     0.85     1.26
    1954       11    30.48     46.7    0.773       NA    33.44     0.83     1.34
    1954       12    30.88     47.4    0.727       NA    34.97     1.28     1.30
    1955        1    31.58     48.2    0.909       NA    35.60     1.39     1.22
    1955        2    31.98     49.1    0.909       NA    36.79     1.29     1.08
    1955        3    32.73     49.7    0.636       NA    36.50     1.35     1.09
    1955        4    33.12     49.9    0.682       NA    37.76     1.43     1.00
    1955        5    33.65     50.2    0.636       NA    37.60     1.43     0.86
    1955        6    33.70     50.3    0.545       NA    39.78     1.64     0.87
    1955        7    33.96     50.8    0.682       NA    42.69     1.68     0.88
    1955        8    33.92     50.6    0.364       NA    42.43     1.96     0.60
    1955        9    34.14     51.0    0.636       NA    44.34     2.18     0.61
    1955       10    34.71     50.9    0.455       NA    42.11     2.24     0.49
    1955       11    34.80     51.2    0.636       NA    44.95     2.35     0.41
    1955       12    34.93     51.3    0.455       NA    45.37     2.48     0.23
    1956        1    35.15     51.1    0.318       NA    44.15     2.45     0.32
    1956        2    34.84     50.7    0.136       NA    44.43     2.50     0.35
    1956        3    34.84     50.9    0.364       NA    47.49     2.50     0.35
    1956        4    35.11     51.0    0.409       NA    48.05     2.62     0.26
    1956        5    34.80     50.2    0.091       NA    46.54     2.75     0.13
    1956        6    34.49     49.9    0.364       NA    46.27     2.71     0.26
    1956        7    33.43     50.3    0.545       NA    48.78     2.75     0.35
    1956        8    34.80     50.4    0.545       NA    48.49     2.73     0.23
    1956        9    35.59     50.2    0.409       NA    46.84     2.95     0.03
    1956       10    35.90     50.4    0.409       NA    46.24     2.96     0.06
    1956       11    35.59     50.6    0.591       NA    45.76     2.88     0.05
    1956       12    36.12     50.4    0.364       NA    46.48     2.94    -0.09
    1957        1    35.99     50.0    0.182       NA    45.43     2.84     0.09
    1957        2    36.34     49.8    0.364       NA    43.47     3.00    -0.04
    1957        3    36.29     49.5    0.318       NA    44.03     2.96    -0.01
    1957        4    35.81     49.1    0.273       NA    45.05     3.00    -0.01
    1957        5    35.68     49.1    0.409       NA    46.78     3.00     0.12
    1957        6    35.77     49.0    0.409       NA    47.55     3.00     0.15
    1957        7    35.99     48.7    0.273       NA    48.51     2.99     0.12
    1957        8    35.99     48.6    0.455       NA    45.84     3.24    -0.08
    1957        9    35.68     47.9    0.227       NA    43.98     3.47    -0.15
    1957       10    35.15     47.5    0.364       NA    41.24     3.50    -0.04
    1957       11    34.31     47.0    0.000       NA    40.35     3.28     0.15
    1957       12    33.65     46.6    0.182       NA    40.33     2.98     0.03
    1958        1    33.04     46.5    0.455       NA    41.12     2.72     0.44
    1958        2    32.33     46.1    0.364       NA    41.26     1.67     1.06
    1958        3    31.93     46.5    0.636       NA    42.11     1.20     1.14
    1958        4    31.41     46.7    0.545       NA    42.34     1.26     1.43
    1958        5    31.71     47.6    0.909       NA    43.70     0.63     1.55
    1958        6    32.55     48.5    0.909       NA    44.74     0.93     1.74
    1958        7    33.03     49.7    0.909       NA    45.98     0.68     1.59
    1958        8    33.70     50.6    0.818       NA    47.70     1.53     1.04
    1958        9    34.00     51.5    0.909       NA    48.96     1.76     0.71
    1958       10    34.40     52.1    0.727       NA    50.95     1.80     0.61
    1958       11    35.41     53.0    0.818       NA    52.50     2.27     0.64
    1958       12    35.46     52.9    0.500       NA    53.49     2.42     0.57
    1959        1    35.99     53.8    0.818       NA    55.62     2.48     0.66
    1959        2    36.69     54.5    0.818       NA    54.77     2.43     0.42
    1959        3    37.22     55.3    0.727       NA    56.15     2.80     0.38
    1959        4    38.01     55.3    0.545       NA    57.10     2.96     0.40
    1959        5    38.58     55.4    0.545       NA    57.96     2.90     0.35
    1959        6    38.63     55.4    0.364       NA    57.46     3.39     0.27
    1959        7    37.70     55.2    0.273       NA    59.74     3.47     0.01
    1959        8    36.43     54.7    0.227       NA    59.40     3.50     0.01
    1959        9    36.38     54.8    0.455       NA    57.05     3.76    -0.32
    1959       10    36.12     54.5    0.318       NA    57.00     3.98    -0.27
    1959       11    36.34     54.0    0.273       NA    57.23     4.00    -0.28
    1959       12    38.58     55.2    0.727       NA    59.06     3.99    -0.45
    1960        1    39.60     55.0    0.273     3.08    58.03     3.99    -0.31
    1960        2    39.25     54.3    0.364     3.08    55.78     3.97    -0.17
    1960        3    38.89     53.4    0.091     3.08    55.02     3.84     0.23
    1960        4    38.58     53.4    0.455     3.08    55.73     3.92     0.24
    1960        5    38.54     53.5    0.545     3.08    55.22     3.85     0.14
    1960        6    38.06     53.5    0.455     3.08    57.26     3.32     0.79
    1960        7    37.92     53.7    0.591     3.08    55.84     3.23     0.70
    1960        8    37.88     53.9    0.364     3.08    56.51     2.98     0.85
    1960        9    37.48     54.0    0.455     3.08    54.81     2.60     0.73
    1960       10    37.44     53.8    0.182     0.88    53.73     2.47     0.85
    1960       11    36.91     53.7    0.455    -0.02    55.47     2.44     0.85
    1960       12    36.21     53.7    0.636     1.10    56.80     1.98     0.98
    1961        1    36.25     54.2    0.545     2.45    59.72     1.45     1.03
    1961        2    36.21     54.6    0.727     1.25    62.17     2.54     0.85
    1961        3    36.43     55.4    0.818     2.36    64.12     2.02     0.86
    1961        4    37.18     56.1    0.909     3.67    65.83     1.49     0.90
    1961        5    37.75     56.8    0.864     4.95    66.50     1.98     0.84
    1961        6    38.28     57.4    0.818     6.60    65.62     1.73     0.82
    1961        7    38.72     57.5    0.727     6.99    65.44     1.17     1.00
    1961        8    39.07     58.4    1.000     7.77    67.79     2.00     0.98
    1961        9    39.03     58.0    0.273     6.82    67.26     1.88     0.92
    1961       10    39.77     58.8    0.818     7.00    68.00     2.26     0.87
    1961       11    40.39     59.4    0.636     6.68    71.08     2.61     0.87
    1961       12    40.74     59.9    0.818     6.54    71.74     2.33     0.88
    1962        1    40.39     59.9    0.545     7.31    69.07     2.15     0.80
    1962        2    41.05     60.6    0.727     8.24    70.22     2.37     0.76
    1962        3    41.27     60.6    0.455     6.42    70.29     2.85     0.87
    1962        4    41.36     60.3    0.455     4.91    68.05     2.78     0.85
    1962        5    41.32     59.7    0.182     3.18    62.99     2.36     0.84
    1962        6    41.23     59.1    0.182     3.24    55.63     2.68     0.88
    1962        7    41.62     59.5    0.682     3.32    56.97     2.71     0.72
    1962        8    41.67     59.6    0.455     2.93    58.52     2.93     0.78
    1962        9    41.93     60.0    0.818     3.93    58.00     2.90     0.92
    1962       10    41.98     60.2    0.636     4.66    56.17     2.90     0.95
    1962       11    42.15     61.2    0.909     4.46    60.04     2.94     0.92
    1962       12    42.15     61.7    0.682     5.24    62.64     2.93     0.85
    1963        1    42.46     62.2    0.773     6.12    65.06     2.92     0.79
    1963        2    42.95     62.8    0.818     6.88    65.92     3.00     0.91
    1963        3    43.21     63.3    0.727     7.42    65.67     2.98     0.90
    1963        4    43.61     63.8    0.636     7.37    68.76     2.90     0.86
    1963        5    44.13     64.3    0.727     8.19    70.14     3.00     0.81
    1963        6    44.27     64.2    0.500     6.72    70.11     2.99     0.79
    1963        7    44.09     64.1    0.455     5.02    69.07     3.02     0.54
    1963        8    44.18     64.3    0.545     3.87    70.98     3.49     0.47
    1963        9    44.62     64.8    0.773     4.24    72.85     3.48     0.51
    1963       10    44.93     65.1    0.727     4.45    73.03     3.50     0.47
    1963       11    45.15     65.4    0.500     4.03    72.62     3.48     0.38
    1963       12    45.06     65.6    0.636     4.80    74.17     3.38     0.32
    1964        1    45.46     66.1    0.727     5.79    76.45     3.48     0.38
    1964        2    45.76     66.6    0.636     6.71    77.39     3.48     0.37
    1964        3    45.76     67.0    0.500     6.18    78.80     3.43     0.31
    1964        4    46.51     67.5    0.727     4.51    79.94     3.47     0.32
    1964        5    46.78     68.0    0.636     3.83    80.72     3.50     0.36
    1964        6    46.91     68.2    0.591     4.44    80.24     3.50     0.34
    1964        7    47.22     69.0    0.727     6.00    83.22     3.42     0.48
    1964        8    47.53     69.5    0.818     7.14    82.00     3.50     0.45
    1964        9    47.70     70.1    0.727     6.78    83.41     3.45     0.36
    1964       10    47.04     70.3    0.636     6.31    84.85     3.36     0.33
    1964       11    48.50     70.8    0.727     6.59    85.44     3.52     0.24
    1964       12    49.07     71.2    0.545     5.02    83.96     3.85     0.16
    1965        1    49.60     71.7    0.727     4.86    86.12     3.90     0.25
    1965        2    49.90     71.7    0.500     4.88    86.75     3.98     0.18
    1965        3    50.57     72.1    0.636     4.83    86.83     4.04     0.15
    1965        4    50.79     72.0    0.591     5.21    87.97     4.09     0.16
    1965        5    51.18     72.6    0.727     4.41    89.28     4.10     0.18
    1965        6    51.58     72.5    0.455     4.90    85.04     4.04     0.22
    1965        7    52.06     72.8    0.409     4.76    84.91     4.09     0.22
    1965        8    52.28     72.7    0.545     4.55    86.49     4.12     0.18
    1965        9    52.42     73.3    0.545     5.39    89.38     4.01     0.09
    1965       10    52.94     74.1    0.818     6.49    91.39     4.08     0.05
    1965       11    53.16     74.9    0.727     7.10    92.15     4.10     0.08
    1965       12    53.83     75.6    0.727     7.31    91.73     4.32    -0.10
    1966        1    54.35     75.8    0.682     7.60    93.32     4.42    -0.27
    1966        2    54.71     76.0    0.636     6.54    92.69     4.60    -0.11
    1966        3    55.45     76.9    0.636     6.19    88.88     4.65    -0.10
    1966        4    55.54     76.7    0.455     5.59    91.60     4.67    -0.15
    1966        5    56.07     75.9    0.182     3.72    86.78     4.90    -0.15
    1966        6    56.34     75.4    0.136     1.13    86.06     5.17    -0.16
    1966        7    56.64     75.1    0.227     0.27    85.84     5.30    -0.15
    1966        8    56.69     74.4    0.364     1.68    80.65     5.53    -0.32
    1966        9    57.22     74.0    0.182     2.51    77.81     5.40    -0.64
    1966       10    57.61     73.3    0.136     2.46    77.13     5.53    -0.57
    1966       11    57.22     73.1    0.182     1.94    80.99     5.76    -0.38
    1966       12    57.35     72.8    0.364     0.39    81.33     5.40    -0.36
    1967        1    57.62     73.3    0.636     0.70    84.45     4.94    -0.17
    1967        2    56.97     72.9    0.545    -2.02    87.36     5.00    -0.08
    1967        3    56.64     72.6    0.409    -2.21    89.42     4.53     0.19
    1967        4    57.18     72.9    0.455    -0.79    90.96     4.05     0.49
    1967        5    56.68     73.6    0.727     1.30    92.59     3.94     0.70
    1967        6    56.67     74.5    0.864     5.42    91.43     3.98     0.54
    1967        7    56.54     75.4    0.773     6.18    93.01     3.79     0.15
    1967        8    57.63     76.5    1.000     6.65    94.49     3.90     0.15
    1967        9    57.53     76.7    0.545     5.68    95.81     3.99     0.06
    1967       10    58.00     76.9    0.545     6.50    95.66     3.88     0.11
    1967       11    58.83     77.2    0.591     5.53    92.66     4.13     0.14
    1967       12    59.47     78.1    0.909     5.15    95.30     4.51    -0.01
    1968        1    59.40     78.0    0.364     5.48    95.04     4.60     0.10
    1968        2    59.61     78.7    0.636     7.42    90.75     4.71     0.15
    1968        3    59.80     79.2    0.636     7.01    89.09     5.05     0.16
    1968        4    59.89     78.5    0.455     4.92    95.67     5.76    -0.07
    1968        5    60.56     79.0    0.500     3.78    97.87     6.11    -0.27
    1968        6    60.78     79.2    0.591     2.02   100.53     6.07    -0.26
    1968        7    60.69     79.6    0.545     0.94   100.30     6.02    -0.15
    1968        8    60.86     79.2    0.364     1.78    98.11     6.03    -0.01
    1968        9    61.09     80.2    0.909     3.58   101.34     5.78     0.02
    1968       10    61.21     81.4    0.864     4.23   103.76     5.91     0.01
    1968       11    62.00     81.9    0.727     4.49   105.40     5.82    -0.05
    1968       12    62.20     82.3    0.591     3.58   106.48     6.02    -0.16
    1969        1    62.57     83.1    0.864     4.20   102.04     6.30    -0.30
    1969        2    62.97     82.9    0.364     4.48   101.40     6.61    -0.22
    1969        3    63.47     82.7    0.318     2.58    99.30     6.79    -0.04
    1969        4    63.23     83.1    0.545     1.02   101.26     7.41    -0.09
    1969        5    62.99     82.9    0.318    -1.12   104.62     8.67    -0.10
    1969        6    63.61     82.3    0.227    -2.47    99.13     8.90    -0.47
    1969        7    63.95     81.8    0.182    -2.90    94.71     8.61    -0.88
    1969        8    64.09     81.7    0.318    -3.22    94.18     9.19    -0.85
    1969        9    64.08     81.9    0.591    -3.54    94.51     9.15    -0.66
    1969       10    64.10     81.1    0.182    -3.16    95.52     9.00    -0.54
    1969       11    63.49     80.5    0.136    -2.99    96.21     8.85    -0.75
    1969       12    63.32     80.2    0.364    -1.40    91.11     8.97    -0.52
    1970        1    62.15     79.0    0.091    -3.45    90.31     8.98    -0.31
    1970        2    62.11     77.9    0.091    -4.11    87.16     8.98    -0.35
    1970        3    62.03     77.2    0.182    -2.87    88.65     7.76     0.10
    1970        4    61.87     76.3    0.091    -2.65    85.95     8.10     0.34
    1970        5    61.80     76.7    0.545    -2.73    76.06     7.94     0.16
    1970        6    61.60     76.9    0.455    -4.14    75.59     7.60     0.29
    1970        7    61.75     76.6    0.636    -4.24    75.72     7.21     0.36
    1970        8    61.64     76.6    0.364    -1.99    77.92     6.61     0.55
    1970        9    61.21     76.7    0.636     0.46    82.58     6.29     0.66
    1970       10    59.99     76.5    0.545     1.65    84.37     6.20     0.90
    1970       11    59.63     76.9    0.545     1.76    84.28     5.60     1.33
    1970       12    60.99     78.7    0.773     2.44    90.05     4.90     1.39
    1971        1    61.46     79.8    0.727     2.23    93.49     4.14     1.67
    1971        2    61.35     80.8    0.727     2.08    97.11     3.72     2.22
    1971        3    61.28     81.8    0.773     2.95    99.60     3.71     2.01
    1971        4    61.62     82.8    0.682     4.09   103.04     4.15     1.53
    1971        5    61.94     83.2    0.500     5.72   101.64     4.63     1.35
    1971        6    62.20     83.3    0.545     6.42    99.72     4.91     0.88
    1971        7    62.02     83.3    0.545     5.23    99.00     5.31     0.69
    1971        8    61.66     83.4    0.455     3.87    97.24     5.56     0.78
    1971        9    62.66     83.9    0.636     3.03    99.40     5.55     0.73
    1971       10    63.13     84.6    0.727     3.37    97.29     5.20     1.02
    1971       11    63.40     85.4    0.682     4.56    92.78     4.91     1.14
    1971       12    64.13     87.0    1.000     5.37    99.17     4.14     1.33
    1972        1    65.55     88.1    0.864     6.84   103.30     3.50     1.67
    1972        2    66.00     89.0    0.818     7.13   105.24     3.29     1.81
    1972        3    66.55     89.9    0.682     5.74   107.69     3.83     1.40
    1972        4    67.61     90.3    0.545     5.55   108.81     4.17     1.23
    1972        5    67.45     90.8    0.545     5.24   107.65     4.27     1.49
    1972        6    67.67     91.4    0.818     5.70   108.01     4.46     1.18
    1972        7    67.57     92.5    0.636     7.28   107.21     4.55     1.15
    1972        8    68.49     93.6    0.818     8.15   111.01     4.80     1.23
    1972        9    69.24     94.5    0.682     8.07   109.39     4.87     1.03
    1972       10    70.23     95.0    0.545     8.78   109.56     5.04     0.96
    1972       11    71.07     95.9    0.818     7.95   115.05     5.06     1.01
    1972       12    71.75     96.8    0.818     8.85   117.50     5.33     0.84
    1973        1    71.81     97.3    0.636     7.83   118.42     5.94     0.57
    1973        2    72.80     97.8    0.545     6.15   114.16     6.58     0.45
    1973        3    72.78     97.8    0.455     3.89   112.42     7.09    -0.14
    1973        4    73.05     97.6    0.364     2.71   110.27     7.12    -0.18
    1973        5    73.36     97.6    0.636     2.32   107.22     7.84    -0.04
    1973        6    73.92     97.5    0.364     1.16   104.75     8.49    -0.41
    1973        7    74.38     97.4    0.545    -2.46   105.83    10.40    -1.26
    1973        8    74.33     96.6    0.273    -3.84   103.80    10.50    -1.42
    1973        9    74.90     97.0    0.455    -3.42   105.61    10.78    -1.22
    1973       10    75.18     96.9    0.455    -3.45   109.84    10.01    -0.61
    1973       11    75.21     97.2    0.545    -4.42   102.03    10.03    -0.84
    1973       12    74.02     95.8    0.273    -3.72    94.78     9.95    -0.53
    1974        1    73.00     95.5    0.455    -1.59    96.11     9.65    -0.43
    1974        2    72.71     94.4    0.182     0.29    93.45     8.97     0.08
    1974        3    73.04     94.9    0.591    -0.62    97.44     9.35    -0.55
    1974        4    72.91     93.1    0.273    -3.24    92.46    10.51    -1.11
    1974        5    73.76     92.7    0.364    -5.98    89.67    11.31    -1.20
    1974        6    73.99     91.1    0.091   -10.83    89.79    11.93    -1.13
    1974        7    73.65     90.0    0.227   -11.02    82.82    12.92    -0.99
    1974        8    73.39     87.7    0.045   -12.70    76.03    12.01    -1.32
    1974        9    73.70     85.3    0.091   -11.06    68.12    11.34    -0.83
    1974       10    73.16     83.4    0.136   -10.32    69.44    10.06    -0.15
    1974       11    71.06     81.1    0.091   -11.73    71.74     9.45     0.02
    1974       12    68.11     79.1    0.182    -6.26    67.07     8.53     0.12
    1975        1    66.34     78.0    0.091    -4.29    72.56     7.13     0.67
    1975        2    65.27     78.0    0.364    -2.71    80.10     6.24     1.41
    1975        3    64.11     78.7    0.455    -0.46    83.78     5.54     1.63
    1975        4    64.71     81.2    0.909     2.15    84.72     5.49     1.33
    1975        5    64.51     83.1    0.909     4.63    90.10     5.22     1.67
    1975        6    65.35     84.3    0.909     5.68    92.40     5.55     1.57
    1975        7    65.73     85.7    0.864     7.37    92.49     6.10     0.95
    1975        8    66.90     86.9    0.682     7.55    85.71     6.14     0.70
    1975        9    67.60     88.2    0.727     6.47    84.67     6.24     0.68
    1975       10    67.93     89.2    0.727     6.67    88.57     5.82     1.19
    1975       11    68.56     90.1    0.818     6.66    90.07     5.22     1.56
    1975       12    69.12     90.6    0.636     6.86    88.70     5.20     1.40
    1976        1    69.87     93.0    0.818     6.76    96.86     4.87     1.93
    1976        2    71.11     94.0    0.545     6.48   100.64     4.77     1.88
    1976        3    70.93     94.4    0.545     6.24   101.08     4.84     1.52
    1976        4    71.21     94.4    0.545     7.46   101.93     4.82     1.64
    1976        5    72.02     95.3    0.636     6.66   101.16     5.29     1.49
    1976        6    72.07     96.1    0.636     6.12   101.77     5.48     1.34
    1976        7    72.50     97.3    0.909     7.72   104.20     5.31     1.63
    1976        8    72.91     97.4    0.455     7.80   103.29     5.29     1.77
    1976        9    73.15     97.7    0.545     8.39   105.45     5.25     1.75
    1976       10    73.40     97.5    0.545     7.31   101.89     5.02     1.91
    1976       11    74.62     98.0    0.545     6.82   101.19     4.95     2.00
    1976       12    75.22     98.8    0.636     7.12   104.66     4.65     1.98
    1977        1    75.51     98.5    0.364     6.89   103.81     4.61     1.92
    1977        2    75.87     99.3    0.545     8.24   100.96     4.68     1.92
    1977        3    76.63     99.7    0.455     8.84   100.57     4.69     1.96
    1977        4    77.67    100.3    0.591     8.61    99.05     4.73     1.93
    1977        5    78.34    100.8    0.591     7.56    98.76     5.35     1.62
    1977        6    78.89    101.3    0.636     7.85    99.29     5.39     1.48
    1977        7    78.92    101.1    0.364     6.55   100.18     5.42     1.39
    1977        8    78.98    101.5    0.636     5.87    97.75     5.90     1.03
    1977        9    79.45    101.9    0.591     5.84    96.23     6.14     0.81
    1977       10    79.42    101.8    0.545     6.10    93.74     6.47     0.56
    1977       11    79.52    102.3    0.773     7.27    94.28     6.51     0.63
    1977       12    79.12    103.1    0.636     8.04    93.82     6.56     0.73
    1978        1    78.82    101.9    0.182     7.45    90.25     6.70     0.68
    1978        2    78.95    102.9    0.545     6.97    88.98     6.78     0.69
    1978        3    79.96    103.0    0.455     7.07    88.82     6.79     0.73
    1978        4    82.04    104.1    0.682     7.60    92.71     6.89     0.70
    1978        5    82.34    104.3    0.455     5.92    97.41     7.36     0.53
    1978        6    83.07    104.5    0.545     5.79    97.66     7.60     0.37
    1978        7    83.31    104.6    0.455     4.18    97.19     7.81     0.25
    1978        8    83.62    105.0    0.682     2.85   103.92     8.04     0.10
    1978        9    84.09    106.1    0.727     3.26   103.86     8.45    -0.22
    1978       10    84.51    107.1    0.636     2.87   100.58     8.96    -0.50
    1978       11    85.16    106.3    0.455     2.45    94.71     9.76    -1.21
    1978       12    85.37    105.4    0.455     2.52    96.10    10.03    -1.29
    1979        1    85.07    105.3    0.364    -0.12    99.71    10.07    -1.31
    1979        2    85.82    105.5    0.500    -0.53    98.23    10.06    -1.14
    1979        3    86.12    106.1    0.818     2.26   100.11    10.09    -1.13
    1979        4    85.23    104.5    0.273     2.69   102.07    10.01    -0.94
    1979        5    86.21    105.1    0.545     3.11    99.73    10.24    -0.87
    1979        6    86.11    104.6    0.318     2.36   101.73    10.29    -0.66
    1979        7    85.60    103.1    0.227     0.43   102.71    10.47    -0.69
    1979        8    85.25    102.7    0.273    -1.50   107.36    10.94    -0.95
    1979        9    85.51    103.0    0.636    -2.53   108.60    11.43    -1.51
    1979       10    86.00    101.6    0.136    -4.00   104.47    13.77    -2.14
    1979       11    85.73    100.8    0.182    -5.08   103.66    13.18    -1.74
    1979       12    85.64    100.6    0.364    -5.45   107.78    13.78    -1.59
    1980        1    85.95    101.2    0.727    -4.63   110.87    13.82    -1.26
    1980        2    86.18    101.5    0.591    -6.22   115.34    14.13    -1.51
    1980        3    86.22     98.0    0.182    -8.24   104.69    17.19    -3.07
    1980        4    84.53     95.3    0.091    -4.01   102.97    17.61    -1.83
    1980        5    82.52     93.6    0.273    -2.31   107.69    10.98     0.79
    1980        6    81.52     94.9    0.636     0.74   114.55     9.47     1.62
    1980        7    81.20     96.9    0.818     3.66   119.83     9.03     1.60
    1980        8    82.37     98.9    0.909     3.88   123.50     9.61     0.86
    1980        9    83.47    101.1    1.000     3.92   126.51    10.87    -0.01
    1980       10    84.03    102.5    0.818     0.88   130.22    12.81    -0.74
    1980       11    85.54    104.1    0.773    -2.01   135.65    15.85    -1.47
    1980       12    85.85    102.8    0.545    -2.08   133.48    18.90    -2.04
    1981        1    85.18    102.8    0.227    -1.91   132.97    19.08    -1.51
    1981        2    85.43    101.5    0.182    -1.28   128.40    15.93    -1.38
    1981        3    85.75    102.0    0.500    -0.51   133.19    14.70    -0.59
    1981        4    84.96    103.5    0.909    -3.34   134.43    15.72    -0.64
    1981        5    85.58    103.7    0.364    -3.26   131.73    18.52    -2.10
    1981        6    86.06    102.7    0.364    -3.86   132.28    19.10    -1.39
    1981        7    87.06    102.0    0.182    -5.27   129.13    19.04    -1.44
    1981        8    86.94    102.1    0.545    -5.43   129.63    17.82    -1.78
    1981        9    86.48    100.0    0.000    -3.12   118.27    15.87    -1.20
    1981       10    85.82     98.5    0.364    -1.32   119.80    15.08    -0.23
    1981       11    84.84     97.8    0.364    -0.73   122.92    13.31     0.98
    1981       12    84.12     97.8    0.545    -0.90   123.79    12.37     0.87
    1982        1    82.44     97.2    0.364     0.72   117.28    13.22     0.27
    1982        2    84.16     98.9    0.545     1.92   114.50    14.78    -0.30
    1982        3    83.74     98.3    0.455     0.06   110.84    14.68    -0.09
    1982        4    83.22     99.1    0.500    -0.59   116.31    14.94    -0.11
    1982        5    82.71     99.4    0.591    -0.75   116.35    14.45     0.28
    1982        6    82.43     98.8    0.273    -0.34   109.70    14.15     0.23
    1982        7    81.96     99.4    0.682     1.15   109.38    12.59     0.71
    1982        8    81.58     98.8    0.364    -0.46   109.65    10.12     1.63
    1982        9    81.00    100.6    0.727    -1.04   122.43    10.31     1.49
    1982       10    80.26    101.8    0.909     0.84   132.66     9.71     1.59
    1982       11    79.99    103.0    0.909     2.20   138.10     9.20     1.39
    1982       12    79.27    104.7    0.909     5.50   139.37     8.95     1.63
    1983        1    80.76    106.8    0.727     7.41   144.27     8.68     1.84
    1983        2    80.73    109.1    0.818     7.75   146.80     8.51     1.80
    1983        3    81.35    111.4    0.909     7.95   151.88     8.77     1.47
    1983        4    82.32    113.5    1.000     8.43   157.71     8.80     1.42
    1983        5    83.24    115.3    0.864     8.89   164.10     8.63     1.48
    1983        6    83.73    116.9    0.909     9.25   166.39     8.98     1.19
    1983        7    85.27    118.2    0.909     8.91   166.96     9.37     1.18
    1983        8    86.50    118.4    0.591     7.48   162.42     9.56     1.32
    1983        9    87.94    119.6    0.818     6.49   167.16     9.45     1.49
    1983       10    88.64    121.1    0.909     7.18   167.65     9.48     1.74
    1983       11    88.80    121.9    0.682     7.63   165.23     9.34     1.75
    1983       12    89.24    122.1    0.591     7.85   164.36     9.47     1.72
    1984        1    91.04    123.1    0.727     8.67   166.39     9.56     1.77
    1984        2    90.92    123.4    0.636     9.19   157.25     9.59     1.80
    1984        3    91.86    123.3    0.545     8.27   157.44     9.91     1.73
    1984        4    92.41    123.3    0.455     7.25   157.60    10.29     1.73
    1984        5    93.03    123.0    0.364     6.52   156.55    10.32     1.75
    1984        6    93.51    121.8    0.227     5.37   153.12    11.06     1.48
    1984        7    93.85    121.3    0.455     3.46   151.08    11.23     1.33
    1984        8    93.98    120.9    0.364     1.65   164.42    11.64     0.90
    1984        9    93.88    120.6    0.455     1.98   166.11    11.30     0.94
    1984       10    93.24    119.4    0.273     2.66   164.82     9.99     1.26
    1984       11    93.28    120.1    0.591     4.94   166.27     9.43     1.75
    1984       12    92.75    120.3    0.318     6.02   164.48     8.38     2.17
    1985        1    93.09    121.4    0.773     6.29   171.61     8.35     2.36
    1985        2    93.76    121.7    0.545     7.06   180.88     8.50     2.22
    1985        3    94.08    122.2    0.455     6.84   179.42     8.58     2.00
    1985        4    94.49    121.8    0.364     6.58   180.62     8.27     2.29
    1985        5    94.68    122.7    0.818     5.57   184.90     7.97     2.39
    1985        6    94.40    123.8    0.636     6.27   188.89     7.53     2.36
    1985        7    94.10    124.1    0.409     6.78   192.54     7.88     2.45
    1985        8    94.55    125.1    0.818     6.79   188.31     7.90     2.28
    1985        9    95.01    125.8    0.591     8.03   184.06     7.92     2.30
    1985       10    94.23    126.0    0.727     6.80   186.18     7.99     2.23
    1985       11    94.65    126.2    0.409     5.13   197.45     8.05     1.90
    1985       12    95.63    127.6    0.727     5.22   207.26     8.27     1.59
    1986        1    96.12    128.4    0.636     5.49   208.19     8.14     1.46
    1986        2    95.46    129.4    0.636     4.43   219.37     7.86     1.09
    1986        3    94.59    130.4    0.545     5.07   232.33     7.48     0.75
    1986        4    94.75    131.5    0.636     5.40   237.97     6.99     0.86
    1986        5    94.78    131.5    0.545     3.94   238.46     6.85     1.06
    1986        6    94.39    132.4    0.545     3.48   245.30     6.92     1.07
    1986        7    94.82    132.8    0.364     4.27   240.18     6.56     1.03
    1986        8    94.97    132.9    0.545     5.70   245.00     6.17     1.24
    1986        9    95.06    132.9    0.455     5.53   238.27     5.89     1.68
    1986       10    95.62    134.4    0.682     4.94   237.36     5.85     1.71
    1986       11    96.20    135.3    0.727     4.00   245.09     6.04     1.45
    1986       12    96.69    137.1    0.909     5.66   248.61     6.91     1.24
    1987        1    96.53    136.8    0.318     5.40   264.51     6.43     1.30
    1987        2    97.64    137.7    0.545     4.98   280.93     6.10     1.29
    1987        3    98.18    138.1    0.636     4.88   292.47     6.13     1.22
    1987        4    98.25    138.6    0.591     4.01   289.32     6.37     1.52
    1987        5    99.17    139.1    0.455     3.51   289.12     6.85     1.61
    1987        6   100.13    140.7    0.864     3.49   301.38     6.73     1.60
    1987        7   100.84    142.4    0.727     3.73   310.09     6.58     1.77
    1987        8   100.98    143.1    0.591     2.53   329.36     6.73     1.73
    1987        9   100.94    143.1    0.455     1.83   318.66     7.22     1.75
    1987       10   102.33    142.2    0.545     2.19   280.16     7.29     1.93
    1987       11   102.24    140.2    0.227     1.31   245.01     6.69     1.90
    1987       12   102.60    139.4    0.273     2.50   240.96     6.77     1.82
    1988        1   103.54    139.4    0.500     2.11   250.48     6.83     1.68
    1988        2   103.47    141.0    0.727     3.33   258.13     6.58     1.57
    1988        3   103.89    141.3    0.591     4.01   265.74     6.58     1.66
    1988        4   104.31    141.4    0.545     3.83   262.61     6.87     1.71
    1988        5   104.82    141.8    0.500     3.84   256.12     7.09     1.69
    1988        6   104.97    143.4    0.636     3.29   270.68     7.51     1.43
    1988        7   106.07    142.5    0.136     1.63   269.05     7.75     1.31
    1988        8   106.40    143.3    0.455     1.04   263.73     8.01     1.09
    1988        9   106.19    143.1    0.455     1.33   267.97     8.19     0.89
    1988       10   106.51    143.3    0.318     2.54   277.70     8.30     0.69
    1988       11   106.93    143.5    0.455     2.56   271.02     8.35     0.48
    1988       12   107.43    144.7    0.545     2.25   276.51     8.76     0.12
    1989        1   107.69    145.9    0.727     2.06   285.41     9.12     0.04
    1989        2   107.58    145.6    0.364     0.23   294.01     9.36    -0.08
    1989        3   107.74    144.7    0.318    -0.76   292.71     9.85    -0.21
    1989        4   108.56    145.7    0.636     0.11   302.25     9.84    -0.18
    1989        5   108.25    144.3    0.182     0.48   313.93     9.81    -0.12
    1989        6   108.45    144.2    0.500     0.96   323.73     9.53    -0.16
    1989        7   107.85    144.0    0.364     0.91   331.92     9.24     0.13
    1989        8   108.20    144.6    0.500     1.03   346.61     8.99    -0.07
    1989        9   108.25    145.3    0.455     1.56   347.33     9.02    -0.03
    1989       10   107.70    144.6    0.318     2.62   347.40     8.84     0.02
    1989       11   108.14    145.0    0.500     3.25   340.22     8.55     0.10
    1989       12   108.58    145.7    0.636     4.64   348.57     8.45     0.12
    1990        1   107.52    145.7    0.364     4.94   339.97     8.23     0.29
    1990        2   108.49    144.3    0.364     4.73   330.48     8.24     0.36
    1990        3   108.93    145.4    0.773     2.94   338.60     8.28     0.24
    1990        4   108.83    145.3    0.364     1.88   338.18     8.26     0.39
    1990        5   109.40    146.3    0.636     1.58   350.25     8.18     0.44
    1990        6   110.07    146.3    0.364     2.06   360.39     8.29     0.38
    1990        7   110.40    146.1    0.318     3.07   360.03     8.15     0.53
    1990        8   110.46    144.3    0.273     3.54   333.99     8.13     0.97
    1990        9   110.59    143.2    0.182     2.71   315.41     8.20     1.13
    1990       10   109.95    141.5    0.273     1.91   307.12     8.11     1.17
    1990       11   108.31    139.9    0.364     1.46   315.29     7.81     1.08
    1990       12   107.20    139.6    0.455     1.52   328.75     7.31     1.03
    1991        1   106.62    138.8    0.364     1.91   325.49     6.91     1.45
    1991        2   105.69    140.4    0.636     2.22   362.26     6.25     1.58
    1991        3   104.95    141.5    0.364     3.00   372.28     6.12     1.71
    1991        4   105.51    141.9    0.591     2.91   379.68     5.91     1.80
    1991        5   106.44    143.0    0.727     3.81   378.27     5.78     1.94
    1991        6   107.34    143.9    0.545     4.05   378.29     5.90     1.92
    1991        7   108.06    145.6    0.773     4.78   380.23     5.82     1.96
    1991        8   107.98    145.6    0.364     5.30   389.40     5.66     2.12
    1991        9   108.38    145.5    0.409     5.22   387.20     5.45     2.08
    1991       10   108.40    145.7    0.455     5.20   386.88     5.21     2.20
    1991       11   108.06    145.5    0.364     5.11   385.87     4.81     2.53
    1991       12   107.36    145.3    0.500     7.10   388.51     4.43     2.71
    1992        1   106.64    146.9    0.636     8.33   416.08     4.03     2.88
    1992        2   107.20    148.0    0.727     8.43   412.56     4.06     3.05
    1992        3   107.58    148.5    0.409     7.15   407.36     3.98     2.91
    1992        4   108.15    149.0    0.455     6.84   407.41     3.73     3.18
    1992        5   108.92    149.9    0.455     6.15   414.81     3.82     3.20
    1992        6   108.49    149.5    0.364     5.55   408.27     3.76     3.09
    1992        7   109.16    149.7    0.500     5.95   415.05     3.25     3.24
    1992        8   108.65    149.4    0.409     6.04   417.93     3.30     3.12
TomDoan
 
Posts: 2720
Joined: Wed Nov 01, 2006 5:36 pm

TVTP forecasting

Postby sareilly » Wed Aug 13, 2008 5:28 am

Does anybody have any ideas on how to use this model to forecast one step ahead (dependant on the time varying variables in the transition probabilities)
Any help would be really appreciated!
Sarah
sareilly
 
Posts: 1
Joined: Thu Jun 05, 2008 6:09 am

Re: Time Varying Markov switching model

Postby GUESMI2010 » Fri Jan 08, 2010 11:25 am

hi,
i used the replication of parts of Filardo, "Business Cycle Phases and Their Transitional Dynamics", JBES, July 1994, vol 12, no 3, pp 299-307. i have this message
## SX21. A END or } Here is Unneeded or Unexpected
## SX11. Identifier %MSPLOGISTIC is Not Recognizable. Incorrect Option Field or Parameter Order?
>>>>theta=%MSPLogistic(<<<<
can you help me?
GUESMI2010
 
Posts: 5
Joined: Wed Jan 06, 2010 4:02 pm

Re: Time Varying Markov switching model

Postby moderator » Wed Jan 13, 2010 11:43 am

This was addressed by email (the user needed to download the current version of markov.src).
moderator
Site Admin
 
Posts: 306
Joined: Thu Oct 19, 2006 4:33 pm


Return to Structural Breaks and Switching Models

Who is online

Users browsing this forum: No registered users and 1 guest