Model Regresi dengan Peubah Lag

Pendahuluan

Latar Belakang

Tekanan udara merupakan tenaga yang bekerja untuk menggerakkan massa udara dalam setiap satuan luas tertentu (Siswanti, 2011). Perubahan tekanan udara akan menyebabkan perubahan pada suhu udara dan curah hujan. Dengan demikian penyebaran curah hujan di seluruh permukaan bumi berhubungan sangat erat dengan sistem tekanan udara. Tekanan udara berkurang dengan bertambahnya ketinggian tempat.Menuut BMKG, Tekanan udara merupakan suatu tenaga yang bekerja untuk menggerakkan massa udara dalam setiap satuan luas tertentu, tenaga yang menggerakkan massa udara tersebut menekan searah gaya gravitasi bumi. Satuan tekanan udara adalah milibar (mb) atau Hecto Pascal (hPa). Tekanan udara dipengaruhi oleh temperature/suhu udara yang terjadi pada suatu tempat dan waktu, apabila temperatur udara tinggi maka volume molekul/partikel udara akan berkembang sehingga tekanan udara menjadi rendah dan berbanding sebaliknya. Tekanan udara juga merupakan salah satu parameter yang berkaitan erat dengan variasi pembentukan arah dan kecepatan angin.

Pada analisis ini, dilakukan analisis regresi pada data rata-rata tekanan udara terhadap rata-rata temperatur udara. Analisis regresi yang digunakan yaitu dengan model koyck, distributed lag, dan autoregressive model.

Data

Beberapa library yang digunakan pada analisis model regresi ini sebagai berikut:

library(dLagM) #bisa otomatis timeseries datanya
## Loading required package: nardl
## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo
## Loading required package: dynlm
## Loading required package: zoo
## 
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric
library(dynlm) #data harus timeseries
library(MLmetrics) #MAPE
## 
## Attaching package: 'MLmetrics'
## The following object is masked from 'package:dLagM':
## 
##     MAPE
## The following object is masked from 'package:base':
## 
##     Recall
library(lmtest)
library(car)
## Loading required package: carData
library(zoo)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following object is masked from 'package:car':
## 
##     recode
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(readxl)
library(readr)

Data yang digunakan merupakan data iklim harian di kota Delhi selama 3 tahun yaitu pada periode tahun 2013-2015yang berasal dari web kaggle: https://www.kaggle.com/datasets/sumanthvrao/daily-climate-time-series-data.

Peubah yang digunakan dalam model regresi adalah date, meanpressure, dan meantemp dengan jumlah amatan sebanyak 1095.

