[R] [FORGED] Fitting arima Models with Exogenous Variables

Rolf Turner r.turner at auckland.ac.nz
Thu Jan 26 20:49:41 CET 2017


I am re-sending this since I have been told that the attachments that I 
made did not get through.  So I am trying again with *.dput attachments.

You will need to read them in using dget().

cheers,

Rolf Turner

-- 
Technical Editor ANZJS
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276

On 24/01/17 10:35, Rolf Turner wrote:
>
> This should have been sent to the R-help mailing list, not to me
> personally.  I am not an expert on this sort of time series modelling
> and cannot thereby provide any useful advice.  My reply to you was of a
> "generic" nature --- when making an enquiry, provide a reproducible
> example!!!
>
> I am cc-ing this email to the R-help list, since someone on that list
> *may* be able to answer your question.  I have (re-) attached the data
> sets that you sent to me.
>
> cheers,
>
> Rolf Turner
>
> On 24/01/17 04:36, Paul Bernal wrote:
>> Hello Rolf,
>>
>> Thank you for your kind reply. I am attaching two datasets, one with the
>> historical data that I used to train the model, and the other one with
>> the exogenous variables.
>>
>> The R code that I used is as follows:
>>
>>> library(forecast)
>>> library(tseries)
>>> library(TSA)
>>> library(stats)
>>> library(stats4)
>>> TrainingDat<-read.csv("Training Data.csv")
>>>
>>> ExogVars<-read.csv("ExogenousVariables5.csv")
>>> #The file ExogVars contains 5 columns, one column for each regressor
>>> Model1<-auto.arima(TrainingDat[,5], xreg=ExogVars)
>>> #In Model1 I was able to incorporate xreg without any trouble
>>> #The problem comes when trying to incorporate newxreg
>>> Model2<-auto.arima(ExoVars[1:5])
>> Error in as.ts(x) : object 'ExoVars' not found
>>>
>>> Model2<-auto.arima(ExogVars[1:5])
>>
>> Error in auto.arima(ExogVars[1:5]) : No suitable ARIMA model found
>>>
>>> Model2<-auto.arima(ExogVars[,1])
>>>
>>> NewXReg<-forecast(Model2, h=12)
>>>
>>> Forec<-forecast(Model1, newxreg=NewXReg)
>> Error in forecast.Arima(Model1, newxreg = NewXReg) :
>>   No regressors provided
>> In addition: Warning message:
>> In forecast.Arima(Model1, newxreg = NewXReg) :
>>   The non-existent newxreg arguments will be ignored.
>>>
>>> Forec<-forecast(Model1, newxreg=NewXReg$mean)
>> Error in forecast.Arima(Model1, newxreg = NewXReg$mean) :
>>   No regressors provided
>> In addition: Warning message:
>> In forecast.Arima(Model1, newxreg = NewXReg$mean) :
>>   The non-existent newxreg arguments will be ignored.
>>
>> I would like to generate the forecasts for all 4 variables included in
>> the Training set, along with all 5 regressors, but it seems like I can
>> only chose one training variable at a time, and one regressor at a time.
>>
>> Please let me know if you can work this out,




