# Word wrap may facilitate visualization

###### Data for 1984-2002 ###################


K.JAN1984 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0) 
K.FEB1984 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.MAR1984 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, 0, 0) 
K.APR1984 = c(0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0) 
K.MAY1984 = c(1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0) 
K.JUN1984 = c(1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) 
K.JUL1984 = c(1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1,
0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1) 
K.AUG1984 = c(1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1,
0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.SEP1984 = c(0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1,
0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.OCT1984 = c(0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1,
1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1984 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.DEC1984 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1985 = c(0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1985 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1985 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0)
K.APR1985 = c(0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0)
K.MAY1985 = c(1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0,
0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0) 
K.JUN1985 = c(0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1,
1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1) 
K.JUL1985 = c(1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1,
0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1) 
K.AUG1985 = c(1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1,
1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1)
K.SEP1985 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1) 
K.OCT1985 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0,
0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1985 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.DEC1985 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1986 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1986 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1986 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1986 = c(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0)
K.MAY1986 = c(1, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0)
K.JUN1986 = c(0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1,
1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1)
K.JUL1986 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1,
1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0)
K.AUG1986 = c(1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1,
1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 1)
K.SEP1986 = c(0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1,
1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.OCT1986 = c(1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0,
0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1986 = c(0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1986 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0)

K.JAN1987 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1987 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)
K.MAR1987 = c(0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1987 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0)
K.MAY1987 = c(1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0)
K.JUN1987 = c(0, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0,
1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0)
K.JUL1987 = c(0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.AUG1987 = c(1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1,
1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.SEP1987 = c(0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1,
1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 1)
K.OCT1987 = c(0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1987 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1987 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)


K.JAN1988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0)
K.MAR1988 = c(0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, 0, 0)
K.APR1988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1)
K.MAY1988 = c(0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1,
1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1)
K.JUN1988 = c(1, 0, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0)
K.JUL1988 = c(0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1)
K.AUG1988 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 0)
K.SEP1988 = c(0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1,
1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0)
K.OCT1988 = c(0, 1, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1,
1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1)
K.DEC1988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0)
K.MAR1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0)
K.APR1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1)
K.MAY1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1)
K.JUN1989 = c(0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1,
1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1)
K.JUL1989 = c(1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1,
1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0)
K.AUG1989 = c(1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1,
1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1)
K.SEP1989 = c(0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1,
1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1)
K.OCT1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0)
K.NOV1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1989 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0)

K.JAN1990 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1990 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0)
K.MAR1990 = c(0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 0)
K.APR1990 = c(1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0)
K.MAY1990 = c(0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1,
1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0)
K.JUN1990 = c(1, 1, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0,
1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0)
K.JUL1990 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0)
K.AUG1990 = c(1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1,
1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 0)
K.SEP1990 = c(0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1,
1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1)
K.OCT1990 = c(1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0, 0,
0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1990 = c(0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1990 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1991 = c(0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1991 = c(0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1991 = c(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1991 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1)
K.MAY1991 = c(0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0,
1, 0, 1, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1)
K.JUN1991 = c(0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 	0)
K.JUL1991 = c(0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0,
0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1)
K.AUG1991 = c(0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1)
K.SEP1991 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0)
K.OCT1991 = c(0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1,
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1)
K.NOV1991 = c(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1991 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0)

K.JAN1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1992 = c(1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0)
K.MAR1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAY1992 = c(1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1)
K.JUN1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0)
K.JUL1992 = c(1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0,
1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1)
K.AUG1992 = c(1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1,
1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0)
K.SEP1992 = c(0, 0, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0)
K.OCT1992 = c(0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0,
0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.DEC1992 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1993 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1993 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0)
K.MAR1993 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0)
K.APR1993 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1,
0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 1)
K.MAY1993 = c(0, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0,
1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1)
K.JUN1993 = c(0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0,
0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0)
K.JUL1993 = c(0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0)
K.AUG1993 = c(1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1,
0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1)
K.SEP1993 = c(1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1,
1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 1)
K.OCT1993 = c(0, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0,
1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1993 = c(1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1993 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1994 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1994 = c(0, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1994 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1994 = c(0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0)
K.MAY1994 = c(1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0,
1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0)
K.JUN1994 = c(0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1,
1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1)
K.JUL1994 = c(0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1)
K.AUG1994 = c(0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 1,
1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0)
K.SEP1994 = c(1, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 0,
0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.OCT1994 = c(0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1994 = c(0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0)
K.DEC1994 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1995 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1995 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1)
K.MAR1995 = c(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1)
K.APR1995 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAY1995 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1,
1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.JUN1995 = c(0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1,
1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0)
K.JUL1995 = c(1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1,
1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1)
K.AUG1995 = c(1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0)
K.SEP1995 = c(1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0, 1)
K.OCT1995 = c(0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0,
0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1)
K.NOV1995 = c(1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0)
K.DEC1995 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1996 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1996 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0)
K.MAR1996 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,0, 0, 0)
K.APR1996 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0)
K.MAY1996 = c(0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0,
0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)
K.JUN1996 = c(0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0,
0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1)
K.JUL1996 = c(1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0,
0, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.AUG1996 = c(1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0,
0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1)
K.SEP1996 = c(0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0,
0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0) 
K.OCT1996 = c(0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0)
K.NOV1996 = c(0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1996 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1997 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1997 = c(1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1997 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1)
K.APR1997 = c(0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 0,
0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0)
K.MAY1997 = c(0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 0,
0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0)
K.JUN1997 = c(0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 1,
1, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.JUL1997 = c(1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1,
1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1)
K.AUG1997 = c(1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1,
1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1)
K.SEP1997 = c(0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1)
K.OCT1997 = c(1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV1997 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1997 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0,
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1998 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1998 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 0)
K.MAR1998 = c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1)
K.APR1998 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0,
0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0)
K.MAY1998 = c(0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0)
K.JUN1998 = c(1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0,
1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0)
K.JUL1998 = c(1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 1)
K.AUG1998 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 1, 0, 1)
K.SEP1998 = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0)
K.OCT1998 = c(0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1,
1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1)
K.NOV1998 = c(1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0)
K.DEC1998 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.APR1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAY1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0,
1, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 0)
K.JUN1999 = c(0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0,
0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0)
K.JUL1999 = c(1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1,
1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1)
K.AUG1999 = c(0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1,
1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1)
K.SEP1999 = c(1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 1, 1,
1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0)
K.OCT1999 = c(1, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0,
0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0)
K.NOV1999 = c(0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC1999 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN2000 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0)
K.FEB2000 = c(0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR2000 = c(0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0,
0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0, 0, 0)
K.APR2000 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1)
K.MAY2000 = c(0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0,
0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 1)
K.JUN2000 = c(1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1,
1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0)
K.JUL2000 = c(1, 0, 0, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0)
K.AUG2000 = c(1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1)
K.SEP2000 = c(1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0) 
K.OCT2000 = c(0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0)
K.NOV2000 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
K.DEC2000 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN2001 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.FEB2001 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0)
K.MAR2001 = c(1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1)
K.APR2001 = c(0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0,
0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAY2001 = c(1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0,
0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0)
K.JUN2001 = c(1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1,
1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1)
K.JUL2001 = c(1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1)
K.AUG2001 = c(1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 0,
1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1)
K.SEP2001 = c(1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0)
K.OCT2001 = c(1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 0,
1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1)
K.NOV2001 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC2001 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.JAN2002 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0)
K.FEB2002 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.MAR2002 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0)
K.APR2002 = c(1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1)
K.MAY2002 = c(0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 0,
0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1)
K.JUN2002 = c(1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1,
1, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1)
K.JUL2002 = c(1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0,
0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1)
K.AUG2002 = c(1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1)
K.SEP2002 = c(1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1,
0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1)
K.OCT2002 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.NOV2002 = c(0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)
K.DEC2002 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)