data <- read.csv("D:/SEMESTER 5/DailyDelhiClimateTrain.csv")
str(data)
## 'data.frame':    1095 obs. of  5 variables:
##  $ date        : chr  "2013-01-01" "2013-01-02" "2013-01-03" "2013-01-04" ...
##  $ meantemp    : num  10 7.4 7.17 8.67 6 ...
##  $ humidity    : num  84.5 92 87 71.3 86.8 ...
##  $ wind_speed  : num  0 2.98 4.63 1.23 3.7 ...
##  $ meanpressure: num  1016 1018 1019 1017 1016 ...
data
##            date  meantemp humidity wind_speed meanpressure
## 1    2013-01-01 10.000000 84.50000  0.0000000    1015.6667
## 2    2013-01-02  7.400000 92.00000  2.9800000    1017.8000
## 3    2013-01-03  7.166667 87.00000  4.6333333    1018.6667
## 4    2013-01-04  8.666667 71.33333  1.2333333    1017.1667
## 5    2013-01-05  6.000000 86.83333  3.7000000    1016.5000
## 6    2013-01-06  7.000000 82.80000  1.4800000    1018.0000
## 7    2013-01-07  7.000000 78.60000  6.3000000    1020.0000
## 8    2013-01-08  8.857143 63.71429  7.1428571    1018.7143
## 9    2013-01-09 14.000000 51.25000 12.5000000    1017.0000
## 10   2013-01-10 11.000000 62.00000  7.4000000    1015.6667
## 11   2013-01-11 15.714286 51.28571 10.5714286    1016.1429
## 12   2013-01-12 14.000000 74.00000 13.2285714    1015.5714
## 13   2013-01-13 15.833333 75.16667  4.6333333    1013.3333
## 14   2013-01-14 12.833333 88.16667  0.6166667    1015.1667
## 15   2013-01-15 14.714286 71.85714  0.5285714    1015.8571
## 16   2013-01-16 13.833333 86.66667  0.0000000    1016.6667
## 17   2013-01-17 16.500000 80.83333  5.2500000    1015.8333
## 18   2013-01-18 13.833333 92.16667  8.9500000    1014.5000
## 19   2013-01-19 12.500000 76.66667  5.8833333    1021.6667
## 20   2013-01-20 11.285714 75.28571  8.4714286    1020.2857
## 21   2013-01-21 11.200000 77.00000  2.2200000    1021.0000
## 22   2013-01-22  9.500000 79.66667  3.0833333    1021.8000
## 23   2013-01-23 14.000000 60.16667  4.0166667    1020.5000
## 24   2013-01-24 13.833333 60.66667  6.1666667    1020.5000
## 25   2013-01-25 12.250000 67.00000  5.5500000    1020.7500
## 26   2013-01-26 12.666667 64.16667  6.8000000    1019.6667
## 27   2013-01-27 12.857143 65.57143  5.5571429    1018.1429
## 28   2013-01-28 14.833333 56.00000  3.7000000    1017.8333
## 29   2013-01-29 14.125000 65.50000  3.2375000    1016.6250
## 30   2013-01-30 14.714286 70.42857  1.0571429    1017.8571
## 31   2013-01-31 16.200000 65.60000  2.9600000    1018.4000
## 32   2013-02-01 16.000000 73.00000  2.2200000    1016.0000
## 33   2013-02-02 16.285714 77.57143  1.3285714    1017.1429
## 34   2013-02-03 18.000000 65.57143  1.8571429    1015.2857
## 35   2013-02-04 17.428571 74.28571 11.1142857    1014.5714
## 36   2013-02-05 16.625000 92.37500  9.7250000    1016.3750
## 37   2013-02-06 16.666667 71.33333  8.6333333    1018.6667
## 38   2013-02-07 15.600000 59.40000 10.7400000    1018.6000
## 39   2013-02-08 14.000000 70.42857  9.2571429    1017.1429
## 40   2013-02-09 15.428571 61.28571  9.2571429    1016.8571
## 41   2013-02-10 15.250000 71.50000  3.4750000    1017.1250
## 42   2013-02-11 15.875000 70.50000  5.3250000    1016.5000
## 43   2013-02-12 15.333333 70.33333  7.4166667    1017.5000
## 44   2013-02-13 16.285714 70.14286  6.0857143    1016.0000
## 45   2013-02-14 17.333333 63.83333  4.3333333    1014.1667
## 46   2013-02-15 19.166667 65.33333 10.1833333    1011.6667
## 47   2013-02-16 14.428571 92.71429  8.4857143    1008.0000
## 48   2013-02-17 13.666667 90.00000  0.0000000    1012.6667
## 49   2013-02-18 15.600000 78.40000 12.2200000    1016.2000
## 50   2013-02-19 15.857143 82.00000  5.8142857    1015.7143
## 51   2013-02-20 17.714286 74.71429  5.8142857    1017.0000
## 52   2013-02-21 20.000000 67.28571  6.6142857    1015.4286
## 53   2013-02-22 20.500000 65.62500 10.8750000    1016.0000
## 54   2013-02-23 17.428571 74.85714  9.2571429    1017.0000
## 55   2013-02-24 16.857143 78.85714  7.4000000    1018.8571
## 56   2013-02-25 16.875000 72.87500  4.6375000    1018.2500
## 57   2013-02-26 17.857143 70.00000 17.5875000    1015.1429
## 58   2013-02-27 20.800000 57.20000  6.6600000    1015.2000
## 59   2013-02-28 19.428571 52.85714 12.9571429    1017.4286
## 60   2013-03-01 17.333333 49.33333 24.0666667    1016.3333
## 61   2013-03-02 19.000000 54.00000 15.7250000    1016.2500
## 62   2013-03-03 19.333333 62.83333  8.6333333    1016.1667
## 63   2013-03-04 17.600000 71.00000  5.5600000    1015.8000
## 64   2013-03-05 20.875000 61.87500  4.1625000    1016.3750
## 65   2013-03-06 20.857143 65.28571  6.8714286    1015.7143
## 66   2013-03-07 23.428571 57.14286  8.7285714    1015.2857
## 67   2013-03-08 24.166667 44.83333  7.1000000    1014.8333
## 68   2013-03-09 25.428571 49.71429  5.2857143    1009.2857
## 69   2013-03-10 23.142857 57.57143  4.2285714    1008.0000
## 70   2013-03-11 24.000000 66.33333  0.9333333    1011.1667
## 71   2013-03-12 23.500000 62.50000  4.9333333    1011.5000
## 72   2013-03-13 21.500000 70.50000  5.5500000    1009.0000
## 73   2013-03-14 22.333333 61.16667 12.0333333    1013.1667
## 74   2013-03-15 24.166667 45.83333  7.7166667    1016.1667
## 75   2013-03-16 20.333333 67.66667  3.7000000    1016.1667
## 76   2013-03-17 22.666667 58.66667  8.9500000    1015.0000
## 77   2013-03-18 23.428571 58.14286 17.4571429    1009.4286
## 78   2013-03-19 22.500000 73.66667 10.4833333    1011.0000
## 79   2013-03-20 29.166667 36.33333  6.8000000    1009.5000
## 80   2013-03-21 23.833333 58.50000 10.5000000    1008.3333
## 81   2013-03-22 25.250000 50.25000  9.7125000    1006.3750
## 82   2013-03-23 27.375000 50.12500  9.9500000    1007.6250
## 83   2013-03-24 27.000000 48.75000 10.8875000    1010.2500
## 84   2013-03-25 23.500000 45.50000 15.9625000    1010.6250
## 85   2013-03-26 24.142857 44.57143 12.9571429    1008.8571
## 86   2013-03-27 21.000000 62.00000  1.8500000    1009.0000
## 87   2013-03-28 22.428571 62.71429  7.4142857    1009.5714
## 88   2013-03-29 21.250000 70.37500  5.5500000    1009.7500
## 89   2013-03-30 23.500000 54.75000 11.1125000    1008.6250
## 90   2013-03-31 23.200000 58.00000  6.6600000    1008.6000
## 91   2013-04-01 25.375000 45.50000  4.4000000    1008.5000
## 92   2013-04-02 25.166667 51.00000  8.6500000    1009.5000
## 93   2013-04-03 26.200000 45.60000  8.1400000    1009.0000
## 94   2013-04-04 24.600000 41.80000 11.1200000    1007.8000
## 95   2013-04-05 25.600000 31.00000 15.5600000    1007.0000
## 96   2013-04-06 25.857143 29.85714 11.9000000    1006.1429
## 97   2013-04-07 29.142857 23.28571 10.3142857    1005.0000
## 98   2013-04-08 28.714286 33.85714  5.3000000    1006.0000
## 99   2013-04-09 30.166667 30.50000  8.6500000    1005.3333
## 100  2013-04-10 30.000000 28.00000  6.1000000    1006.7143
## 101  2013-04-11 30.000000 24.20000  7.7800000    1006.4000
## 102  2013-04-12 28.857143 32.57143  6.3428571    1007.5714
## 103  2013-04-13 30.200000 29.20000 10.0000000    1008.4000
## 104  2013-04-14 28.250000 39.37500  6.4875000    1007.0000
## 105  2013-04-15 28.250000 41.37500  6.2500000    1003.8750
## 106  2013-04-16 32.125000 24.62500 10.4250000    1000.0000
## 107  2013-04-17 29.200000 24.20000  6.6600000    1002.2000
## 108  2013-04-18 30.285714 30.28571  4.7571429    1002.8571
## 109  2013-04-19 28.285714 31.28571  3.9714286    1002.5714
## 110  2013-04-20 30.625000 29.00000  7.6375000    1003.3750
## 111  2013-04-21 27.666667 38.66667 13.8833333    1006.8333
## 112  2013-04-22 27.375000 45.37500  7.6500000    1010.0000
## 113  2013-04-23 28.625000 44.12500  4.6250000    1009.8750
## 114  2013-04-24 30.285714 41.71429  2.1142857    1008.5714
## 115  2013-04-25 31.142857 38.28571  3.7000000    1007.7143
## 116  2013-04-26 29.875000 45.87500  6.7125000    1008.3750
## 117  2013-04-27 31.142857 31.42857 13.4857143    1007.0000
## 118  2013-04-28 30.571429 28.00000 12.9714286    1005.4286
## 119  2013-04-29 32.125000 26.37500  7.8750000    1004.8750
## 120  2013-04-30 31.142857 32.00000  7.9285714    1004.8571
## 121  2013-05-01 31.857143 15.85714 12.6857143    1002.8333
## 122  2013-05-02 29.833333 22.16667 11.4333333    1001.2000
## 123  2013-05-03 28.571429 31.57143  9.0000000     999.4286
## 124  2013-05-04 32.857143 31.42857  2.9142857     999.0000
## 125  2013-05-05 32.625000 31.12500  3.0125000    1000.6250
## 126  2013-05-06 32.750000 39.25000  3.7000000    1001.6250
## 127  2013-05-07 32.875000 33.25000  7.1750000    1002.0000
## 128  2013-05-08 34.500000 23.00000  9.2500000    1001.1667
## 129  2013-05-09 34.285714 26.00000  8.9857143    1000.8571
## 130  2013-05-10 34.000000 27.71429  9.5285714    1000.1429
## 131  2013-05-11 30.750000 30.37500 14.8125000     999.8750
## 132  2013-05-12 29.857143 40.14286  5.8285714    1003.7143
## 133  2013-05-13 31.714286 34.00000  3.9714286    1003.5714
## 134  2013-05-14 32.285714 34.28571  4.7714286    1001.5714
## 135  2013-05-15 33.000000 33.00000  3.2375000     999.8571
## 136  2013-05-16 33.000000 34.75000  7.1750000     998.5000
## 137  2013-05-17 32.833333 28.16667  9.8666667    1000.3333
## 138  2013-05-18 31.400000 42.20000 12.9666667     999.0000
## 139  2013-05-19 35.333333 22.33333 15.7500000     999.6667
## 140  2013-05-20 36.400000 24.20000  7.4000000     998.4000
## 141  2013-05-21 36.000000 19.00000 11.3714286     998.6667
## 142  2013-05-22 36.750000 22.12500 17.5875000     998.6250
## 143  2013-05-23 37.500000 23.33333 13.5666667     997.1667
## 144  2013-05-24 38.428571 27.42857 11.3857143     996.4286
## 145  2013-05-25 38.714286 22.42857 10.3142857     998.1429
## 146  2013-05-26 37.800000 21.20000 10.7400000     998.2000
## 147  2013-05-27 35.857143 31.57143  8.7285714     999.1429
## 148  2013-05-28 35.333333 29.00000  8.3333333    1000.1667
## 149  2013-05-29 34.142857 27.85714  5.5571429     999.4286
## 150  2013-05-30 32.200000 28.20000  5.5600000     999.6000
## 151  2013-05-31 33.625000 40.12500 10.6375000     998.7143
## 152  2013-06-01 32.000000 54.00000 13.4375000     998.7500
## 153  2013-06-02 32.400000 56.60000 14.2000000    1001.0000
## 154  2013-06-03 35.600000 47.00000 12.2400000     999.6000
## 155  2013-06-04 35.857143 42.57143  5.0285714     999.8571
## 156  2013-06-05 37.166667 36.50000  9.8666667     998.0000
## 157  2013-06-06 31.285714 61.71429 10.0428571     997.2857
## 158  2013-06-07 34.000000 49.25000  7.4125000     998.0000
## 159  2013-06-08 34.200000 57.40000 13.3400000     997.0000
## 160  2013-06-09 36.166667 40.00000 11.4166667     995.1667
## 161  2013-06-10 36.625000 42.75000  8.1000000     995.1250
## 162  2013-06-11 30.166667 61.66667 16.3500000    1000.0000
## 163  2013-06-12 34.142857 54.14286  7.4142857     998.0000
## 164  2013-06-13 29.833333 68.66667  9.8666667     998.6667
## 165  2013-06-14 30.142857 68.14286  7.9428571     997.4286
## 166  2013-06-15 30.714286 65.00000 13.2285714     997.8571
## 167  2013-06-16 27.000000 88.85714  8.4857143     996.4286
## 168  2013-06-17 26.875000 87.37500  8.3375000     994.7500
## 169  2013-06-18 28.400000 83.60000  5.1800000     996.4000
## 170  2013-06-19 29.857143 72.85714  4.2285714     999.0000
## 171  2013-06-20 33.000000 52.71429  4.7714286     997.1667
## 172  2013-06-21 34.833333 47.33333  5.8666667     995.3333
## 173  2013-06-22 35.600000 39.40000 10.7400000     996.2000
## 174  2013-06-23 35.166667 53.00000  8.3500000     995.6667
## 175  2013-06-24 33.142857 55.00000 13.2285714     996.2857
## 176  2013-06-25 30.571429 70.42857 11.1142857     999.1667
## 177  2013-06-26 30.666667 71.50000  6.4833333     999.0000
## 178  2013-06-27 31.428571 61.85714  7.9571429     996.7143
## 179  2013-06-28 31.500000 64.25000  8.8125000     996.0000
## 180  2013-06-29 33.250000 57.00000 13.8750000     995.2500
## 181  2013-06-30 32.833333 52.16667 10.5000000     997.1667
## 182  2013-07-01 33.857143 54.00000  8.7285714     996.4286
## 183  2013-07-02 33.142857 58.28571  9.8000000     996.0000
## 184  2013-07-03 31.571429 70.00000  6.3571429     997.5714
## 185  2013-07-04 32.375000 67.50000  9.2500000     997.0000
## 186  2013-07-05 32.800000 62.20000 12.9600000     997.8000
## 187  2013-07-06 31.000000 73.14286 11.3714286    1000.5714
## 188  2013-07-07 31.166667 77.33333  7.1000000    1000.1667
## 189  2013-07-08 29.833333 75.16667  6.4833333     997.5000
## 190  2013-07-09 29.000000 78.83333  8.0166667     996.8333
## 191  2013-07-10 29.750000 75.37500  6.4875000     995.7500
## 192  2013-07-11 32.500000 65.50000  7.1166667     997.3333
## 193  2013-07-12 28.875000 86.62500  5.5625000    1000.0000
## 194  2013-07-13 31.750000 70.25000 12.0625000     998.8750
## 195  2013-07-14 30.750000 71.62500  4.8125000     997.3750
## 196  2013-07-15 31.750000 73.25000 13.5375000     998.3750
## 197  2013-07-16 29.875000 83.87500  9.0000000     999.0000
## 198  2013-07-17 31.125000 78.37500  8.6500000     998.0000
## 199  2013-07-18 31.750000 76.00000  5.5500000     995.0000
## 200  2013-07-19 30.500000 73.50000  9.2666667     995.3333
## 201  2013-07-20 26.833333 90.66667  4.9333333     999.1667
## 202  2013-07-21 28.200000 88.00000  2.2400000     996.4000
## 203  2013-07-22 29.500000 80.33333  4.3333333     996.8333
## 204  2013-07-23 31.857143 70.71429  3.9714286     998.0000
## 205  2013-07-24 29.714286 84.14286  3.9714286     997.1429
## 206  2013-07-25 28.333333 84.50000  8.0166667     995.1667
## 207  2013-07-26 30.000000 73.85714  5.5714286     995.4286
## 208  2013-07-27 30.142857 79.00000  4.2285714     997.2857
## 209  2013-07-28 32.285714 67.42857  6.0857143     998.1429
## 210  2013-07-29 31.000000 77.42857  4.7714286     997.1429
## 211  2013-07-30 30.500000 75.66667  8.0166667     996.5000
## 212  2013-07-31 28.833333 78.50000  9.8666667     996.6667
## 213  2013-08-01 30.000000 73.50000  5.5500000     998.0000
## 214  2013-08-02 31.571429 66.42857 11.1142857     999.4286
## 215  2013-08-03 32.500000 62.33333 10.8166667    1000.8333
## 216  2013-08-04 32.500000 62.50000  5.1000000     993.2500
## 217  2013-08-05 29.500000 77.00000  6.7833333    1001.0000
## 218  2013-08-06 27.166667 90.66667  4.0166667    1002.1667
## 219  2013-08-07 28.375000 85.25000  3.2500000    1001.8750
## 220  2013-08-08 28.500000 84.00000  6.1666667    1000.5000
## 221  2013-08-09 27.166667 87.50000  3.7166667    1001.0000
## 222  2013-08-10 29.428571 79.85714  6.1000000    1000.3333
## 223  2013-08-11 30.000000 76.83333  6.8000000    1001.1667
## 224  2013-08-12 29.500000 81.50000  8.3500000    1001.0000
## 225  2013-08-13 30.750000 72.00000 11.1250000     999.2500
## 226  2013-08-14 29.666667 74.33333  9.8666667     999.8333
## 227  2013-08-15 27.714286 85.00000  8.2000000    1001.7143
## 228  2013-08-16 26.600000 92.00000  4.4600000    1000.8000
## 229  2013-08-17 27.428571 91.00000 14.2857143     999.4286
## 230  2013-08-18 28.333333 83.83333  1.2333333    1000.1667
## 231  2013-08-19 28.166667 84.66667  5.5500000    1000.5000
## 232  2013-08-20 27.375000 88.50000  5.3250000     999.0000
## 233  2013-08-21 28.333333 81.16667 10.3142857     999.2857
## 234  2013-08-22 29.166667 74.83333  6.4833333    1001.1667
## 235  2013-08-23 31.285714 73.14286  2.6428571    1001.4286
## 236  2013-08-24 32.000000 71.25000 18.5250000    1000.3333
## 237  2013-08-25 29.000000 90.75000  0.9250000     998.7500
## 238  2013-08-26 32.000000 69.87500  6.4750000     997.8571
## 239  2013-08-27 32.142857 57.28571 17.0714286     999.8571
## 240  2013-08-28 30.142857 60.71429 11.6428571    1003.7143
## 241  2013-08-29 31.500000 57.00000  9.2750000    1003.2500
## 242  2013-08-30 29.166667 73.00000  3.1000000    1003.5000
## 243  2013-08-31 29.000000 71.75000  4.6500000    1003.2500
## 244  2013-09-01 30.000000 65.20000 11.4800000    1003.2000
## 245  2013-09-02 31.857143 50.85714 15.3428571    1003.7143
## 246  2013-09-03 32.750000 45.50000 14.8000000    1001.5000
## 247  2013-09-04 31.714286 48.00000 14.5428571    1002.1667
## 248  2013-09-05 31.666667 46.00000 13.5666667    1001.5000
## 249  2013-09-06 30.571429 55.71429 10.0428571    1003.4286
## 250  2013-09-07 30.833333 51.83333  6.7833333    1004.1667
## 251  2013-09-08 28.000000 71.75000  2.3250000    1007.2500
## 252  2013-09-09 31.000000 57.50000 42.2200000    1007.0000
## 253  2013-09-10 29.666667 66.33333  3.7000000    1005.3333
## 254  2013-09-11 30.857143 61.14286  3.9714286    1005.5714
## 255  2013-09-12 31.166667 60.33333  2.7833333    1008.1667
## 256  2013-09-13 31.142857 67.14286  3.1714286    1007.1429
## 257  2013-09-14 30.500000 67.00000  2.7833333    1003.8333
## 258  2013-09-15 31.625000 56.37500  3.9375000    1001.3750
## 259  2013-09-16 29.666667 60.00000  7.1000000    1002.3333
## 260  2013-09-17 29.250000 56.50000  3.4750000    1002.8750
## 261  2013-09-18 29.142857 58.57143  5.5571429    1003.1429
## 262  2013-09-19 29.800000 61.60000  2.9600000    1001.8000
## 263  2013-09-20 28.666667 70.00000  9.5833333    1003.0000
## 264  2013-09-21 25.200000 89.00000  5.1800000    1003.0000
## 265  2013-09-22 28.333333 79.66667  3.0833333    1001.6667
## 266  2013-09-23 30.285714 68.00000  4.5000000    1001.0000
## 267  2013-09-24 30.750000 65.75000  5.7875000    1002.7500
## 268  2013-09-25 28.571429 75.00000  6.6142857    1005.4286
## 269  2013-09-26 28.200000 78.80000  4.0800000    1006.6000
## 270  2013-09-27 29.000000 71.25000  7.4125000    1005.6250
## 271  2013-09-28 30.000000 67.33333  1.2333333    1005.3333
## 272  2013-09-29 28.142857 79.85714  5.5571429    1003.8571
## 273  2013-09-30 26.857143 84.57143  4.2285714    1004.7143
## 274  2013-10-01 28.285714 68.71429  6.3714286    1007.7143
## 275  2013-10-02 29.200000 61.40000  7.4200000    1010.0000
## 276  2013-10-03 28.600000 59.00000  7.4000000    1008.6000
## 277  2013-10-04 24.833333 84.33333  6.7833333    1007.5000
## 278  2013-10-05 28.500000 69.00000  2.0875000    1004.1250
## 279  2013-10-06 29.750000 66.50000  1.8500000    1005.2500
## 280  2013-10-07 29.714286 66.57143  1.5857143    1008.8571
## 281  2013-10-08 29.200000 69.40000  0.0000000    1009.8000
## 282  2013-10-09 30.833333 65.50000  5.5666667    1010.3333
## 283  2013-10-10 29.428571 66.28571  3.9714286    1010.8571
## 284  2013-10-11 23.857143 89.85714  4.2428571    1010.5714
## 285  2013-10-12 26.142857 82.57143  2.9142857    1008.2857
## 286  2013-10-13 27.166667 84.33333  0.6166667    1007.6667
## 287  2013-10-14 27.714286 74.85714  5.5571429    1007.0000
## 288  2013-10-15 25.857143 76.42857  3.7000000    1009.1429
## 289  2013-10-16 26.428571 75.00000  0.0000000    1009.0000
## 290  2013-10-17 26.857143 76.71429  1.0571429    1010.4286
## 291  2013-10-18 26.333333 73.16667  0.0000000    1010.5000
## 292  2013-10-19 24.571429 69.57143  2.1142857    1011.5714
## 293  2013-10-20 24.333333 67.83333  1.5500000    1010.6667
## 294  2013-10-21 24.875000 63.62500  1.3875000    1009.8750
## 295  2013-10-22 27.800000 53.60000  0.0000000    1011.0000
## 296  2013-10-23 25.000000 67.50000  1.1625000    1013.6250
## 297  2013-10-24 23.857143 73.14286  1.0571429    1013.2857
## 298  2013-10-25 22.800000 75.80000  1.8600000    1011.0000
## 299  2013-10-26 23.000000 78.20000  0.7400000    1011.2000
## 300  2013-10-27 22.875000 66.87500  2.5500000    1010.6250
## 301  2013-10-28 22.666667 63.16667  4.9333333    1010.0000
## 302  2013-10-29 23.000000 62.71429  3.9714286    1011.1429
## 303  2013-10-30 22.857143 67.42857  2.3857143    1012.8571
## 304  2013-10-31 23.666667 58.83333  3.0833333    1012.8333
## 305  2013-11-01 23.428571 69.14286  1.0571429    1013.2857
## 306  2013-11-02 21.625000 65.00000  4.6375000    1014.5000
## 307  2013-11-03 20.625000 57.87500  5.5625000    1016.0000
## 308  2013-11-04 21.166667 54.16667  2.4666667    1017.5000
## 309  2013-11-05 18.833333 74.33333  0.0000000    1017.5000
## 310  2013-11-06 20.571429 66.71429  3.1857143    1015.8571
## 311  2013-11-07 20.142857 73.14286  4.5000000    1015.6667
## 312  2013-11-08 18.000000 82.20000 19.9125000    1015.6250
## 313  2013-11-09 19.857143 73.57143  0.5285714    1016.7143
## 314  2013-11-10 16.833333 72.00000  2.4666667    1016.0000
## 315  2013-11-11 18.857143 61.85714  5.8142857    1013.4286
## 316  2013-11-12 16.571429 71.71429  0.8000000    1013.8571
## 317  2013-11-13 17.750000 65.00000  3.0125000    1014.2500
## 318  2013-11-14 17.625000 63.37500  3.4750000    1015.2500
## 319  2013-11-15 17.000000 67.00000  0.9250000    1016.5000
## 320  2013-11-16 16.625000 69.62500  0.0000000    1018.0000
## 321  2013-11-17 15.571429 67.42857  2.3857143    1017.7143
## 322  2013-11-18 18.500000 55.12500  7.4125000    1016.8750
## 323  2013-11-19 17.875000 63.50000  7.1750000    1015.5000
## 324  2013-11-20 18.250000 58.87500  2.5500000    1012.7500
## 325  2013-11-21 17.875000 66.25000  2.0875000    1013.1250
## 326  2013-11-22 17.625000 74.87500  1.1625000    1015.7500
## 327  2013-11-23 18.142857 77.85714  0.8000000    1014.5714
## 328  2013-11-24 19.125000 58.87500  5.0875000    1014.5000
## 329  2013-11-25 21.250000 49.87500  9.5000000    1012.8750
## 330  2013-11-26 21.250000 51.37500  9.7500000    1011.7500
## 331  2013-11-27 19.125000 68.50000  2.7875000    1014.7500
## 332  2013-11-28 18.625000 71.87500  1.8500000    1017.7500
## 333  2013-11-29 17.750000 66.62500  4.8625000    1016.8750
## 334  2013-11-30 17.875000 69.62500 11.8125000    1015.8750
## 335  2013-12-01 18.000000 65.00000  1.5857143    1016.1429
## 336  2013-12-02 17.250000 73.25000  1.8500000    1015.8750
## 337  2013-12-03 17.500000 73.50000  1.3875000    1016.7500
## 338  2013-12-04 17.142857 71.85714  4.2285714    1017.0000
## 339  2013-12-05 17.142857 69.00000  6.8857143    1015.2857
## 340  2013-12-06 16.125000 72.37500  3.7125000    1014.0000
## 341  2013-12-07 16.000000 80.50000  1.1625000    1014.2500
## 342  2013-12-08 16.500000 83.62500  0.0000000    1014.0000
## 343  2013-12-09 16.125000 79.25000  2.3125000    1013.5000
## 344  2013-12-10 15.500000 81.75000  2.5500000    1013.7500
## 345  2013-12-11 16.625000 71.87500  6.4750000    1014.1250
## 346  2013-12-12 17.750000 62.37500 16.9000000    1012.5000
## 347  2013-12-13 17.142857 67.71429 11.6285714    1012.0000
## 348  2013-12-14 16.142857 79.00000  0.8000000    1012.5714
## 349  2013-12-15 15.500000 86.25000  2.5500000    1015.1250
## 350  2013-12-16 15.250000 89.50000  1.1625000    1015.6250
## 351  2013-12-17 14.750000 92.87500  0.9250000    1014.2500
## 352  2013-12-18 14.875000 93.00000  2.0875000    1012.6250
## 353  2013-12-19 16.125000 89.12500  3.2500000    1012.8750
## 354  2013-12-20 15.375000 91.37500  4.4125000    1016.0000
## 355  2013-12-21 14.750000 94.00000  0.4625000    1017.0000
## 356  2013-12-22 15.250000 94.00000  1.3875000    1018.3750
## 357  2013-12-23 14.250000 88.25000  9.0250000    1020.5000
## 358  2013-12-24 13.500000 84.75000  5.7875000    1021.3750
## 359  2013-12-25 13.666667 63.16667 34.4875000    1020.6250
## 360  2013-12-26 12.125000 67.62500  4.4000000    1019.2500
## 361  2013-12-27 11.875000 79.87500  1.1625000    1018.6250
## 362  2013-12-28 10.875000 70.00000  5.3250000    1019.2500
## 363  2013-12-29 10.571429 69.42857  5.3250000    1018.5000
## 364  2013-12-30 12.375000 79.50000  6.4750000    1018.1250
## 365  2013-12-31 14.500000 89.37500  4.8625000    1020.5000
## 366  2014-01-01 13.375000 89.62500  7.6500000    1021.0000
## 367  2014-01-02 11.000000 78.37500  8.1000000    1020.2500
## 368  2014-01-03 12.500000 74.87500  5.3250000    1017.7500
## 369  2014-01-04 12.875000 88.12500  1.1625000    1016.2500
## 370  2014-01-05 12.375000 89.00000  0.4625000    1014.5000
## 371  2014-01-06 11.428571 86.28571 27.3375000    1017.6250
## 372  2014-01-07 12.142857 72.28571 20.6000000    1017.7500
## 373  2014-01-08 11.875000 76.87500  1.8500000    1017.8750
## 374  2014-01-09 12.833333 83.00000  1.2333333    1017.0000
## 375  2014-01-10 12.375000 85.37500  0.0000000    1016.6250
## 376  2014-01-11 12.125000 82.62500  6.7250000    1019.6250
## 377  2014-01-12 12.875000 80.87500  3.2500000    1020.7500
## 378  2014-01-13 14.000000 72.37500  1.1625000    1019.3750
## 379  2014-01-14 14.875000 76.75000  2.7875000    1017.0000
## 380  2014-01-15 12.000000 93.50000  0.4625000    1019.7500
## 381  2014-01-16 12.285714 90.57143  1.5857143    1019.5714
## 382  2014-01-17 12.000000 95.75000  1.1625000    1017.6250
## 383  2014-01-18 12.500000 94.25000  1.8500000    1020.2500
## 384  2014-01-19 14.500000 84.50000  3.4750000    1023.0000
## 385  2014-01-20 14.625000 83.50000  4.4125000    1021.3750
## 386  2014-01-21 13.571429 96.85714  7.1714286    1020.5714
## 387  2014-01-22 15.250000 96.12500  6.7125000    1019.8750
## 388  2014-01-23 15.625000 90.12500  6.7125000    1020.2500
## 389  2014-01-24 13.875000 88.87500  8.1125000    1019.2500
## 390  2014-01-25 11.750000 95.62500  8.3375000    1019.7500
## 391  2014-01-26 14.375000 85.75000  8.8125000    1019.1250
## 392  2014-01-27 16.625000 72.37500  9.9500000    1018.7500
## 393  2014-01-28 16.625000 69.87500  7.4125000    1016.7500
## 394  2014-01-29 14.875000 81.87500  4.6250000    1018.6250
## 395  2014-01-30 14.500000 83.25000  0.9250000    1019.3750
## 396  2014-01-31 14.750000 82.50000  3.0125000    1019.2500
## 397  2014-02-01 14.000000 79.37500  2.5500000    1019.7500
## 398  2014-02-02 15.250000 76.12500  1.8500000    1016.8750
## 399  2014-02-03 18.428571 63.85714 12.2875000    1014.1250
## 400  2014-02-04 18.857143 65.42857 13.4250000    1011.0000
## 401  2014-02-05 17.000000 72.75000  1.6250000    1011.5000
## 402  2014-02-06 18.500000 65.37500  6.0125000    1010.8750
## 403  2014-02-07 19.625000 67.25000 15.0500000    1009.2500
## 404  2014-02-08 13.625000 85.37500  8.1125000    1011.3750
## 405  2014-02-09 12.875000 63.37500  9.2625000    1012.0000
## 406  2014-02-10 12.625000 61.75000  3.7000000    1012.5000
## 407  2014-02-11 13.375000 64.62500  3.7000000    1013.5000
## 408  2014-02-12 13.250000 67.00000  9.2625000    1013.5000
## 409  2014-02-13 15.250000 64.50000  4.8750000    1013.6250
## 410  2014-02-14 13.000000 95.50000  5.3375000    1011.8750
## 411  2014-02-15 12.375000 89.50000  3.9375000    1012.8750
## 412  2014-02-16 13.500000 79.00000  8.7875000    1016.0000
## 413  2014-02-17 15.125000 66.12500 12.7250000    1017.7500
## 414  2014-02-18 15.250000 60.12500  8.8000000    1017.3750
## 415  2014-02-19 15.125000 75.37500  1.8500000    1016.6250
## 416  2014-02-20 16.125000 78.62500  1.8625000    1014.7500
## 417  2014-02-21 17.125000 72.25000  8.7875000    1014.2500
## 418  2014-02-22 16.625000 85.62500  6.2500000    1017.3750
## 419  2014-02-23 16.500000 80.00000  7.8875000    1019.3750
## 420  2014-02-24 16.375000 69.37500  5.3375000    1018.6250
## 421  2014-02-25 18.125000 69.00000  0.9250000    1017.2500
## 422  2014-02-26 18.750000 71.12500  2.0875000    1014.2500
## 423  2014-02-27 18.625000 76.62500  7.1750000    1013.1250
## 424  2014-02-28 15.625000 88.37500  5.8000000    1013.1250
## 425  2014-03-01 16.142857 87.28571  4.2428571    1015.2857
## 426  2014-03-02 16.125000 79.00000  5.7875000    1017.1250
## 427  2014-03-03 17.875000 65.62500  3.7000000    1016.2500
## 428  2014-03-04 18.750000 66.37500  3.7125000    1016.5000
## 429  2014-03-05 18.875000 67.00000  5.3250000    1016.1250
## 430  2014-03-06 18.375000 64.00000  4.6375000    1016.3750
## 431  2014-03-07 18.500000 62.87500  6.7125000    1017.1250
## 432  2014-03-08 19.875000 63.87500  9.9625000    1017.0000
## 433  2014-03-09 22.000000 56.62500  7.6500000    1014.5000
## 434  2014-03-10 20.000000 63.87500 10.8875000    1014.3750
## 435  2014-03-11 20.714286 70.00000 15.6000000    1013.7143
## 436  2014-03-12 20.000000 68.87500  6.4875000    1016.2500
## 437  2014-03-13 19.000000 60.62500  9.9500000    1016.2500
## 438  2014-03-14 19.625000 61.25000  5.8000000    1015.0000
## 439  2014-03-15 21.625000 62.62500  2.7875000    1016.3750
## 440  2014-03-16 23.125000 58.50000  6.2500000    1015.3750
## 441  2014-03-17 26.250000 52.12500  6.7125000    1011.2500
## 442  2014-03-18 24.857143 53.42857  9.8000000    1012.8571
## 443  2014-03-19 21.750000 57.37500 10.1875000    1013.2500
## 444  2014-03-20 21.500000 55.00000  7.4000000    1010.7500
## 445  2014-03-21 22.500000 64.75000  5.3375000    1007.1250
## 446  2014-03-22 25.000000 59.50000  3.0125000    1007.0000
## 447  2014-03-23 22.750000 61.87500  9.0250000    1010.6250
## 448  2014-03-24 22.375000 73.12500  9.4875000    1011.6250
## 449  2014-03-25 24.625000 60.25000  2.7750000    1010.5000
## 450  2014-03-26 25.125000 58.87500  2.7750000    1010.6250
## 451  2014-03-27 23.500000 63.25000  4.8875000    1011.7500
## 452  2014-03-28 24.428571 67.57143  4.5000000    1011.1429
## 453  2014-03-29 25.875000 55.62500  8.3500000    1011.2500
## 454  2014-03-30 24.125000 46.75000 13.9000000    1009.8750
## 455  2014-03-31 24.125000 44.12500 13.2000000    1008.2500
## 456  2014-04-01 24.875000 45.75000  9.7125000    1006.1250
## 457  2014-04-02 26.750000 43.37500  7.1750000    1005.7500
## 458  2014-04-03 25.000000 48.00000 10.4250000    1010.3750
## 459  2014-04-04 25.857143 47.85714 11.9142857    1009.8571
## 460  2014-04-05 27.375000 44.00000  9.2625000    1008.8750
## 461  2014-04-06 30.125000 45.25000  6.7125000    1007.5000
## 462  2014-04-07 27.250000 47.62500 12.9625000    1009.6250
## 463  2014-04-08 25.750000 41.87500 12.7250000    1010.3750
## 464  2014-04-09 25.750000 35.25000 12.7375000    1010.2500
## 465  2014-04-10 26.375000 29.25000 12.5000000    1010.8750
## 466  2014-04-11 28.142857 29.00000 10.0571429    1008.5714
## 467  2014-04-12 29.125000 39.87500  5.1000000    1007.1250
## 468  2014-04-13 26.875000 40.37500 15.5125000    1009.6250
## 469  2014-04-14 27.625000 46.62500  5.3500000    1010.6250
## 470  2014-04-15 29.000000 48.12500  2.5500000    1009.1250
## 471  2014-04-16 27.125000 53.62500  4.4000000    1008.2500
## 472  2014-04-17 26.375000 49.62500  9.2500000    1007.0000
## 473  2014-04-18 24.250000 64.37500 14.3375000    1011.6250
## 474  2014-04-19 25.625000 57.87500  5.1000000    1011.3750
## 475  2014-04-20 27.000000 48.00000  4.6375000    1010.5000
## 476  2014-04-21 28.000000 33.12500  8.8125000    1010.0000
## 477  2014-04-22 29.375000 29.37500 12.5000000    1011.7500
## 478  2014-04-23 29.125000 37.00000  6.0250000    1009.0000
## 479  2014-04-24 29.375000 35.87500  7.8750000    1005.8750
## 480  2014-04-25 30.625000 34.37500  8.1125000    1006.1250
## 481  2014-04-26 31.500000 32.00000  8.1125000    1006.3750
## 482  2014-04-27 30.625000 30.00000  8.3500000    1007.1250
## 483  2014-04-28 31.500000 27.87500 10.1875000    1005.6250
## 484  2014-04-29 32.125000 28.87500  8.3250000    1003.8750
## 485  2014-04-30 33.250000 25.50000  5.3375000    1003.7500
## 486  2014-05-01 34.875000 22.50000  9.7375000    1001.5000
## 487  2014-05-02 34.375000 24.75000  6.4875000    1001.8750
## 488  2014-05-03 31.750000 33.37500  6.4750000    1003.6250
## 489  2014-05-04 31.875000 34.37500  5.8000000    1003.0000
## 490  2014-05-05 29.750000 53.25000  9.0375000    1005.7500
## 491  2014-05-06 29.571429 51.14286  6.8857143    1006.8571
## 492  2014-05-07 31.750000 48.62500  1.8500000    1004.2500
## 493  2014-05-08 31.375000 43.00000 10.6625000    1002.6250
## 494  2014-05-09 30.500000 39.62500  6.4875000    1000.4286
## 495  2014-05-10 30.000000 37.62500  7.8875000    1000.7500
## 496  2014-05-11 29.625000 40.37500  7.6500000    1004.0000
## 497  2014-05-12 28.000000 53.37500 10.8875000    1004.8750
## 498  2014-05-13 25.250000 65.75000  9.9625000    1008.3750
## 499  2014-05-14 26.875000 55.25000  8.5750000    1006.5000
## 500  2014-05-15 29.500000 39.87500 11.1250000    1004.2500
## 501  2014-05-16 31.125000 40.25000  9.4875000    1003.5000
## 502  2014-05-17 29.875000 47.62500  7.1750000    1007.0000
## 503  2014-05-18 31.500000 35.37500 13.9000000    1007.3750
## 504  2014-05-19 32.375000 29.50000 11.5750000    1006.6250
## 505  2014-05-20 31.125000 34.00000  7.6375000    1006.6250
## 506  2014-05-21 33.000000 32.50000  9.7125000    1002.1250
## 507  2014-05-22 33.250000 33.37500 11.1125000     999.5000
## 508  2014-05-23 34.625000 32.75000 19.4500000    1000.0000
## 509  2014-05-24 31.500000 50.50000  6.7250000    1004.1250
## 510  2014-05-25 31.500000 43.25000 16.4375000    1005.5000
## 511  2014-05-26 30.125000 54.12500  4.1750000    1007.1250
## 512  2014-05-27 31.875000 49.75000  2.5500000    1004.7500
## 513  2014-05-28 33.625000 44.50000  3.9500000    1001.2500
## 514  2014-05-29 35.250000 40.12500  4.4000000     999.3750
## 515  2014-05-30 33.125000 44.50000 11.5875000     999.5714
## 516  2014-05-31 34.625000 39.00000  6.4875000     999.6250
## 517  2014-06-01 32.750000 34.25000 16.2000000    1003.2500
## 518  2014-06-02 32.750000 41.37500  4.6375000    1003.2500
## 519  2014-06-03 34.875000 33.62500  6.4875000    1000.0000
## 520  2014-06-04 35.250000 33.12500  7.1875000     998.1250
## 521  2014-06-05 36.875000 26.87500  8.5625000     996.1250
## 522  2014-06-06 37.500000 27.62500 10.2000000     994.8750
## 523  2014-06-07 38.500000 22.75000 13.1875000     994.0000
## 524  2014-06-08 37.625000 26.25000  6.7250000     993.1250
## 525  2014-06-09 37.875000 32.12500  6.0250000     993.1250
## 526  2014-06-10 37.250000 20.87500 12.0375000     994.7500
## 527  2014-06-11 37.625000 23.12500 14.8125000     995.6250
## 528  2014-06-12 32.875000 39.25000 24.3250000     997.0000
## 529  2014-06-13 30.250000 55.12500 16.6625000     997.6250
## 530  2014-06-14 30.500000 66.87500  6.4875000     997.3750
## 531  2014-06-15 33.875000 52.00000  8.1000000     997.1250
## 532  2014-06-16 36.375000 41.37500 10.8875000     995.6250
## 533  2014-06-17 36.875000 42.12500 15.9625000     994.1250
## 534  2014-06-18 33.875000 58.00000 10.4250000     994.2500
## 535  2014-06-19 35.750000 47.87500 10.8875000     993.2500
## 536  2014-06-20 38.000000 36.12500 17.8250000     991.3750
## 537  2014-06-21 36.875000 40.12500  9.7375000     993.8750
## 538  2014-06-22 34.500000 46.75000  4.4000000     996.3750
## 539  2014-06-23 32.125000 49.75000  6.7125000     996.2500
## 540  2014-06-24 33.000000 46.25000 10.1875000     995.5000
## 541  2014-06-25 32.250000 53.25000 13.2000000     997.6250
## 542  2014-06-26 32.500000 60.37500  3.2375000    1001.0000
## 543  2014-06-27 34.750000 48.50000  9.4875000    1001.6250
## 544  2014-06-28 34.875000 49.00000  5.7875000     999.3750
## 545  2014-06-29 33.000000 57.75000  9.7125000     997.7500
## 546  2014-06-30 32.000000 67.75000  3.2500000     997.5000
## 547  2014-07-01 31.375000 65.12500  6.2500000    1000.8750
## 548  2014-07-02 29.500000 71.37500  5.7875000    1001.1250
## 549  2014-07-03 28.875000 75.12500  5.0875000     998.2857
## 550  2014-07-04 31.388889 69.00000  4.3333333     998.7778
## 551  2014-07-05 34.705882 52.05882  9.2029412    1000.7059
## 552  2014-07-06 30.625000 61.12500  3.7125000    1000.5000
## 553  2014-07-07 34.000000 48.75000  5.5625000     998.2500
## 554  2014-07-08 35.000000 39.50000  9.7250000     996.8750
## 555  2014-07-09 36.000000 37.12500 13.6500000     996.5000
## 556  2014-07-10 36.125000 38.75000 11.3375000     996.8750
## 557  2014-07-11 36.625000 38.62500  8.7875000     996.1250
## 558  2014-07-12 35.875000 51.37500  5.3375000     996.0000
## 559  2014-07-13 31.625000 69.75000  6.0250000     998.0000
## 560  2014-07-14 32.375000 67.87500  6.7250000     997.1250
## 561  2014-07-15 33.125000 62.62500  8.8250000     996.0000
## 562  2014-07-16 32.625000 63.00000 10.4500000     997.2500
## 563  2014-07-17 30.750000 74.25000  6.7250000     996.1250
## 564  2014-07-18 27.250000 83.87500  2.5500000     994.8750
## 565  2014-07-19 29.875000 71.00000  6.2500000     995.1250
## 566  2014-07-20 30.875000 68.87500  6.2625000     996.7500
## 567  2014-07-21 32.125000 65.50000  4.6500000     996.0000
## 568  2014-07-22 31.625000 69.75000  7.1750000     997.0000
## 569  2014-07-23 30.625000 72.87500  8.1125000     997.7500
## 570  2014-07-24 30.500000 63.62500 12.0375000     998.7143
## 571  2014-07-25 30.750000 63.25000  9.9625000    1000.7500
## 572  2014-07-26 32.500000 61.50000  7.1750000    1001.8750
## 573  2014-07-27 31.714286 65.42857  6.2500000    1001.8750
## 574  2014-07-28 30.000000 75.75000 18.0500000     999.3750
## 575  2014-07-29 30.250000 69.87500  5.7875000     999.2500
## 576  2014-07-30 31.875000 65.75000  4.1750000     999.5714
## 577  2014-07-31 32.500000 65.75000  7.4250000     999.6250
## 578  2014-08-01 33.250000 59.50000  6.5000000     998.0000
## 579  2014-08-02 30.125000 77.75000  0.9250000     997.1250
## 580  2014-08-03 31.500000 66.25000  5.7875000     995.3750
## 581  2014-08-04 30.500000 69.75000  5.1000000     996.1250
## 582  2014-08-05 29.750000 77.50000  1.8500000     997.5000
## 583  2014-08-06 32.125000 65.00000  5.5625000     996.8750
## 584  2014-08-07 31.125000 65.37500  6.9500000     997.0000
## 585  2014-08-08 31.000000 67.87500  6.0125000     998.3750
## 586  2014-08-09 29.235294 80.47059  6.6529412    1001.6471
## 587  2014-08-10 28.500000 89.90625  5.2250000    1002.1562
## 588  2014-08-11 31.125000 72.12500  3.0125000     999.8750
## 589  2014-08-12 32.125000 57.50000  8.3375000     998.7500
## 590  2014-08-13 32.125000 57.71429 20.8250000     998.2500
## 591  2014-08-14 30.375000 57.50000 18.5125000     997.3750
## 592  2014-08-15 31.125000 54.87500 14.5875000     998.8750
## 593  2014-08-16 30.428571 59.85714  8.3500000    1002.5000
## 594  2014-08-17 31.375000 53.87500  7.1875000    1003.2500
## 595  2014-08-18 31.875000 56.62500  5.5625000    1005.3750
## 596  2014-08-19 33.000000 52.42857  6.8857143    1004.8571
## 597  2014-08-20 33.250000 47.12500  8.7875000    1002.3750
## 598  2014-08-21 32.500000 43.00000 10.8875000    1001.2500
## 599  2014-08-22 33.000000 47.12500  8.5625000    1002.8571
## 600  2014-08-23 33.125000 45.75000  9.2625000    1003.2500
## 601  2014-08-24 33.125000 46.25000 11.3500000    1002.8750
## 602  2014-08-25 33.875000 46.62500  7.8750000    1002.3750
## 603  2014-08-26 34.125000 40.25000  9.9625000    1001.3750
## 604  2014-08-27 33.857143 44.71429  8.1125000    1001.1250
## 605  2014-08-28 29.625000 67.75000  6.7250000    1003.7500
## 606  2014-08-29 30.375000 69.25000  3.7000000    1003.5000
## 607  2014-08-30 28.625000 78.12500  4.1625000    1002.6250
## 608  2014-08-31 27.625000 86.12500  4.4125000    1001.2500
## 609  2014-09-01 29.428571 65.85714 15.5125000    1001.1250
## 610  2014-09-02 29.250000 71.00000 10.4250000    1001.8750
## 611  2014-09-03 28.625000 77.87500 15.0375000    1001.4286
## 612  2014-09-04 27.875000 80.75000  9.2750000    1000.3750
## 613  2014-09-05 26.500000 85.75000  8.3375000    1000.1429
## 614  2014-09-06 28.500000 76.87500  2.0875000    1002.3750
## 615  2014-09-07 31.000000 67.62500  1.4000000    1003.2500
## 616  2014-09-08 31.250000 66.37500  3.4750000    1003.1250
## 617  2014-09-09 31.000000 70.12500  1.8500000    1002.3750
## 618  2014-09-10 28.750000 77.50000  4.8750000    1003.2500
## 619  2014-09-11 26.571429 85.71429 17.8250000    1005.5714
## 620  2014-09-12 28.500000 73.00000  1.6250000    1005.1250
## 621  2014-09-13 29.375000 68.25000  4.1625000    1004.0000
## 622  2014-09-14 29.125000 65.50000  6.7125000    1004.2500
## 623  2014-09-15 29.625000 63.12500  6.7125000    1004.2500
## 624  2014-09-16 30.625000 65.12500  6.2500000    1002.2500
## 625  2014-09-17 30.625000 61.87500  7.4000000    1002.1250
## 626  2014-09-18 30.500000 62.37500  6.7125000    1003.8750
## 627  2014-09-19 31.000000 56.62500  4.1750000    1003.5000
## 628  2014-09-20 31.125000 55.75000  4.4000000    1002.8750
## 629  2014-09-21 30.500000 54.87500  5.5625000    1003.6250
## 630  2014-09-22 30.625000 50.00000 13.4375000    1004.5000
## 631  2014-09-23 31.000000 47.00000 30.6857143    1004.2857
## 632  2014-09-24 29.625000 53.00000  8.1000000    1006.2500
## 633  2014-09-25 30.000000 52.50000  9.7250000    1007.7500
## 634  2014-09-26 30.000000 52.50000  5.7875000    1008.1250
## 635  2014-09-27 30.500000 46.87500  8.8000000    1009.2500
## 636  2014-09-28 30.250000 48.87500  7.1750000    1009.6250
## 637  2014-09-29 30.500000 48.75000  7.4125000    1009.1250
## 638  2014-09-30 30.625000 50.12500  6.7125000    1009.6250
## 639  2014-10-01 30.428571 52.85714  3.1714286    1010.4286
## 640  2014-10-02 30.375000 53.37500  2.7875000    1009.1250
## 641  2014-10-03 30.500000 52.62500  1.8500000    1008.0000
## 642  2014-10-04 31.000000 53.50000  0.0000000    1008.5000
## 643  2014-10-05 30.625000 56.12500  0.0000000    1010.3750
## 644  2014-10-06 30.625000 56.75000  0.0000000    1010.1250
## 645  2014-10-07 31.000000 59.37500  1.1625000    1008.3750
## 646  2014-10-08 28.750000 52.75000  3.7125000    1009.8750
## 647  2014-10-09 27.625000 49.87500  3.2375000    1010.1250
## 648  2014-10-10 27.125000 57.37500  3.2375000    1008.7500
## 649  2014-10-11 25.714286 57.14286  6.9375000    1008.2500
## 650  2014-10-12 26.000000 58.62500  0.4625000    1008.3750
## 651  2014-10-13 27.625000 59.75000  8.8000000    1009.8750
## 652  2014-10-14 24.500000 67.62500  5.3250000    1010.8750
## 653  2014-10-15 24.500000 68.00000  3.2625000    1013.3750
## 654  2014-10-16 23.000000 58.87500  4.8625000    1014.0000
## 655  2014-10-17 25.000000 47.66667 25.0125000    1013.6250
## 656  2014-10-18 24.000000 56.14286 21.3000000    1012.8750
## 657  2014-10-19 23.750000 57.75000  4.1625000    1012.6250
## 658  2014-10-20 24.750000 55.62500  4.4000000    1013.1250
## 659  2014-10-21 25.625000 56.12500  2.0875000    1014.6250
## 660  2014-10-22 25.750000 64.87500  0.4625000    1014.5000
## 661  2014-10-23 25.250000 66.00000  0.9250000    1012.7500
## 662  2014-10-24 25.500000 67.50000  1.8500000    1011.7500
## 663  2014-10-25 25.250000 67.12500  0.9250000    1012.3750
## 664  2014-10-26 25.125000 65.50000  2.0875000    1014.8750
## 665  2014-10-27 24.875000 62.87500  1.8625000    1013.7500
## 666  2014-10-28 25.125000 64.75000  2.5625000    1012.0000
## 667  2014-10-29 24.750000 65.75000  1.6250000    1012.1250
## 668  2014-10-30 23.875000 61.37500  1.8500000    1014.5000
## 669  2014-10-31 23.375000 58.25000  3.7125000    1013.6250
## 670  2014-11-01 22.500000 50.87500  5.3250000    1014.1250
## 671  2014-11-02 22.875000 49.12500  7.1875000    1014.5000
## 672  2014-11-03 23.375000 52.75000  4.1625000    1013.5000
## 673  2014-11-04 23.625000 56.62500  3.9375000    1012.1250
## 674  2014-11-05 22.750000 56.25000  4.6250000    1012.3750
## 675  2014-11-06 23.000000 56.25000  6.2500000    1011.6250
## 676  2014-11-07 24.000000 66.12500  2.3125000    1012.8750
## 677  2014-11-08 23.875000 66.25000  1.6250000    1014.6250
## 678  2014-11-09 24.125000 56.37500  4.1625000    1015.1250
## 679  2014-11-10 21.375000 52.12500  3.4750000    1014.0000
## 680  2014-11-11 19.250000 55.12500  2.3250000    1014.7500
## 681  2014-11-12 19.625000 52.00000  2.3125000    1015.8750
## 682  2014-11-13 18.500000 50.50000  3.2375000    1016.5000
## 683  2014-11-14 18.375000 43.00000  7.4125000    1016.3750
## 684  2014-11-15 18.625000 44.75000  2.7750000    1014.8750
## 685  2014-11-16 18.250000 47.62500  2.5500000    1015.3750
## 686  2014-11-17 17.750000 51.75000  2.7750000    1016.1250
## 687  2014-11-18 17.625000 53.00000  2.0875000    1017.8750
## 688  2014-11-19 17.500000 58.12500  1.1625000    1016.5000
## 689  2014-11-20 20.375000 54.62500  3.0250000    1013.3750
## 690  2014-11-21 18.250000 62.12500  3.2375000    1012.7500
## 691  2014-11-22 19.875000 58.87500  1.6250000    1011.3750
## 692  2014-11-23 17.125000 57.37500  1.3875000    1014.6250
## 693  2014-11-24 16.875000 52.25000  4.1750000    1016.0000
## 694  2014-11-25 16.500000 51.50000  2.5500000    1014.8750
## 695  2014-11-26 19.625000 51.00000  5.5500000    1009.8750
## 696  2014-11-27 18.125000 63.87500  0.0000000    1011.7500
## 697  2014-11-28 19.375000 67.37500  0.0000000    1012.5000
## 698  2014-11-29 20.250000 67.62500  0.0000000    1012.3750
## 699  2014-11-30 19.750000 54.75000  2.5500000    1014.0000
## 700  2014-12-01 20.000000 42.62500  8.8000000    1013.7500
## 701  2014-12-02 21.375000 47.62500  7.8750000    1012.7500
## 702  2014-12-03 25.000000 66.57143 11.3714286    1007.0000
## 703  2014-12-04 21.250000 68.87500  4.8625000    1010.3750
## 704  2014-12-05 17.750000 59.87500  2.5500000    1017.1250
## 705  2014-12-06 18.250000 48.87500  4.8625000    1016.3750
## 706  2014-12-07 17.250000 48.37500  4.8625000    1014.7500
## 707  2014-12-08 17.000000 49.87500  2.7750000    1015.0000
## 708  2014-12-09 15.285714 67.57143  0.0000000    1016.8571
## 709  2014-12-10 21.500000 65.50000  4.8625000    1010.5000
## 710  2014-12-11 19.625000 70.50000  3.2500000    1011.6250
## 711  2014-12-12 17.500000 65.62500  0.4625000    1014.3750
## 712  2014-12-13 16.000000 71.87500  8.8000000    1016.1250
## 713  2014-12-14 16.875000 86.62500  3.9500000    1014.5000
## 714  2014-12-15 16.625000 75.75000  2.3250000    1016.1250
## 715  2014-12-16 17.750000 72.37500  7.4250000    1017.3750
## 716  2014-12-17 12.875000 72.12500  6.4875000    1020.2500
## 717  2014-12-18 11.750000 79.75000  3.0125000    1021.3750
## 718  2014-12-19 11.750000 81.75000  0.4625000    1022.1250
## 719  2014-12-20 11.250000 77.00000  5.8000000    1018.7500
## 720  2014-12-21  9.500000 79.12500  0.9250000    1018.0000
## 721  2014-12-22  9.375000 85.00000  1.3875000    1017.1250
## 722  2014-12-23  9.625000 78.62500  1.6250000    1017.8750
## 723  2014-12-24  9.875000 83.00000  3.4750000    1021.2500
## 724  2014-12-25  9.250000 81.12500  8.1000000    1020.1250
## 725  2014-12-26 10.750000 77.00000  7.6375000    1017.7500
## 726  2014-12-27 10.375000 69.00000  2.7750000    1018.6250
## 727  2014-12-28  9.000000 86.00000  0.7000000    1019.7500
## 728  2014-12-29 11.125000 72.62500  1.3875000    1017.2500
## 729  2014-12-30 11.625000 70.62500  2.5500000    1014.6250
## 730  2014-12-31 12.375000 67.12500  2.7875000    1016.8750
## 731  2015-01-01 14.750000 72.00000  0.9250000    1017.5000
## 732  2015-01-02 14.875000 96.62500  3.0125000    1017.8750
## 733  2015-01-03 15.125000 92.00000  0.9250000    1017.3750
## 734  2015-01-04 14.125000 78.75000  9.5125000    1019.6250
## 735  2015-01-05 14.000000 69.37500 15.0500000    1016.0000
## 736  2015-01-06 12.000000 79.87500  4.4000000    1015.5000
## 737  2015-01-07  9.625000 86.00000  3.9375000    1016.6250
## 738  2015-01-08 10.000000 87.87500  3.4875000    1019.1250
## 739  2015-01-09 10.625000 80.62500  1.6250000    1020.1250
## 740  2015-01-10 11.125000 73.12500  4.4125000    1018.8750
## 741  2015-01-11 11.000000 80.75000  4.8625000    1019.1250
## 742  2015-01-12 10.625000 79.50000  0.9250000    1019.5000
## 743  2015-01-13 12.250000 77.87500  4.8750000    1017.1250
## 744  2015-01-14 12.000000 82.25000  5.1000000    1017.5000
## 745  2015-01-15 12.375000 79.12500  0.9250000    1019.3750
## 746  2015-01-16 13.000000 78.12500  4.6375000    1019.2500
## 747  2015-01-17 13.500000 76.62500  1.8625000    1019.2500
## 748  2015-01-18 13.000000 73.87500  2.5500000    1018.7500
## 749  2015-01-19 13.750000 78.00000  1.8500000    1018.8750
## 750  2015-01-20 13.375000 74.37500  7.1875000    1020.0000
## 751  2015-01-21 13.250000 73.37500  3.7000000    1021.0000
## 752  2015-01-22 12.625000 94.12500  3.2375000    1020.6250
## 753  2015-01-23 14.625000 85.00000  1.1625000    1020.0000
## 754  2015-01-24 12.375000 86.25000  9.7375000    1020.7500
## 755  2015-01-25 13.250000 77.75000  4.1625000    1020.3750
## 756  2015-01-26 13.750000 88.12500  0.9250000    1019.8750
## 757  2015-01-27 13.250000 74.25000  5.3375000    1019.3333
## 758  2015-01-28 11.714286 67.00000 27.7750000    1018.2500
## 759  2015-01-29 11.500000 67.50000  3.9375000    1019.7500
## 760  2015-01-30 12.750000 56.12500 12.0375000    1020.3750
## 761  2015-01-31 13.750000 57.62500  9.0250000    1017.8750
## 762  2015-02-01 15.000000 65.25000  5.1000000    1017.7500
## 763  2015-02-02 17.375000 63.87500 11.8125000    1017.5000
## 764  2015-02-03 16.750000 71.25000  7.1750000    1020.5000
## 765  2015-02-04 14.625000 81.12500  9.2500000    1022.0000
## 766  2015-02-05 14.625000 62.87500  7.4000000    1019.2500
## 767  2015-02-06 15.125000 62.62500  3.0125000    1019.0000
## 768  2015-02-07 15.875000 67.62500  1.8500000    1017.0000
## 769  2015-02-08 16.625000 63.87500  4.4125000    1017.3750
## 770  2015-02-09 16.875000 63.50000  5.5500000    1019.7500
## 771  2015-02-10 15.750000 64.00000  4.6250000    1019.7500
## 772  2015-02-11 16.000000 66.62500  1.8500000    1017.1250
## 773  2015-02-12 15.500000 69.62500  6.0250000    1016.1250
## 774  2015-02-13 17.250000 59.00000  6.7250000    1012.6250
## 775  2015-02-14 19.375000 63.12500  4.1625000    1009.6250
## 776  2015-02-15 18.875000 70.12500  1.3875000    1012.6250
## 777  2015-02-16 20.625000 71.25000 10.4250000    1013.6250
## 778  2015-02-17 21.142857 72.57143  5.0285714    1011.8571
## 779  2015-02-18 21.250000 72.87500  5.1000000    1011.7500
## 780  2015-02-19 22.750000 73.50000  7.6500000    1012.0000
## 781  2015-02-20 22.750000 67.12500  7.1875000    1013.8750
## 782  2015-02-21 21.875000 63.50000 11.3375000    1014.1250
## 783  2015-02-22 22.142857 58.85714 14.8125000    1013.3750
## 784  2015-02-23 21.500000 74.37500  3.0250000    1014.2500
## 785  2015-02-24 22.875000 72.37500  6.0375000    1010.6250
## 786  2015-02-25 23.125000 58.62500  6.5000000    1008.5000
## 787  2015-02-26 20.000000 34.37500  6.7125000    1009.8750
## 788  2015-02-27 19.250000 43.12500  5.8000000    1013.0000
## 789  2015-02-28 21.250000 45.50000  6.9500000    1013.0000
## 790  2015-03-01 17.375000 92.75000 14.8125000    1011.1250
## 791  2015-03-02 17.500000 89.87500  6.7250000    1012.5000
## 792  2015-03-03 16.250000 82.50000  5.5625000    1013.7500
## 793  2015-03-04 16.750000 69.25000  6.7125000    1014.0000
## 794  2015-03-05 18.750000 64.00000  3.7000000    1015.7500
## 795  2015-03-06 20.000000 62.37500  5.0875000    1016.3750
## 796  2015-03-07 18.375000 80.37500  4.1750000    1013.7143
## 797  2015-03-08 19.500000 74.12500  7.8875000    1015.3750
## 798  2015-03-09 18.125000 54.25000 12.2625000    1021.0000
## 799  2015-03-10 17.750000 68.50000  8.5625000    1019.1250
## 800  2015-03-11 18.250000 69.87500  6.0250000    1018.1250
## 801  2015-03-12 20.500000 63.37500  5.3250000    1017.0000
## 802  2015-03-13 21.375000 65.25000  4.1750000    1014.7500
## 803  2015-03-14 21.500000 64.62500  7.8750000    1015.1250
## 804  2015-03-15 17.625000 89.12500 11.5750000    1015.5000
## 805  2015-03-16 19.857143 81.85714  8.5625000    1017.3750
## 806  2015-03-17 20.875000 70.37500 13.2000000    1015.8750
## 807  2015-03-18 20.250000 69.62500 12.2625000    1013.5000
## 808  2015-03-19 21.125000 66.00000  5.3375000    1014.8750
## 809  2015-03-20 22.500000 65.75000  9.0250000    1014.8750
## 810  2015-03-21 24.875000 60.50000  5.7875000    1014.8750
## 811  2015-03-22 25.125000 64.00000 10.2000000    1012.0000
## 812  2015-03-23 24.750000 65.75000  7.6375000    1010.0000
## 813  2015-03-24 25.750000 65.50000  1.3875000    1010.2500
## 814  2015-03-25 26.625000 65.62500  2.5500000    1011.8571
## 815  2015-03-26 28.500000 54.37500  9.2625000    1013.0000
## 816  2015-03-27 27.375000 51.37500  7.6500000    1012.0000
## 817  2015-03-28 27.000000 53.62500  3.7125000    1009.6250
## 818  2015-03-29 25.000000 52.75000 13.1875000    1008.6250
## 819  2015-03-30 23.500000 66.25000 14.8125000    1010.3750
## 820  2015-03-31 25.250000 64.87500  3.7000000    1010.1429
## 821  2015-04-01 26.428571 59.71429  6.7125000    1007.5000
## 822  2015-04-02 28.714286 45.28571 10.1750000    1004.6250
## 823  2015-04-03 24.500000 71.12500  6.9625000    1004.7500
## 824  2015-04-04 21.000000 77.50000  1.3875000    1005.2500
## 825  2015-04-05 23.571429 67.71429 17.6000000    1007.7500
## 826  2015-04-06 25.000000 60.12500  2.7750000    1007.7500
## 827  2015-04-07 23.750000 62.12500  9.7375000    1010.2500
## 828  2015-04-08 25.625000 50.37500  6.7125000    1012.0000
## 829  2015-04-09 26.000000 52.87500  9.0375000    1010.5000
## 830  2015-04-10 27.125000 52.25000  6.2500000    1008.5000
## 831  2015-04-11 28.250000 44.87500  8.1000000    1008.3750
## 832  2015-04-12 24.375000 63.00000  5.3250000    1012.3750
## 833  2015-04-13 24.750000 70.87500  6.9500000    1013.7500
## 834  2015-04-14 26.250000 62.75000  3.0125000    1012.7500
## 835  2015-04-15 27.625000 59.25000  5.7875000    1012.8750
## 836  2015-04-16 25.250000 58.25000  8.1125000    1013.8750
## 837  2015-04-17 28.125000 48.75000  8.3375000    1007.1250
## 838  2015-04-18 30.125000 44.50000  5.5625000    1002.3750
## 839  2015-04-19 32.125000 41.00000  4.1750000    1003.8750
## 840  2015-04-20 32.875000 37.87500  7.1875000    1005.7500
## 841  2015-04-21 31.750000 28.50000 11.1125000    1005.6250
## 842  2015-04-22 30.750000 30.25000 11.3500000    1003.1250
## 843  2015-04-23 30.750000 33.00000  8.1000000    1004.0000
## 844  2015-04-24 31.125000 39.50000  2.5500000    1004.3750
## 845  2015-04-25 30.500000 43.12500 12.9750000    1005.8750
## 846  2015-04-26 28.750000 51.12500 15.7375000    1008.0000
## 847  2015-04-27 30.125000 53.00000 17.5875000    1006.7500
## 848  2015-04-28 31.000000 48.25000 12.5125000    1004.5000
## 849  2015-04-29 31.250000 47.25000  9.0250000    1005.0000
## 850  2015-04-30 32.625000 44.75000 15.0500000    1004.0000
## 851  2015-05-01 31.250000 43.50000  5.3375000    1005.8750
## 852  2015-05-02 30.125000 40.25000  4.4000000    1006.2500
## 853  2015-05-03 31.000000 39.25000  3.0125000    1006.6250
## 854  2015-05-04 30.875000 39.25000  3.9375000    1006.6250
## 855  2015-05-05 32.125000 31.12500  7.6500000    1003.8750
## 856  2015-05-06 34.125000 22.00000 11.3625000    1001.8750
## 857  2015-05-07 35.428571 13.42857 15.3428571    1003.5714
## 858  2015-05-08 34.125000 22.37500 12.2750000    1003.2500
## 859  2015-05-09 35.125000 27.00000  8.3375000    1001.6250
## 860  2015-05-10 35.125000 25.00000  6.9625000    1001.3750
## 861  2015-05-11 34.000000 32.62500 10.8875000    1003.6250
## 862  2015-05-12 33.125000 46.00000  6.0375000    1004.8750
## 863  2015-05-13 27.750000 61.00000  6.7125000    1006.1250
## 864  2015-05-14 29.625000 49.62500  4.4125000    1007.2500
## 865  2015-05-15 28.500000 52.12500  7.1750000    1006.5000
## 866  2015-05-16 30.000000 51.62500  4.1625000    1003.2500
## 867  2015-05-17 33.250000 40.37500  8.3500000    1001.1250
## 868  2015-05-18 35.375000 35.75000  5.5625000    1000.0000
## 869  2015-05-19 31.500000 44.00000 14.8125000    1002.1250
## 870  2015-05-20 33.125000 31.62500 14.6000000    1002.0000
## 871  2015-05-21 35.625000 24.62500 11.3375000    1000.3750
## 872  2015-05-22 36.250000 19.50000 12.5125000     998.1250
## 873  2015-05-23 37.375000 21.75000 12.7250000     996.1250
## 874  2015-05-24 36.000000 29.25000  5.1000000     997.3750
## 875  2015-05-25 37.750000 22.25000 10.8750000     996.8750
## 876  2015-05-26 35.250000 26.62500 17.1250000    1001.5000
## 877  2015-05-27 33.750000 20.50000  6.7125000    1001.2500
## 878  2015-05-28 34.750000 29.50000  8.3375000    1000.6250
## 879  2015-05-29 35.125000 27.25000  8.1250000    1001.6250
## 880  2015-05-30 32.250000 38.75000 12.2625000    1003.8750
## 881  2015-05-31 34.250000 25.00000 14.5875000    1001.7500
## 882  2015-06-01 29.375000 41.75000 12.5125000    1003.2500
## 883  2015-06-02 31.625000 39.25000  8.5750000    1003.5000
## 884  2015-06-03 28.625000 50.50000  6.7250000    1004.6250
## 885  2015-06-04 31.875000 48.37500  4.8625000    1004.0000
## 886  2015-06-05 32.625000 39.87500  6.0125000    1003.5000
## 887  2015-06-06 33.375000 39.75000  7.8625000    1001.3750
## 888  2015-06-07 35.000000 36.25000  7.1750000     999.2500
## 889  2015-06-08 36.500000 33.00000  6.2500000     998.8750
## 890  2015-06-09 37.625000 23.87500 18.7500000     998.8750
## 891  2015-06-10 36.625000 32.28571 18.0500000     999.6250
## 892  2015-06-11 35.625000 34.75000 18.2875000    1000.0000
## 893  2015-06-12 35.875000 34.62500 12.7375000     998.5000
## 894  2015-06-13 31.125000 48.12500 10.1875000     998.8571
## 895  2015-06-14 27.000000 77.62500  3.9500000    1000.2857
## 896  2015-06-15 31.000000 60.75000  3.0125000     999.5714
## 897  2015-06-16 32.375000 52.12500  2.7875000     999.1429
## 898  2015-06-17 33.875000 49.75000  3.2500000     998.0000
## 899  2015-06-18 35.125000 45.12500  3.0125000     997.6250
## 900  2015-06-19 35.875000 46.75000  8.5750000     996.2500
## 901  2015-06-20 33.125000 48.50000 12.9750000     997.8750
## 902  2015-06-21 30.375000 63.62500  8.1125000     998.1250
## 903  2015-06-22 32.125000 59.50000  4.4000000     995.6250
## 904  2015-06-23 31.500000 61.62500  6.0125000     995.0000
## 905  2015-06-24 37.250000 63.87500  6.5000000     995.3750
## 906  2015-06-25 26.625000 82.75000  6.2625000     994.7500
## 907  2015-06-26 31.125000 54.25000 11.5625000     995.2500
## 908  2015-06-27 32.625000 48.25000 14.1125000     997.1250
## 909  2015-06-28 34.250000 47.50000  9.0375000    1000.2857
## 910  2015-06-29 31.125000 56.62500 11.5625000    1003.6250
## 911  2015-06-30 30.875000 58.62500  9.9625000    1002.6250
## 912  2015-07-01 32.000000 53.25000  4.1750000    1001.8750
## 913  2015-07-02 33.750000 45.87500  6.9500000    1001.1250
## 914  2015-07-03 33.000000 49.00000 14.3500000     999.6250
## 915  2015-07-04 34.125000 51.50000 16.9000000     997.5000
## 916  2015-07-05 35.714286 47.28571 14.0142857     995.4286
## 917  2015-07-06 27.875000 83.25000  9.2750000     997.8750
## 918  2015-07-07 28.125000 86.12500  1.8500000     998.1250
## 919  2015-07-08 31.250000 72.00000  2.7750000     998.6250
## 920  2015-07-09 29.625000 77.25000  2.5500000     998.6250
## 921  2015-07-10 26.625000 92.12500  6.0250000     998.7500
## 922  2015-07-11 25.500000 98.00000  3.2375000     999.3333
## 923  2015-07-12 25.500000 94.75000  7.4250000     998.6250
## 924  2015-07-13 29.875000 73.12500  4.8625000    1000.6250
## 925  2015-07-14 30.875000 67.87500 10.8750000     997.8571
## 926  2015-07-15 33.375000 52.87500 17.1375000     993.7500
## 927  2015-07-16 31.750000 61.37500 13.4375000     994.7500
## 928  2015-07-17 30.750000 75.62500  4.4000000     997.2500
## 929  2015-07-18 29.500000 78.87500  4.1750000     998.0000
## 930  2015-07-19 30.125000 70.00000  0.7000000     997.6250
## 931  2015-07-20 29.125000 78.00000  4.8625000     999.1250
## 932  2015-07-21 31.375000 69.25000  8.1250000    1000.0000
## 933  2015-07-22 32.250000 64.12500  7.8750000    1001.1250
## 934  2015-07-23 32.125000 64.87500  8.5750000    1001.7500
## 935  2015-07-24 32.000000 64.75000  6.7125000    1001.0000
## 936  2015-07-25 30.857143 71.00000 19.9125000     998.5714
## 937  2015-07-26 30.000000 71.62500 16.8875000     998.8750
## 938  2015-07-27 29.625000 69.00000 11.3375000    1000.2857
## 939  2015-07-28 29.500000 61.00000 11.5875000    1000.1250
## 940  2015-07-29 28.500000 68.12500 10.8750000     998.0000
## 941  2015-07-30 29.000000 67.75000  9.9625000     998.8750
## 942  2015-07-31 29.000000 71.00000  4.1750000     999.7500
## 943  2015-08-01 28.250000 71.75000  0.9250000    1002.7500
## 944  2015-08-02 29.625000 69.50000  2.5500000    1003.6250
## 945  2015-08-03 31.500000 58.62500  5.3250000    1001.8750
## 946  2015-08-04 29.428571 77.78571  4.7714286    1000.0714
## 947  2015-08-05 29.593750 79.65625  8.2218750    1000.2812
## 948  2015-08-06 29.500000 79.12500  3.0250000     998.0000
## 949  2015-08-07 29.250000 83.00000  2.3125000    1000.0000
## 950  2015-08-08 28.875000 84.87500  1.6250000    1002.0000
## 951  2015-08-09 27.375000 90.12500  3.0125000    1000.2500
## 952  2015-08-10 30.375000 77.00000  0.0000000    1000.0000
## 953  2015-08-11 28.375000 86.62500  0.4625000    1003.0000
## 954  2015-08-12 30.750000 74.25000  0.0000000    1003.1250
## 955  2015-08-13 32.000000 69.50000  0.4625000    1001.8750
## 956  2015-08-14 30.875000 72.62500  3.4750000    1002.0000
## 957  2015-08-15 28.500000 85.50000  0.0000000    1002.7500
## 958  2015-08-16 28.625000 82.62500  1.8500000    1002.8750
## 959  2015-08-17 30.750000 72.62500  0.7000000    1001.3750
## 960  2015-08-18 29.125000 72.87500  4.8750000    1002.5714
## 961  2015-08-19 31.250000 61.25000 33.3250000    1002.3750
## 962  2015-08-20 30.750000 58.62500 14.6000000    1001.3750
## 963  2015-08-21 29.375000 72.00000  8.5750000    1002.0000
## 964  2015-08-22 30.250000 65.37500  5.8000000    1000.8750
## 965  2015-08-23 29.875000 69.37500  5.5625000    1000.8571
## 966  2015-08-24 30.750000 65.62500  7.1875000    1001.7500
## 967  2015-08-25 32.000000 57.28571 10.0714286    1001.0000
## 968  2015-08-26 32.000000 61.00000  4.8625000    1001.1250
## 969  2015-08-27 32.500000 61.00000  5.3250000    1002.0000
## 970  2015-08-28 32.000000 63.50000  2.3125000    1002.2500
## 971  2015-08-29 32.250000 58.00000  5.8125000    1001.3750
## 972  2015-08-30 32.625000 56.50000  9.2750000    1001.3750
## 973  2015-08-31 31.857143 54.42857 13.4857143    1001.8333
## 974  2015-09-01 32.125000 47.62500 15.2750000    1000.8571
## 975  2015-09-02 31.875000 45.62500 12.7375000    1001.7500
## 976  2015-09-03 31.625000 49.12500 10.1875000    1005.5000
## 977  2015-09-04 32.000000 42.62500 11.3500000    1007.8750
## 978  2015-09-05 30.875000 47.25000  7.6375000    1007.1250
## 979  2015-09-06 31.250000 41.80000  7.6375000    1006.2500
## 980  2015-09-07 31.500000 46.37500  9.4875000    1005.7500
## 981  2015-09-08 31.625000 47.00000 11.8125000    1005.5000
## 982  2015-09-09 31.000000 47.87500 10.1875000    1005.6250
## 983  2015-09-10 31.250000 46.50000  8.5750000    1006.7500
## 984  2015-09-11 30.750000 49.75000  6.0125000    1006.2500
## 985  2015-09-12 32.142857 54.42857  3.1714286    1006.0000
## 986  2015-09-13 32.125000 55.25000  1.3875000    1005.8750
## 987  2015-09-14 31.875000 54.37500  2.5500000    1003.6250
## 988  2015-09-15 33.000000 55.00000  0.4625000    1003.2500
## 989  2015-09-16 32.375000 55.00000  1.4000000    1003.6250
## 990  2015-09-17 31.375000 58.37500  1.1625000    1004.3750
## 991  2015-09-18 31.000000 67.62500  7.4125000    1005.1250
## 992  2015-09-19 30.625000 74.25000  2.0875000    1003.6250
## 993  2015-09-20 29.142857 69.14286  5.5571429    1002.8571
## 994  2015-09-21 30.125000 64.50000  5.8000000    1002.5000
## 995  2015-09-22 30.250000 64.37500 10.1875000    1002.5000
## 996  2015-09-23 28.625000 61.25000 13.6750000    1006.3750
## 997  2015-09-24 28.250000 52.37500 11.5750000    1004.6250
## 998  2015-09-25 28.500000 49.37500  8.5750000    1003.7500
## 999  2015-09-26 29.250000 53.50000  6.2500000    1005.2500
## 1000 2015-09-27 28.714286 58.14286  2.9142857    1007.0000
## 1001 2015-09-28 28.250000 58.50000  3.9375000    1009.0000
## 1002 2015-09-29 29.000000 54.57143  9.5285714    1010.5714
## 1003 2015-09-30 29.000000 46.87500  8.3375000    1011.3750
## 1004 2015-10-01 28.125000 51.25000  5.5750000    1012.0000
## 1005 2015-10-02 28.285714 53.71429 14.5875000    1010.7500
## 1006 2015-10-03 30.000000 52.33333 18.7857143    1008.7143
## 1007 2015-10-04 29.000000 54.37500  3.0250000    1008.1250
## 1008 2015-10-05 28.875000 59.37500  1.3875000    1009.2500
## 1009 2015-10-06 28.375000 61.12500  4.6375000    1009.0000
## 1010 2015-10-07 28.625000 65.75000  6.4750000    1009.1250
## 1011 2015-10-08 28.375000 65.12500  4.4125000    1007.5000
## 1012 2015-10-09 28.500000 59.75000  6.0250000    1004.2500
## 1013 2015-10-10 28.500000 54.62500  6.7125000    1005.0000
## 1014 2015-10-11 28.571429 51.57143  6.8857143    1008.7143
## 1015 2015-10-12 29.250000 54.12500  3.2500000    1010.0000
## 1016 2015-10-13 29.500000 62.87500  3.7250000    1011.0000
## 1017 2015-10-14 28.125000 63.62500  4.8750000    1013.2500
## 1018 2015-10-15 27.500000 66.00000  0.4625000    1013.5000
## 1019 2015-10-16 27.875000 62.62500  0.9250000    1013.5000
## 1020 2015-10-17 28.000000 63.75000  3.7125000    1012.2500
## 1021 2015-10-18 26.857143 71.71429 11.1125000    1013.1250
## 1022 2015-10-19 28.142857 66.28571  3.4571429    1011.7143
## 1023 2015-10-20 27.714286 64.85714  2.9142857    1009.8571
## 1024 2015-10-21 26.500000 62.12500  2.0875000    1011.8750
## 1025 2015-10-22 24.750000 56.00000  4.4000000    1012.2500
## 1026 2015-10-23 24.000000 51.62500  4.8625000    1011.8750
## 1027 2015-10-24 24.875000 53.25000  1.4000000    1011.6250
## 1028 2015-10-25 26.375000 47.12500  5.5750000    1012.7500
## 1029 2015-10-26 25.625000 52.12500  5.5750000    1014.8750
## 1030 2015-10-27 24.750000 54.25000  3.4750000    1013.8750
## 1031 2015-10-28 20.875000 57.25000  7.8750000    1015.8750
## 1032 2015-10-29 22.125000 63.50000  1.8500000    1016.8750
## 1033 2015-10-30 22.375000 67.25000  1.6250000    1016.3750
## 1034 2015-10-31 22.375000 69.37500  1.1625000    1016.1250
## 1035 2015-11-01 22.875000 70.50000  0.9250000    1016.2500
## 1036 2015-11-02 23.375000 71.37500  0.4625000    1015.8750
## 1037 2015-11-03 24.000000 67.87500  0.4750000    1014.2500
## 1038 2015-11-04 23.750000 69.00000  0.9250000    1012.2500
## 1039 2015-11-05 22.125000 75.25000  5.1000000    1012.3750
## 1040 2015-11-06 22.875000 74.37500  2.3250000    1012.7500
## 1041 2015-11-07 21.875000 74.87500  1.1625000    1014.6250
## 1042 2015-11-08 21.875000 70.25000  1.8625000    1015.8750
## 1043 2015-11-09 21.750000 73.62500  0.7000000    1015.0000
## 1044 2015-11-10 22.750000 72.62500  0.7000000    1013.5000
## 1045 2015-11-11 21.875000 59.50000  4.6375000    1014.3750
## 1046 2015-11-12 21.875000 60.50000  5.7875000    1015.3750
## 1047 2015-11-13 22.375000 59.37500  5.3250000    1014.8750
## 1048 2015-11-14 22.375000 59.62500  3.4750000    1012.0000
## 1049 2015-11-15 22.500000 62.00000  0.0000000    1012.0000
## 1050 2015-11-16 21.875000 50.37500  2.7875000    1012.6250
## 1051 2015-11-17 19.875000 48.37500  5.1000000    1013.5000
## 1052 2015-11-18 18.500000 57.62500  4.4000000    1016.3750
## 1053 2015-11-19 18.500000 64.75000  0.9250000    1017.5000
## 1054 2015-11-20 18.500000 66.00000  0.4625000    1016.1250
## 1055 2015-11-21 18.750000 67.62500  1.3875000    1015.5000
## 1056 2015-11-22 18.625000 64.87500  0.4625000    1015.1250
## 1057 2015-11-23 17.000000 73.50000  0.0000000    1015.8750
## 1058 2015-11-24 18.500000 63.00000  1.6250000    1016.0000
## 1059 2015-11-25 19.375000 56.37500  5.7875000    1014.7500
## 1060 2015-11-26 20.000000 59.50000  0.9250000    1012.7500
## 1061 2015-11-27 18.500000 66.25000  2.5500000    1011.2500
## 1062 2015-11-28 18.750000 82.12500  1.1625000    1015.3750
## 1063 2015-11-29 18.125000 81.25000  1.1625000    1018.0000
## 1064 2015-11-30 19.500000 82.37500  0.9250000    1016.7500
## 1065 2015-12-01 19.250000 78.75000  1.6250000    1014.7500
## 1066 2015-12-02 18.125000 72.12500  2.0875000    1016.1250
## 1067 2015-12-03 17.625000 70.25000  7.1875000    1016.0000
## 1068 2015-12-04 17.625000 67.50000  4.6375000    1015.0000
## 1069 2015-12-05 17.500000 72.37500  0.9250000    1015.8750
## 1070 2015-12-06 17.125000 75.87500  1.4000000    1016.5000
## 1071 2015-12-07 17.625000 76.87500  2.3125000    1017.2500
## 1072 2015-12-08 18.428571 71.28571  1.0571429    1015.7143
## 1073 2015-12-09 18.375000 78.12500  3.4750000    1015.2500
## 1074 2015-12-10 19.000000 77.25000  2.0875000    1011.7500
## 1075 2015-12-11 17.250000 77.87500  5.1000000    1011.3750
## 1076 2015-12-12 13.250000 91.00000  4.4000000    1014.2500
## 1077 2015-12-13 12.875000 67.00000  3.4750000    1016.5000
## 1078 2015-12-14 13.250000 63.62500  7.4000000    1017.0000
## 1079 2015-12-15 13.125000 63.87500  5.3250000    1016.5000
## 1080 2015-12-16 13.000000 65.50000  1.1625000    1016.3750
## 1081 2015-12-17 12.250000 73.62500  0.9250000    1016.8750
## 1082 2015-12-18 12.875000 71.25000  4.4000000    1018.6250
## 1083 2015-12-19 12.375000 71.87500  3.4875000    1019.0000
## 1084 2015-12-20 11.750000 70.25000  1.1625000    1018.5000
## 1085 2015-12-21 12.625000 74.75000  1.3875000    1018.8750
## 1086 2015-12-22 12.000000 79.87500  0.7000000    1020.8750
## 1087 2015-12-23 12.500000 81.37500  3.4875000    1020.6250
## 1088 2015-12-24 11.250000 68.62500  6.0250000    1021.2500
## 1089 2015-12-25 11.500000 66.12500  5.5625000    1022.0000
## 1090 2015-12-26 12.750000 64.50000  6.4875000    1021.1250
## 1091 2015-12-27 15.375000 63.25000  7.8875000    1020.6250
## 1092 2015-12-28 17.125000 58.12500 10.8875000    1020.8750
## 1093 2015-12-29 16.375000 65.00000  7.4125000    1018.1250
## 1094 2015-12-30 15.500000 71.75000  2.1000000    1017.5000
## 1095 2015-12-31 15.000000 71.37500  2.0875000    1020.5000

