[R] concatenating two successive time series

Jędrzej Bojanowski jelus at wp.pl
Mon Mar 31 16:48:31 CEST 2008


Dear Helpers,

I am looking for methods and tools to compare and then to concatenate 
two successive time series. They are both in the same frequency and they 
describe one phenomena. There is no time gap between them. The problem 
is that the method of measurements has changed between both time series 
and they are no statistically the same. I would like to merge them to 
receive one homogeneous time series.

Here is an example of time series ts1 (1982-1999) and ts2 (2000-2006). 
Zeros are NA's.

> ts1 <- structure(c(100, 48, 31, 10, 16, 25, 56, 81, 51, 45, 73, 112,
103, 60, 30, 21, 10, 26, 60, 55, 40, 46, 78, 125, 55, 21, 22,
11, 13, 20, 38, 25, 17, 27, 49, 68, 40, 29, 44, 10, 21, 28, 58,
71, 50, 68, 88, 38, 71, 37, 51, 11, 30, 30, 52, 59, 40, 47, 91,
95, 62, 31, 44, 12, 27, 35, 36, 69, 32, 50, 81, 88, 72, 47, 28,
20, 9, 15, 34, 38, 10, 18, 56, 124, 79, 26, 29, 9, 16, 21, 71,
74, 52, 49, 80, 84, 52, 36, 37, 12, 21, 23, 70, 75, 69, 69, 81,
66, 53, 26, 18, 12, 12, 11, 45, 60, 40, 39, 64, 80, 72, 30, 23,
9, 18, 22, 53, 77, 45, 31, 54, 59, 26, 16, 14, 4, 4, 12, 20,
21, 7, 13, 28, 56, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 50, 36,
17, 18, 12, 28, 99, 90, 72, 69, 98, 78, 56, 31, 48, 11, 15, 31,
70, 79, 67, 80, 107, 76, 43, 15, 29, 13, 17, 30, 70, 86, 52,
35, 45, 93, 80, 49, 29, 16, 12, 23, 68, 80, 38, 52, 86, 94, 67,
30, 14, 10, 14, 21, 35, 49, 33, 20, 39, 77), .Tsp = c(1982, 
1999.91666666667,
12), class = "ts")

> ts2 <- structure(c(0, 0, 0, 220, 72, 45, 45, 56, 60, 80, 81, 99, 97,
92, 220, 141, 85, 40, 74, 115, 88, 97, 127, 107, 90, 86, 185,
193, 100, 43, 45, 73, 74, 95, 87, 107, 77, 63, 163, 123, 106,
35, 102, 122, 103, 96, 104, 140, 135, 87, 215, 173, 56, 38, 80,
90, 113, 72, 86, 95, 88, 93, 198, 128, 40, 28, 76, 133, 80, 89,
72, 108, 102, 86, 199, 204, 84, 60, 49, 103, 121, 98, 108, 165
), .Tsp = c(2000, 2006.91666666667, 12), class = "ts")

When we plot them together the difference is visible:

> ts12 <- ts(c(ts1,ts2),start=1982,frequency=12)
> plot(ts12)

but
> plot(decompose(ts12))
confirm that both describe one phenomena.

I have dozens of such pairs and I would like to find a method to 
recalculate one ts to the other to make one homogeneous ts.

Do you have any suggestions how to statistically concatenate those two 
ts or zoo (I tried this package) objects?   Thanks in advance for the 
answer and the precious time you might waste.

Regards,

Jedrzej

----------------------------------------------------
PS. Włącz radio w Twoim komputerze!
Muzyka idealna do pracy i do zabawy! - Kliknij:
http://klik.wp.pl/?adr=http%3A%2F%2Fcorto.www.wp.pl%2Fas%2Fradiokosci.html&sid=291



More information about the R-help mailing list