K.rain1984<-c(K.JAN1984,K.FEB1984,K.MAR1984,K.APR1984,K.MAY1984,K.JUN1984,K.JUL1984,K.AUG1984,K.SEP1984,K.OCT1984,K.NOV1984,K.DEC1984)
K.rain1985<-c(K.JAN1985,K.FEB1985,K.MAR1985,K.APR1985,K.MAY1985,K.JUN1985,K.JUL1985,K.AUG1985,K.SEP1985,K.OCT1985,K.NOV1985,K.DEC1985)
K.rain1986<-c(K.JAN1986,K.FEB1986,K.MAR1986,K.APR1986,K.MAY1986,K.JUN1986,K.JUL1986,K.AUG1986,K.SEP1986,K.OCT1986,K.NOV1986,K.DEC1986)
K.rain1987<-c(K.JAN1987,K.FEB1987,K.MAR1987,K.APR1987,K.MAY1987,K.JUN1987,K.JUL1987,K.AUG1987,K.SEP1987,K.OCT1987,K.NOV1987,K.DEC1987)
K.rain1988<-c(K.JAN1988,K.FEB1988,K.MAR1988,K.APR1988,K.MAY1988,K.JUN1988,K.JUL1988,K.AUG1988,K.SEP1988,K.OCT1988,K.NOV1988,K.DEC1988)
K.rain1989<-c(K.JAN1989,K.FEB1989,K.MAR1989,K.APR1989,K.MAY1989,K.JUN1989,K.JUL1989,K.AUG1989,K.SEP1989,K.OCT1989,K.NOV1989,K.DEC1989)
K.rain1990<-c(K.JAN1990,K.FEB1990,K.MAR1990,K.APR1990,K.MAY1990,K.JUN1990,K.JUL1990,K.AUG1990,K.SEP1990,K.OCT1990,K.NOV1990,K.DEC1990)
K.rain1991<-c(K.JAN1991,K.FEB1991,K.MAR1991,K.APR1991,K.MAY1991,K.JUN1991,K.JUL1991,K.AUG1991,K.SEP1991,K.OCT1991,K.NOV1991,K.DEC1991)
K.rain1992<-c(K.JAN1992,K.FEB1992,K.MAR1992,K.APR1992,K.MAY1992,K.JUN1992,K.JUL1992,K.AUG1992,K.SEP1992,K.OCT1992,K.NOV1992,K.DEC1992)
K.rain1993<-c(K.JAN1993,K.FEB1993,K.MAR1993,K.APR1993,K.MAY1993,K.JUN1993,K.JUL1993,K.AUG1993,K.SEP1993,K.OCT1993,K.NOV1993,K.DEC1993)
K.rain1994<-c(K.JAN1994,K.FEB1994,K.MAR1994,K.APR1994,K.MAY1994,K.JUN1994,K.JUL1994,K.AUG1994,K.SEP1994,K.OCT1994,K.NOV1994,K.DEC1994)
K.rain1995<-c(K.JAN1995,K.FEB1995,K.MAR1995,K.APR1995,K.MAY1995,K.JUN1995,K.JUL1995,K.AUG1995,K.SEP1995,K.OCT1995,K.NOV1995,K.DEC1995)
K.rain1996<-c(K.JAN1996,K.FEB1996,K.MAR1996,K.APR1996,K.MAY1996,K.JUN1996,K.JUL1996,K.AUG1996,K.SEP1996,K.OCT1996,K.NOV1996,K.DEC1996)
K.rain1997<-c(K.JAN1997,K.FEB1997,K.MAR1997,K.APR1997,K.MAY1997,K.JUN1997,K.JUL1997,K.AUG1997,K.SEP1997,K.OCT1997,K.NOV1997,K.DEC1997)
K.rain1998<-c(K.JAN1998,K.FEB1998,K.MAR1998,K.APR1998,K.MAY1998,K.JUN1998,K.JUL1998,K.AUG1998,K.SEP1998,K.OCT1998,K.NOV1998,K.DEC1998)
K.rain1999<-c(K.JAN1999,K.FEB1999,K.MAR1999,K.APR1999,K.MAY1999,K.JUN1999,K.JUL1999,K.AUG1999,K.SEP1999,K.OCT1999,K.NOV1999,K.DEC1999)
K.rain2000<-c(K.JAN2000,K.FEB2000,K.MAR2000,K.APR2000,K.MAY2000,K.JUN2000,K.JUL2000,K.AUG2000,K.SEP2000,K.OCT2000,K.NOV2000,K.DEC2000)
K.rain2001<-c(K.JAN2001,K.FEB2001,K.MAR2001,K.APR2001,K.MAY2001,K.JUN2001,K.JUL2001,K.AUG2001,K.SEP2001,K.OCT2001,K.NOV2001,K.DEC2001)
K.rain2002<-c(K.JAN2002,K.FEB2002,K.MAR2002,K.APR2002,K.MAY2002,K.JUN2002,K.JUL2002,K.AUG2002,K.SEP2002,K.OCT2002,K.NOV2002,K.DEC2002)