Peubah yang dijadikan sebagi peubah bebas (xt) adalah meantemp (rata-rata temperatur udara) dan peubah responnya (yt) adalah meanpressure (rata-rata tekanan udara)

t <- data$date
t <- tibble::num(1:1095)
xt <- data$meantemp
yt <- data$meanpressure
xt <- as.vector(xt)
yt <- as.vector(yt)
data1 <- cbind(t, xt, yt)
data1 <- as.data.frame(data1)
str(data1)
## 'data.frame':    1095 obs. of  3 variables:
##  $ t : num  1 2 3 4 5 6 7 8 9 10 ...
##  $ xt: num  10 7.4 7.17 8.67 6 ...
##  $ yt: num  1016 1018 1019 1017 1016 ...
data1
##         t        xt        yt
## 1       1 10.000000 1015.6667
## 2       2  7.400000 1017.8000
## 3       3  7.166667 1018.6667
## 4       4  8.666667 1017.1667
## 5       5  6.000000 1016.5000
## 6       6  7.000000 1018.0000
## 7       7  7.000000 1020.0000
## 8       8  8.857143 1018.7143
## 9       9 14.000000 1017.0000
## 10     10 11.000000 1015.6667
## 11     11 15.714286 1016.1429
## 12     12 14.000000 1015.5714
## 13     13 15.833333 1013.3333
## 14     14 12.833333 1015.1667
## 15     15 14.714286 1015.8571
## 16     16 13.833333 1016.6667
## 17     17 16.500000 1015.8333
## 18     18 13.833333 1014.5000
## 19     19 12.500000 1021.6667
## 20     20 11.285714 1020.2857
## 21     21 11.200000 1021.0000
## 22     22  9.500000 1021.8000
## 23     23 14.000000 1020.5000
## 24     24 13.833333 1020.5000
## 25     25 12.250000 1020.7500
## 26     26 12.666667 1019.6667
## 27     27 12.857143 1018.1429
## 28     28 14.833333 1017.8333
## 29     29 14.125000 1016.6250
## 30     30 14.714286 1017.8571
## 31     31 16.200000 1018.4000
## 32     32 16.000000 1016.0000
## 33     33 16.285714 1017.1429
## 34     34 18.000000 1015.2857
## 35     35 17.428571 1014.5714
## 36     36 16.625000 1016.3750
## 37     37 16.666667 1018.6667
## 38     38 15.600000 1018.6000
## 39     39 14.000000 1017.1429
## 40     40 15.428571 1016.8571
## 41     41 15.250000 1017.1250
## 42     42 15.875000 1016.5000
## 43     43 15.333333 1017.5000
## 44     44 16.285714 1016.0000
## 45     45 17.333333 1014.1667
## 46     46 19.166667 1011.6667
## 47     47 14.428571 1008.0000
## 48     48 13.666667 1012.6667
## 49     49 15.600000 1016.2000
## 50     50 15.857143 1015.7143
## 51     51 17.714286 1017.0000
## 52     52 20.000000 1015.4286
## 53     53 20.500000 1016.0000
## 54     54 17.428571 1017.0000
## 55     55 16.857143 1018.8571
## 56     56 16.875000 1018.2500
## 57     57 17.857143 1015.1429
## 58     58 20.800000 1015.2000
## 59     59 19.428571 1017.4286
## 60     60 17.333333 1016.3333
## 61     61 19.000000 1016.2500
## 62     62 19.333333 1016.1667
## 63     63 17.600000 1015.8000
## 64     64 20.875000 1016.3750
## 65     65 20.857143 1015.7143
## 66     66 23.428571 1015.2857
## 67     67 24.166667 1014.8333
## 68     68 25.428571 1009.2857
## 69     69 23.142857 1008.0000
## 70     70 24.000000 1011.1667
## 71     71 23.500000 1011.5000
## 72     72 21.500000 1009.0000
## 73     73 22.333333 1013.1667
## 74     74 24.166667 1016.1667
## 75     75 20.333333 1016.1667
## 76     76 22.666667 1015.0000
## 77     77 23.428571 1009.4286
## 78     78 22.500000 1011.0000
## 79     79 29.166667 1009.5000
## 80     80 23.833333 1008.3333
## 81     81 25.250000 1006.3750
## 82     82 27.375000 1007.6250
## 83     83 27.000000 1010.2500
## 84     84 23.500000 1010.6250
## 85     85 24.142857 1008.8571
## 86     86 21.000000 1009.0000
## 87     87 22.428571 1009.5714
## 88     88 21.250000 1009.7500
## 89     89 23.500000 1008.6250
## 90     90 23.200000 1008.6000
## 91     91 25.375000 1008.5000
## 92     92 25.166667 1009.5000
## 93     93 26.200000 1009.0000
## 94     94 24.600000 1007.8000
## 95     95 25.600000 1007.0000
## 96     96 25.857143 1006.1429
## 97     97 29.142857 1005.0000
## 98     98 28.714286 1006.0000
## 99     99 30.166667 1005.3333
## 100   100 30.000000 1006.7143
## 101   101 30.000000 1006.4000
## 102   102 28.857143 1007.5714
## 103   103 30.200000 1008.4000
## 104   104 28.250000 1007.0000
## 105   105 28.250000 1003.8750
## 106   106 32.125000 1000.0000
## 107   107 29.200000 1002.2000
## 108   108 30.285714 1002.8571
## 109   109 28.285714 1002.5714
## 110   110 30.625000 1003.3750
## 111   111 27.666667 1006.8333
## 112   112 27.375000 1010.0000
## 113   113 28.625000 1009.8750
## 114   114 30.285714 1008.5714
## 115   115 31.142857 1007.7143
## 116   116 29.875000 1008.3750
## 117   117 31.142857 1007.0000
## 118   118 30.571429 1005.4286
## 119   119 32.125000 1004.8750
## 120   120 31.142857 1004.8571
## 121   121 31.857143 1002.8333
## 122   122 29.833333 1001.2000
## 123   123 28.571429  999.4286
## 124   124 32.857143  999.0000
## 125   125 32.625000 1000.6250
## 126   126 32.750000 1001.6250
## 127   127 32.875000 1002.0000
## 128   128 34.500000 1001.1667
## 129   129 34.285714 1000.8571
## 130   130 34.000000 1000.1429
## 131   131 30.750000  999.8750
## 132   132 29.857143 1003.7143
## 133   133 31.714286 1003.5714
## 134   134 32.285714 1001.5714
## 135   135 33.000000  999.8571
## 136   136 33.000000  998.5000
## 137   137 32.833333 1000.3333
## 138   138 31.400000  999.0000
## 139   139 35.333333  999.6667
## 140   140 36.400000  998.4000
## 141   141 36.000000  998.6667
## 142   142 36.750000  998.6250
## 143   143 37.500000  997.1667
## 144   144 38.428571  996.4286
## 145   145 38.714286  998.1429
## 146   146 37.800000  998.2000
## 147   147 35.857143  999.1429
## 148   148 35.333333 1000.1667
## 149   149 34.142857  999.4286
## 150   150 32.200000  999.6000
## 151   151 33.625000  998.7143
## 152   152 32.000000  998.7500
## 153   153 32.400000 1001.0000
## 154   154 35.600000  999.6000
## 155   155 35.857143  999.8571
## 156   156 37.166667  998.0000
## 157   157 31.285714  997.2857
## 158   158 34.000000  998.0000
## 159   159 34.200000  997.0000
## 160   160 36.166667  995.1667
## 161   161 36.625000  995.1250
## 162   162 30.166667 1000.0000
## 163   163 34.142857  998.0000
## 164   164 29.833333  998.6667
## 165   165 30.142857  997.4286
## 166   166 30.714286  997.8571
## 167   167 27.000000  996.4286
## 168   168 26.875000  994.7500
## 169   169 28.400000  996.4000
## 170   170 29.857143  999.0000
## 171   171 33.000000  997.1667
## 172   172 34.833333  995.3333
## 173   173 35.600000  996.2000
## 174   174 35.166667  995.6667
## 175   175 33.142857  996.2857
## 176   176 30.571429  999.1667
## 177   177 30.666667  999.0000
## 178   178 31.428571  996.7143
## 179   179 31.500000  996.0000
## 180   180 33.250000  995.2500
## 181   181 32.833333  997.1667
## 182   182 33.857143  996.4286
## 183   183 33.142857  996.0000
## 184   184 31.571429  997.5714
## 185   185 32.375000  997.0000
## 186   186 32.800000  997.8000
## 187   187 31.000000 1000.5714
## 188   188 31.166667 1000.1667
## 189   189 29.833333  997.5000
## 190   190 29.000000  996.8333
## 191   191 29.750000  995.7500
## 192   192 32.500000  997.3333
## 193   193 28.875000 1000.0000
## 194   194 31.750000  998.8750
## 195   195 30.750000  997.3750
## 196   196 31.750000  998.3750
## 197   197 29.875000  999.0000
## 198   198 31.125000  998.0000
## 199   199 31.750000  995.0000
## 200   200 30.500000  995.3333
## 201   201 26.833333  999.1667
## 202   202 28.200000  996.4000
## 203   203 29.500000  996.8333
## 204   204 31.857143  998.0000
## 205   205 29.714286  997.1429
## 206   206 28.333333  995.1667
## 207   207 30.000000  995.4286
## 208   208 30.142857  997.2857
## 209   209 32.285714  998.1429
## 210   210 31.000000  997.1429
## 211   211 30.500000  996.5000
## 212   212 28.833333  996.6667
## 213   213 30.000000  998.0000
## 214   214 31.571429  999.4286
## 215   215 32.500000 1000.8333
## 216   216 32.500000  993.2500
## 217   217 29.500000 1001.0000
## 218   218 27.166667 1002.1667
## 219   219 28.375000 1001.8750
## 220   220 28.500000 1000.5000
## 221   221 27.166667 1001.0000
## 222   222 29.428571 1000.3333
## 223   223 30.000000 1001.1667
## 224   224 29.500000 1001.0000
## 225   225 30.750000  999.2500
## 226   226 29.666667  999.8333
## 227   227 27.714286 1001.7143
## 228   228 26.600000 1000.8000
## 229   229 27.428571  999.4286
## 230   230 28.333333 1000.1667
## 231   231 28.166667 1000.5000
## 232   232 27.375000  999.0000
## 233   233 28.333333  999.2857
## 234   234 29.166667 1001.1667
## 235   235 31.285714 1001.4286
## 236   236 32.000000 1000.3333
## 237   237 29.000000  998.7500
## 238   238 32.000000  997.8571
## 239   239 32.142857  999.8571
## 240   240 30.142857 1003.7143
## 241   241 31.500000 1003.2500
## 242   242 29.166667 1003.5000
## 243   243 29.000000 1003.2500
## 244   244 30.000000 1003.2000
## 245   245 31.857143 1003.7143
## 246   246 32.750000 1001.5000
## 247   247 31.714286 1002.1667
## 248   248 31.666667 1001.5000
## 249   249 30.571429 1003.4286
## 250   250 30.833333 1004.1667
## 251   251 28.000000 1007.2500
## 252   252 31.000000 1007.0000
## 253   253 29.666667 1005.3333
## 254   254 30.857143 1005.5714
## 255   255 31.166667 1008.1667
## 256   256 31.142857 1007.1429
## 257   257 30.500000 1003.8333
## 258   258 31.625000 1001.3750
## 259   259 29.666667 1002.3333
## 260   260 29.250000 1002.8750
## 261   261 29.142857 1003.1429
## 262   262 29.800000 1001.8000
## 263   263 28.666667 1003.0000
## 264   264 25.200000 1003.0000
## 265   265 28.333333 1001.6667
## 266   266 30.285714 1001.0000
## 267   267 30.750000 1002.7500
## 268   268 28.571429 1005.4286
## 269   269 28.200000 1006.6000
## 270   270 29.000000 1005.6250
## 271   271 30.000000 1005.3333
## 272   272 28.142857 1003.8571
## 273   273 26.857143 1004.7143
## 274   274 28.285714 1007.7143
## 275   275 29.200000 1010.0000
## 276   276 28.600000 1008.6000
## 277   277 24.833333 1007.5000
## 278   278 28.500000 1004.1250
## 279   279 29.750000 1005.2500
## 280   280 29.714286 1008.8571
## 281   281 29.200000 1009.8000
## 282   282 30.833333 1010.3333
## 283   283 29.428571 1010.8571
## 284   284 23.857143 1010.5714
## 285   285 26.142857 1008.2857
## 286   286 27.166667 1007.6667
## 287   287 27.714286 1007.0000
## 288   288 25.857143 1009.1429
## 289   289 26.428571 1009.0000
## 290   290 26.857143 1010.4286
## 291   291 26.333333 1010.5000
## 292   292 24.571429 1011.5714
## 293   293 24.333333 1010.6667
## 294   294 24.875000 1009.8750
## 295   295 27.800000 1011.0000
## 296   296 25.000000 1013.6250
## 297   297 23.857143 1013.2857
## 298   298 22.800000 1011.0000
## 299   299 23.000000 1011.2000
## 300   300 22.875000 1010.6250
## 301   301 22.666667 1010.0000
## 302   302 23.000000 1011.1429
## 303   303 22.857143 1012.8571
## 304   304 23.666667 1012.8333
## 305   305 23.428571 1013.2857
## 306   306 21.625000 1014.5000
## 307   307 20.625000 1016.0000
## 308   308 21.166667 1017.5000
## 309   309 18.833333 1017.5000
## 310   310 20.571429 1015.8571
## 311   311 20.142857 1015.6667
## 312   312 18.000000 1015.6250
## 313   313 19.857143 1016.7143
## 314   314 16.833333 1016.0000
## 315   315 18.857143 1013.4286
## 316   316 16.571429 1013.8571
## 317   317 17.750000 1014.2500
## 318   318 17.625000 1015.2500
## 319   319 17.000000 1016.5000
## 320   320 16.625000 1018.0000
## 321   321 15.571429 1017.7143
## 322   322 18.500000 1016.8750
## 323   323 17.875000 1015.5000
## 324   324 18.250000 1012.7500
## 325   325 17.875000 1013.1250
## 326   326 17.625000 1015.7500
## 327   327 18.142857 1014.5714
## 328   328 19.125000 1014.5000
## 329   329 21.250000 1012.8750
## 330   330 21.250000 1011.7500
## 331   331 19.125000 1014.7500
## 332   332 18.625000 1017.7500
## 333   333 17.750000 1016.8750
## 334   334 17.875000 1015.8750
## 335   335 18.000000 1016.1429
## 336   336 17.250000 1015.8750
## 337   337 17.500000 1016.7500
## 338   338 17.142857 1017.0000
## 339   339 17.142857 1015.2857
## 340   340 16.125000 1014.0000
## 341   341 16.000000 1014.2500
## 342   342 16.500000 1014.0000
## 343   343 16.125000 1013.5000
## 344   344 15.500000 1013.7500
## 345   345 16.625000 1014.1250
## 346   346 17.750000 1012.5000
## 347   347 17.142857 1012.0000
## 348   348 16.142857 1012.5714
## 349   349 15.500000 1015.1250
## 350   350 15.250000 1015.6250
## 351   351 14.750000 1014.2500
## 352   352 14.875000 1012.6250
## 353   353 16.125000 1012.8750
## 354   354 15.375000 1016.0000
## 355   355 14.750000 1017.0000
## 356   356 15.250000 1018.3750
## 357   357 14.250000 1020.5000
## 358   358 13.500000 1021.3750
## 359   359 13.666667 1020.6250
## 360   360 12.125000 1019.2500
## 361   361 11.875000 1018.6250
## 362   362 10.875000 1019.2500
## 363   363 10.571429 1018.5000
## 364   364 12.375000 1018.1250
## 365   365 14.500000 1020.5000
## 366   366 13.375000 1021.0000
## 367   367 11.000000 1020.2500
## 368   368 12.500000 1017.7500
## 369   369 12.875000 1016.2500
## 370   370 12.375000 1014.5000
## 371   371 11.428571 1017.6250
## 372   372 12.142857 1017.7500
## 373   373 11.875000 1017.8750
## 374   374 12.833333 1017.0000
## 375   375 12.375000 1016.6250
## 376   376 12.125000 1019.6250
## 377   377 12.875000 1020.7500
## 378   378 14.000000 1019.3750
## 379   379 14.875000 1017.0000
## 380   380 12.000000 1019.7500
## 381   381 12.285714 1019.5714
## 382   382 12.000000 1017.6250
## 383   383 12.500000 1020.2500
## 384   384 14.500000 1023.0000
## 385   385 14.625000 1021.3750
## 386   386 13.571429 1020.5714
## 387   387 15.250000 1019.8750
## 388   388 15.625000 1020.2500
## 389   389 13.875000 1019.2500
## 390   390 11.750000 1019.7500
## 391   391 14.375000 1019.1250
## 392   392 16.625000 1018.7500
## 393   393 16.625000 1016.7500
## 394   394 14.875000 1018.6250
## 395   395 14.500000 1019.3750
## 396   396 14.750000 1019.2500
## 397   397 14.000000 1019.7500
## 398   398 15.250000 1016.8750
## 399   399 18.428571 1014.1250
## 400   400 18.857143 1011.0000
## 401   401 17.000000 1011.5000
## 402   402 18.500000 1010.8750
## 403   403 19.625000 1009.2500
## 404   404 13.625000 1011.3750
## 405   405 12.875000 1012.0000
## 406   406 12.625000 1012.5000
## 407   407 13.375000 1013.5000
## 408   408 13.250000 1013.5000
## 409   409 15.250000 1013.6250
## 410   410 13.000000 1011.8750
## 411   411 12.375000 1012.8750
## 412   412 13.500000 1016.0000
## 413   413 15.125000 1017.7500
## 414   414 15.250000 1017.3750
## 415   415 15.125000 1016.6250
## 416   416 16.125000 1014.7500
## 417   417 17.125000 1014.2500
## 418   418 16.625000 1017.3750
## 419   419 16.500000 1019.3750
## 420   420 16.375000 1018.6250
## 421   421 18.125000 1017.2500
## 422   422 18.750000 1014.2500
## 423   423 18.625000 1013.1250
## 424   424 15.625000 1013.1250
## 425   425 16.142857 1015.2857
## 426   426 16.125000 1017.1250
## 427   427 17.875000 1016.2500
## 428   428 18.750000 1016.5000
## 429   429 18.875000 1016.1250
## 430   430 18.375000 1016.3750
## 431   431 18.500000 1017.1250
## 432   432 19.875000 1017.0000
## 433   433 22.000000 1014.5000
## 434   434 20.000000 1014.3750
## 435   435 20.714286 1013.7143
## 436   436 20.000000 1016.2500
## 437   437 19.000000 1016.2500
## 438   438 19.625000 1015.0000
## 439   439 21.625000 1016.3750
## 440   440 23.125000 1015.3750
## 441   441 26.250000 1011.2500
## 442   442 24.857143 1012.8571
## 443   443 21.750000 1013.2500
## 444   444 21.500000 1010.7500
## 445   445 22.500000 1007.1250
## 446   446 25.000000 1007.0000
## 447   447 22.750000 1010.6250
## 448   448 22.375000 1011.6250
## 449   449 24.625000 1010.5000
## 450   450 25.125000 1010.6250
## 451   451 23.500000 1011.7500
## 452   452 24.428571 1011.1429
## 453   453 25.875000 1011.2500
## 454   454 24.125000 1009.8750
## 455   455 24.125000 1008.2500
## 456   456 24.875000 1006.1250
## 457   457 26.750000 1005.7500
## 458   458 25.000000 1010.3750
## 459   459 25.857143 1009.8571
## 460   460 27.375000 1008.8750
## 461   461 30.125000 1007.5000
## 462   462 27.250000 1009.6250
## 463   463 25.750000 1010.3750
## 464   464 25.750000 1010.2500
## 465   465 26.375000 1010.8750
## 466   466 28.142857 1008.5714
## 467   467 29.125000 1007.1250
## 468   468 26.875000 1009.6250
## 469   469 27.625000 1010.6250
## 470   470 29.000000 1009.1250
## 471   471 27.125000 1008.2500
## 472   472 26.375000 1007.0000
## 473   473 24.250000 1011.6250
## 474   474 25.625000 1011.3750
## 475   475 27.000000 1010.5000
## 476   476 28.000000 1010.0000
## 477   477 29.375000 1011.7500
## 478   478 29.125000 1009.0000
## 479   479 29.375000 1005.8750
## 480   480 30.625000 1006.1250
## 481   481 31.500000 1006.3750
## 482   482 30.625000 1007.1250
## 483   483 31.500000 1005.6250
## 484   484 32.125000 1003.8750
## 485   485 33.250000 1003.7500
## 486   486 34.875000 1001.5000
## 487   487 34.375000 1001.8750
## 488   488 31.750000 1003.6250
## 489   489 31.875000 1003.0000
## 490   490 29.750000 1005.7500
## 491   491 29.571429 1006.8571
## 492   492 31.750000 1004.2500
## 493   493 31.375000 1002.6250
## 494   494 30.500000 1000.4286
## 495   495 30.000000 1000.7500
## 496   496 29.625000 1004.0000
## 497   497 28.000000 1004.8750
## 498   498 25.250000 1008.3750
## 499   499 26.875000 1006.5000
## 500   500 29.500000 1004.2500
## 501   501 31.125000 1003.5000
## 502   502 29.875000 1007.0000
## 503   503 31.500000 1007.3750
## 504   504 32.375000 1006.6250
## 505   505 31.125000 1006.6250
## 506   506 33.000000 1002.1250
## 507   507 33.250000  999.5000
## 508   508 34.625000 1000.0000
## 509   509 31.500000 1004.1250
## 510   510 31.500000 1005.5000
## 511   511 30.125000 1007.1250
## 512   512 31.875000 1004.7500
## 513   513 33.625000 1001.2500
## 514   514 35.250000  999.3750
## 515   515 33.125000  999.5714
## 516   516 34.625000  999.6250
## 517   517 32.750000 1003.2500
## 518   518 32.750000 1003.2500
## 519   519 34.875000 1000.0000
## 520   520 35.250000  998.1250
## 521   521 36.875000  996.1250
## 522   522 37.500000  994.8750
## 523   523 38.500000  994.0000
## 524   524 37.625000  993.1250
## 525   525 37.875000  993.1250
## 526   526 37.250000  994.7500
## 527   527 37.625000  995.6250
## 528   528 32.875000  997.0000
## 529   529 30.250000  997.6250
## 530   530 30.500000  997.3750
## 531   531 33.875000  997.1250
## 532   532 36.375000  995.6250
## 533   533 36.875000  994.1250
## 534   534 33.875000  994.2500
## 535   535 35.750000  993.2500
## 536   536 38.000000  991.3750
## 537   537 36.875000  993.8750
## 538   538 34.500000  996.3750
## 539   539 32.125000  996.2500
## 540   540 33.000000  995.5000
## 541   541 32.250000  997.6250
## 542   542 32.500000 1001.0000
## 543   543 34.750000 1001.6250
## 544   544 34.875000  999.3750
## 545   545 33.000000  997.7500
## 546   546 32.000000  997.5000
## 547   547 31.375000 1000.8750
## 548   548 29.500000 1001.1250
## 549   549 28.875000  998.2857
## 550   550 31.388889  998.7778
## 551   551 34.705882 1000.7059
## 552   552 30.625000 1000.5000
## 553   553 34.000000  998.2500
## 554   554 35.000000  996.8750
## 555   555 36.000000  996.5000
## 556   556 36.125000  996.8750
## 557   557 36.625000  996.1250
## 558   558 35.875000  996.0000
## 559   559 31.625000  998.0000
## 560   560 32.375000  997.1250
## 561   561 33.125000  996.0000
## 562   562 32.625000  997.2500
## 563   563 30.750000  996.1250
## 564   564 27.250000  994.8750
## 565   565 29.875000  995.1250
## 566   566 30.875000  996.7500
## 567   567 32.125000  996.0000
## 568   568 31.625000  997.0000
## 569   569 30.625000  997.7500
## 570   570 30.500000  998.7143
## 571   571 30.750000 1000.7500
## 572   572 32.500000 1001.8750
## 573   573 31.714286 1001.8750
## 574   574 30.000000  999.3750
## 575   575 30.250000  999.2500
## 576   576 31.875000  999.5714
## 577   577 32.500000  999.6250
## 578   578 33.250000  998.0000
## 579   579 30.125000  997.1250
## 580   580 31.500000  995.3750
## 581   581 30.500000  996.1250
## 582   582 29.750000  997.5000
## 583   583 32.125000  996.8750
## 584   584 31.125000  997.0000
## 585   585 31.000000  998.3750
## 586   586 29.235294 1001.6471
## 587   587 28.500000 1002.1562
## 588   588 31.125000  999.8750
## 589   589 32.125000  998.7500
## 590   590 32.125000  998.2500
## 591   591 30.375000  997.3750
## 592   592 31.125000  998.8750
## 593   593 30.428571 1002.5000
## 594   594 31.375000 1003.2500
## 595   595 31.875000 1005.3750
## 596   596 33.000000 1004.8571
## 597   597 33.250000 1002.3750
## 598   598 32.500000 1001.2500
## 599   599 33.000000 1002.8571
## 600   600 33.125000 1003.2500
## 601   601 33.125000 1002.8750
## 602   602 33.875000 1002.3750
## 603   603 34.125000 1001.3750
## 604   604 33.857143 1001.1250
## 605   605 29.625000 1003.7500
## 606   606 30.375000 1003.5000
## 607   607 28.625000 1002.6250
## 608   608 27.625000 1001.2500
## 609   609 29.428571 1001.1250
## 610   610 29.250000 1001.8750
## 611   611 28.625000 1001.4286
## 612   612 27.875000 1000.3750
## 613   613 26.500000 1000.1429
## 614   614 28.500000 1002.3750
## 615   615 31.000000 1003.2500
## 616   616 31.250000 1003.1250
## 617   617 31.000000 1002.3750
## 618   618 28.750000 1003.2500
## 619   619 26.571429 1005.5714
## 620   620 28.500000 1005.1250
## 621   621 29.375000 1004.0000
## 622   622 29.125000 1004.2500
## 623   623 29.625000 1004.2500
## 624   624 30.625000 1002.2500
## 625   625 30.625000 1002.1250
## 626   626 30.500000 1003.8750
## 627   627 31.000000 1003.5000
## 628   628 31.125000 1002.8750
## 629   629 30.500000 1003.6250
## 630   630 30.625000 1004.5000
## 631   631 31.000000 1004.2857
## 632   632 29.625000 1006.2500
## 633   633 30.000000 1007.7500
## 634   634 30.000000 1008.1250
## 635   635 30.500000 1009.2500
## 636   636 30.250000 1009.6250
## 637   637 30.500000 1009.1250
## 638   638 30.625000 1009.6250
## 639   639 30.428571 1010.4286
## 640   640 30.375000 1009.1250
## 641   641 30.500000 1008.0000
## 642   642 31.000000 1008.5000
## 643   643 30.625000 1010.3750
## 644   644 30.625000 1010.1250
## 645   645 31.000000 1008.3750
## 646   646 28.750000 1009.8750
## 647   647 27.625000 1010.1250
## 648   648 27.125000 1008.7500
## 649   649 25.714286 1008.2500
## 650   650 26.000000 1008.3750
## 651   651 27.625000 1009.8750
## 652   652 24.500000 1010.8750
## 653   653 24.500000 1013.3750
## 654   654 23.000000 1014.0000
## 655   655 25.000000 1013.6250
## 656   656 24.000000 1012.8750
## 657   657 23.750000 1012.6250
## 658   658 24.750000 1013.1250
## 659   659 25.625000 1014.6250
## 660   660 25.750000 1014.5000
## 661   661 25.250000 1012.7500
## 662   662 25.500000 1011.7500
## 663   663 25.250000 1012.3750
## 664   664 25.125000 1014.8750
## 665   665 24.875000 1013.7500
## 666   666 25.125000 1012.0000
## 667   667 24.750000 1012.1250
## 668   668 23.875000 1014.5000
## 669   669 23.375000 1013.6250
## 670   670 22.500000 1014.1250
## 671   671 22.875000 1014.5000
## 672   672 23.375000 1013.5000
## 673   673 23.625000 1012.1250
## 674   674 22.750000 1012.3750
## 675   675 23.000000 1011.6250
## 676   676 24.000000 1012.8750
## 677   677 23.875000 1014.6250
## 678   678 24.125000 1015.1250
## 679   679 21.375000 1014.0000
## 680   680 19.250000 1014.7500
## 681   681 19.625000 1015.8750
## 682   682 18.500000 1016.5000
## 683   683 18.375000 1016.3750
## 684   684 18.625000 1014.8750
## 685   685 18.250000 1015.3750
## 686   686 17.750000 1016.1250
## 687   687 17.625000 1017.8750
## 688   688 17.500000 1016.5000
## 689   689 20.375000 1013.3750
## 690   690 18.250000 1012.7500
## 691   691 19.875000 1011.3750
## 692   692 17.125000 1014.6250
## 693   693 16.875000 1016.0000
## 694   694 16.500000 1014.8750
## 695   695 19.625000 1009.8750
## 696   696 18.125000 1011.7500
## 697   697 19.375000 1012.5000
## 698   698 20.250000 1012.3750
## 699   699 19.750000 1014.0000
## 700   700 20.000000 1013.7500
## 701   701 21.375000 1012.7500
## 702   702 25.000000 1007.0000
## 703   703 21.250000 1010.3750
## 704   704 17.750000 1017.1250
## 705   705 18.250000 1016.3750
## 706   706 17.250000 1014.7500
## 707   707 17.000000 1015.0000
## 708   708 15.285714 1016.8571
## 709   709 21.500000 1010.5000
## 710   710 19.625000 1011.6250
## 711   711 17.500000 1014.3750
## 712   712 16.000000 1016.1250
## 713   713 16.875000 1014.5000
## 714   714 16.625000 1016.1250
## 715   715 17.750000 1017.3750
## 716   716 12.875000 1020.2500
## 717   717 11.750000 1021.3750
## 718   718 11.750000 1022.1250
## 719   719 11.250000 1018.7500
## 720   720  9.500000 1018.0000
## 721   721  9.375000 1017.1250
## 722   722  9.625000 1017.8750
## 723   723  9.875000 1021.2500
## 724   724  9.250000 1020.1250
## 725   725 10.750000 1017.7500
## 726   726 10.375000 1018.6250
## 727   727  9.000000 1019.7500
## 728   728 11.125000 1017.2500
## 729   729 11.625000 1014.6250
## 730   730 12.375000 1016.8750
## 731   731 14.750000 1017.5000
## 732   732 14.875000 1017.8750
## 733   733 15.125000 1017.3750
## 734   734 14.125000 1019.6250
## 735   735 14.000000 1016.0000
## 736   736 12.000000 1015.5000
## 737   737  9.625000 1016.6250
## 738   738 10.000000 1019.1250
## 739   739 10.625000 1020.1250
## 740   740 11.125000 1018.8750
## 741   741 11.000000 1019.1250
## 742   742 10.625000 1019.5000
## 743   743 12.250000 1017.1250
## 744   744 12.000000 1017.5000
## 745   745 12.375000 1019.3750
## 746   746 13.000000 1019.2500
## 747   747 13.500000 1019.2500
## 748   748 13.000000 1018.7500
## 749   749 13.750000 1018.8750
## 750   750 13.375000 1020.0000
## 751   751 13.250000 1021.0000
## 752   752 12.625000 1020.6250
## 753   753 14.625000 1020.0000
## 754   754 12.375000 1020.7500
## 755   755 13.250000 1020.3750
## 756   756 13.750000 1019.8750
## 757   757 13.250000 1019.3333
## 758   758 11.714286 1018.2500
## 759   759 11.500000 1019.7500
## 760   760 12.750000 1020.3750
## 761   761 13.750000 1017.8750
## 762   762 15.000000 1017.7500
## 763   763 17.375000 1017.5000
## 764   764 16.750000 1020.5000
## 765   765 14.625000 1022.0000
## 766   766 14.625000 1019.2500
## 767   767 15.125000 1019.0000
## 768   768 15.875000 1017.0000
## 769   769 16.625000 1017.3750
## 770   770 16.875000 1019.7500
## 771   771 15.750000 1019.7500
## 772   772 16.000000 1017.1250
## 773   773 15.500000 1016.1250
## 774   774 17.250000 1012.6250
## 775   775 19.375000 1009.6250
## 776   776 18.875000 1012.6250
## 777   777 20.625000 1013.6250
## 778   778 21.142857 1011.8571
## 779   779 21.250000 1011.7500
## 780   780 22.750000 1012.0000
## 781   781 22.750000 1013.8750
## 782   782 21.875000 1014.1250
## 783   783 22.142857 1013.3750
## 784   784 21.500000 1014.2500
## 785   785 22.875000 1010.6250
## 786   786 23.125000 1008.5000
## 787   787 20.000000 1009.8750
## 788   788 19.250000 1013.0000
## 789   789 21.250000 1013.0000
## 790   790 17.375000 1011.1250
## 791   791 17.500000 1012.5000
## 792   792 16.250000 1013.7500
## 793   793 16.750000 1014.0000
## 794   794 18.750000 1015.7500
## 795   795 20.000000 1016.3750
## 796   796 18.375000 1013.7143
## 797   797 19.500000 1015.3750
## 798   798 18.125000 1021.0000
## 799   799 17.750000 1019.1250
## 800   800 18.250000 1018.1250
## 801   801 20.500000 1017.0000
## 802   802 21.375000 1014.7500
## 803   803 21.500000 1015.1250
## 804   804 17.625000 1015.5000
## 805   805 19.857143 1017.3750
## 806   806 20.875000 1015.8750
## 807   807 20.250000 1013.5000
## 808   808 21.125000 1014.8750
## 809   809 22.500000 1014.8750
## 810   810 24.875000 1014.8750
## 811   811 25.125000 1012.0000
## 812   812 24.750000 1010.0000
## 813   813 25.750000 1010.2500
## 814   814 26.625000 1011.8571
## 815   815 28.500000 1013.0000
## 816   816 27.375000 1012.0000
## 817   817 27.000000 1009.6250
## 818   818 25.000000 1008.6250
## 819   819 23.500000 1010.3750
## 820   820 25.250000 1010.1429
## 821   821 26.428571 1007.5000
## 822   822 28.714286 1004.6250
## 823   823 24.500000 1004.7500
## 824   824 21.000000 1005.2500
## 825   825 23.571429 1007.7500
## 826   826 25.000000 1007.7500
## 827   827 23.750000 1010.2500
## 828   828 25.625000 1012.0000
## 829   829 26.000000 1010.5000
## 830   830 27.125000 1008.5000
## 831   831 28.250000 1008.3750
## 832   832 24.375000 1012.3750
## 833   833 24.750000 1013.7500
## 834   834 26.250000 1012.7500
## 835   835 27.625000 1012.8750
## 836   836 25.250000 1013.8750
## 837   837 28.125000 1007.1250
## 838   838 30.125000 1002.3750
## 839   839 32.125000 1003.8750
## 840   840 32.875000 1005.7500
## 841   841 31.750000 1005.6250
## 842   842 30.750000 1003.1250
## 843   843 30.750000 1004.0000
## 844   844 31.125000 1004.3750
## 845   845 30.500000 1005.8750
## 846   846 28.750000 1008.0000
## 847   847 30.125000 1006.7500
## 848   848 31.000000 1004.5000
## 849   849 31.250000 1005.0000
## 850   850 32.625000 1004.0000
## 851   851 31.250000 1005.8750
## 852   852 30.125000 1006.2500
## 853   853 31.000000 1006.6250
## 854   854 30.875000 1006.6250
## 855   855 32.125000 1003.8750
## 856   856 34.125000 1001.8750
## 857   857 35.428571 1003.5714
## 858   858 34.125000 1003.2500
## 859   859 35.125000 1001.6250
## 860   860 35.125000 1001.3750
## 861   861 34.000000 1003.6250
## 862   862 33.125000 1004.8750
## 863   863 27.750000 1006.1250
## 864   864 29.625000 1007.2500
## 865   865 28.500000 1006.5000
## 866   866 30.000000 1003.2500
## 867   867 33.250000 1001.1250
## 868   868 35.375000 1000.0000
## 869   869 31.500000 1002.1250
## 870   870 33.125000 1002.0000
## 871   871 35.625000 1000.3750
## 872   872 36.250000  998.1250
## 873   873 37.375000  996.1250
## 874   874 36.000000  997.3750
## 875   875 37.750000  996.8750
## 876   876 35.250000 1001.5000
## 877   877 33.750000 1001.2500
## 878   878 34.750000 1000.6250
## 879   879 35.125000 1001.6250
## 880   880 32.250000 1003.8750
## 881   881 34.250000 1001.7500
## 882   882 29.375000 1003.2500
## 883   883 31.625000 1003.5000
## 884   884 28.625000 1004.6250
## 885   885 31.875000 1004.0000
## 886   886 32.625000 1003.5000
## 887   887 33.375000 1001.3750
## 888   888 35.000000  999.2500
## 889   889 36.500000  998.8750
## 890   890 37.625000  998.8750
## 891   891 36.625000  999.6250
## 892   892 35.625000 1000.0000
## 893   893 35.875000  998.5000
## 894   894 31.125000  998.8571
## 895   895 27.000000 1000.2857
## 896   896 31.000000  999.5714
## 897   897 32.375000  999.1429
## 898   898 33.875000  998.0000
## 899   899 35.125000  997.6250
## 900   900 35.875000  996.2500
## 901   901 33.125000  997.8750
## 902   902 30.375000  998.1250
## 903   903 32.125000  995.6250
## 904   904 31.500000  995.0000
## 905   905 37.250000  995.3750
## 906   906 26.625000  994.7500
## 907   907 31.125000  995.2500
## 908   908 32.625000  997.1250
## 909   909 34.250000 1000.2857
## 910   910 31.125000 1003.6250
## 911   911 30.875000 1002.6250
## 912   912 32.000000 1001.8750
## 913   913 33.750000 1001.1250
## 914   914 33.000000  999.6250
## 915   915 34.125000  997.5000
## 916   916 35.714286  995.4286
## 917   917 27.875000  997.8750
## 918   918 28.125000  998.1250
## 919   919 31.250000  998.6250
## 920   920 29.625000  998.6250
## 921   921 26.625000  998.7500
## 922   922 25.500000  999.3333
## 923   923 25.500000  998.6250
## 924   924 29.875000 1000.6250
## 925   925 30.875000  997.8571
## 926   926 33.375000  993.7500
## 927   927 31.750000  994.7500
## 928   928 30.750000  997.2500
## 929   929 29.500000  998.0000
## 930   930 30.125000  997.6250
## 931   931 29.125000  999.1250
## 932   932 31.375000 1000.0000
## 933   933 32.250000 1001.1250
## 934   934 32.125000 1001.7500
## 935   935 32.000000 1001.0000
## 936   936 30.857143  998.5714
## 937   937 30.000000  998.8750
## 938   938 29.625000 1000.2857
## 939   939 29.500000 1000.1250
## 940   940 28.500000  998.0000
## 941   941 29.000000  998.8750
## 942   942 29.000000  999.7500
## 943   943 28.250000 1002.7500
## 944   944 29.625000 1003.6250
## 945   945 31.500000 1001.8750
## 946   946 29.428571 1000.0714
## 947   947 29.593750 1000.2812
## 948   948 29.500000  998.0000
## 949   949 29.250000 1000.0000
## 950   950 28.875000 1002.0000
## 951   951 27.375000 1000.2500
## 952   952 30.375000 1000.0000
## 953   953 28.375000 1003.0000
## 954   954 30.750000 1003.1250
## 955   955 32.000000 1001.8750
## 956   956 30.875000 1002.0000
## 957   957 28.500000 1002.7500
## 958   958 28.625000 1002.8750
## 959   959 30.750000 1001.3750
## 960   960 29.125000 1002.5714
## 961   961 31.250000 1002.3750
## 962   962 30.750000 1001.3750
## 963   963 29.375000 1002.0000
## 964   964 30.250000 1000.8750
## 965   965 29.875000 1000.8571
## 966   966 30.750000 1001.7500
## 967   967 32.000000 1001.0000
## 968   968 32.000000 1001.1250
## 969   969 32.500000 1002.0000
## 970   970 32.000000 1002.2500
## 971   971 32.250000 1001.3750
## 972   972 32.625000 1001.3750
## 973   973 31.857143 1001.8333
## 974   974 32.125000 1000.8571
## 975   975 31.875000 1001.7500
## 976   976 31.625000 1005.5000
## 977   977 32.000000 1007.8750
## 978   978 30.875000 1007.1250
## 979   979 31.250000 1006.2500
## 980   980 31.500000 1005.7500
## 981   981 31.625000 1005.5000
## 982   982 31.000000 1005.6250
## 983   983 31.250000 1006.7500
## 984   984 30.750000 1006.2500
## 985   985 32.142857 1006.0000
## 986   986 32.125000 1005.8750
## 987   987 31.875000 1003.6250
## 988   988 33.000000 1003.2500
## 989   989 32.375000 1003.6250
## 990   990 31.375000 1004.3750
## 991   991 31.000000 1005.1250
## 992   992 30.625000 1003.6250
## 993   993 29.142857 1002.8571
## 994   994 30.125000 1002.5000
## 995   995 30.250000 1002.5000
## 996   996 28.625000 1006.3750
## 997   997 28.250000 1004.6250
## 998   998 28.500000 1003.7500
## 999   999 29.250000 1005.2500
## 1000 1000 28.714286 1007.0000
## 1001 1001 28.250000 1009.0000
## 1002 1002 29.000000 1010.5714
## 1003 1003 29.000000 1011.3750
## 1004 1004 28.125000 1012.0000
## 1005 1005 28.285714 1010.7500
## 1006 1006 30.000000 1008.7143
## 1007 1007 29.000000 1008.1250
## 1008 1008 28.875000 1009.2500
## 1009 1009 28.375000 1009.0000
## 1010 1010 28.625000 1009.1250
## 1011 1011 28.375000 1007.5000
## 1012 1012 28.500000 1004.2500
## 1013 1013 28.500000 1005.0000
## 1014 1014 28.571429 1008.7143
## 1015 1015 29.250000 1010.0000
## 1016 1016 29.500000 1011.0000
## 1017 1017 28.125000 1013.2500
## 1018 1018 27.500000 1013.5000
## 1019 1019 27.875000 1013.5000
## 1020 1020 28.000000 1012.2500
## 1021 1021 26.857143 1013.1250
## 1022 1022 28.142857 1011.7143
## 1023 1023 27.714286 1009.8571
## 1024 1024 26.500000 1011.8750
## 1025 1025 24.750000 1012.2500
## 1026 1026 24.000000 1011.8750
## 1027 1027 24.875000 1011.6250
## 1028 1028 26.375000 1012.7500
## 1029 1029 25.625000 1014.8750
## 1030 1030 24.750000 1013.8750
## 1031 1031 20.875000 1015.8750
## 1032 1032 22.125000 1016.8750
## 1033 1033 22.375000 1016.3750
## 1034 1034 22.375000 1016.1250
## 1035 1035 22.875000 1016.2500
## 1036 1036 23.375000 1015.8750
## 1037 1037 24.000000 1014.2500
## 1038 1038 23.750000 1012.2500
## 1039 1039 22.125000 1012.3750
## 1040 1040 22.875000 1012.7500
## 1041 1041 21.875000 1014.6250
## 1042 1042 21.875000 1015.8750
## 1043 1043 21.750000 1015.0000
## 1044 1044 22.750000 1013.5000
## 1045 1045 21.875000 1014.3750
## 1046 1046 21.875000 1015.3750
## 1047 1047 22.375000 1014.8750
## 1048 1048 22.375000 1012.0000
## 1049 1049 22.500000 1012.0000
## 1050 1050 21.875000 1012.6250
## 1051 1051 19.875000 1013.5000
## 1052 1052 18.500000 1016.3750
## 1053 1053 18.500000 1017.5000
## 1054 1054 18.500000 1016.1250
## 1055 1055 18.750000 1015.5000
## 1056 1056 18.625000 1015.1250
## 1057 1057 17.000000 1015.8750
## 1058 1058 18.500000 1016.0000
## 1059 1059 19.375000 1014.7500
## 1060 1060 20.000000 1012.7500
## 1061 1061 18.500000 1011.2500
## 1062 1062 18.750000 1015.3750
## 1063 1063 18.125000 1018.0000
## 1064 1064 19.500000 1016.7500
## 1065 1065 19.250000 1014.7500
## 1066 1066 18.125000 1016.1250
## 1067 1067 17.625000 1016.0000
## 1068 1068 17.625000 1015.0000
## 1069 1069 17.500000 1015.8750
## 1070 1070 17.125000 1016.5000
## 1071 1071 17.625000 1017.2500
## 1072 1072 18.428571 1015.7143
## 1073 1073 18.375000 1015.2500
## 1074 1074 19.000000 1011.7500
## 1075 1075 17.250000 1011.3750
## 1076 1076 13.250000 1014.2500
## 1077 1077 12.875000 1016.5000
## 1078 1078 13.250000 1017.0000
## 1079 1079 13.125000 1016.5000
## 1080 1080 13.000000 1016.3750
## 1081 1081 12.250000 1016.8750
## 1082 1082 12.875000 1018.6250
## 1083 1083 12.375000 1019.0000
## 1084 1084 11.750000 1018.5000
## 1085 1085 12.625000 1018.8750
## 1086 1086 12.000000 1020.8750
## 1087 1087 12.500000 1020.6250
## 1088 1088 11.250000 1021.2500
## 1089 1089 11.500000 1022.0000
## 1090 1090 12.750000 1021.1250
## 1091 1091 15.375000 1020.6250
## 1092 1092 17.125000 1020.8750
## 1093 1093 16.375000 1018.1250
## 1094 1094 15.500000 1017.5000
## 1095 1095 15.000000 1020.5000