-------------- next part --------------
structure(list(date = c("1994/Oct", "1994/Nov", "1994/Dec", "1995/Jan", 
"1995/Feb", "1995/Mar", "1995/Apr", "1995/May", "1995/Jun", "1995/Jul", 
"1995/Aug", "1995/Sep", "1995/Oct", "1995/Nov", "1995/Dec", "1996/Jan", 
"1996/Feb", "1996/Mar", "1996/Apr", "1996/May", "1996/Jun", "1996/Jul", 
"1996/Aug", "1996/Sep", "1996/Oct", "1996/Nov", "1996/Dec", "1997/Jan", 
"1997/Feb", "1997/Mar", "1997/Apr", "1997/May", "1997/Jun", "1997/Jul", 
"1997/Aug", "1997/Sep", "1997/Oct", "1997/Nov", "1997/Dec", "1998/Jan", 
"1998/Feb", "1998/Mar", "1998/Apr", "1998/May", "1998/Jun", "1998/Jul", 
"1998/Aug", "1998/Sep", "1998/Oct", "1998/Nov", "1998/Dec", "1999/Jan", 
"1999/Feb", "1999/Mar", "1999/Apr", "1999/May", "1999/Jun", "1999/Jul", 
"1999/Aug", "1999/Sep", "1999/Oct", "1999/Nov", "1999/Dec", "2000/Jan", 
"2000/Feb", "2000/Mar", "2000/Apr", "2000/May", "2000/Jun", "2000/Jul", 
"2000/Aug", "2000/Sep", "2000/Oct", "2000/Nov", "2000/Dec", "2001/Jan", 
"2001/Feb", "2001/Mar", "2001/Apr", "2001/May", "2001/Jun", "2001/Jul", 
"2001/Aug", "2001/Sep", "2001/Oct", "2001/Nov", "2001/Dec", "2002/Jan", 
"2002/Feb", "2002/Mar", "2002/Apr", "2002/May", "2002/Jun", "2002/Jul", 
"2002/Aug", "2002/Sep", "2002/Oct", "2002/Nov", "2002/Dec", "2003/Jan", 
"2003/Feb", "2003/Mar", "2003/Apr", "2003/May", "2003/Jun", "2003/Jul", 
"2003/Aug", "2003/Sep", "2003/Oct", "2003/Nov", "2003/Dec", "2004/Jan", 
"2004/Feb", "2004/Mar", "2004/Apr", "2004/May", "2004/Jun", "2004/Jul", 
"2004/Aug", "2004/Sep", "2004/Oct", "2004/Nov", "2004/Dec", "2005/Jan", 
"2005/Feb", "2005/Mar", "2005/Apr", "2005/May", "2005/Jun", "2005/Jul", 
"2005/Aug", "2005/Sep", "2005/Oct", "2005/Nov", "2005/Dec", "2006/Jan", 
"2006/Feb", "2006/Mar", "2006/Apr", "2006/May", "2006/Jun", "2006/Jul", 
"2006/Aug", "2006/Sep", "2006/Oct", "2006/Nov", "2006/Dec", "2007/Jan", 
"2007/Feb", "2007/Mar", "2007/Apr", "2007/May", "2007/Jun", "2007/Jul", 
"2007/Aug", "2007/Sep", "2007/Oct", "2007/Nov", "2007/Dec", "2008/Jan", 
"2008/Feb", "2008/Mar", "2008/Apr", "2008/May", "2008/Jun", "2008/Jul", 
"2008/Aug", "2008/Sep", "2008/Oct", "2008/Nov", "2008/Dec", "2009/Jan", 
"2009/Feb", "2009/Mar", "2009/Apr", "2009/May", "2009/Jun", "2009/Jul", 
"2009/Aug", "2009/Sep", "2009/Oct", "2009/Nov", "2009/Dec", "2010/Jan", 
"2010/Feb", "2010/Mar", "2010/Apr", "2010/May", "2010/Jun", "2010/Jul", 
"2010/Aug", "2010/Sep", "2010/Oct", "2010/Nov", "2010/Dec", "2011/Jan", 
"2011/Feb", "2011/Mar", "2011/Apr", "2011/May", "2011/Jun", "2011/Jul", 
"2011/Aug", "2011/Sep", "2011/Oct", "2011/Nov", "2011/Dec", "2012/Jan", 
"2012/Feb", "2012/Mar", "2012/Apr", "2012/May", "2012/Jun", "2012/Jul", 
"2012/Aug", "2012/Sep", "2012/Oct", "2012/Nov", "2012/Dec", "2013/Jan", 
"2013/Feb", "2013/Mar", "2013/Apr", "2013/May", "2013/Jun", "2013/Jul", 
"2013/Aug", "2013/Sep", "2013/Oct", "2013/Nov", "2013/Dec", "2014/Jan", 
"2014/Feb", "2014/Mar", "2014/Apr", "2014/May", "2014/Jun", "2014/Jul", 
"2014/Aug", "2014/Sep", "2014/Oct", "2014/Nov", "2014/Dec", "2015/Jan", 
"2015/Feb", "2015/Mar", "2015/Apr", "2015/May", "2015/Jun", "2015/Jul", 
"2015/Aug", "2015/Sep", "2015/Oct", "2015/Nov", "2015/Dec"), 
    Transits.Class.A = c(22L, 18L, 19L, 18L, 18L, 19L, 21L, 22L, 
    17L, 17L, 17L, 16L, 16L, 17L, 18L, 15L, 17L, 19L, 19L, 19L, 
    13L, 11L, 11L, 6L, 7L, 6L, 11L, 7L, 9L, 8L, 7L, 11L, 14L, 
    17L, 16L, 12L, 14L, 12L, 13L, 11L, 11L, 15L, 11L, 11L, 12L, 
    17L, 13L, 11L, 11L, 12L, 11L, 11L, 10L, 4L, 5L, 4L, 5L, 7L, 
    5L, 2L, 5L, 6L, 9L, 8L, 11L, 4L, 2L, 2L, 2L, 3L, 3L, 6L, 
    3L, 4L, 1L, 7L, 6L, 7L, 10L, 4L, 4L, 8L, 12L, 12L, 10L, 10L, 
    7L, 12L, 7L, 7L, 8L, 12L, 12L, 13L, 13L, 9L, 10L, 10L, 6L, 
    7L, 11L, 17L, 16L, 14L, 13L, 14L, 13L, 16L, 14L, 12L, 13L, 
    12L, 11L, 12L, 14L, 14L, 15L, 14L, 12L, 12L, 15L, 11L, 15L, 
    13L, 11L, 14L, 14L, 14L, 15L, 12L, 14L, 10L, 5L, 6L, 4L, 
    3L, 3L, 5L, 9L, 6L, 5L, 5L, 2L, 3L, 5L, 4L, 3L, 6L, 3L, 7L, 
    4L, 3L, 2L, 1L, 3L, 5L, 3L, 3L, 6L, 7L, 10L, 10L, 13L, 10L, 
    6L, 7L, 4L, 5L, 6L, 4L, 5L, 3L, 6L, 9L, 3L, 3L, 6L, 7L, 8L, 
    6L, 6L, 6L, 7L, 8L, 4L, 8L, 8L, 6L, 5L, 5L, 7L, 4L, 5L, 6L, 
    5L, 8L, 9L, 12L, 10L, 10L, 7L, 4L, 4L, 4L, 7L, 6L, 7L, 6L, 
    7L, 4L, 5L, 7L, 2L, 4L, 4L, 4L, 3L, 2L, 3L, 3L, 1L, 5L, 2L, 
    3L, 3L, 2L, 2L, 3L, 1L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 4L, 
    2L, 2L, 2L, 2L, 2L, 4L, 5L, 5L, 5L, 8L, 7L, 5L, 3L, 3L, 2L, 
    2L, 5L), Transits.Class.B = c(15L, 15L, 14L, 14L, 14L, 15L, 
    16L, 15L, 18L, 21L, 18L, 17L, 21L, 18L, 19L, 19L, 18L, 21L, 
    22L, 25L, 23L, 28L, 25L, 25L, 25L, 21L, 28L, 21L, 22L, 23L, 
    20L, 24L, 22L, 26L, 24L, 24L, 28L, 31L, 34L, 36L, 33L, 29L, 
    34L, 33L, 33L, 31L, 33L, 29L, 32L, 26L, 24L, 30L, 27L, 31L, 
    23L, 23L, 21L, 21L, 18L, 21L, 19L, 19L, 16L, 14L, 22L, 24L, 
    28L, 26L, 29L, 23L, 28L, 21L, 17L, 17L, 14L, 12L, 14L, 10L, 
    12L, 13L, 14L, 15L, 21L, 10L, 15L, 12L, 15L, 16L, 17L, 14L, 
    14L, 12L, 7L, 8L, 8L, 11L, 9L, 12L, 11L, 6L, 8L, 11L, 7L, 
    8L, 14L, 9L, 11L, 8L, 8L, 6L, 3L, 4L, 4L, 5L, 4L, 3L, 3L, 
    5L, 8L, 6L, 8L, 6L, 7L, 14L, 14L, 15L, 13L, 16L, 18L, 20L, 
    21L, 18L, 19L, 24L, 27L, 26L, 18L, 28L, 21L, 22L, 25L, 29L, 
    26L, 22L, 22L, 14L, 18L, 18L, 13L, 15L, 16L, 19L, 9L, 14L, 
    11L, 9L, 14L, 16L, 11L, 13L, 11L, 16L, 16L, 20L, 17L, 19L, 
    18L, 24L, 19L, 19L, 15L, 21L, 21L, 24L, 28L, 14L, 20L, 15L, 
    24L, 18L, 24L, 24L, 25L, 21L, 24L, 27L, 28L, 27L, 26L, 29L, 
    25L, 24L, 25L, 28L, 23L, 30L, 23L, 29L, 26L, 31L, 33L, 40L, 
    36L, 39L, 36L, 34L, 32L, 32L, 30L, 27L, 32L, 30L, 31L, 26L, 
    31L, 29L, 26L, 28L, 24L, 22L, 11L, 15L, 17L, 12L, 13L, 14L, 
    9L, 12L, 11L, 17L, 9L, 8L, 7L, 11L, 8L, 15L, 8L, 10L, 13L, 
    12L, 13L, 12L, 12L, 19L, 21L, 18L, 20L, 21L, 20L, 22L, 19L, 
    21L, 20L, 21L, 18L), Transits.Class.C = c(15L, 12L, 15L, 
    14L, 13L, 17L, 15L, 12L, 18L, 15L, 11L, 16L, 13L, 14L, 16L, 
    14L, 12L, 15L, 13L, 12L, 14L, 13L, 12L, 14L, 11L, 15L, 13L, 
    14L, 14L, 14L, 15L, 13L, 17L, 14L, 15L, 13L, 18L, 16L, 19L, 
    19L, 14L, 18L, 20L, 13L, 19L, 17L, 23L, 21L, 23L, 19L, 25L, 
    24L, 21L, 23L, 24L, 28L, 27L, 29L, 28L, 23L, 27L, 26L, 26L, 
    24L, 20L, 23L, 21L, 28L, 24L, 29L, 28L, 32L, 26L, 36L, 25L, 
    28L, 24L, 26L, 28L, 29L, 24L, 24L, 28L, 21L, 23L, 23L, 24L, 
    21L, 27L, 23L, 22L, 27L, 26L, 24L, 28L, 24L, 27L, 27L, 34L, 
    26L, 27L, 23L, 27L, 23L, 30L, 35L, 34L, 31L, 33L, 31L, 34L, 
    37L, 29L, 37L, 29L, 32L, 26L, 33L, 28L, 27L, 29L, 26L, 29L, 
    30L, 29L, 34L, 38L, 31L, 37L, 33L, 34L, 32L, 33L, 27L, 28L, 
    30L, 32L, 35L, 32L, 37L, 36L, 41L, 43L, 39L, 41L, 40L, 49L, 
    56L, 55L, 64L, 73L, 68L, 64L, 58L, 69L, 51L, 64L, 57L, 67L, 
    70L, 65L, 65L, 60L, 66L, 56L, 48L, 58L, 56L, 64L, 58L, 53L, 
    55L, 56L, 64L, 53L, 69L, 60L, 59L, 51L, 54L, 56L, 49L, 48L, 
    50L, 45L, 50L, 56L, 54L, 47L, 52L, 46L, 46L, 45L, 48L, 41L, 
    43L, 38L, 45L, 46L, 41L, 40L, 37L, 30L, 31L, 42L, 34L, 37L, 
    39L, 44L, 40L, 44L, 43L, 26L, 32L, 24L, 25L, 30L, 27L, 28L, 
    26L, 25L, 27L, 32L, 29L, 23L, 25L, 25L, 22L, 22L, 19L, 21L, 
    18L, 17L, 17L, 15L, 19L, 15L, 18L, 13L, 14L, 14L, 14L, 15L, 
    14L, 12L, 18L, 15L, 16L, 15L, 18L, 15L, 15L, 14L, 13L, 19L
    ), Transits.Class.D = c(58L, 55L, 53L, 60L, 54L, 60L, 62L, 
    63L, 60L, 67L, 59L, 62L, 64L, 62L, 65L, 67L, 61L, 65L, 64L, 
    64L, 68L, 65L, 73L, 66L, 71L, 65L, 69L, 67L, 61L, 71L, 63L, 
    71L, 67L, 66L, 71L, 68L, 67L, 69L, 66L, 72L, 62L, 70L, 69L, 
    79L, 80L, 87L, 86L, 85L, 89L, 81L, 81L, 81L, 65L, 75L, 71L, 
    82L, 80L, 86L, 81L, 80L, 83L, 78L, 87L, 88L, 79L, 87L, 79L, 
    96L, 95L, 98L, 102L, 93L, 108L, 97L, 98L, 102L, 93L, 101L, 
    96L, 106L, 98L, 104L, 107L, 105L, 109L, 107L, 107L, 117L, 
    106L, 112L, 122L, 117L, 134L, 143L, 141L, 141L, 138L, 136L, 
    151L, 152L, 129L, 146L, 137L, 150L, 147L, 160L, 161L, 157L, 
    165L, 162L, 171L, 151L, 152L, 154L, 158L, 168L, 167L, 160L, 
    176L, 161L, 172L, 163L, 168L, 175L, 166L, 179L, 176L, 200L, 
    182L, 195L, 200L, 191L, 207L, 204L, 207L, 215L, 183L, 202L, 
    201L, 211L, 210L, 225L, 232L, 237L, 250L, 236L, 238L, 228L, 
    198L, 217L, 209L, 228L, 207L, 238L, 222L, 239L, 227L, 222L, 
    211L, 204L, 195L, 189L, 198L, 217L, 197L, 222L, 227L, 222L, 
    219L, 206L, 213L, 211L, 178L, 187L, 195L, 200L, 198L, 186L, 
    188L, 181L, 180L, 164L, 162L, 158L, 135L, 164L, 145L, 180L, 
    179L, 190L, 199L, 195L, 207L, 182L, 201L, 189L, 181L, 186L, 
    189L, 205L, 186L, 202L, 196L, 196L, 197L, 197L, 188L, 191L, 
    172L, 198L, 202L, 208L, 216L, 239L, 233L, 225L, 228L, 204L, 
    210L, 218L, 197L, 204L, 210L, 218L, 210L, 215L, 219L, 216L, 
    219L, 194L, 201L, 206L, 190L, 214L, 207L, 225L, 211L, 220L, 
    214L, 224L, 222L, 202L, 211L, 217L, 188L, 226L, 205L, 218L, 
    219L, 241L, 244L, 226L, 238L, 214L, 218L)), .Names = c("date", 
"Transits.Class.A", "Transits.Class.B", "Transits.Class.C", "Transits.Class.D"
), class = "data.frame", row.names = c(NA, -255L))
-------------- next part --------------
structure(list(MGOFujairah = c(183.75, 185, 185, 185, 185, 185, 
185.25, 186.25, 188.5, 188.63, 184.13, 183.3, 182.5, 182.5, 193, 
205.25, 206.63, 211.9, 211.25, 214.5, 209.5, 205, 205.5, 220, 
231.5, 236.25, 247.5, 230.6, 220, 220, 220, 216.2, 210.5, 207.5, 
212.62, 212.74, 215, 214.75, 201.25, 175.3, 172.25, 172.63, 177.75, 
172.7, 165.25, 152, 146.75, 153.38, 147.8, 142.25, 132.5, 133.75, 
119.5, 129.13, 140.4, 142.75, 144, 153.2, 178.75, 190.75, 201, 
204.25, 201.9, 209.38, 226.25, 269.3, 254.38, 255, 253, 254.5, 
281.25, 323.6, 327.25, 305.63, 279.4, 274, 256, 252, 260, 258.75, 
254.7, 252.75, 251.7, 258.75, 242.75, 228.7, 215.38, 211.25, 
207, 206.3, 219.63, 222.1, 221.88, 223.75, 227, 245.25, 253.75, 
249, 244.38, 260.5, 302.5, 319.38, 298.75, 275.5, 255.63, 245.63, 
249, 249.63, 254.2, 259.38, 270.63, 307.1, 317.5, 311.38, 314.2, 
338.13, 335.63, 349.4, 390.63, 398.13, 434.5, 441.25, 431, 424.38, 
425, 470, 511.5, 505, 510.63, 517, 523.75, 552.8, 557.5, 546.25, 
534, 533.75, 541.25, 548.5, 624.38, 646.88, 646.9, 680.88, 692.63, 
672.5, 634.38, 613.75, 590, 581.88, 586, 592.5, 615, 618.75, 
633.9, 672.5, 698, 707.5, 717.5, 807.5, 812.94, 833.13, 843, 
923.75, 1010, 1204.5, 1274.38, 1353.25, 1296, 1105, 928.5, 778.13, 
680.63, 653.5, 543.75, 477.5, 503.13, 545.3, 614.38, 610, 644.38, 
626.88, 634, 646.88, 645, 658, 660.63, 668.13, 697.9, 720, 720, 
727, 731.88, 732, 739, 752.5, 791.4, 870, 959.63, 990, 1028.2, 
1036.25, 1035, 1065.6, 1078.75, 1061, 1047.5, 1050, 1043.1, 1060.5, 
1055.88, 1058.8, 1043.38, 1055, 1024, 1017.25, 1028, 1035, 1028.75, 
1026.5, 1016.88, 999.38, 1007.5, 1016.5, 1016.25, 1001.5, 995, 
991.88, 1002.5, 999.38, 1000.63, 1000.5, 991.25, 986.1, 981.25, 
981.25, 979.25, 978.5, 975.13, 972.13, 979.8, 991.25, 939.5, 
887.25, 873.25, 779.3, 775.75, 746.88, 730, 732.5, 744.38, 715, 
680, 634.38, 599.8, 612.75, 556.75), MGORotterdam = c(150.5, 
150, 141.8, 141, 140, 144, 156, 155.25, 153, 145, 148, 154.3, 
146.75, 149.5, 165.2, 168.38, 176.13, 189.3, 194.63, 169.2, 163.13, 
174.88, 182.6, 218.75, 234.6, 218.75, 217.5, 219.5, 186.25, 167.13, 
159.88, 169.2, 164.38, 162.88, 186.54, 163.5, 178, 176.75, 162, 
144.3, 139.5, 129.5, 132.63, 130.2, 118.25, 114.8, 104.25, 114.5, 
117.5, 106.63, 98.9, 100.38, 97.5, 110.63, 125.9, 122.88, 120.25, 
145.3, 165.13, 178.5, 183.8, 194.5, 209, 217, 225.38, 230.7, 
216.25, 224.5, 236.5, 244.25, 263.75, 315.9, 298, 297.5, 265.4, 
255, 230.5, 223, 227.5, 225, 227.7, 219.38, 223.3, 228.38, 203.13, 
182.3, 168.88, 170.13, 166.38, 191.5, 207.38, 207.3, 199.75, 
210.25, 220.7, 239.25, 241.88, 215.2, 248.88, 270.3, 309.5, 304.13, 
237.5, 226.9, 231.88, 234.63, 245, 233, 245, 262.38, 272.88, 
281.9, 263.38, 297.13, 299.9, 326.88, 321.63, 347.6, 381.75, 
391.13, 475, 445.75, 398, 391.63, 413.75, 493.13, 494.7, 448.75, 
508.88, 532, 567.5, 620.5, 593.38, 515.13, 502, 535.63, 521.88, 
549.5, 610.88, 611.25, 619, 630, 639.38, 575, 530.13, 524.75, 
535.5, 484.38, 506.88, 537.6, 574.25, 579.38, 602.1, 637.13, 
609.7, 673.63, 696.88, 814.5, 815, 810, 838.5, 961.88, 1023.13, 
1178.8, 1225, 1234.5, 1040.5, 927.5, 749.9, 577.88, 453.13, 450.3, 
404.5, 407.25, 452.75, 477.8, 575, 536.9, 613.75, 574.88, 609.5, 
635, 608.75, 640.44, 620, 668.75, 719.5, 665.38, 660, 648.5, 
661.25, 666.5, 728.8, 741.25, 772.4, 813.88, 898.75, 973.75, 
1016.5, 951.25, 955, 964.9, 946.13, 955, 936, 981.63, 941.2, 
968, 993.63, 1023, 1007.75, 961.25, 848.3, 891, 959.6, 975.63, 
980.63, 929.5, 925.63, 940, 990, 905.5, 851.25, 838.4, 861.25, 
897.5, 920.7, 926.75, 912.84, 893.7, 907.5, 891.5, 885, 866.38, 
874.25, 879.2, 880.75, 863.75, 847.5, 815.75, 743.5, 693.5, 558.25, 
479.5, 551.13, 537.5, 530.5, 580.3, 559.75, 494.1, 433, 442.75, 
424.5, 401.88, 319.88), X380CSTFujairah = c(85.5, 92.25, 88.4, 
96.25, 98.38, 102.3, 104.38, 103.63, 94.1, 83.13, 84, 83.9, 87.75, 
91.88, 105.8, 109.88, 99.38, 103.2, 107.75, 96.3, 88.13, 92.88, 
97.5, 111.88, 115.2, 116.13, 115, 105.8, 94, 96.5, 95.63, 96.4, 
97.38, 96.63, 117.55, 107.75, 110.1, 104.75, 89.75, 67.4, 59.38, 
70.63, 77.75, 75.9, 62.25, 62.4, 59.38, 76.88, 82.7, 73.5, 62.4, 
62, 55.13, 60.25, 74.7, 83.13, 83.63, 96, 115.13, 127.88, 138, 
143.25, 142.9, 135.88, 139.63, 170.4, 153.88, 165.13, 166.6, 
147.63, 147.25, 164.8, 177.13, 163.63, 131.5, 126.5, 126.63, 
129.9, 135.88, 146.38, 132, 127.25, 137.3, 149.13, 126.25, 106.1, 
112.5, 109.63, 119.38, 129.9, 144.25, 151.1, 149.13, 153, 160.7, 
172.25, 162.13, 144, 155.38, 176.8, 200.13, 176.63, 148.88, 159.3, 
162.88, 173.88, 166.8, 156.75, 159.4, 158.25, 156.25, 162.3, 
169.63, 163.13, 172.9, 187.75, 179.5, 179.6, 185.13, 178, 193.1, 
177.75, 171.8, 184.63, 196.88, 210.63, 247.6, 256.88, 257.75, 
264.9, 275, 312.6, 301.38, 293.25, 277.6, 297, 319.25, 330.4, 
345.75, 342.88, 323.8, 337, 326.5, 282.9, 282.5, 272, 270.6, 
264.25, 302, 314.8, 343.63, 344, 349.5, 379.25, 381.3, 404.38, 
428.75, 500.6, 472.5, 484.13, 476.6, 507.38, 531.38, 582.35, 
628.5, 715.25, 703.6, 595, 418.7, 242.63, 226.75, 254.6, 258.75, 
243.88, 284.13, 339.8, 406.75, 406.8, 454.38, 456.75, 441.5, 
466.75, 459.25, 481.8, 465.25, 472.25, 481.2, 460.38, 443.25, 
442.2, 453.5, 445.5, 473.1, 492.63, 507, 543.5, 631.5, 638.5, 
672.7, 651, 659, 674.5, 671.5, 661.9, 659.75, 690.5, 680.8, 725.75, 
731.88, 745.2, 729.13, 683.63, 604.4, 620.38, 668.2, 665.25, 
638.5, 604.8, 605, 619.5, 653.5, 637.4, 617.25, 608.8, 614.5, 
593.13, 605.8, 608.25, 616.75, 618.9, 607.75, 617.4, 605.5, 601.75, 
596, 596.5, 608.63, 602.75, 602.5, 590.63, 505.9, 454.88, 358.75, 
295.3, 354.75, 335.63, 337.88, 379.8, 348.63, 311.8, 250.25, 
238.88, 245.7, 226.25, 178.75), X380CSTRotterdam = c(89.13, 102.88, 
96.1, 108.5, 103.38, 106.2, 105.19, 107.88, 93.6, 81.38, 83.25, 
87.5, 85.88, 85.5, 102.8, 103.63, 95.5, 106.6, 115.38, 99.6, 
88.63, 88.13, 97.4, 114.5, 122.5, 116.5, 120.38, 110.2, 96.38, 
89.75, 86.13, 84.6, 86.63, 89.38, 108.69, 97.5, 104.6, 107.5, 
88.69, 73.5, 69, 67.25, 79.75, 73.2, 67.38, 66.7, 62.63, 66.25, 
69.5, 60.5, 55.8, 64.5, 57.75, 61.88, 72.8, 71.88, 80.75, 95.7, 
112.88, 117.75, 127.7, 130, 127.3, 127.5, 133.88, 147.3, 125.5, 
126.5, 144.3, 132.88, 132.5, 154.5, 159.88, 149.75, 126.7, 117.9, 
121.5, 119.8, 116.75, 122, 121.6, 118, 125.4, 130.13, 110.63, 
101.9, 103.75, 104.5, 103.25, 120.5, 136.13, 140.8, 135, 142.75, 
145.9, 161.38, 154, 125.7, 134.38, 173.5, 175.13, 148.25, 126.38, 
137.9, 147.5, 170.5, 159.7, 147.75, 151.9, 154.75, 141, 141.4, 
139, 144.5, 152.6, 170.63, 159.5, 162.9, 167, 161.75, 173.9, 
146.5, 143.5, 157.5, 171.25, 203, 231.9, 230.13, 232.25, 248.9, 
261.75, 288.4, 270.5, 256.38, 255.8, 282.88, 294.75, 300.3, 320.13, 
324.63, 301.4, 316.5, 311, 280.5, 266.13, 262.38, 255.9, 229.13, 
251.5, 272.7, 312, 325.88, 325.9, 359.88, 353.7, 374, 412.5, 
476.1, 447.5, 447.75, 436.9, 477.38, 494.88, 542.9, 593.76, 679.5, 
635.6, 544.13, 398, 217.63, 194.5, 225.7, 239.13, 244.5, 276.25, 
327.5, 384.25, 382.8, 429.38, 412.38, 423.1, 462, 438.75, 457.2, 
445.5, 453.25, 467, 437.38, 424.38, 423.7, 439.63, 431.5, 458.6, 
476, 488.6, 514.75, 575.25, 606.38, 641.6, 621.75, 631.25, 647.5, 
634.38, 640.5, 633, 645.13, 623.8, 682.38, 694.63, 712.6, 697, 
651.5, 572.4, 597.13, 640.5, 640, 616.75, 588, 582.75, 608.75, 
634.5, 605.6, 584, 580.2, 580, 596.75, 601.6, 597.5, 588.5, 575.9, 
584.25, 570.4, 580, 574.5, 578.75, 576.4, 590.25, 575, 563.4, 
547.5, 481.6, 420, 327.88, 250.1, 304.63, 301.88, 307.25, 343.7, 
328.5, 292.7, 233.75, 226.63, 222.5, 201, 157.13), X380CSTHouston = c(91, 
94.25, 89.2, 90, 96.75, 97.05, 100.19, 106.75, 101.5, 85.5, 83.5, 
84.7, 86.5, 86.75, 99.2, 99.75, 99.13, 102.7, 107.38, 101.4, 
96.75, 93.5, 96.1, 111.38, 130.3, 119.75, 112, 106.3, 90.5, 90.75, 
93.38, 96.4, 94.88, 95.38, 107.55, 99.63, 109, 108.63, 86.88, 
74.4, 72.75, 59, 78.13, 72.1, 67.88, 68.4, 62.63, 65.13, 68.9, 
63.38, 53.3, 56.13, 52, 63, 82.2, 81, 81.88, 92.2, 114.5, 119.13, 
125.5, 126.5, 125.8, 128, 135.63, 141.1, 123.75, 132, 145.6, 
133.25, 128.88, 148.5, 152.88, 142, 120.1, 116.4, 129.5, 120.9, 
104.13, 114.38, 111.1, 113.5, 117, 128.75, 105.75, 93.8, 98.38, 
99.5, 97.75, 117.9, 139.25, 142, 139, 141.75, 148.6, 158.88, 
153, 128.4, 142.5, 195, 186.25, 160.25, 134.88, 139.9, 152, 167.75, 
172.1, 149.13, 157.7, 158, 150, 148.9, 153.5, 151.13, 160.1, 
181.13, 167.5, 164.5, 170.38, 171.13, 213.3, 152.88, 173.4, 182.63, 
177, 200.25, 247, 255.13, 253.63, 253.5, 264.63, 305.3, 294.63, 
269.13, 276.4, 300, 310.25, 309.4, 330.68, 334.13, 309.4, 326.13, 
331.38, 284.9, 269, 263.38, 267.5, 247.75, 267.25, 272.7, 305.75, 
334.75, 348.3, 359.88, 366.4, 376.38, 403.25, 482, 457, 458.5, 
454.6, 486.13, 498.13, 568.4, 629.5, 721.25, 661.5, 592, 420.5, 
242, 229.5, 256.9, 265.13, 247.38, 275.88, 328.4, 389, 388, 422.38, 
416.25, 425.2, 466.5, 447.75, 458.5, 450.13, 449.75, 464, 439.25, 
424.75, 427.1, 432.88, 433.38, 456.1, 468.25, 487.9, 511.38, 
572.75, 629.75, 656.3, 633.88, 651.75, 653.4, 631.25, 635.1, 
647.25, 657.25, 628.1, 669, 700.25, 722.1, 709.38, 661.5, 577.9, 
587, 644.3, 647.5, 625, 603.9, 611.13, 624.25, 645.63, 612, 591, 
584.2, 581.63, 584.75, 597.9, 620.75, 600.38, 589.5, 600, 591.7, 
598, 592.25, 591.75, 598.5, 601, 584.13, 574.1, 558.13, 492.7, 
424.25, 340, 277.4, 326.88, 312.38, 315.63, 345.82, 338.5, 296.1, 
235, 224.88, 221.1, 205.63, 152.75), ContainerTCRIndex = c(101L, 
103L, 102L, 104L, 105L, 103L, 106L, 106L, 110L, 109L, 110L, 111L, 
114L, 107L, 109L, 109L, 110L, 110L, 110L, 110L, 110L, 109L, 108L, 
104L, 103L, 101L, 100L, 97L, 98L, 96L, 91L, 89L, 89L, 89L, 86L, 
85L, 84L, 83L, 84L, 83L, 81L, 80L, 80L, 78L, 75L, 72L, 69L, 66L, 
65L, 65L, 63L, 58L, 58L, 59L, 61L, 61L, 64L, 67L, 72L, 74L, 74L, 
71L, 68L, 71L, 78L, 81L, 86L, 90L, 91L, 91L, 90L, 90L, 90L, 85L, 
83L, 81L, 82L, 83L, 82L, 79L, 76L, 72L, 66L, 59L, 55L, 49L, 47L, 
48L, 49L, 51L, 56L, 57L, 57L, 59L, 64L, 62L, 63L, 62L, 59L, 62L, 
71L, 75L, 79L, 85L, 88L, 93L, 95L, 97L, 95L, 95L, 94L, 104L, 
113L, 120L, 124L, 126L, 128L, 131L, 136L, 142L, 146L, 154L, 163L, 
169L, 171L, 172L, 172L, 170L, 164L, 159L, 150L, 142L, 124L, 114L, 
116L, 111L, 107L, 107L, 111L, 113L, 112L, 112L, 108L, 105L, 99L, 
94L, 91L, 96L, 100L, 102L, 103L, 105L, 107L, 109L, 112L, 116L, 
116L, 114L, 114L, 113L, 113L, 113L, 111L, 108L, 101L, 97L, 92L, 
87L, 65L, 58L, 47L, 45L, 40L, 37L, 35L, 35L, 34L, 34L, 34L, 33L, 
33L, 32L, 32L, 32L, 34L, 36L, 41L, 47L, 58L, 61L, 64L, 64L, 61L, 
57L, 59L, 66L, 71L, 76L, 75L, 75L, 73L, 65L, 61L, 55L, 48L, 46L, 
42L, 40L, 41L, 42L, 43L, 45L, 45L, 44L, 44L, 44L, 43L, 42L, 42L, 
42L, 43L, 43L, 45L, 46L, 47L, 47L, 48L, 47L, 48L, 47L, 47L, 47L, 
47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 48L, 50L, 
52L, 59L, 63L, 62L, 60L, 55L, 51L, 46L, 45L, 43L), USAInterestRates = c(0.08, 
0.08, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 
0.09, 0.09, 0.09, 0.09, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 
0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.09, 0.09, 0.09, 0.09, 
0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 0.09, 
0.09, 0.09, 0.09, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 
0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.09, 0.09, 0.09, 0.09, 0.09, 
0.09, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.09, 0.09, 0.08, 0.08, 
0.07, 0.07, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.05, 0.05, 0.05, 
0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.04, 0.04, 0.04, 0.04, 
0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 
0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.04, 0.05, 0.05, 0.05, 0.05, 
0.05, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.07, 0.07, 0.07, 
0.07, 0.07, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 
0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 
0.08, 0.08, 0.07, 0.07, 0.06, 0.06, 0.05, 0.05, 0.05, 0.05, 0.05, 
0.05, 0.05, 0.04, 0.04, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 
0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03
)), .Names = c("MGOFujairah", "MGORotterdam", "X380CSTFujairah", 
"X380CSTRotterdam", "X380CSTHouston", "ContainerTCRIndex", "USAInterestRates"
), class = "data.frame", row.names = c(NA, -255L))


More information about the R-help mailing list