K.JAN<-c(K.JAN1984,K.JAN1985,K.JAN1986,K.JAN1987,K.JAN1988,K.JAN1989,K.JAN1990,K.JAN1991,K.JAN1992,K.JAN1993,K.JAN1994,K.JAN1995,K.JAN1996,K.JAN1997,K.JAN1998,K.JAN1999,K.JAN2000,K.JAN2001,K.JAN2002)
K.FEB<-c(K.FEB1984,K.FEB1985,K.FEB1986,K.FEB1987,K.FEB1988,K.FEB1989,K.FEB1990,K.FEB1991,K.FEB1992,K.FEB1993,K.FEB1994,K.FEB1995,K.FEB1996,K.FEB1997,K.FEB1998,K.FEB1999,K.FEB2000,K.FEB2001,K.FEB2002)
K.MAR<-c(K.MAR1984,K.MAR1985,K.MAR1986,K.MAR1987,K.MAR1988,K.MAR1989,K.MAR1990,K.MAR1991,K.MAR1992,K.MAR1993,K.MAR1994,K.MAR1995,K.MAR1996,K.MAR1997,K.MAR1998,K.MAR1999,K.MAR2000,K.MAR2001,K.MAR2002)
K.APR<-c(K.APR1984,K.APR1985,K.APR1986,K.APR1987,K.APR1988,K.APR1989,K.APR1990,K.APR1991,K.APR1992,K.APR1993,K.APR1994,K.APR1995,K.APR1996,K.APR1997,K.APR1998,K.APR1999,K.APR2000,K.APR2001,K.APR2002)
K.MAY<-c(K.MAY1984,K.MAY1985,K.MAY1986,K.MAY1987,K.MAY1988,K.MAY1989,K.MAY1990,K.MAY1991,K.MAY1992,K.MAY1993,K.MAY1994,K.MAY1995,K.MAY1996,K.MAY1997,K.MAY1998,K.MAY1999,K.MAY2000,K.MAY2001,K.MAY2002)
K.JUN<-c(K.JUN1984,K.JUN1985,K.JUN1986,K.JUN1987,K.JUN1988,K.JUN1989,K.JUN1990,K.JUN1991,K.JUN1992,K.JUN1993,K.JUN1994,K.JUN1995,K.JUN1996,K.JUN1997,K.JUN1998,K.JUN1999,K.JUN2000,K.JUN2001,K.JUN2002)
K.JUL<-c(K.JUL1984,K.JUL1985,K.JUL1986,K.JUL1987,K.JUL1988,K.JUL1989,K.JUL1990,K.JUL1991,K.JUL1992,K.JUL1993,K.JUL1994,K.JUL1995,K.JUL1996,K.JUL1997,K.JUL1998,K.JUL1999,K.JUL2000,K.JUL2001,K.JUL2002)
K.AUG<-c(K.AUG1984,K.AUG1985,K.AUG1986,K.AUG1987,K.AUG1988,K.AUG1989,K.AUG1990,K.AUG1991,K.AUG1992,K.AUG1993,K.AUG1994,K.AUG1995,K.AUG1996,K.AUG1997,K.AUG1998,K.AUG1999,K.AUG2000,K.AUG2001,K.AUG2002)
K.SEP<-c(K.SEP1984,K.SEP1985,K.SEP1986,K.SEP1987,K.SEP1988,K.SEP1989,K.SEP1990,K.SEP1991,K.SEP1992,K.SEP1993,K.SEP1994,K.SEP1995,K.SEP1996,K.SEP1997,K.SEP1998,K.SEP1999,K.SEP2000,K.SEP2001,K.SEP2002)
K.OCT<-c(K.OCT1984,K.OCT1985,K.OCT1986,K.OCT1987,K.OCT1988,K.OCT1989,K.OCT1990,K.OCT1991,K.OCT1992,K.OCT1993,K.OCT1994,K.OCT1995,K.OCT1996,K.OCT1997,K.OCT1998,K.OCT1999,K.OCT2000,K.OCT2001,K.OCT2002)
K.NOV<-c(K.NOV1984,K.NOV1985,K.NOV1986,K.NOV1987,K.NOV1988,K.NOV1989,K.NOV1990,K.NOV1991,K.NOV1992,K.NOV1993,K.NOV1994,K.NOV1995,K.NOV1996,K.NOV1997,K.NOV1998,K.NOV1999,K.NOV2000,K.NOV2001,K.NOV2002)
K.DEC<-c(K.DEC1984,K.DEC1985,K.DEC1986,K.DEC1987,K.DEC1988,K.DEC1989,K.DEC1990,K.DEC1991,K.DEC1992,K.DEC1993,K.DEC1994,K.DEC1995,K.DEC1996,K.DEC1997,K.DEC1998,K.DEC1999,K.DEC2000,K.DEC2001,K.DEC2002)