Kemudian, data dibagi menjadi 2 yaitu data training dan data testing. Data training digunakan untuk menentukan model yang sesuai dengan keseluruhan data, sedangkan data testing digunakan untuk menguji performa dari model yang telah didapatkan. Data dibagi dengan perbandingan 80 : 20, yaitu dengan data training sebanyak 876 amatan (amatan 1 sampai 876) dan data testing sebanyak 219 amatan (amatan 877 sampai 1095).

#SPLIT DATA
train<-data1[1:876,]
test<-data1[877:1095,]
train <- as.data.frame(train)
test <- as.data.frame(test)
#data time series
train.ts<-ts(train)
test.ts<-ts(test)
data.ts<-ts(data1)

Eksplorasi Data

Plot Time Series

data.ts1 <- ts(yt)
plot(data.ts1, main = " Time Series Plot of Daily Delhi Climate", xlab = "Periode", ylab ="Mean Pressure")
points(data.ts1)

Berdasarkan plot time series di atas, dapat dilihat bahwa peubah mean pressure atau rata-rata tekanan udara memiliki pola musiman.

#Time Series Plot
data.ts<-ts(xt)
plot(data.ts, main = "Time Series Plot of Daily Delhi Climate", xlab = "Period", ylab="Mean Temperature")
points(data.ts)