###### program to calculate probabilities for alternate renewal process  using Green (1964) paper ###################
n<-length(K.JUL)
I<-table(K.JUL[-n])
C<- c(I[2],I[1])
J<-table(K.JUL[-n],K.JUL[-1])
TCM<-matrix(c(J[2,2],J[2,1],J[1,2],J[1,1]),2,2,byrow=T)
# Its not usual TCM, but the row-columns are in 1-0 format, not 0-1
K<-table(K.JUL[-((n-1):n)],K.JUL[-c(1,n)],K.JUL[-c(1,2)])
A<-matrix(c(K[2,2,2],K[2,1,2],K[1,2,2],K[1,1,2],K[2,2,1],K[2,1,1],K[1,2,1],K[1,1,1]),4,2,byrow=T)
L<-
table(K.JUL[-((n-2):n)],K.JUL[-c(1,n,n-1)],K.JUL[-c(1,2,n)],K.JUL[-c(1,2,3)])
B<-matrix(c(L[2,2,2,2],L[2,1,2,2],L[1,2,2,2],L[1,1,2,2],L[2,2,1,2],L[2,1,1,2],L[1,2,1,2],L[1,1,1,2],L[2,2,2,1],L[2,1,2,1],L[1,2,2,1],L[1,1,2,1],L[2,2,1,1],L[2,1,1,1],L[1,2,1,1],L[1,1,1,1]),8,2,byrow=T)
Total<-
c(C[1]+C[2],TCM[1,1]+TCM[1,2],TCM[2,1]+TCM[2,2],A[1,1]+A[1,2],A[2,1]+A[2,2],A[3,1]+A[3,2],A[4,1]+A[4,2],B[1,1]+B[1,2],B[2,1]+B[2,2],B[3,1]+B[3,2],B[4,1]+B[4,2],B[5,1]+B[5,2],B[6,1]+B[6,2],B[7,1]+B[7,2],B[8,1]+B[8,2])
RFWD<-
c(C[1]/(C[1]+C[2]),TCM[1,1]/(TCM[1,1]+TCM[1,2]),TCM[2,1]/(TCM[2,1]+TCM[2,2]),A[1,1]/(A[1,1]+A[1,2]),A[2,1]/(A[2,1]+A[2,2]),A[3,1]/(A[3,1]+A[3,2]),A[4,1]/(A[4,1]+A[4,2]),B[1,1]/(B[1,1]+B[1,2]),B[2,1]/(B[2,1]+B[2,2]),B[3,1]/(B[3,1]+B[3,2]),B[4,1]/(B[4,1]+B[4,2]),B[5,1]/(B[5,1]+B[5,2]),B[6,1]/(B[6,1]+B[6,2]),B[7,1]/(B[7,1]+B[7,2]),B[8,1]/(B[8,1]+B[8,2]))
actual<-rbind(C,TCM,A,B)
cbind(actual,Total,RFWD)
TPM<-matrix(c(TCM[1,1]/(TCM[1,1]+TCM[1,2]),TCM[1,2]/(TCM[1,1]+TCM[1,2]),TCM[2,1]/(TCM[2,1]+TCM[2,2]),TCM[2,2]/(TCM[2,1]+TCM[2,2])),2,byrow=T)
Alpha<--log(1-TPM[2,1]);
Beta<-Alpha*(1-TPM[1,1])/(TPM[2,1]*(TPM[1,1]-TPM[2,1]));
d<-Beta/(Alpha+Beta);a<-exp(-Alpha)*d*(1-exp(-Alpha-Beta));
b<-exp(-Alpha)*(1-exp(-Beta));

Q0<-c(C[1]/(C[1]+C[2]))
Q1<-(1-b)+(b-a)/Q0;
Q2<-(1-b)+(b-a)/Q1;
Q3<-(1-b)+(b-a)/Q2;
Q4<-(1-b)+(b-a)/Q3;
Q5<-(1-b)+(b-a)/Q4;
Q6<-(1-b)+(b-a)/Q5;
Q7<-(1-b)+(b-a)/Q6;
Q8<-(1-b)+(b-a)/Q7;
Q9<-(1-b)+(b-a)/Q8;
Q10<-(1-b)+(b-a)/Q9;
Q11<-(1-b)+(b-a)/Q10;
Q12<-(1-b)+(b-a)/Q11;
Q13<-(1-b)+(b-a)/Q12;
Q14<-(1-b)+(b-a)/Q13;
Q15<-(1-b)+(b-a)/Q14;
Q16<-(1-b)+(b-a)/Q15;
Q17<-(1-b)+(b-a)/Q16;
Q18<-(1-b)+(b-a)/Q17;
Q19<-(1-b)+(b-a)/Q18;
Q20<-(1-b)+(b-a)/Q19;
Q21<-(1-b)+(b-a)/Q20;
Q22<-(1-b)+(b-a)/Q21;
Q23<-(1-b)+(b-a)/Q22;
Q24<-(1-b)+(b-a)/Q23;
Q25<-(1-b)+(b-a)/Q24;
Q26<-(1-b)+(b-a)/Q25;
Q27<-(1-b)+(b-a)/Q26;
Q28<-(1-b)+(b-a)/Q27;
Q29<-(1-b)+(b-a)/Q28;
Q30<-(1-b)+(b-a)/Q29; 

P0<-TPM[2,1];
P1<--b+a/(1-Q1);
P2<--b+a/(1-Q2);
P3<--b+a/(1-Q3);
P4<--b+a/(1-Q4);
P5<--b+a/(1-Q5);
P6<--b+a/(1-Q6);
P7<--b+a/(1-Q7);
P8<--b+a/(1-Q8);
P9<--b+a/(1-Q9);
P10<--b+a/(1-Q10);
P11<--b+a/(1-Q11);
P12<--b+a/(1-Q12);
P13<--b+a/(1-Q13);
P14<--b+a/(1-Q14);
P15<--b+a/(1-Q15);
P16<--b+a/(1-Q16);
P17<--b+a/(1-Q17);
P18<--b+a/(1-Q18);
P19<--b+a/(1-Q19);
P20<--b+a/(1-Q20);
P21<--b+a/(1-Q21);
P22<--b+a/(1-Q22);
P23<--b+a/(1-Q23);
P24<--b+a/(1-Q24);
P25<--b+a/(1-Q25);
P26<--b+a/(1-Q26);
P27<--b+a/(1-Q27);
P28<--b+a/(1-Q28);
P29<--b+a/(1-Q29);
P30<--b+a/(1-Q30);

green.prob<-c(Q0,Q1,P0,Q2,P1,P0,P0,Q3,P2,P1,P1,P0,P0,P0,P0);
gn.prob<-c(NA,Q1,P0,Q1,Q1,P0,P0,Q1,Q1,Q1,Q1,P0,P0,P0,P0);
cbind(actual,Total,RFWD,green.prob,gn.prob)
W0<-1
DW0<-C[2]/(C[1]+C[2])
W1<-W0-DW0
DW1<-a*W0-b*DW0
W2<-W1-DW1
DW2<-a*W1-b*DW1
W3<-W2-DW2
DW3<-a*W2-b*DW2
W4<-W3-DW3
DW4<-a*W3-b*DW3
W5<-W4-DW4
DW5<-a*W4-b*DW4
W6<-W5-DW5
DW6<-a*W5-b*DW5
W7<-W6-DW6
DW7<-a*W6-b*DW6
W8<-W7-DW7
DW8<-a*W7-b*DW7
W9<-W8-DW8
DW9<-a*W8-b*DW8
W10<-W9-DW9
DW10<-a*W9-b*DW9
W11<-W10-DW10
DW11<-a*W10-b*DW10
W12<-W11-DW11
DW12<-a*W11-b*DW11
W13<-W12-DW12
DW13<-a*W12-b*DW12
W14<-W13-DW13
DW14<-a*W13-b*DW13
W15<-W14-DW14
DW15<-a*W14-b*DW14
W16<-W15-DW15
DW16<-a*W15-b*DW15
W17<-W16-DW16
DW17<-a*W16-b*DW16
W18<-W17-DW17
DW18<-a*W17-b*DW17
W19<-W18-DW18
DW19<-a*W18-b*DW18
W20<-W19-DW19
DW20<-a*W19-b*DW19
W21<-W20-DW20
DW21<-a*W20-b*DW20
W22<-W21-DW21
DW22<-a*W21-b*DW21
W23<-W22-DW22
DW23<-a*W22-b*DW22
W24<-W23-DW23
DW24<-a*W23-b*DW23
W25<-W24-DW24
DW25<-a*W24-b*DW24
W26<-W25-DW25
DW26<-a*W25-b*DW25
W27<-W26-DW26
DW27<-a*W26-b*DW26
W28<-W27-DW27
DW28<-a*W27-b*DW27
W29<-W28-DW28
DW29<-a*W28-b*DW28
W30<-W29-DW29
DW30<-a*W29-b*DW29
W31<-W30-DW30
DW31<-a*W30-b*DW30


DWn<-c(DW1,
DW2,DW3,DW4,DW5,DW6,DW7,DW8,DW9,DW10,DW11,DW12,DW13,DW14,DW15,DW16,DW17,DW18,DW19,DW20,DW21,DW22,DW23,DW24,DW25,DW26,DW27,DW28,DW29,DW30)
Wn<-c(W1,
W2,W3,W4,W5,W6,W7,W8,W9,W10,W11,W12,W13,W14,W15,W16,W17,W18,W19,W20,W21,W22,W23,W24,W25,W26,W27,W28,W29,W30)
cbind(DWn,Wn)

W11D.DW<-(DW1-DW2)/DW1
W21D.DW<-(DW2-DW3)/DW1
W31D.DW<-(DW3-DW4)/DW1
W41D.DW<-(DW4-DW5)/DW1
W51D.DW<-(DW5-DW6)/DW1
W61D.DW<-(DW6-DW7)/DW1
W71D.DW<-(DW7-DW8)/DW1
W81D.DW<-(DW8-DW9)/DW1
W91D.DW<-(DW9-DW10)/DW1
W101D.DW<-(DW10-DW11)/DW1
W111D.DW<-(DW11-DW12)/DW1
W121D.DW<-(DW12-DW13)/DW1
W131D.DW<-(DW13-DW14)/DW1
W141D.DW<-(DW14-DW15)/DW1
W151D.DW<-(DW15-DW16)/DW1
W161D.DW<-(DW16-DW17)/DW1
W171D.DW<-(DW17-DW18)/DW1
W181D.DW<-(DW18-DW19)/DW1
W191D.DW<-(DW19-DW20)/DW1
W201D.DW<-(DW20-DW21)/DW1
W211D.DW<-(DW21-DW22)/DW1
W221D.DW<-(DW22-DW23)/DW1
W231D.DW<-(DW23-DW24)/DW1
W241D.DW<-(DW24-DW25)/DW1
W251D.DW<-(DW25-DW26)/DW1
W261D.DW<-(DW26-DW27)/DW1
W271D.DW<-(DW27-DW28)/DW1
W281D.DW<-(DW28-DW29)/DW1
W291D.DW<-(DW29-DW30)/DW1
W301D.DW<-(DW30-DW31)/DW1
n<-1:30;Dn1W.WD<-(exp(-Alpha))^(n-1)*(1-(exp(-Alpha)))
Wn1D.DW<-c(W11D.DW,W21D.DW,W31D.DW,W41D.DW,W51D.DW,W61D.DW,W71D.DW,W81D.DW,W91D.DW,W101D.DW,W111D.DW,W121D.DW,W131D.DW,W141D.DW,W151D.DW,W161D.DW,W171D.DW,W181D.DW,W191D.DW,W201D.DW,W211D.DW,W221D.DW,W231D.DW,W241D.DW,W251D.DW,W261D.DW,W271D.DW,W281D.DW,W291D.DW,W301D.DW)
cbind(Wn1D.DW,Dn1W.WD)