Berdasarkan plot time series di atas, dapat dilihat bahwa peubah mean temperature atau rata-rata temperatur udara memiliki pola musiman.

Korelasi Peubah Unique Visits (X) dan Page Loads (Y)

corrplot::corrplot(corr=cor(data1), method="number", type="upper")

cor(xt, yt)
## [1] -0.8899641

Scatter Plot Peubah Unique Visits (X) dan Page Loads (Y)

plot(xt, yt, pch = 20, col = "blue", main = "Scatter Plot Rata-Rata Tekanan Udara vs Rata-Rata Temperatur Udara")

Berdasarkan corrplot dan scatter plot di atas, dapat dilihat bahwa hubungan peubah Mean Temperature (xt) dan peubah Mean Pressure (yt) memiliki hubungan yang negatif kuat yaitu dengan nilai korelasi sebesar -0.8899641, artinya dengan kenaikan temperatur udara, tekanan udara cenderung menurun.

Model KOYCK

Metode Koyck didasarkan asumsi bahwa semakin jauh jarak lag peubah independen dari periode sekarang maka semakin kecil pengaruh peubah lag terhadap peubah dependen. Koyck mengusulkan suatu metode untuk menduga model dinamis distributed lag dengan mengasumsikan bahwa semua koefisien 𝛽 mempunyai tanda sama.

Model Koyck merupakan jenis paling umum dari model infinite distributed lag dan juga dikenal sebagai geometric lag.

Pemodelan Koyck

#MODEL KOYCK
model.koyck <- dLagM :: koyckDlm(x=train$xt , y =train$yt)
summary(model.koyck)
## 
## Call:
## "Y ~ (Intercept) + Y.1 + X.t"
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -7.65185 -1.04223 -0.02695  1.00772  6.61856 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 120.48015   17.79467   6.771 2.35e-11 ***
## Y.1           0.88306    0.01725  51.199  < 2e-16 ***
## X.t          -0.10389    0.01786  -5.816 8.47e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.649 on 872 degrees of freedom
## Multiple R-Squared: 0.9535,  Adjusted R-squared: 0.9534 
## Wald test:  8907 on 2 and 872 DF,  p-value: < 2.2e-16 
## 
## Diagnostic tests:
## NULL
## 
##                             alpha       beta       phi
## Geometric coefficients:  1030.292 -0.1038871 0.8830621

AIC dan BIC

AIC (Akaike’s Information Criteria) dan BIC (Bayesian Information Criteria) digunakan untuk mengukur kebaikan model. BIC memiliki perhitungan yang lebih kompleks. Kebaikan model dapat dilihat dari, semakin kecil nilai AIC/BIC, semakin baik suatu model.