DRY<-Wn1D.DW;WET<-Dn1W.WD
CYCLE.n<-c(NA,
WET[1]*DRY[1],
WET[1]*DRY[2]+WET[2]*DRY[1],
WET[1]*DRY[3]+WET[2]*DRY[2]+WET[3]*DRY[1],
WET[1]*DRY[4]+WET[2]*DRY[3]+WET[3]*DRY[2]+WET[4]*DRY[1],
WET[1]*DRY[5]+WET[2]*DRY[4]+WET[3]*DRY[3]+WET[4]*DRY[2]+WET[5]*DRY[1],
WET[1]*DRY[6]+WET[2]*DRY[5]+WET[3]*DRY[4]+WET[4]*DRY[3]+WET[5]*DRY[2]+WET[6]*DRY[1],
WET[1]*DRY[7]+WET[2]*DRY[6]+WET[3]*DRY[5]+WET[4]*DRY[4]+WET[5]*DRY[3]+WET[6]*DRY[2]+WET[7]*DRY[1],
WET[1]*DRY[8]+WET[2]*DRY[7]+WET[3]*DRY[6]+WET[4]*DRY[5]+WET[5]*DRY[4]+WET[6]*DRY[3]+WET[7]*DRY[2]+WET[8]*DRY[1],
WET[1]*DRY[9]+WET[2]*DRY[8]+WET[3]*DRY[7]+WET[4]*DRY[6]+WET[5]*DRY[5]+WET[6]*DRY[4]+WET[7]*DRY[3]+WET[8]*DRY[2]+WET[9]*DRY[1],
WET[1]*DRY[10]+WET[2]*DRY[9]+WET[3]*DRY[8]+WET[4]*DRY[7]+WET[5]*DRY[6]+WET[6]*DRY[5]+WET[7]*DRY[4]+WET[8]*DRY[3]+WET[9]*DRY[2]+WET[10]*DRY[1],
WET[1]*DRY[11]+WET[2]*DRY[10]+WET[3]*DRY[9]+WET[4]*DRY[8]+WET[5]*DRY[7]+WET[6]*DRY[6]+WET[7]*DRY[5]+WET[8]*DRY[4]+WET[9]*DRY[3]+WET[10]*DRY[2]+WET[11]*DRY[1],
WET[1]*DRY[12]+WET[2]*DRY[11]+WET[3]*DRY[10]+WET[4]*DRY[9]+WET[5]*DRY[8]+WET[6]*DRY[7]+WET[7]*DRY[6]+WET[8]*DRY[5]+WET[9]*DRY[4]+WET[10]*DRY[3]+WET[11]*DRY[2]+WET[12]*DRY[1],
WET[1]*DRY[13]+WET[2]*DRY[12]+WET[3]*DRY[11]+WET[4]*DRY[10]+WET[5]*DRY[9]+WET[6]*DRY[8]+WET[7]*DRY[7]+WET[8]*DRY[6]+WET[9]*DRY[5]+WET[10]*DRY[4]+WET[11]*DRY[3]+WET[12]*DRY[2]+WET[13]*DRY[1],
WET[1]*DRY[14]+WET[2]*DRY[13]+WET[3]*DRY[12]+WET[4]*DRY[11]+WET[5]*DRY[10]+WET[6]*DRY[9]+WET[7]*DRY[8]+WET[8]*DRY[7]+WET[9]*DRY[6]+WET[10]*DRY[5]+WET[11]*DRY[4]+WET[12]*DRY[3]+WET[13]*DRY[2]+WET[14]*DRY[1],
WET[1]*DRY[15]+WET[2]*DRY[14]+WET[3]*DRY[13]+WET[4]*DRY[12]+WET[5]*DRY[11]+WET[6]*DRY[10]+WET[7]*DRY[9]+WET[8]*DRY[8]+WET[9]*DRY[7]+WET[10]*DRY[6]+WET[11]*DRY[5]+WET[12]*DRY[4]+WET[13]*DRY[3]+WET[14]*DRY[2]+WET[15]*DRY[1],
WET[1]*DRY[16]+WET[2]*DRY[15]+WET[3]*DRY[14]+WET[4]*DRY[13]+WET[5]*DRY[12]+WET[6]*DRY[11]+WET[7]*DRY[10]+WET[8]*DRY[9]+WET[9]*DRY[8]+WET[10]*DRY[7]+WET[11]*DRY[6]+WET[12]*DRY[5]+WET[13]*DRY[4]+WET[14]*DRY[3]+WET[15]*DRY[2]+WET[16]*DRY[1],
WET[1]*DRY[17]+WET[2]*DRY[16]+WET[3]*DRY[15]+WET[4]*DRY[14]+WET[5]*DRY[13]+WET[6]*DRY[12]+WET[7]*DRY[11]+WET[8]*DRY[10]+WET[9]*DRY[9]+WET[10]*DRY[8]+WET[11]*DRY[7]+WET[12]*DRY[6]+WET[13]*DRY[5]+WET[14]*DRY[4]+WET[15]*DRY[3]+WET[16]*DRY[2]+WET[17]*DRY[1],
WET[1]*DRY[18]+WET[2]*DRY[17]+WET[3]*DRY[16]+WET[4]*DRY[15]+WET[5]*DRY[14]+WET[6]*DRY[13]+WET[7]*DRY[12]+WET[8]*DRY[11]+WET[9]*DRY[10]+WET[10]*DRY[9]+WET[11]*DRY[8]+WET[12]*DRY[7]+WET[13]*DRY[6]+WET[14]*DRY[5]+WET[15]*DRY[4]+WET[16]*DRY[3]+WET[17]*DRY[2]+WET[18]*DRY[1],
WET[1]*DRY[19]+WET[2]*DRY[18]+WET[3]*DRY[17]+WET[4]*DRY[16]+WET[5]*DRY[15]+WET[6]*DRY[15]+WET[7]*DRY[13]+WET[8]*DRY[12]+WET[9]*DRY[11]+WET[10]*DRY[10]+WET[11]*DRY[9]+WET[12]*DRY[8]+WET[13]*DRY[7]+WET[14]*DRY[6]+WET[15]*DRY[5]+WET[16]*DRY[4]+WET[17]*DRY[3]+WET[18]*DRY[2]+WET[19]*DRY[1],
WET[1]*DRY[20]+WET[2]*DRY[19]+WET[3]*DRY[18]+WET[4]*DRY[17]+WET[5]*DRY[16]+WET[6]*DRY[15]+WET[7]*DRY[14]+WET[8]*DRY[13]+WET[9]*DRY[12]+WET[10]*DRY[11]+WET[11]*DRY[10]+WET[12]*DRY[9]+WET[13]*DRY[8]+WET[14]*DRY[7]+WET[15]*DRY[6]+WET[16]*DRY[5]+WET[17]*DRY[4]+WET[18]*DRY[3]+WET[19]*DRY[2]+WET[20]*DRY[1],
WET[1]*DRY[21]+WET[2]*DRY[20]+WET[3]*DRY[19]+WET[4]*DRY[18]+WET[5]*DRY[17]+WET[6]*DRY[16]+WET[7]*DRY[15]+WET[8]*DRY[14]+WET[9]*DRY[13]+WET[10]*DRY[12]+WET[11]*DRY[11]+WET[12]*DRY[10]+WET[13]*DRY[9]+WET[14]*DRY[8]+WET[15]*DRY[7]+WET[16]*DRY[6]+WET[17]*DRY[5]+WET[18]*DRY[4]+WET[19]*DRY[3]+WET[20]*DRY[2]+WET[21]*DRY[1],
WET[1]*DRY[22]+WET[2]*DRY[21]+WET[3]*DRY[20]+WET[4]*DRY[19]+WET[5]*DRY[18]+WET[6]*DRY[17]+WET[7]*DRY[16]+WET[8]*DRY[15]+WET[9]*DRY[14]+WET[10]*DRY[13]+WET[11]*DRY[12]+WET[12]*DRY[11]+WET[13]*DRY[10]+WET[14]*DRY[9]+WET[15]*DRY[8]+WET[16]*DRY[7]+WET[17]*DRY[6]+WET[18]*DRY[5]+WET[19]*DRY[4]+WET[20]*DRY[3]+WET[21]*DRY[2]+WET[22]*DRY[1],
WET[1]*DRY[23]+WET[2]*DRY[22]+WET[3]*DRY[21]+WET[4]*DRY[20]+WET[5]*DRY[19]+WET[6]*DRY[18]+WET[7]*DRY[17]+WET[8]*DRY[16]+WET[9]*DRY[15]+WET[10]*DRY[14]+WET[11]*DRY[13]+WET[12]*DRY[12]+WET[13]*DRY[11]+WET[14]*DRY[10]+WET[15]*DRY[9]+WET[16]*DRY[8]+WET[17]*DRY[7]+WET[18]*DRY[6]+WET[19]*DRY[5]+WET[20]*DRY[4]+WET[21]*DRY[3]+WET[22]*DRY[2]+WET[23]*DRY[1],
WET[1]*DRY[24]+WET[2]*DRY[23]+WET[3]*DRY[22]+WET[4]*DRY[21]+WET[5]*DRY[20]+WET[6]*DRY[19]+WET[7]*DRY[18]+WET[8]*DRY[17]+WET[9]*DRY[16]+WET[10]*DRY[15]+WET[11]*DRY[14]+WET[12]*DRY[13]+WET[13]*DRY[12]+WET[14]*DRY[11]+WET[15]*DRY[10]+WET[16]*DRY[9]+WET[17]*DRY[8]+WET[18]*DRY[7]+WET[19]*DRY[6]+WET[20]*DRY[5]+WET[21]*DRY[4]+WET[22]*DRY[3]+WET[23]*DRY[2]+WET[24]*DRY[1],
WET[1]*DRY[25]+WET[2]*DRY[24]+WET[3]*DRY[23]+WET[4]*DRY[22]+WET[5]*DRY[21]+WET[6]*DRY[20]+WET[7]*DRY[19]+WET[8]*DRY[18]+WET[9]*DRY[17]+WET[10]*DRY[16]+WET[11]*DRY[15]+WET[12]*DRY[14]+WET[13]*DRY[13]+WET[14]*DRY[12]+WET[15]*DRY[11]+WET[16]*DRY[10]+WET[17]*DRY[9]+WET[18]*DRY[8]+WET[19]*DRY[7]+WET[20]*DRY[6]+WET[21]*DRY[5]+WET[22]*DRY[4]+WET[23]*DRY[3]+WET[24]*DRY[2]+WET[25]*DRY[1],
WET[1]*DRY[26]+WET[2]*DRY[25]+WET[3]*DRY[24]+WET[4]*DRY[23]+WET[5]*DRY[22]+WET[6]*DRY[21]+WET[7]*DRY[20]+WET[8]*DRY[19]+WET[9]*DRY[18]+WET[10]*DRY[17]+WET[11]*DRY[16]+WET[12]*DRY[15]+WET[13]*DRY[14]+WET[14]*DRY[13]+WET[15]*DRY[12]+WET[16]*DRY[11]+WET[17]*DRY[10]+WET[18]*DRY[9]+WET[19]*DRY[8]+WET[20]*DRY[7]+WET[21]*DRY[6]+WET[22]*DRY[5]+WET[23]*DRY[4]+WET[24]*DRY[3]+WET[25]*DRY[2]+WET[26]*DRY[1],
WET[1]*DRY[27]+WET[2]*DRY[26]+WET[3]*DRY[25]+WET[4]*DRY[24]+WET[5]*DRY[23]+WET[6]*DRY[22]+WET[7]*DRY[21]+WET[8]*DRY[20]+WET[9]*DRY[19]+WET[10]*DRY[18]+WET[11]*DRY[17]+WET[12]*DRY[16]+WET[13]*DRY[15]+WET[14]*DRY[14]+WET[15]*DRY[13]+WET[16]*DRY[12]+WET[17]*DRY[11]+WET[18]*DRY[10]+WET[19]*DRY[9]+WET[20]*DRY[8]+WET[21]*DRY[7]+WET[22]*DRY[6]+WET[23]*DRY[5]+WET[24]*DRY[4]+WET[25]*DRY[3]+WET[26]*DRY[2]+WET[27]*DRY[1],
WET[1]*DRY[28]+WET[2]*DRY[27]+WET[3]*DRY[26]+WET[4]*DRY[25]+WET[5]*DRY[24]+WET[6]*DRY[23]+WET[7]*DRY[22]+WET[8]*DRY[21]+WET[9]*DRY[20]+WET[10]*DRY[19]+WET[11]*DRY[18]+WET[12]*DRY[17]+WET[13]*DRY[16]+WET[14]*DRY[15]+WET[15]*DRY[14]+WET[16]*DRY[13]+WET[17]*DRY[12]+WET[18]*DRY[11]+WET[19]*DRY[10]+WET[20]*DRY[9]+WET[21]*DRY[8]+WET[22]*DRY[7]+WET[23]*DRY[6]+WET[24]*DRY[5]+WET[25]*DRY[4]+WET[26]*DRY[3]+WET[27]*DRY[2]+WET[28]*DRY[1],
WET[1]*DRY[29]+WET[2]*DRY[28]+WET[3]*DRY[27]+WET[4]*DRY[26]+WET[5]*DRY[25]+WET[6]*DRY[24]+WET[7]*DRY[23]+WET[8]*DRY[22]+WET[9]*DRY[21]+WET[10]*DRY[20]+WET[11]*DRY[19]+WET[12]*DRY[18]+WET[13]*DRY[17]+WET[14]*DRY[16]+WET[15]*DRY[15]+WET[16]*DRY[14]+WET[17]*DRY[13]+WET[18]*DRY[12]+WET[19]*DRY[11]+WET[20]*DRY[10]+WET[21]*DRY[9]+WET[22]*DRY[8]+WET[23]*DRY[7]+WET[24]*DRY[6]+WET[25]*DRY[5]+WET[26]*DRY[4]+WET[27]*DRY[3]+WET[28]*DRY[2]+WET[29]*DRY[1]
)
cbind(1:30,WET,DRY,CYCLE.n)