AIC(model.koyck)
## [1] 3363.836
BIC(model.koyck)
## [1] 3382.933

Forecasting

#ramalan
(fore.koyck <- forecast(model = model.koyck, x=test$xt, h=nrow(test)))
## $forecasts
##   [1] 1001.3607 1001.1338 1000.8945 1000.9818 1000.8512 1001.2422 1001.3538
##   [8] 1001.7641 1001.7887 1001.7325 1001.6049 1001.3235 1000.9192 1000.4452
##  [15] 1000.1306  999.9566  999.7770 1000.1119 1000.8362 1001.0602 1001.1152
##  [22] 1001.0079 1000.7833 1000.5070 1000.5488 1000.8713 1000.9744 1001.1303
##  [29] 1000.6706 1001.3685 1001.5173 1001.4928 1001.3024 1001.4589 1001.6231
##  [36] 1001.6512 1001.4942 1001.4335 1001.2630 1000.9474 1001.4830 1001.9301
##  [43] 1002.0002 1002.2309 1002.7464 1003.3184 1003.8235 1003.8151 1003.7037
##  [50] 1003.3457 1003.1983 1003.1721 1003.2788 1003.3081 1003.4378 1003.3186
##  [57] 1003.1225 1002.9623 1002.8338 1002.8390 1002.9327 1003.0544 1003.1749
##  [64] 1003.3851 1003.5188 1003.6369 1003.8191 1003.8372 1003.6583 1003.7156
##  [71] 1003.7489 1003.7882 1003.8488 1003.9413 1004.1788 1004.0768 1004.1946
##  [78] 1004.0519 1003.7960 1003.6868 1003.8372 1003.9570 1003.8421 1003.9094
##  [85] 1003.7480 1003.6575 1003.7204 1003.6851 1003.6928 1003.6088 1003.4047
##  [92] 1003.2244 1003.0133 1002.8789 1002.7342 1002.5674 1002.4999 1002.4125
##  [99] 1002.3613 1002.3420 1002.2861 1002.3535 1002.3741 1002.3664 1002.3465
## [106] 1002.3939 1002.4098 1002.4758 1002.3893 1002.3148 1002.2750 1002.1230
## [113] 1002.0537 1002.0964 1002.1730 1002.2797 1002.5278 1002.6449 1002.7353
## [120] 1002.9840 1003.2425 1003.4449 1003.5456 1003.6903 1003.8662 1003.9437
## [127] 1004.0121 1004.1634 1004.2803 1004.2055 1004.2433 1004.2896 1004.3825
## [134] 1004.4385 1004.5140 1004.5677 1004.6150 1004.6495 1004.6094 1004.5480
## [141] 1004.6366 1004.7798 1004.8673 1004.9316 1005.1071 1005.1285 1005.1919
## [148] 1005.3741 1005.7167 1006.0972 1006.3423 1006.4030 1006.5344 1006.7414
## [155] 1007.3267 1007.7137 1008.0295 1008.3084 1008.5027 1008.6223 1008.6631
## [162] 1008.7250 1008.9485 1009.0680 1009.2773 1009.4622 1009.6385 1009.6903
## [169] 1009.8269 1009.9475 1010.0021 1010.0503 1010.0799 1010.1709 1010.4591
## [176] 1010.8564 1011.2073 1011.5171 1011.7647 1011.9964 1012.3698 1012.5436
## [183] 1012.6063 1012.5967 1012.7441 1012.8482 1013.0051 1013.0008 1013.0230
## [190] 1013.1595 1013.3319 1013.4842 1013.6316 1013.8008 1013.8983 1013.9008
## [197] 1013.9087 1013.8507 1013.9812 1014.5121 1015.0198 1015.4292 1015.8037
## [204] 1016.1475 1016.5289 1016.8008 1017.0928 1017.4156 1017.6098 1017.8462
## [211] 1018.0030 1018.2714 1018.4823 1018.5388 1018.3159 1017.9373 1017.6809
## [218] 1017.5454 1017.4777
## 
## $call
## forecast.koyckDlm(model = model.koyck, x = test$xt, h = nrow(test))
## 
## attr(,"class")
## [1] "forecast.koyckDlm" "dLagM"

Akurasi

mape.koyck <- dLagM::GoF(model.koyck)["MAPE"]

Regression with Distributed Lag

Pemodelan Regression with Distributed Lag

model.dlm = dLagM::dlm(x = train$xt,y = train$yt , q = 2)
summary(model.dlm)
## 
## Call:
## lm(formula = model.formula, data = design)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.743  -2.220   0.230   2.417   7.405 
## 
## Coefficients:
##               Estimate Std. Error  t value Pr(>|t|)    
## (Intercept) 1031.20661    0.38456 2681.539  < 2e-16 ***
## x.t           -0.53001    0.06855   -7.731 2.94e-14 ***
## x.1           -0.14692    0.08855   -1.659  0.09745 .  
## x.2           -0.24428    0.06847   -3.568  0.00038 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.311 on 870 degrees of freedom
## Multiple R-squared:  0.8127, Adjusted R-squared:  0.812 
## F-statistic:  1258 on 3 and 870 DF,  p-value: < 2.2e-16
## 
## AIC and BIC values for the model:
##        AIC      BIC
## 1 4579.104 4602.969

AIC dan BIC

AIC(model.dlm)
## [1] 4579.104
BIC(model.dlm)
## [1] 4602.969

Akurasi

#akurasi data training
mape.dlm <- GoF(model.dlm)["MAPE"]

Lag Optimum

#penentuan lag optimum 
finiteDLMauto(formula = yt ~ xt,
              data = data.frame(train), q.min = 1, q.max = 5 ,
              model.type = "dlm", error.type = "AIC", trace = TRUE) ##q max lag maksimum
##   q - k    MASE      AIC      BIC   GMRAE   MBRAE R.Adj.Sq Ljung-Box
## 5     5 1.97631 4534.450 4572.607 3.20129 0.79117  0.81857         0
## 4     4 1.98502 4549.724 4583.120 3.16957 0.55588  0.81627         0
## 3     3 2.00707 4564.009 4592.641 3.24695 5.51654  0.81419         0
## 2     2 2.01415 4579.104 4602.969 3.19890 0.63890  0.81202         0
## 1     1 2.01895 4596.848 4615.945 3.13677 0.73277  0.80920         0

Pemodelan Lag Optimum

#model dlm dengan lag optimum
model.dlm2 = dLagM::dlm(x = train$xt,y = train$yt , q = 5) #terdapat lag yang tidak signifikan sehingga dapat dikurangi jumlah lagnya 
summary(model.dlm2)
## 
## Call:
## lm(formula = model.formula, data = design)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -9.7981 -2.2534  0.2189  2.3736  7.4943 
## 
## Coefficients:
##               Estimate Std. Error  t value Pr(>|t|)    
## (Intercept) 1031.59603    0.38341 2690.574  < 2e-16 ***
## x.t           -0.48326    0.06856   -7.049  3.7e-12 ***
## x.1           -0.11436    0.08753   -1.306  0.19175    
## x.2           -0.06154    0.08733   -0.705  0.48122    
## x.3           -0.05114    0.08737   -0.585  0.55851    
## x.4           -0.04800    0.08761   -0.548  0.58394    
## x.5           -0.17893    0.06836   -2.617  0.00902 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.251 on 864 degrees of freedom
## Multiple R-squared:  0.8198, Adjusted R-squared:  0.8186 
## F-statistic: 655.2 on 6 and 864 DF,  p-value: < 2.2e-16
## 
## AIC and BIC values for the model:
##       AIC      BIC
## 1 4534.45 4572.607

AIC dan BIC

AIC(model.dlm2)
## [1] 4534.45
BIC(model.dlm2)
## [1] 4572.607

Akurasi

#akurasi data training
mape.dlm2 <- GoF(model.dlm2)["MAPE"]

Model Autoregressive / Dynamic Regression

Apabila peubah dependen dipengaruhi oleh peubah independen pada waktu sekarang, serta dipengaruhii juga oleh peubah dependen itu sendiri pada satu waktu yang lalu maka model tersebut disebut autoregressive.

Pemodelan Model Autoregressive

#MODEL AUTOREGRESSIVE 
#library(dLagM)
model.ardl = ardlDlm(x = train$xt, y = train$yt, p = 1 , q = 1) #p:lag x, q:lag y
#model untuk p=1, q=1: yt=b0+b1yt-1+b2xt+b3xt-1
#model untuk p=2, q=3: yt=b0+b1yt-1+b2yt-2+b3xt+b4xt-1+b5xt-2

summary(model.ardl)
## 
## Time series regression with "ts" data:
## Start = 2, End = 876
## 
## Call:
## dynlm(formula = as.formula(model.text), data = data, start = 1)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -7.612 -1.092  0.026  1.013  6.117 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 134.59797   16.40353   8.205 8.22e-16 ***
## X.t          -0.37987    0.03260 -11.654  < 2e-16 ***
## X.1           0.25680    0.03428   7.490 1.69e-13 ***
## Y.1           0.86954    0.01591  54.653  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.586 on 871 degrees of freedom
## Multiple R-squared:  0.957,  Adjusted R-squared:  0.9569 
## F-statistic:  6465 on 3 and 871 DF,  p-value: < 2.2e-16

AIC dan BIC

AIC(model.ardl)
## [1] 3296.63
BIC(model.ardl)
## [1] 3320.501

Ramalan

#ramalan
(fore.ardl <- forecast(model = model.ardl, x=test$xt, h=219))
## $forecasts
##   [1] 1001.6729 1001.0581 1000.6379 1001.4609 1000.6786 1002.3637 1001.7225
##   [8] 1002.8823 1001.8858 1001.5690 1001.2012 1000.4567  999.6568  998.9192
##  [15]  998.9465  999.0933  998.8692 1000.5429 1002.3455 1001.3342 1000.9596
##  [22] 1000.4172  999.8559  999.4039 1000.2482 1001.3207 1000.8824 1001.1880
##  [29]  999.1091 1002.8140 1001.5978 1001.1260 1000.4836 1001.5294 1001.7313
##  [36] 1001.4153 1000.7646 1000.9331 1000.4597  999.7332 1002.4875 1002.7745
##  [43] 1001.9011 1002.5614 1003.8579 1004.6423 1005.0354 1003.7153 1003.3110
##  [50] 1002.2666 1002.6177 1002.8856 1003.3365 1003.1703 1003.5661 1002.7987
##  [57] 1002.3769 1002.2822 1002.2153 1002.5592 1002.8903 1003.1006 1003.2347
##  [64] 1003.6990 1003.6560 1003.7470 1004.1110 1003.7127 1003.0071 1003.6620
##  [71] 1003.6367 1003.6928 1003.8124 1003.9947 1004.6267 1003.6515 1004.3336
##  [78] 1003.5109 1002.9306 1003.1744 1003.9997 1004.0599 1003.3372 1003.8717
##  [85] 1003.1119 1003.1869 1003.6461 1003.3599 1003.4781 1003.1522 1002.6188
##  [92] 1002.4759 1002.1617 1002.2068 1002.0227 1001.7843 1001.9650 1001.8233
##  [99] 1001.8637 1001.9297 1001.7804 1002.1742 1002.0853 1002.0093 1001.9600
## [106] 1002.1866 1002.1282 1002.3316 1001.8509 1001.7973 1001.8412 1001.3878
## [113] 1001.5198 1001.8540 1002.0302 1002.2296 1002.8698 1002.6727 1002.7060
## [120] 1003.3844 1003.6995 1003.7822 1003.6333 1003.9001 1004.1708 1004.0020
## [127] 1004.0479 1004.4202 1004.4581 1003.8812 1004.1996 1004.2672 1004.4838
## [134] 1004.4488 1004.5775 1004.5777 1004.6100 1004.6110 1004.3724 1004.2442
## [141] 1004.7193 1005.0167 1004.9723 1004.9826 1005.4578 1005.0890 1005.2614
## [148] 1005.7625 1006.5512 1007.0725 1007.0008 1006.5933 1006.9091 1007.3235
## [155] 1008.9311 1008.8591 1009.0225 1009.2288 1009.2182 1009.1475 1008.9770
## [162] 1009.0842 1009.7305 1009.5903 1010.0408 1010.1758 1010.3406 1010.0720
## [169] 1010.4276 1010.5121 1010.3957 1010.4228 1010.3989 1010.6477 1011.4632
## [176] 1012.1811 1012.4523 1012.6880 1012.7981 1013.0054 1013.7709 1013.4495
## [183] 1013.2227 1013.0129 1013.5607 1013.5569 1013.8552 1013.4318 1013.5116
## [190] 1013.9442 1014.2214 1014.3341 1014.4795 1014.7163 1014.6360 1014.3893
## [197] 1014.4015 1014.1610 1014.7770 1016.3829 1016.8944 1017.1005 1017.4235
## [204] 1017.7198 1018.2302 1018.2439 1018.6064 1019.0305 1018.9065 1019.2607
## [211] 1019.2183 1019.7847 1019.8612 1019.5171 1018.5417 1017.7029 1017.7078
## [218] 1017.8518 1017.9423
## 
## $call
## forecast.ardlDlm(model = model.ardl, x = test$xt, h = 219)
## 
## attr(,"class")
## [1] "forecast.ardlDlm" "dLagM"

AKurasi

mape.ardl <- GoF(model.ardl)["MAPE"]
#PEMODELAN DLM dan ARDL dengan library dynlm
#library(dynlm)

#sama dengan model dlm p=1
cons_lm1 <- dynlm(yt ~ xt+L(xt),data = train.ts)

#sama dengan model ardl p=0 q=1
cons_lm2 <- dynlm(yt ~ xt+L(yt),data = train.ts)

#sama dengan ardl p=1 q=1
cons_lm3 <- dynlm(yt ~ xt+L(xt)+L(yt),data = train.ts)

#sama dengan dlm p=2
cons_lm4 <- dynlm(yt ~ xt+L(xt)+L(xt,2),data = train.ts)
summary(cons_lm1)
## 
## Time series regression with "ts" data:
## Start = 2, End = 876
## 
## Call:
## dynlm(formula = yt ~ xt + L(xt), data = train.ts)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -11.6727  -2.2683   0.2323   2.3889   7.4706 
## 
## Coefficients:
##               Estimate Std. Error  t value Pr(>|t|)    
## (Intercept) 1031.04429    0.38506 2677.622  < 2e-16 ***
## xt            -0.55871    0.06822   -8.190 9.23e-16 ***
## L(xt)         -0.35579    0.06815   -5.221 2.23e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.337 on 872 degrees of freedom
## Multiple R-squared:  0.8096, Adjusted R-squared:  0.8092 
## F-statistic:  1854 on 2 and 872 DF,  p-value: < 2.2e-16
summary(cons_lm2)
## 
## Time series regression with "ts" data:
## Start = 2, End = 876
## 
## Call:
## dynlm(formula = yt ~ xt + L(yt), data = train.ts)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -7.5832 -1.0472 -0.0479  1.0470  6.2044 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 174.88826   15.97860   10.95   <2e-16 ***
## xt           -0.16389    0.01567  -10.46   <2e-16 ***
## L(yt)         0.83058    0.01550   53.57   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.636 on 872 degrees of freedom
## Multiple R-squared:  0.9543, Adjusted R-squared:  0.9541 
## F-statistic:  9095 on 2 and 872 DF,  p-value: < 2.2e-16
summary(cons_lm3)
## 
## Time series regression with "ts" data:
## Start = 2, End = 876
## 
## Call:
## dynlm(formula = yt ~ xt + L(xt) + L(yt), data = train.ts)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -7.612 -1.092  0.026  1.013  6.117 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 134.59797   16.40353   8.205 8.22e-16 ***
## xt           -0.37987    0.03260 -11.654  < 2e-16 ***
## L(xt)         0.25680    0.03428   7.490 1.69e-13 ***
## L(yt)         0.86954    0.01591  54.653  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.586 on 871 degrees of freedom
## Multiple R-squared:  0.957,  Adjusted R-squared:  0.9569 
## F-statistic:  6465 on 3 and 871 DF,  p-value: < 2.2e-16
summary(cons_lm4)
## 
## Time series regression with "ts" data:
## Start = 3, End = 876
## 
## Call:
## dynlm(formula = yt ~ xt + L(xt) + L(xt, 2), data = train.ts)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.743  -2.220   0.230   2.417   7.405 
## 
## Coefficients:
##               Estimate Std. Error  t value Pr(>|t|)    
## (Intercept) 1031.20661    0.38456 2681.539  < 2e-16 ***
## xt            -0.53001    0.06855   -7.731 2.94e-14 ***
## L(xt)         -0.14692    0.08855   -1.659  0.09745 .  
## L(xt, 2)      -0.24428    0.06847   -3.568  0.00038 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.311 on 870 degrees of freedom
## Multiple R-squared:  0.8127, Adjusted R-squared:  0.812 
## F-statistic:  1258 on 3 and 870 DF,  p-value: < 2.2e-16

SSE

deviance(cons_lm1)
## [1] 9708.382
deviance(cons_lm2)
## [1] 2333.016
deviance(cons_lm3)
## [1] 2191.837
deviance(cons_lm4)
## [1] 9537.731

Berdasarkan output di atas, model yang memiliki nilai SSE terendah yaitu cons_lm3, model tersebut merupakan model yang sama seperti ardl dengan p=1 q=1 # Diagnostik Model

Uji Non Autokorelasi

Hipotesis

H0: Tidak terdapat autokorelasi

H1: Terdapat autokorelasi

if(require("lmtest")) encomptest(cons_lm1, cons_lm2)
## Encompassing test
## 
## Model 1: yt ~ xt + L(xt)
## Model 2: yt ~ xt + L(yt)
## Model E: yt ~ xt + L(xt) + L(yt)
##           Res.Df Df        F    Pr(>F)    
## M1 vs. ME    871 -1 2986.951 < 2.2e-16 ***
## M2 vs. ME    871 -1   56.102 1.688e-13 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#durbin watson
dwtest(cons_lm1)
## 
##  Durbin-Watson test
## 
## data:  cons_lm1
## DW = 0.26376, p-value < 2.2e-16
## alternative hypothesis: true autocorrelation is greater than 0
dwtest(cons_lm2)
## 
##  Durbin-Watson test
## 
## data:  cons_lm2
## DW = 1.6421, p-value = 4.156e-08
## alternative hypothesis: true autocorrelation is greater than 0
dwtest(cons_lm3)
## 
##  Durbin-Watson test
## 
## data:  cons_lm3
## DW = 1.721, p-value = 1.404e-05
## alternative hypothesis: true autocorrelation is greater than 0
dwtest(cons_lm4)
## 
##  Durbin-Watson test
## 
## data:  cons_lm4
## DW = 0.2637, p-value < 2.2e-16
## alternative hypothesis: true autocorrelation is greater than 0

Berdasarkan output di atas, pengujian non-autokorelasi dengan Durbin-Watson test semua model memiliki p-value < alpha, sehingga tolak H0. Artinya semua model terdapat autokorelasi.

Uji Heterogenitas

Hipotesis

H0: Tidak terdapat Heterogenitas

H1: Terdapat Heterogenitas

bptest(cons_lm1)
## 
##  studentized Breusch-Pagan test
## 
## data:  cons_lm1
## BP = 15.011, df = 2, p-value = 0.00055
bptest(cons_lm2)
## 
##  studentized Breusch-Pagan test
## 
## data:  cons_lm2
## BP = 3.0001, df = 2, p-value = 0.2231
bptest(cons_lm3)
## 
##  studentized Breusch-Pagan test
## 
## data:  cons_lm3
## BP = 8.0233, df = 3, p-value = 0.04553
bptest(cons_lm4)
## 
##  studentized Breusch-Pagan test
## 
## data:  cons_lm4
## BP = 20.289, df = 3, p-value = 0.0001478

Berdasarkan output di atas, pengujian heterogenitas dengan Breusch-Pagan test hanya satu model yang memiliki p-value > alpha yaitu model cons_lm2 yang artinya tidak terdapat heterogenitas, dan model lainnya memiliki memiliki p-value < alpha yang artinya terdapat heterogenitas.

Normalitas

Hipotesis

H0: Sisaan berdistribusi normal

H1: Sisaan tidak berdistribusi normal

shapiro.test(residuals(cons_lm1))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(cons_lm1)
## W = 0.98982, p-value = 9.312e-06
shapiro.test(residuals(cons_lm2))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(cons_lm2)
## W = 0.99251, p-value = 0.0002126
shapiro.test(residuals(cons_lm3))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(cons_lm3)
## W = 0.99397, p-value = 0.001393
shapiro.test(residuals(cons_lm4))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(cons_lm4)
## W = 0.99029, p-value = 1.585e-05

Berdasarkan output di atas, pengujian normalitas dengan Shapiro-Wilk normality test semua model memiliki p-value < alpha, sehingga tolak H0. Artinya semua model memiliki sisaan yang tidak berdistribusi normal.

Perbandingan Keakuratan Ramalan

#PERBANDINGAN
akurasi <- matrix(c(mape.koyck, mape.dlm, mape.dlm2, mape.ardl))
row.names(akurasi)<- c("Koyck","DLM 1","DLM 2","Autoregressive")
colnames(akurasi) <- c("MAPE")
akurasi
##                MAPE       
## Koyck          0.001257333
## DLM 1          0.002665095
## DLM 2          0.002615808
## Autoregressive 0.001222942

Simpulan

Berdasarkan perbangingan keakuratan ramalan, model yang terbaik yaitu model autoregressive karena model tersebut memiliki nilai MAPE terkecil.

Daftar Pustaka

https://www.kaggle.com/datasets/sumanthvrao/daily-climate-time-series-data?select=DailyDelhiClimateTrain.csv

Siswanti, K.Y. 2011. Model Fungsi Transfer Multivariat dan Aplikasinya untuk Meramalkan Curah Hujan Di Kota Yogyakarta. Skripsi. Yogyakarta : FMIPA Universitas Negeri Yogyakarta.