# for july month:
kod6<-c(35, 24, 13, 9, 4, 4, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1)
kow6<-c(31, 15, 11, 7, 8, 2, 4, 2, 3, 2, 3, 1, 2, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)
kodw6<-c(10, 14, 10, 11, 11, 7, 8, 4, 1, 4, 4, 2, 1, 1, 2, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1)
kowd6<-c(11, 14, 13, 13, 8, 10, 4, 1, 6, 2, 5, 1, 3, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0)
kEC6<-CYCLE.n[-c(1,29,30)]*sum(kodw6)

Chi.Sq43<-sum((c(kodw6[1:7], sum(kodw6[8:9]), sum(kodw6[10:11]), sum(kodw6[12:27]))-c(kEC6[1:7], sum(kEC6[8:9]), sum(kEC6[10:11]), sum(kEC6[12:27])))^2/c(kEC6[1:7], sum(kEC6[8:9]), sum(kEC6[10:11]), sum(kEC6[12:27])))
pchisq(Chi.Sq43, df=9, ncp=0, lower.tail = FALSE, log.p = FALSE)
kEC6<-CYCLE.n[-c(1,29,30)]*sum(kowD6)
Chi.Sq44<-sum((c(kowd6[1:6], sum(kowd6[7:8]), sum(kowd6[9]), sum(kowd6[10:11]), sum(kowd6[12:27]))-c(kEC6[1:6], sum(kEC6[7:8]), sum(kEC6[9]), sum(kEC6[10:11]), sum(kEC6[12:27])))^2/c(kEC6[1:6], sum(kEC6[7:8]), sum(kEC6[9]), sum(kEC6[10:11]), sum(kEC6[12:27])))
pchisq(Chi.Sq44, df=9, ncp=0, lower.tail = FALSE, log.p = FALSE)