R version 2.6.1 (2007-11-26)
Copyright (C) 2007 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(209
+ ,71
+ ,3
+ ,283
+ ,3810
+ ,216
+ ,71
+ ,3
+ ,469
+ ,3558
+ ,216
+ ,71
+ ,1
+ ,536
+ ,3397
+ ,223
+ ,72
+ ,4
+ ,781
+ ,3407
+ ,214
+ ,72
+ ,3
+ ,1056
+ ,3646
+ ,214
+ ,72
+ ,4
+ ,1032
+ ,3614
+ ,214
+ ,72
+ ,2
+ ,942
+ ,3621
+ ,214
+ ,72
+ ,4
+ ,918
+ ,3607
+ ,215
+ ,72
+ ,-1
+ ,903
+ ,3642
+ ,220
+ ,72
+ ,-3
+ ,781
+ ,3600
+ ,205
+ ,72
+ ,-4
+ ,763
+ ,3758
+ ,217
+ ,72
+ ,-2
+ ,719
+ ,3510
+ ,220
+ ,72
+ ,-1
+ ,768
+ ,3576
+ ,220
+ ,72
+ ,-1
+ ,768
+ ,3579
+ ,219
+ ,72
+ ,0
+ ,762
+ ,3599
+ ,223
+ ,72
+ ,1
+ ,734
+ ,3552
+ ,224
+ ,72
+ ,3
+ ,715
+ ,3520
+ ,204
+ ,72
+ ,3
+ ,719
+ ,3732
+ ,206
+ ,72
+ ,5
+ ,711
+ ,3556
+ ,204
+ ,72
+ ,6
+ ,725
+ ,3596
+ ,205
+ ,72
+ ,6
+ ,719
+ ,3587
+ ,210
+ ,72
+ ,4
+ ,724
+ ,3560
+ ,210
+ ,72
+ ,5
+ ,760
+ ,3594
+ ,209
+ ,72
+ ,7
+ ,767
+ ,3599
+ ,221
+ ,72
+ ,7
+ ,791
+ ,3539
+ ,223
+ ,72
+ ,7
+ ,758
+ ,3563
+ ,223
+ ,72
+ ,10
+ ,768
+ ,3579
+ ,223
+ ,72
+ ,6
+ ,762
+ ,3586
+ ,225
+ ,72
+ ,2
+ ,757
+ ,3561
+ ,214
+ ,72
+ ,1
+ ,747
+ ,3650
+ ,201
+ ,72
+ ,1
+ ,744
+ ,3672
+ ,194
+ ,72
+ ,3
+ ,715
+ ,3675
+ ,210
+ ,72
+ ,8
+ ,703
+ ,3543
+ ,211
+ ,72
+ ,6
+ ,796
+ ,3615
+ ,210
+ ,72
+ ,7
+ ,831
+ ,3736
+ ,229
+ ,72
+ ,3
+ ,827
+ ,3549
+ ,226
+ ,72
+ ,-1
+ ,795
+ ,3600
+ ,226
+ ,72
+ ,-1
+ ,852
+ ,3674
+ ,226
+ ,72
+ ,-3
+ ,882
+ ,3660
+ ,226
+ ,66
+ ,-6
+ ,894
+ ,3696
+ ,226
+ ,66
+ ,-4
+ ,905
+ ,3712
+ ,226
+ ,66
+ ,-1
+ ,922
+ ,3719
+ ,224
+ ,66
+ ,2
+ ,1012
+ ,3747
+ ,224
+ ,66
+ ,4
+ ,991
+ ,3695
+ ,224
+ ,66
+ ,3
+ ,918
+ ,3680
+ ,224
+ ,72
+ ,5
+ ,901
+ ,3638
+ ,226
+ ,72
+ ,4
+ ,901
+ ,3583
+ ,225
+ ,72
+ ,8
+ ,901
+ ,3614
+ ,226
+ ,72
+ ,2
+ ,901
+ ,3610
+ ,183
+ ,72
+ ,5
+ ,901
+ ,3929
+ ,208
+ ,72
+ ,3
+ ,901
+ ,3676
+ ,197
+ ,72
+ ,3
+ ,901
+ ,3809
+ ,213
+ ,72
+ ,7
+ ,901
+ ,3635
+ ,214
+ ,72
+ ,4
+ ,901
+ ,3633
+ ,213
+ ,72
+ ,5
+ ,901
+ ,3633
+ ,213
+ ,72
+ ,6
+ ,901
+ ,3630
+ ,213
+ ,72
+ ,7
+ ,901
+ ,3634
+ ,178
+ ,72
+ ,6
+ ,901
+ ,3943
+ ,200
+ ,72
+ ,7
+ ,901
+ ,3620
+ ,198
+ ,72
+ ,5
+ ,901
+ ,3657
+ ,197
+ ,72
+ ,5
+ ,901
+ ,3669
+ ,198
+ ,72
+ ,4
+ ,901
+ ,3653
+ ,198
+ ,72
+ ,4
+ ,901
+ ,3658
+ ,200
+ ,72
+ ,4
+ ,901
+ ,3527
+ ,230
+ ,72
+ ,4
+ ,901
+ ,3432
+ ,230
+ ,72
+ ,5
+ ,901
+ ,3468
+ ,231
+ ,72
+ ,4
+ ,901
+ ,3572
+ ,231
+ ,72
+ ,3
+ ,901
+ ,3577
+ ,231
+ ,72
+ ,2
+ ,901
+ ,3584
+ ,230
+ ,72
+ ,2
+ ,901
+ ,3592
+ ,231
+ ,72
+ ,1
+ ,901
+ ,3584
+ ,231
+ ,72
+ ,3
+ ,901
+ ,3585
+ ,229
+ ,72
+ ,3
+ ,901
+ ,3593
+ ,222
+ ,72
+ ,6
+ ,766
+ ,3547
+ ,218
+ ,72
+ ,8
+ ,766
+ ,3569
+ ,218
+ ,72
+ ,4
+ ,766
+ ,3570
+ ,218
+ ,72
+ ,5
+ ,766
+ ,3573
+ ,235
+ ,72
+ ,5
+ ,766
+ ,3459
+ ,218
+ ,72
+ ,8
+ ,766
+ ,3667
+ ,218
+ ,72
+ ,11
+ ,766
+ ,3567
+ ,233
+ ,72
+ ,8
+ ,766
+ ,3482
+ ,232
+ ,72
+ ,8
+ ,766
+ ,3535
+ ,232
+ ,72
+ ,8
+ ,766
+ ,3533
+ ,232
+ ,72
+ ,8
+ ,766
+ ,3532
+ ,233
+ ,72
+ ,9
+ ,766
+ ,3523
+ ,232
+ ,72
+ ,10
+ ,766
+ ,3535
+ ,228
+ ,72
+ ,9
+ ,766
+ ,3589
+ ,210
+ ,72
+ ,6
+ ,766
+ ,3673
+ ,212
+ ,72
+ ,7
+ ,766
+ ,3562
+ ,212
+ ,72
+ ,6
+ ,766
+ ,3588
+ ,211
+ ,72
+ ,11
+ ,766
+ ,3580
+ ,224
+ ,72
+ ,4
+ ,766
+ ,3476
+ ,220
+ ,72
+ ,3
+ ,766
+ ,3601
+ ,218
+ ,72
+ ,3
+ ,766
+ ,3593
+ ,218
+ ,72
+ ,5
+ ,766
+ ,3601
+ ,220
+ ,72
+ ,5
+ ,766
+ ,3574
+ ,211
+ ,72
+ ,7
+ ,766
+ ,3755
+ ,220
+ ,72
+ ,9
+ ,766
+ ,3460
+ ,219
+ ,72
+ ,8
+ ,766
+ ,3471
+ ,220
+ ,72
+ ,9
+ ,766
+ ,3449
+ ,223
+ ,72
+ ,11
+ ,766
+ ,3551
+ ,223
+ ,72
+ ,8
+ ,766
+ ,3553
+ ,213
+ ,72
+ ,10
+ ,766
+ ,3649
+ ,214
+ ,72
+ ,11
+ ,766
+ ,3598
+ ,214
+ ,72
+ ,12
+ ,701
+ ,3573
+ ,214
+ ,72
+ ,11
+ ,701
+ ,3573
+ ,215
+ ,72
+ ,11
+ ,701
+ ,3568
+ ,227
+ ,72
+ ,7
+ ,701
+ ,3439
+ ,224
+ ,72
+ ,8
+ ,701
+ ,3501
+ ,224
+ ,72
+ ,11
+ ,701
+ ,3505
+ ,224
+ ,72
+ ,10
+ ,701
+ ,3508
+ ,224
+ ,72
+ ,9
+ ,701
+ ,3508
+ ,224
+ ,72
+ ,11
+ ,701
+ ,3506
+ ,218
+ ,72
+ ,12
+ ,701
+ ,3609
+ ,217
+ ,72
+ ,13
+ ,701
+ ,3564
+ ,231
+ ,72
+ ,14
+ ,701
+ ,3446
+ ,221
+ ,72
+ ,16
+ ,701
+ ,3540
+ ,220
+ ,72
+ ,16
+ ,701
+ ,3521
+ ,221
+ ,72
+ ,14
+ ,701
+ ,3522
+ ,220
+ ,72
+ ,15
+ ,701
+ ,3534
+ ,223
+ ,72
+ ,15
+ ,701
+ ,3524
+ ,225
+ ,72
+ ,11
+ ,701
+ ,3482
+ ,231
+ ,72
+ ,8
+ ,701
+ ,3469
+ ,229
+ ,72
+ ,9
+ ,701
+ ,3544
+ ,227
+ ,72
+ ,10
+ ,701
+ ,3528
+ ,227
+ ,72
+ ,10
+ ,701
+ ,3516
+ ,219
+ ,72
+ ,13
+ ,701
+ ,3571
+ ,218
+ ,72
+ ,15
+ ,701
+ ,3541
+ ,226
+ ,72
+ ,15
+ ,701
+ ,3472
+ ,227
+ ,72
+ ,15
+ ,701
+ ,3496
+ ,227
+ ,72
+ ,17
+ ,701
+ ,3486
+ ,227
+ ,72
+ ,18
+ ,701
+ ,3502
+ ,227
+ ,72
+ ,13
+ ,701
+ ,3525
+ ,227
+ ,72
+ ,11
+ ,701
+ ,3521
+ ,228
+ ,72
+ ,8
+ ,813
+ ,3551
+ ,183
+ ,72
+ ,9
+ ,813
+ ,3845
+ ,162
+ ,72
+ ,12
+ ,813
+ ,3835
+ ,169
+ ,72
+ ,16
+ ,813
+ ,3775
+ ,169
+ ,72
+ ,13
+ ,813
+ ,3795
+ ,168
+ ,72
+ ,13
+ ,813
+ ,3807
+ ,180
+ ,72
+ ,18
+ ,813
+ ,3677
+ ,181
+ ,72
+ ,15
+ ,813
+ ,3747
+ ,183
+ ,72
+ ,11
+ ,813
+ ,3713
+ ,184
+ ,72
+ ,15
+ ,813
+ ,3705
+ ,182
+ ,72
+ ,19
+ ,813
+ ,3735
+ ,182
+ ,72
+ ,14
+ ,813
+ ,3726
+ ,182
+ ,72
+ ,19
+ ,813
+ ,3732
+ ,189
+ ,67
+ ,25
+ ,813
+ ,3649
+ ,189
+ ,67
+ ,19
+ ,813
+ ,3681
+ ,213
+ ,67
+ ,12
+ ,813
+ ,3573
+ ,216
+ ,67
+ ,12
+ ,813
+ ,3582
+ ,215
+ ,67
+ ,13
+ ,813
+ ,3595
+ ,215
+ ,67
+ ,18
+ ,813
+ ,3598
+ ,215
+ ,67
+ ,21
+ ,813
+ ,3596
+ ,215
+ ,67
+ ,22
+ ,813
+ ,3597
+ ,214
+ ,67
+ ,22
+ ,813
+ ,3609
+ ,214
+ ,67
+ ,23
+ ,813
+ ,3614
+ ,214
+ ,67
+ ,23
+ ,813
+ ,3621
+ ,214
+ ,67
+ ,25
+ ,813
+ ,3613
+ ,214
+ ,67
+ ,25
+ ,813
+ ,3619
+ ,214
+ ,67
+ ,24
+ ,813
+ ,3619
+ ,214
+ ,67
+ ,18
+ ,813
+ ,3603
+ ,214
+ ,67
+ ,17
+ ,813
+ ,3620
+ ,214
+ ,67
+ ,14
+ ,813
+ ,3609
+ ,214
+ ,67
+ ,16
+ ,813
+ ,3591
+ ,227
+ ,72
+ ,18
+ ,817
+ ,3515
+ ,228
+ ,72
+ ,18
+ ,815
+ ,3558
+ ,228
+ ,72
+ ,17
+ ,823
+ ,3556
+ ,227
+ ,72
+ ,21
+ ,855
+ ,3572
+ ,224
+ ,72
+ ,20
+ ,852
+ ,3585
+ ,233
+ ,72
+ ,20
+ ,864
+ ,3516
+ ,234
+ ,72
+ ,20
+ ,870
+ ,3541
+ ,234
+ ,72
+ ,22
+ ,864
+ ,3536
+ ,234
+ ,72
+ ,23
+ ,813
+ ,3526
+ ,234
+ ,72
+ ,22
+ ,833
+ ,3527
+ ,233
+ ,73
+ ,23
+ ,821
+ ,3537
+ ,225
+ ,72
+ ,24
+ ,765
+ ,3567
+ ,205
+ ,72
+ ,21
+ ,726
+ ,3646
+ ,218
+ ,72
+ ,22
+ ,758
+ ,3535
+ ,219
+ ,72
+ ,23
+ ,802
+ ,3558
+ ,219
+ ,72
+ ,22
+ ,958
+ ,3599
+ ,219
+ ,72
+ ,20
+ ,966
+ ,3597
+ ,218
+ ,72
+ ,21
+ ,982
+ ,3609
+ ,225
+ ,72
+ ,22
+ ,910
+ ,3535
+ ,227
+ ,72
+ ,25
+ ,968
+ ,3552
+ ,222
+ ,72
+ ,22
+ ,1000
+ ,3582
+ ,222
+ ,72
+ ,22
+ ,933
+ ,3535
+ ,221
+ ,72
+ ,20
+ ,977
+ ,3557
+ ,221
+ ,72
+ ,21
+ ,978
+ ,3561
+ ,221
+ ,72
+ ,23
+ ,992
+ ,3566
+ ,221
+ ,72
+ ,17
+ ,963
+ ,3559
+ ,221
+ ,72
+ ,18
+ ,763
+ ,3511
+ ,221
+ ,72
+ ,16
+ ,781
+ ,3525
+ ,220
+ ,72
+ ,18
+ ,770
+ ,3530
+ ,221
+ ,72
+ ,18
+ ,779
+ ,3527
+ ,221
+ ,72
+ ,24
+ ,776
+ ,3525
+ ,221
+ ,72
+ ,21
+ ,814
+ ,3537
+ ,221
+ ,72
+ ,22
+ ,846
+ ,3579
+ ,228
+ ,72
+ ,16
+ ,933
+ ,3520
+ ,230
+ ,72
+ ,15
+ ,877
+ ,3496
+ ,230
+ ,72
+ ,17
+ ,879
+ ,3504
+ ,231
+ ,72
+ ,19
+ ,786
+ ,3480
+ ,231
+ ,72
+ ,21
+ ,786
+ ,3475
+ ,219
+ ,72
+ ,20
+ ,771
+ ,3564
+ ,218
+ ,70
+ ,15
+ ,771
+ ,3555
+ ,217
+ ,70
+ ,17
+ ,745
+ ,3553
+ ,228
+ ,70
+ ,16
+ ,752
+ ,3498
+ ,219
+ ,70
+ ,19
+ ,745
+ ,3593
+ ,217
+ ,70
+ ,23
+ ,755
+ ,3556
+ ,217
+ ,70
+ ,26
+ ,753
+ ,3561
+ ,217
+ ,70
+ ,27
+ ,872
+ ,3609
+ ,217
+ ,70
+ ,27
+ ,802
+ ,3575
+ ,227
+ ,70
+ ,30
+ ,718
+ ,3473
+ ,228
+ ,70
+ ,31
+ ,763
+ ,3523
+ ,227
+ ,70
+ ,29
+ ,757
+ ,3542
+ ,227
+ ,70
+ ,28
+ ,815
+ ,3557
+ ,227
+ ,70
+ ,25
+ ,816
+ ,3564
+ ,227
+ ,70
+ ,23
+ ,761
+ ,3543
+ ,227
+ ,70
+ ,21
+ ,869
+ ,3572
+ ,227
+ ,70
+ ,20
+ ,862
+ ,3546
+ ,227
+ ,70
+ ,21
+ ,826
+ ,3503
+ ,227
+ ,70
+ ,23
+ ,797
+ ,3499
+ ,227
+ ,70
+ ,23
+ ,900
+ ,3510
+ ,228
+ ,70
+ ,22
+ ,831
+ ,3508
+ ,227
+ ,70
+ ,17
+ ,769
+ ,3487
+ ,215
+ ,70
+ ,17
+ ,775
+ ,3579
+ ,218
+ ,67
+ ,18
+ ,867
+ ,3562
+ ,218
+ ,67
+ ,20
+ ,873
+ ,3574
+ ,217
+ ,67
+ ,21
+ ,958
+ ,3609
+ ,217
+ ,67
+ ,20
+ ,975
+ ,3593
+ ,218
+ ,67
+ ,20
+ ,987
+ ,3555
+ ,201
+ ,67
+ ,23
+ ,939
+ ,3747
+ ,220
+ ,67
+ ,23
+ ,942
+ ,3521
+ ,210
+ ,66
+ ,25
+ ,771
+ ,3705
+ ,223
+ ,65
+ ,23
+ ,833
+ ,3508
+ ,211
+ ,65
+ ,22
+ ,899
+ ,3705
+ ,211
+ ,65
+ ,22
+ ,767
+ ,3645
+ ,212
+ ,65
+ ,19
+ ,781
+ ,3640
+ ,221
+ ,63
+ ,22
+ ,772
+ ,3566
+ ,222
+ ,63
+ ,22
+ ,780
+ ,3619
+ ,224
+ ,63
+ ,25
+ ,760
+ ,3585
+ ,224
+ ,63
+ ,27
+ ,790
+ ,3602
+ ,237
+ ,63
+ ,22
+ ,777
+ ,3510
+ ,236
+ ,63
+ ,22
+ ,788
+ ,3565
+ ,236
+ ,63
+ ,19
+ ,803
+ ,3570
+ ,237
+ ,63
+ ,20
+ ,821
+ ,3568
+ ,232
+ ,63
+ ,20
+ ,804
+ ,3588
+ ,229
+ ,63
+ ,22
+ ,802
+ ,3597
+ ,229
+ ,63
+ ,25
+ ,767
+ ,3567
+ ,228
+ ,63
+ ,24
+ ,803
+ ,3635
+ ,228
+ ,63
+ ,23
+ ,829
+ ,3666
+ ,227
+ ,63
+ ,23
+ ,761
+ ,3630
+ ,225
+ ,63
+ ,22
+ ,771
+ ,3617
+ ,226
+ ,63
+ ,18
+ ,790
+ ,3588
+ ,226
+ ,63
+ ,13
+ ,775
+ ,3581
+ ,237
+ ,63
+ ,16
+ ,796
+ ,3546
+ ,234
+ ,63
+ ,17
+ ,802
+ ,3605
+ ,233
+ ,63
+ ,16
+ ,808
+ ,3602
+ ,233
+ ,63
+ ,16
+ ,784
+ ,3562
+ ,234
+ ,63
+ ,16
+ ,789
+ ,3561
+ ,234
+ ,63
+ ,12
+ ,793
+ ,3571
+ ,229
+ ,61
+ ,9
+ ,777
+ ,3622
+ ,229
+ ,61
+ ,10
+ ,791
+ ,3573
+ ,233
+ ,61
+ ,11
+ ,790
+ ,3557
+ ,235
+ ,61
+ ,15
+ ,802
+ ,3572
+ ,235
+ ,61
+ ,15
+ ,787
+ ,3566
+ ,235
+ ,61
+ ,13
+ ,786
+ ,3564
+ ,235
+ ,61
+ ,12
+ ,844
+ ,3582
+ ,235
+ ,61
+ ,12
+ ,768
+ ,3567
+ ,235
+ ,61
+ ,14
+ ,804
+ ,3574
+ ,235
+ ,61
+ ,14
+ ,774
+ ,3567
+ ,235
+ ,61
+ ,13
+ ,821
+ ,3573
+ ,235
+ ,61
+ ,13
+ ,790
+ ,3569
+ ,230
+ ,61
+ ,17
+ ,764
+ ,3603
+ ,229
+ ,61
+ ,16
+ ,781
+ ,3601
+ ,230
+ ,61
+ ,12
+ ,798
+ ,3592
+ ,230
+ ,61
+ ,9
+ ,772
+ ,3584
+ ,240
+ ,61
+ ,11
+ ,780
+ ,3514
+ ,242
+ ,61
+ ,12
+ ,815
+ ,3552
+ ,242
+ ,61
+ ,14
+ ,845
+ ,3563
+ ,242
+ ,61
+ ,13
+ ,820
+ ,3555
+ ,236
+ ,61
+ ,14
+ ,794
+ ,3576
+ ,236
+ ,61
+ ,14
+ ,774
+ ,3566
+ ,236
+ ,61
+ ,15
+ ,880
+ ,3604
+ ,230
+ ,61
+ ,12
+ ,761
+ ,3590
+ ,230
+ ,61
+ ,12
+ ,957
+ ,3694
+ ,229
+ ,61
+ ,14
+ ,1029
+ ,3717
+ ,217
+ ,61
+ ,14
+ ,977
+ ,3744
+ ,217
+ ,61
+ ,16
+ ,967
+ ,3693
+ ,217
+ ,61
+ ,16
+ ,891
+ ,3661
+ ,218
+ ,61
+ ,15
+ ,1029
+ ,3707
+ ,234
+ ,61
+ ,14
+ ,1029
+ ,3601
+ ,234
+ ,61
+ ,11
+ ,1052
+ ,3639
+ ,234
+ ,61
+ ,9
+ ,993
+ ,3610
+ ,229
+ ,61
+ ,9
+ ,967
+ ,3604
+ ,229
+ ,61
+ ,9
+ ,1008
+ ,3591
+ ,221
+ ,61
+ ,10
+ ,1118
+ ,3770
+ ,220
+ ,61
+ ,9
+ ,936
+ ,3705
+ ,220
+ ,61
+ ,8
+ ,881
+ ,3676
+ ,220
+ ,61
+ ,7
+ ,898
+ ,3682
+ ,220
+ ,61
+ ,5
+ ,901
+ ,3683
+ ,228
+ ,61
+ ,6
+ ,934
+ ,3668
+ ,229
+ ,61
+ ,6
+ ,956
+ ,3672
+ ,230
+ ,61
+ ,7
+ ,944
+ ,3648
+ ,231
+ ,61
+ ,8
+ ,958
+ ,3640
+ ,232
+ ,61
+ ,4
+ ,914
+ ,3627
+ ,231
+ ,61
+ ,4
+ ,970
+ ,3634
+ ,231
+ ,61
+ ,4
+ ,952
+ ,3616
+ ,232
+ ,61
+ ,6
+ ,1007
+ ,3627
+ ,232
+ ,61
+ ,6
+ ,1024
+ ,3610
+ ,232
+ ,61
+ ,10
+ ,1008
+ ,3593
+ ,232
+ ,61
+ ,8
+ ,1033
+ ,3601
+ ,231
+ ,61
+ ,10
+ ,1024
+ ,3603
+ ,241
+ ,65
+ ,9
+ ,1031
+ ,3462
+ ,232
+ ,65
+ ,9
+ ,1039
+ ,3607
+ ,232
+ ,65
+ ,12
+ ,994
+ ,3590
+ ,227
+ ,65
+ ,10
+ ,983
+ ,3603
+ ,227
+ ,65
+ ,7
+ ,941
+ ,3550
+ ,227
+ ,65
+ ,5
+ ,957
+ ,3555
+ ,227
+ ,65
+ ,7
+ ,974
+ ,3558
+ ,227
+ ,65
+ ,10
+ ,1014
+ ,3578
+ ,227
+ ,65
+ ,9
+ ,972
+ ,3560
+ ,223
+ ,65
+ ,3
+ ,1027
+ ,3658
+ ,222
+ ,65
+ ,-1
+ ,1013
+ ,3592
+ ,222
+ ,65
+ ,3
+ ,1007
+ ,3617
+ ,222
+ ,65
+ ,5
+ ,1017
+ ,3637
+ ,222
+ ,65
+ ,4
+ ,995
+ ,3631
+ ,222
+ ,65
+ ,2
+ ,1039
+ ,3652
+ ,222
+ ,65
+ ,0
+ ,992
+ ,3636
+ ,221
+ ,61
+ ,0
+ ,927
+ ,3640
+ ,220
+ ,61
+ ,1
+ ,881
+ ,3643
+ ,193
+ ,61
+ ,0
+ ,1023
+ ,3978
+ ,193
+ ,58
+ ,1
+ ,991
+ ,3822
+ ,208
+ ,58
+ ,-2
+ ,978
+ ,3723
+ ,208
+ ,58
+ ,-2
+ ,985
+ ,3765
+ ,208
+ ,58
+ ,1
+ ,1063
+ ,3774
+ ,209
+ ,58
+ ,2
+ ,996
+ ,3762
+ ,208
+ ,58
+ ,3
+ ,1017
+ ,3782
+ ,208
+ ,58
+ ,4
+ ,1020
+ ,3794
+ ,201
+ ,58
+ ,8
+ ,1008
+ ,3850
+ ,224
+ ,58
+ ,9
+ ,1014
+ ,3648
+ ,229
+ ,58
+ ,7
+ ,849
+ ,3586
+ ,229
+ ,58
+ ,5
+ ,1053
+ ,3695
+ ,195
+ ,58
+ ,2
+ ,1011
+ ,4000
+ ,194
+ ,58
+ ,1
+ ,917
+ ,3814
+ ,194
+ ,58
+ ,-2
+ ,903
+ ,3815
+ ,194
+ ,58
+ ,-2
+ ,908
+ ,3800
+ ,188
+ ,58
+ ,-2
+ ,982
+ ,3949
+ ,188
+ ,58
+ ,1
+ ,1007
+ ,3874
+ ,189
+ ,58
+ ,3
+ ,997
+ ,3858
+ ,188
+ ,58
+ ,5
+ ,1009
+ ,3877
+ ,183
+ ,58
+ ,4
+ ,969
+ ,3919
+ ,178
+ ,58
+ ,3
+ ,951
+ ,3923
+ ,177
+ ,58
+ ,3
+ ,954
+ ,3905
+ ,185
+ ,58
+ ,3
+ ,949
+ ,3810
+ ,185
+ ,58
+ ,3
+ ,945
+ ,3839
+ ,185
+ ,58
+ ,2
+ ,958
+ ,3859
+ ,185
+ ,58
+ ,4
+ ,986
+ ,3855
+ ,169
+ ,58
+ ,4
+ ,1043
+ ,4143
+ ,193
+ ,58
+ ,4
+ ,1044
+ ,3701
+ ,196
+ ,58
+ ,2
+ ,1054
+ ,3800
+ ,208
+ ,58
+ ,-2
+ ,1053
+ ,3748
+ ,209
+ ,58
+ ,0
+ ,1033
+ ,3763
+ ,209
+ ,58
+ ,-2
+ ,1051
+ ,3767
+ ,229
+ ,58
+ ,-4
+ ,1034
+ ,3645
+ ,228
+ ,58
+ ,-2
+ ,1087
+ ,3681
+ ,228
+ ,58
+ ,-3
+ ,1091
+ ,3669
+ ,241
+ ,58
+ ,-4
+ ,1125
+ ,3619
+ ,240
+ ,58
+ ,-5
+ ,1161
+ ,3697
+ ,240
+ ,58
+ ,-7
+ ,1098
+ ,3666
+ ,242
+ ,58
+ ,-7
+ ,1079
+ ,3616
+ ,243
+ ,58
+ ,-5
+ ,1034
+ ,3597
+ ,239
+ ,58
+ ,-5
+ ,1094
+ ,3665
+ ,238
+ ,58
+ ,-6
+ ,1066
+ ,3622
+ ,239
+ ,58
+ ,-5
+ ,1022
+ ,3600
+ ,239
+ ,58
+ ,-5
+ ,1037
+ ,3632
+ ,239
+ ,58
+ ,-8
+ ,1027
+ ,3652
+ ,239
+ ,58
+ ,-6
+ ,993
+ ,3637
+ ,239
+ ,58
+ ,-3
+ ,1031
+ ,3639
+ ,239
+ ,58
+ ,-3
+ ,986
+ ,3652
+ ,235
+ ,58
+ ,-3
+ ,988
+ ,3640
+ ,235
+ ,58
+ ,-3
+ ,1022
+ ,3606
+ ,236
+ ,58
+ ,-3
+ ,1017
+ ,3593
+ ,236
+ ,58
+ ,-2
+ ,983
+ ,3583
+ ,236
+ ,58
+ ,0
+ ,1086
+ ,3621)
+ ,dim=c(5
+ ,385)
+ ,dimnames=list(c('Pull'
+ ,'Cullet'
+ ,'Temp'
+ ,'Boost'
+ ,'SEC')
+ ,1:385))
> y <- array(NA,dim=c(5,385),dimnames=list(c('Pull','Cullet','Temp','Boost','SEC'),1:385))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '5'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
SEC Pull Cullet Temp Boost
1 3810 209 71 3 283
2 3558 216 71 3 469
3 3397 216 71 1 536
4 3407 223 72 4 781
5 3646 214 72 3 1056
6 3614 214 72 4 1032
7 3621 214 72 2 942
8 3607 214 72 4 918
9 3642 215 72 -1 903
10 3600 220 72 -3 781
11 3758 205 72 -4 763
12 3510 217 72 -2 719
13 3576 220 72 -1 768
14 3579 220 72 -1 768
15 3599 219 72 0 762
16 3552 223 72 1 734
17 3520 224 72 3 715
18 3732 204 72 3 719
19 3556 206 72 5 711
20 3596 204 72 6 725
21 3587 205 72 6 719
22 3560 210 72 4 724
23 3594 210 72 5 760
24 3599 209 72 7 767
25 3539 221 72 7 791
26 3563 223 72 7 758
27 3579 223 72 10 768
28 3586 223 72 6 762
29 3561 225 72 2 757
30 3650 214 72 1 747
31 3672 201 72 1 744
32 3675 194 72 3 715
33 3543 210 72 8 703
34 3615 211 72 6 796
35 3736 210 72 7 831
36 3549 229 72 3 827
37 3600 226 72 -1 795
38 3674 226 72 -1 852
39 3660 226 72 -3 882
40 3696 226 66 -6 894
41 3712 226 66 -4 905
42 3719 226 66 -1 922
43 3747 224 66 2 1012
44 3695 224 66 4 991
45 3680 224 66 3 918
46 3638 224 72 5 901
47 3583 226 72 4 901
48 3614 225 72 8 901
49 3610 226 72 2 901
50 3929 183 72 5 901
51 3676 208 72 3 901
52 3809 197 72 3 901
53 3635 213 72 7 901
54 3633 214 72 4 901
55 3633 213 72 5 901
56 3630 213 72 6 901
57 3634 213 72 7 901
58 3943 178 72 6 901
59 3620 200 72 7 901
60 3657 198 72 5 901
61 3669 197 72 5 901
62 3653 198 72 4 901
63 3658 198 72 4 901
64 3527 200 72 4 901
65 3432 230 72 4 901
66 3468 230 72 5 901
67 3572 231 72 4 901
68 3577 231 72 3 901
69 3584 231 72 2 901
70 3592 230 72 2 901
71 3584 231 72 1 901
72 3585 231 72 3 901
73 3593 229 72 3 901
74 3547 222 72 6 766
75 3569 218 72 8 766
76 3570 218 72 4 766
77 3573 218 72 5 766
78 3459 235 72 5 766
79 3667 218 72 8 766
80 3567 218 72 11 766
81 3482 233 72 8 766
82 3535 232 72 8 766
83 3533 232 72 8 766
84 3532 232 72 8 766
85 3523 233 72 9 766
86 3535 232 72 10 766
87 3589 228 72 9 766
88 3673 210 72 6 766
89 3562 212 72 7 766
90 3588 212 72 6 766
91 3580 211 72 11 766
92 3476 224 72 4 766
93 3601 220 72 3 766
94 3593 218 72 3 766
95 3601 218 72 5 766
96 3574 220 72 5 766
97 3755 211 72 7 766
98 3460 220 72 9 766
99 3471 219 72 8 766
100 3449 220 72 9 766
101 3551 223 72 11 766
102 3553 223 72 8 766
103 3649 213 72 10 766
104 3598 214 72 11 766
105 3573 214 72 12 701
106 3573 214 72 11 701
107 3568 215 72 11 701
108 3439 227 72 7 701
109 3501 224 72 8 701
110 3505 224 72 11 701
111 3508 224 72 10 701
112 3508 224 72 9 701
113 3506 224 72 11 701
114 3609 218 72 12 701
115 3564 217 72 13 701
116 3446 231 72 14 701
117 3540 221 72 16 701
118 3521 220 72 16 701
119 3522 221 72 14 701
120 3534 220 72 15 701
121 3524 223 72 15 701
122 3482 225 72 11 701
123 3469 231 72 8 701
124 3544 229 72 9 701
125 3528 227 72 10 701
126 3516 227 72 10 701
127 3571 219 72 13 701
128 3541 218 72 15 701
129 3472 226 72 15 701
130 3496 227 72 15 701
131 3486 227 72 17 701
132 3502 227 72 18 701
133 3525 227 72 13 701
134 3521 227 72 11 701
135 3551 228 72 8 813
136 3845 183 72 9 813
137 3835 162 72 12 813
138 3775 169 72 16 813
139 3795 169 72 13 813
140 3807 168 72 13 813
141 3677 180 72 18 813
142 3747 181 72 15 813
143 3713 183 72 11 813
144 3705 184 72 15 813
145 3735 182 72 19 813
146 3726 182 72 14 813
147 3732 182 72 19 813
148 3649 189 67 25 813
149 3681 189 67 19 813
150 3573 213 67 12 813
151 3582 216 67 12 813
152 3595 215 67 13 813
153 3598 215 67 18 813
154 3596 215 67 21 813
155 3597 215 67 22 813
156 3609 214 67 22 813
157 3614 214 67 23 813
158 3621 214 67 23 813
159 3613 214 67 25 813
160 3619 214 67 25 813
161 3619 214 67 24 813
162 3603 214 67 18 813
163 3620 214 67 17 813
164 3609 214 67 14 813
165 3591 214 67 16 813
166 3515 227 72 18 817
167 3558 228 72 18 815
168 3556 228 72 17 823
169 3572 227 72 21 855
170 3585 224 72 20 852
171 3516 233 72 20 864
172 3541 234 72 20 870
173 3536 234 72 22 864
174 3526 234 72 23 813
175 3527 234 72 22 833
176 3537 233 73 23 821
177 3567 225 72 24 765
178 3646 205 72 21 726
179 3535 218 72 22 758
180 3558 219 72 23 802
181 3599 219 72 22 958
182 3597 219 72 20 966
183 3609 218 72 21 982
184 3535 225 72 22 910
185 3552 227 72 25 968
186 3582 222 72 22 1000
187 3535 222 72 22 933
188 3557 221 72 20 977
189 3561 221 72 21 978
190 3566 221 72 23 992
191 3559 221 72 17 963
192 3511 221 72 18 763
193 3525 221 72 16 781
194 3530 220 72 18 770
195 3527 221 72 18 779
196 3525 221 72 24 776
197 3537 221 72 21 814
198 3579 221 72 22 846
199 3520 228 72 16 933
200 3496 230 72 15 877
201 3504 230 72 17 879
202 3480 231 72 19 786
203 3475 231 72 21 786
204 3564 219 72 20 771
205 3555 218 70 15 771
206 3553 217 70 17 745
207 3498 228 70 16 752
208 3593 219 70 19 745
209 3556 217 70 23 755
210 3561 217 70 26 753
211 3609 217 70 27 872
212 3575 217 70 27 802
213 3473 227 70 30 718
214 3523 228 70 31 763
215 3542 227 70 29 757
216 3557 227 70 28 815
217 3564 227 70 25 816
218 3543 227 70 23 761
219 3572 227 70 21 869
220 3546 227 70 20 862
221 3503 227 70 21 826
222 3499 227 70 23 797
223 3510 227 70 23 900
224 3508 228 70 22 831
225 3487 227 70 17 769
226 3579 215 70 17 775
227 3562 218 67 18 867
228 3574 218 67 20 873
229 3609 217 67 21 958
230 3593 217 67 20 975
231 3555 218 67 20 987
232 3747 201 67 23 939
233 3521 220 67 23 942
234 3705 210 66 25 771
235 3508 223 65 23 833
236 3705 211 65 22 899
237 3645 211 65 22 767
238 3640 212 65 19 781
239 3566 221 63 22 772
240 3619 222 63 22 780
241 3585 224 63 25 760
242 3602 224 63 27 790
243 3510 237 63 22 777
244 3565 236 63 22 788
245 3570 236 63 19 803
246 3568 237 63 20 821
247 3588 232 63 20 804
248 3597 229 63 22 802
249 3567 229 63 25 767
250 3635 228 63 24 803
251 3666 228 63 23 829
252 3630 227 63 23 761
253 3617 225 63 22 771
254 3588 226 63 18 790
255 3581 226 63 13 775
256 3546 237 63 16 796
257 3605 234 63 17 802
258 3602 233 63 16 808
259 3562 233 63 16 784
260 3561 234 63 16 789
261 3571 234 63 12 793
262 3622 229 61 9 777
263 3573 229 61 10 791
264 3557 233 61 11 790
265 3572 235 61 15 802
266 3566 235 61 15 787
267 3564 235 61 13 786
268 3582 235 61 12 844
269 3567 235 61 12 768
270 3574 235 61 14 804
271 3567 235 61 14 774
272 3573 235 61 13 821
273 3569 235 61 13 790
274 3603 230 61 17 764
275 3601 229 61 16 781
276 3592 230 61 12 798
277 3584 230 61 9 772
278 3514 240 61 11 780
279 3552 242 61 12 815
280 3563 242 61 14 845
281 3555 242 61 13 820
282 3576 236 61 14 794
283 3566 236 61 14 774
284 3604 236 61 15 880
285 3590 230 61 12 761
286 3694 230 61 12 957
287 3717 229 61 14 1029
288 3744 217 61 14 977
289 3693 217 61 16 967
290 3661 217 61 16 891
291 3707 218 61 15 1029
292 3601 234 61 14 1029
293 3639 234 61 11 1052
294 3610 234 61 9 993
295 3604 229 61 9 967
296 3591 229 61 9 1008
297 3770 221 61 10 1118
298 3705 220 61 9 936
299 3676 220 61 8 881
300 3682 220 61 7 898
301 3683 220 61 5 901
302 3668 228 61 6 934
303 3672 229 61 6 956
304 3648 230 61 7 944
305 3640 231 61 8 958
306 3627 232 61 4 914
307 3634 231 61 4 970
308 3616 231 61 4 952
309 3627 232 61 6 1007
310 3610 232 61 6 1024
311 3593 232 61 10 1008
312 3601 232 61 8 1033
313 3603 231 61 10 1024
314 3462 241 65 9 1031
315 3607 232 65 9 1039
316 3590 232 65 12 994
317 3603 227 65 10 983
318 3550 227 65 7 941
319 3555 227 65 5 957
320 3558 227 65 7 974
321 3578 227 65 10 1014
322 3560 227 65 9 972
323 3658 223 65 3 1027
324 3592 222 65 -1 1013
325 3617 222 65 3 1007
326 3637 222 65 5 1017
327 3631 222 65 4 995
328 3652 222 65 2 1039
329 3636 222 65 0 992
330 3640 221 61 0 927
331 3643 220 61 1 881
332 3978 193 61 0 1023
333 3822 193 58 1 991
334 3723 208 58 -2 978
335 3765 208 58 -2 985
336 3774 208 58 1 1063
337 3762 209 58 2 996
338 3782 208 58 3 1017
339 3794 208 58 4 1020
340 3850 201 58 8 1008
341 3648 224 58 9 1014
342 3586 229 58 7 849
343 3695 229 58 5 1053
344 4000 195 58 2 1011
345 3814 194 58 1 917
346 3815 194 58 -2 903
347 3800 194 58 -2 908
348 3949 188 58 -2 982
349 3874 188 58 1 1007
350 3858 189 58 3 997
351 3877 188 58 5 1009
352 3919 183 58 4 969
353 3923 178 58 3 951
354 3905 177 58 3 954
355 3810 185 58 3 949
356 3839 185 58 3 945
357 3859 185 58 2 958
358 3855 185 58 4 986
359 4143 169 58 4 1043
360 3701 193 58 4 1044
361 3800 196 58 2 1054
362 3748 208 58 -2 1053
363 3763 209 58 0 1033
364 3767 209 58 -2 1051
365 3645 229 58 -4 1034
366 3681 228 58 -2 1087
367 3669 228 58 -3 1091
368 3619 241 58 -4 1125
369 3697 240 58 -5 1161
370 3666 240 58 -7 1098
371 3616 242 58 -7 1079
372 3597 243 58 -5 1034
373 3665 239 58 -5 1094
374 3622 238 58 -6 1066
375 3600 239 58 -5 1022
376 3632 239 58 -5 1037
377 3652 239 58 -8 1027
378 3637 239 58 -6 993
379 3639 239 58 -3 1031
380 3652 239 58 -3 986
381 3640 235 58 -3 988
382 3606 235 58 -3 1022
383 3593 236 58 -3 1017
384 3583 236 58 -2 983
385 3621 236 58 0 1086
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Pull Cullet Temp Boost
5062.2146 -4.8805 -7.3104 -1.6016 0.1527
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-172.225 -25.229 -3.003 21.317 248.441
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5062.2146 67.8036 74.660 < 2e-16 ***
Pull -4.8805 0.1651 -29.553 < 2e-16 ***
Cullet -7.3104 0.5560 -13.149 < 2e-16 ***
Temp -1.6016 0.3079 -5.201 3.25e-07 ***
Boost 0.1527 0.0267 5.718 2.18e-08 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 47.64 on 380 degrees of freedom
Multiple R-Squared: 0.7839, Adjusted R-squared: 0.7816
F-statistic: 344.6 on 4 and 380 DF, p-value: < 2.2e-16
> postscript(file="/var/www/html/rcomp/tmp/1c4su1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/25nxz1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/3wn5w1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/4wx1t1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/506m21199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 385
Frequency = 1
1 2 3 4 5 6
248.4408433 2.2072778 -172.2249231 -153.3507988 -1.8611243 -28.5953986
7 8 9 10 11 12
-11.0582485 -18.1909017 15.9715218 13.7964627 99.7359660 -79.7777035
13 14 15 16 17 18
-5.0155718 -2.0155718 15.6216118 -5.9801239 -26.9956686 86.7843752
19 20 21 22 23 24
-75.0300872 -45.3267877 -48.5302979 -55.0945772 -24.9891131 -22.7350316
25 26 27 28 29 30
-27.8335730 10.9654981 30.2436511 31.7531929 10.8709908 46.1109803
31 32 33 34 35 36
5.1229665 -18.4095768 -62.4820013 -3.0031853 109.3744863 9.3074946
37 38 39 40 41 42
44.1450922 109.4428437 87.6594708 73.1600533 90.6839142 99.8933701
43 44 45 46 47 48
109.1969139 63.6062476 58.1496114 65.8107628 18.9700690 51.4960897
49 50 51 52 53 54
42.7668268 156.7117550 22.5201030 101.8350033 12.3289054 10.4045057
55 56 57 58 59 60
7.1256632 5.7272843 11.3289054 147.9110581 -66.1171215 -42.0812909
61 62 63 64 65 66
-34.9617545 -47.6829120 -42.6829120 -163.9219848 -112.5080766 -74.9064555
67 68 69 70 71 72
32.3723870 35.7707659 41.1691448 44.2886812 39.5675237 43.7707659
73 74 75 76 77 78
42.0098387 -12.7379548 -7.0565671 -12.4630515 -7.8614304 -38.8935491
79 80 81 82 83 84
90.9434329 -4.2517038 -20.8496130 27.2699234 25.2699234 24.2699234
85 86 87 88 89 90
21.7520081 30.4731656 63.3496901 54.6964819 -44.9409698 -20.5425909
91 92 93 94 95 96
-25.4149490 -77.1802698 26.6962546 8.9353274 20.1385696 2.8994968
97 98 99 100 101 102
143.1785666 -104.6940188 -100.1761035 -115.6940188 4.1506143 1.3457510
103 104 105 106 107 108
51.7443571 7.2264418 -6.2483204 -7.8499415 -7.9694779 -84.8103991
109 110 111 112 113 114
-35.8501688 -27.0453055 -25.6469266 -27.2485477 -26.0453055 49.2735340
115 116 117 118 119 120
0.9946915 -47.0771970 1.3214092 -22.5590544 -19.8818330 -11.1606755
121 122 123 124 125 126
-6.5192847 -45.1648419 -33.6869236 33.1537703 8.9944642 -3.0055358
127 128 129 130 131 132
17.7556187 -13.9216027 -43.8778939 -14.9974303 -21.7941881 -4.1925670
133 134 135 136 137 138
10.7993275 3.5960853 16.5725308 92.5532897 -15.1315828 -34.5618531
139 140 141 142 143 144
-19.3667164 -12.2471800 -75.6735113 -5.5979110 -36.2434681 -32.9565201
145 146 147 148 149 150
-6.3109630 -23.3190685 -9.3109630 -85.0900760 -62.6998026 -64.7800238
151 152 153 154 155 156
-41.1386330 -31.4174755 -20.4093700 -17.6045067 -15.0028856 -7.8833492
157 158 159 160 161 162
-1.2817281 5.7182719 0.9215141 6.9215141 5.3198930 -20.2898336
163 164 165 166 167 168
-4.8914547 -20.6963180 -35.4930758 -8.9024059 39.2833997 34.4604104
169 170 171 172 173 174
47.1009584 44.3159596 17.4080797 46.3725172 45.4917855 44.8796289
175 176 177 178 179 180
41.2245873 57.0882141 50.8852868 33.4253214 -17.4125035 5.3520561
181 182 183 184 185 186
20.9337551 14.5091447 20.7875657 -6.4552541 16.2556169 12.1629628
187 188 189 190 191 192
-24.6080785 -17.4093094 -11.9603593 -5.8945115 -18.0767784 -33.9409522
193 194 195 196 197 198
-25.8922728 -20.8901130 -20.3836886 -12.3159489 -10.9223112 27.7938371
199 200 201 202 203 204
-19.9350235 -27.2261399 -16.3282398 -18.0461286 -19.8428864 11.2799946
205 206 207 208 209 210
-25.2294085 -24.9371832 -28.9224018 28.0269862 -13.8541669 -3.7439615
211 212 213 214 215 216
27.6898076 4.3767793 -31.1893552 18.4225334 30.2548537 34.7983132
217 218 219 220 221 222
36.8407788 21.0344430 30.3427301 3.8098062 -32.0924158 -28.4617139
223 224 225 226 227 228
-33.1868295 -21.3736862 -45.7966518 -13.2782412 -50.0122145 -35.7249985
229 230 231 232 233 234
-16.9808781 -37.1779067 -72.1294954 49.0356959 -84.6935087 72.5014258
235 236 237 238 239 240
-81.0318101 45.7247179 5.8772932 -1.1845008 -39.7022598 16.9568356
241 242 243 244 245 246
0.5760466 16.1991580 -18.3781973 30.0619579 27.9670292 29.7010354
247 248 249 250 251 252
27.8941248 25.7613183 5.9096675 61.9314258 87.3603581 56.8615242
253 254 255 256 257 258
30.9722656 -2.4545046 -15.1725447 5.1113267 50.1555308 39.7574199
259 260 261 262 263 264
3.4215246 6.5386330 9.5214645 19.1361857 -30.3995876 -25.1234410
265 266 267 268 269 270
4.2119183 0.5019836 -4.5485875 2.9948719 -0.4021302 4.3049551
271 272 273 274 275 276
1.8850859 -0.8920734 -0.1592716 19.8143414 8.7368493 -4.3845790
277 278 279 280 281 282
-13.2199956 -32.4334856 11.5855769 21.2086883 15.4238428 12.7121290
283 284 285 286 287 288
5.7655495 29.1840419 -0.7357510 73.3407280 83.6711928 60.0445228
289 290 291 292 293 294
13.7744753 -6.6225268 21.5877142 -7.9264892 21.7572139 -1.4384378
295 296 297 298 299 300
-27.8713091 -47.1308212 77.6332783 33.9373202 11.7326055 13.5355770
301 302 303 304 305 306
10.8743217 31.4815078 37.0032088 21.3173458 17.6620362 9.8535405
307 308 309 310 311 312
3.4234995 -11.8284220 -1.1416226 -20.7370301 -28.8878093 -27.9078271
313 314 315 316 317 318
-26.2110093 -91.8350236 8.0194358 2.6944952 -10.2316838 -61.6243640
319 320 321 322 323 324
-62.2703426 -58.6625078 -39.9644855 -53.1539236 7.3175890 -67.8319646
325 326 327 328 329 330
-35.5094541 -13.8329221 -18.0757807 -6.9965480 -19.0242520 -39.2227669
331 332 333 334 335 336
-32.4787422 147.4678338 -23.9763233 -52.5895091 -11.6582063 -9.7616830
337 338 339 340 341 342
-5.0506396 8.4644264 21.6080344 51.6833258 -37.3804163 -52.9906213
343 344 345 346 347 348
21.6612474 162.3328045 -15.7982038 -17.4656727 -33.2290278 75.1905346
349 350 351 352 353 354
1.1786223 -5.2109616 10.2797647 32.3826665 13.1268059 -10.2116708
355 356 357 358 359 360
-65.4046068 -35.7939227 -19.3802672 -24.4518137 176.7585202 -148.2630243
361 362 363 364 365 366
-39.3515859 -39.0398360 -12.9027097 -14.8540304 -39.8525931 -13.6213788
367 368 369 370 371 372
-27.8336840 -21.1800931 44.8416653 20.2566977 -17.0816256 -21.1277237
373 374 375 376 377 378
18.1901604 -27.0171356 -35.8175258 -6.1075912 10.6142558 4.0083128
379 380 381 382 383 384
5.0116772 24.8818733 -6.9453232 -46.1361380 -53.4923193 -56.6998833
385
-31.2217567
> postscript(file="/var/www/html/rcomp/tmp/6kmav1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 385
Frequency = 1
lag(myerror, k = 1) myerror
0 248.4408433 NA
1 2.2072778 248.4408433
2 -172.2249231 2.2072778
3 -153.3507988 -172.2249231
4 -1.8611243 -153.3507988
5 -28.5953986 -1.8611243
6 -11.0582485 -28.5953986
7 -18.1909017 -11.0582485
8 15.9715218 -18.1909017
9 13.7964627 15.9715218
10 99.7359660 13.7964627
11 -79.7777035 99.7359660
12 -5.0155718 -79.7777035
13 -2.0155718 -5.0155718
14 15.6216118 -2.0155718
15 -5.9801239 15.6216118
16 -26.9956686 -5.9801239
17 86.7843752 -26.9956686
18 -75.0300872 86.7843752
19 -45.3267877 -75.0300872
20 -48.5302979 -45.3267877
21 -55.0945772 -48.5302979
22 -24.9891131 -55.0945772
23 -22.7350316 -24.9891131
24 -27.8335730 -22.7350316
25 10.9654981 -27.8335730
26 30.2436511 10.9654981
27 31.7531929 30.2436511
28 10.8709908 31.7531929
29 46.1109803 10.8709908
30 5.1229665 46.1109803
31 -18.4095768 5.1229665
32 -62.4820013 -18.4095768
33 -3.0031853 -62.4820013
34 109.3744863 -3.0031853
35 9.3074946 109.3744863
36 44.1450922 9.3074946
37 109.4428437 44.1450922
38 87.6594708 109.4428437
39 73.1600533 87.6594708
40 90.6839142 73.1600533
41 99.8933701 90.6839142
42 109.1969139 99.8933701
43 63.6062476 109.1969139
44 58.1496114 63.6062476
45 65.8107628 58.1496114
46 18.9700690 65.8107628
47 51.4960897 18.9700690
48 42.7668268 51.4960897
49 156.7117550 42.7668268
50 22.5201030 156.7117550
51 101.8350033 22.5201030
52 12.3289054 101.8350033
53 10.4045057 12.3289054
54 7.1256632 10.4045057
55 5.7272843 7.1256632
56 11.3289054 5.7272843
57 147.9110581 11.3289054
58 -66.1171215 147.9110581
59 -42.0812909 -66.1171215
60 -34.9617545 -42.0812909
61 -47.6829120 -34.9617545
62 -42.6829120 -47.6829120
63 -163.9219848 -42.6829120
64 -112.5080766 -163.9219848
65 -74.9064555 -112.5080766
66 32.3723870 -74.9064555
67 35.7707659 32.3723870
68 41.1691448 35.7707659
69 44.2886812 41.1691448
70 39.5675237 44.2886812
71 43.7707659 39.5675237
72 42.0098387 43.7707659
73 -12.7379548 42.0098387
74 -7.0565671 -12.7379548
75 -12.4630515 -7.0565671
76 -7.8614304 -12.4630515
77 -38.8935491 -7.8614304
78 90.9434329 -38.8935491
79 -4.2517038 90.9434329
80 -20.8496130 -4.2517038
81 27.2699234 -20.8496130
82 25.2699234 27.2699234
83 24.2699234 25.2699234
84 21.7520081 24.2699234
85 30.4731656 21.7520081
86 63.3496901 30.4731656
87 54.6964819 63.3496901
88 -44.9409698 54.6964819
89 -20.5425909 -44.9409698
90 -25.4149490 -20.5425909
91 -77.1802698 -25.4149490
92 26.6962546 -77.1802698
93 8.9353274 26.6962546
94 20.1385696 8.9353274
95 2.8994968 20.1385696
96 143.1785666 2.8994968
97 -104.6940188 143.1785666
98 -100.1761035 -104.6940188
99 -115.6940188 -100.1761035
100 4.1506143 -115.6940188
101 1.3457510 4.1506143
102 51.7443571 1.3457510
103 7.2264418 51.7443571
104 -6.2483204 7.2264418
105 -7.8499415 -6.2483204
106 -7.9694779 -7.8499415
107 -84.8103991 -7.9694779
108 -35.8501688 -84.8103991
109 -27.0453055 -35.8501688
110 -25.6469266 -27.0453055
111 -27.2485477 -25.6469266
112 -26.0453055 -27.2485477
113 49.2735340 -26.0453055
114 0.9946915 49.2735340
115 -47.0771970 0.9946915
116 1.3214092 -47.0771970
117 -22.5590544 1.3214092
118 -19.8818330 -22.5590544
119 -11.1606755 -19.8818330
120 -6.5192847 -11.1606755
121 -45.1648419 -6.5192847
122 -33.6869236 -45.1648419
123 33.1537703 -33.6869236
124 8.9944642 33.1537703
125 -3.0055358 8.9944642
126 17.7556187 -3.0055358
127 -13.9216027 17.7556187
128 -43.8778939 -13.9216027
129 -14.9974303 -43.8778939
130 -21.7941881 -14.9974303
131 -4.1925670 -21.7941881
132 10.7993275 -4.1925670
133 3.5960853 10.7993275
134 16.5725308 3.5960853
135 92.5532897 16.5725308
136 -15.1315828 92.5532897
137 -34.5618531 -15.1315828
138 -19.3667164 -34.5618531
139 -12.2471800 -19.3667164
140 -75.6735113 -12.2471800
141 -5.5979110 -75.6735113
142 -36.2434681 -5.5979110
143 -32.9565201 -36.2434681
144 -6.3109630 -32.9565201
145 -23.3190685 -6.3109630
146 -9.3109630 -23.3190685
147 -85.0900760 -9.3109630
148 -62.6998026 -85.0900760
149 -64.7800238 -62.6998026
150 -41.1386330 -64.7800238
151 -31.4174755 -41.1386330
152 -20.4093700 -31.4174755
153 -17.6045067 -20.4093700
154 -15.0028856 -17.6045067
155 -7.8833492 -15.0028856
156 -1.2817281 -7.8833492
157 5.7182719 -1.2817281
158 0.9215141 5.7182719
159 6.9215141 0.9215141
160 5.3198930 6.9215141
161 -20.2898336 5.3198930
162 -4.8914547 -20.2898336
163 -20.6963180 -4.8914547
164 -35.4930758 -20.6963180
165 -8.9024059 -35.4930758
166 39.2833997 -8.9024059
167 34.4604104 39.2833997
168 47.1009584 34.4604104
169 44.3159596 47.1009584
170 17.4080797 44.3159596
171 46.3725172 17.4080797
172 45.4917855 46.3725172
173 44.8796289 45.4917855
174 41.2245873 44.8796289
175 57.0882141 41.2245873
176 50.8852868 57.0882141
177 33.4253214 50.8852868
178 -17.4125035 33.4253214
179 5.3520561 -17.4125035
180 20.9337551 5.3520561
181 14.5091447 20.9337551
182 20.7875657 14.5091447
183 -6.4552541 20.7875657
184 16.2556169 -6.4552541
185 12.1629628 16.2556169
186 -24.6080785 12.1629628
187 -17.4093094 -24.6080785
188 -11.9603593 -17.4093094
189 -5.8945115 -11.9603593
190 -18.0767784 -5.8945115
191 -33.9409522 -18.0767784
192 -25.8922728 -33.9409522
193 -20.8901130 -25.8922728
194 -20.3836886 -20.8901130
195 -12.3159489 -20.3836886
196 -10.9223112 -12.3159489
197 27.7938371 -10.9223112
198 -19.9350235 27.7938371
199 -27.2261399 -19.9350235
200 -16.3282398 -27.2261399
201 -18.0461286 -16.3282398
202 -19.8428864 -18.0461286
203 11.2799946 -19.8428864
204 -25.2294085 11.2799946
205 -24.9371832 -25.2294085
206 -28.9224018 -24.9371832
207 28.0269862 -28.9224018
208 -13.8541669 28.0269862
209 -3.7439615 -13.8541669
210 27.6898076 -3.7439615
211 4.3767793 27.6898076
212 -31.1893552 4.3767793
213 18.4225334 -31.1893552
214 30.2548537 18.4225334
215 34.7983132 30.2548537
216 36.8407788 34.7983132
217 21.0344430 36.8407788
218 30.3427301 21.0344430
219 3.8098062 30.3427301
220 -32.0924158 3.8098062
221 -28.4617139 -32.0924158
222 -33.1868295 -28.4617139
223 -21.3736862 -33.1868295
224 -45.7966518 -21.3736862
225 -13.2782412 -45.7966518
226 -50.0122145 -13.2782412
227 -35.7249985 -50.0122145
228 -16.9808781 -35.7249985
229 -37.1779067 -16.9808781
230 -72.1294954 -37.1779067
231 49.0356959 -72.1294954
232 -84.6935087 49.0356959
233 72.5014258 -84.6935087
234 -81.0318101 72.5014258
235 45.7247179 -81.0318101
236 5.8772932 45.7247179
237 -1.1845008 5.8772932
238 -39.7022598 -1.1845008
239 16.9568356 -39.7022598
240 0.5760466 16.9568356
241 16.1991580 0.5760466
242 -18.3781973 16.1991580
243 30.0619579 -18.3781973
244 27.9670292 30.0619579
245 29.7010354 27.9670292
246 27.8941248 29.7010354
247 25.7613183 27.8941248
248 5.9096675 25.7613183
249 61.9314258 5.9096675
250 87.3603581 61.9314258
251 56.8615242 87.3603581
252 30.9722656 56.8615242
253 -2.4545046 30.9722656
254 -15.1725447 -2.4545046
255 5.1113267 -15.1725447
256 50.1555308 5.1113267
257 39.7574199 50.1555308
258 3.4215246 39.7574199
259 6.5386330 3.4215246
260 9.5214645 6.5386330
261 19.1361857 9.5214645
262 -30.3995876 19.1361857
263 -25.1234410 -30.3995876
264 4.2119183 -25.1234410
265 0.5019836 4.2119183
266 -4.5485875 0.5019836
267 2.9948719 -4.5485875
268 -0.4021302 2.9948719
269 4.3049551 -0.4021302
270 1.8850859 4.3049551
271 -0.8920734 1.8850859
272 -0.1592716 -0.8920734
273 19.8143414 -0.1592716
274 8.7368493 19.8143414
275 -4.3845790 8.7368493
276 -13.2199956 -4.3845790
277 -32.4334856 -13.2199956
278 11.5855769 -32.4334856
279 21.2086883 11.5855769
280 15.4238428 21.2086883
281 12.7121290 15.4238428
282 5.7655495 12.7121290
283 29.1840419 5.7655495
284 -0.7357510 29.1840419
285 73.3407280 -0.7357510
286 83.6711928 73.3407280
287 60.0445228 83.6711928
288 13.7744753 60.0445228
289 -6.6225268 13.7744753
290 21.5877142 -6.6225268
291 -7.9264892 21.5877142
292 21.7572139 -7.9264892
293 -1.4384378 21.7572139
294 -27.8713091 -1.4384378
295 -47.1308212 -27.8713091
296 77.6332783 -47.1308212
297 33.9373202 77.6332783
298 11.7326055 33.9373202
299 13.5355770 11.7326055
300 10.8743217 13.5355770
301 31.4815078 10.8743217
302 37.0032088 31.4815078
303 21.3173458 37.0032088
304 17.6620362 21.3173458
305 9.8535405 17.6620362
306 3.4234995 9.8535405
307 -11.8284220 3.4234995
308 -1.1416226 -11.8284220
309 -20.7370301 -1.1416226
310 -28.8878093 -20.7370301
311 -27.9078271 -28.8878093
312 -26.2110093 -27.9078271
313 -91.8350236 -26.2110093
314 8.0194358 -91.8350236
315 2.6944952 8.0194358
316 -10.2316838 2.6944952
317 -61.6243640 -10.2316838
318 -62.2703426 -61.6243640
319 -58.6625078 -62.2703426
320 -39.9644855 -58.6625078
321 -53.1539236 -39.9644855
322 7.3175890 -53.1539236
323 -67.8319646 7.3175890
324 -35.5094541 -67.8319646
325 -13.8329221 -35.5094541
326 -18.0757807 -13.8329221
327 -6.9965480 -18.0757807
328 -19.0242520 -6.9965480
329 -39.2227669 -19.0242520
330 -32.4787422 -39.2227669
331 147.4678338 -32.4787422
332 -23.9763233 147.4678338
333 -52.5895091 -23.9763233
334 -11.6582063 -52.5895091
335 -9.7616830 -11.6582063
336 -5.0506396 -9.7616830
337 8.4644264 -5.0506396
338 21.6080344 8.4644264
339 51.6833258 21.6080344
340 -37.3804163 51.6833258
341 -52.9906213 -37.3804163
342 21.6612474 -52.9906213
343 162.3328045 21.6612474
344 -15.7982038 162.3328045
345 -17.4656727 -15.7982038
346 -33.2290278 -17.4656727
347 75.1905346 -33.2290278
348 1.1786223 75.1905346
349 -5.2109616 1.1786223
350 10.2797647 -5.2109616
351 32.3826665 10.2797647
352 13.1268059 32.3826665
353 -10.2116708 13.1268059
354 -65.4046068 -10.2116708
355 -35.7939227 -65.4046068
356 -19.3802672 -35.7939227
357 -24.4518137 -19.3802672
358 176.7585202 -24.4518137
359 -148.2630243 176.7585202
360 -39.3515859 -148.2630243
361 -39.0398360 -39.3515859
362 -12.9027097 -39.0398360
363 -14.8540304 -12.9027097
364 -39.8525931 -14.8540304
365 -13.6213788 -39.8525931
366 -27.8336840 -13.6213788
367 -21.1800931 -27.8336840
368 44.8416653 -21.1800931
369 20.2566977 44.8416653
370 -17.0816256 20.2566977
371 -21.1277237 -17.0816256
372 18.1901604 -21.1277237
373 -27.0171356 18.1901604
374 -35.8175258 -27.0171356
375 -6.1075912 -35.8175258
376 10.6142558 -6.1075912
377 4.0083128 10.6142558
378 5.0116772 4.0083128
379 24.8818733 5.0116772
380 -6.9453232 24.8818733
381 -46.1361380 -6.9453232
382 -53.4923193 -46.1361380
383 -56.6998833 -53.4923193
384 -31.2217567 -56.6998833
385 NA -31.2217567
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 2.2072778 248.4408433
[2,] -172.2249231 2.2072778
[3,] -153.3507988 -172.2249231
[4,] -1.8611243 -153.3507988
[5,] -28.5953986 -1.8611243
[6,] -11.0582485 -28.5953986
[7,] -18.1909017 -11.0582485
[8,] 15.9715218 -18.1909017
[9,] 13.7964627 15.9715218
[10,] 99.7359660 13.7964627
[11,] -79.7777035 99.7359660
[12,] -5.0155718 -79.7777035
[13,] -2.0155718 -5.0155718
[14,] 15.6216118 -2.0155718
[15,] -5.9801239 15.6216118
[16,] -26.9956686 -5.9801239
[17,] 86.7843752 -26.9956686
[18,] -75.0300872 86.7843752
[19,] -45.3267877 -75.0300872
[20,] -48.5302979 -45.3267877
[21,] -55.0945772 -48.5302979
[22,] -24.9891131 -55.0945772
[23,] -22.7350316 -24.9891131
[24,] -27.8335730 -22.7350316
[25,] 10.9654981 -27.8335730
[26,] 30.2436511 10.9654981
[27,] 31.7531929 30.2436511
[28,] 10.8709908 31.7531929
[29,] 46.1109803 10.8709908
[30,] 5.1229665 46.1109803
[31,] -18.4095768 5.1229665
[32,] -62.4820013 -18.4095768
[33,] -3.0031853 -62.4820013
[34,] 109.3744863 -3.0031853
[35,] 9.3074946 109.3744863
[36,] 44.1450922 9.3074946
[37,] 109.4428437 44.1450922
[38,] 87.6594708 109.4428437
[39,] 73.1600533 87.6594708
[40,] 90.6839142 73.1600533
[41,] 99.8933701 90.6839142
[42,] 109.1969139 99.8933701
[43,] 63.6062476 109.1969139
[44,] 58.1496114 63.6062476
[45,] 65.8107628 58.1496114
[46,] 18.9700690 65.8107628
[47,] 51.4960897 18.9700690
[48,] 42.7668268 51.4960897
[49,] 156.7117550 42.7668268
[50,] 22.5201030 156.7117550
[51,] 101.8350033 22.5201030
[52,] 12.3289054 101.8350033
[53,] 10.4045057 12.3289054
[54,] 7.1256632 10.4045057
[55,] 5.7272843 7.1256632
[56,] 11.3289054 5.7272843
[57,] 147.9110581 11.3289054
[58,] -66.1171215 147.9110581
[59,] -42.0812909 -66.1171215
[60,] -34.9617545 -42.0812909
[61,] -47.6829120 -34.9617545
[62,] -42.6829120 -47.6829120
[63,] -163.9219848 -42.6829120
[64,] -112.5080766 -163.9219848
[65,] -74.9064555 -112.5080766
[66,] 32.3723870 -74.9064555
[67,] 35.7707659 32.3723870
[68,] 41.1691448 35.7707659
[69,] 44.2886812 41.1691448
[70,] 39.5675237 44.2886812
[71,] 43.7707659 39.5675237
[72,] 42.0098387 43.7707659
[73,] -12.7379548 42.0098387
[74,] -7.0565671 -12.7379548
[75,] -12.4630515 -7.0565671
[76,] -7.8614304 -12.4630515
[77,] -38.8935491 -7.8614304
[78,] 90.9434329 -38.8935491
[79,] -4.2517038 90.9434329
[80,] -20.8496130 -4.2517038
[81,] 27.2699234 -20.8496130
[82,] 25.2699234 27.2699234
[83,] 24.2699234 25.2699234
[84,] 21.7520081 24.2699234
[85,] 30.4731656 21.7520081
[86,] 63.3496901 30.4731656
[87,] 54.6964819 63.3496901
[88,] -44.9409698 54.6964819
[89,] -20.5425909 -44.9409698
[90,] -25.4149490 -20.5425909
[91,] -77.1802698 -25.4149490
[92,] 26.6962546 -77.1802698
[93,] 8.9353274 26.6962546
[94,] 20.1385696 8.9353274
[95,] 2.8994968 20.1385696
[96,] 143.1785666 2.8994968
[97,] -104.6940188 143.1785666
[98,] -100.1761035 -104.6940188
[99,] -115.6940188 -100.1761035
[100,] 4.1506143 -115.6940188
[101,] 1.3457510 4.1506143
[102,] 51.7443571 1.3457510
[103,] 7.2264418 51.7443571
[104,] -6.2483204 7.2264418
[105,] -7.8499415 -6.2483204
[106,] -7.9694779 -7.8499415
[107,] -84.8103991 -7.9694779
[108,] -35.8501688 -84.8103991
[109,] -27.0453055 -35.8501688
[110,] -25.6469266 -27.0453055
[111,] -27.2485477 -25.6469266
[112,] -26.0453055 -27.2485477
[113,] 49.2735340 -26.0453055
[114,] 0.9946915 49.2735340
[115,] -47.0771970 0.9946915
[116,] 1.3214092 -47.0771970
[117,] -22.5590544 1.3214092
[118,] -19.8818330 -22.5590544
[119,] -11.1606755 -19.8818330
[120,] -6.5192847 -11.1606755
[121,] -45.1648419 -6.5192847
[122,] -33.6869236 -45.1648419
[123,] 33.1537703 -33.6869236
[124,] 8.9944642 33.1537703
[125,] -3.0055358 8.9944642
[126,] 17.7556187 -3.0055358
[127,] -13.9216027 17.7556187
[128,] -43.8778939 -13.9216027
[129,] -14.9974303 -43.8778939
[130,] -21.7941881 -14.9974303
[131,] -4.1925670 -21.7941881
[132,] 10.7993275 -4.1925670
[133,] 3.5960853 10.7993275
[134,] 16.5725308 3.5960853
[135,] 92.5532897 16.5725308
[136,] -15.1315828 92.5532897
[137,] -34.5618531 -15.1315828
[138,] -19.3667164 -34.5618531
[139,] -12.2471800 -19.3667164
[140,] -75.6735113 -12.2471800
[141,] -5.5979110 -75.6735113
[142,] -36.2434681 -5.5979110
[143,] -32.9565201 -36.2434681
[144,] -6.3109630 -32.9565201
[145,] -23.3190685 -6.3109630
[146,] -9.3109630 -23.3190685
[147,] -85.0900760 -9.3109630
[148,] -62.6998026 -85.0900760
[149,] -64.7800238 -62.6998026
[150,] -41.1386330 -64.7800238
[151,] -31.4174755 -41.1386330
[152,] -20.4093700 -31.4174755
[153,] -17.6045067 -20.4093700
[154,] -15.0028856 -17.6045067
[155,] -7.8833492 -15.0028856
[156,] -1.2817281 -7.8833492
[157,] 5.7182719 -1.2817281
[158,] 0.9215141 5.7182719
[159,] 6.9215141 0.9215141
[160,] 5.3198930 6.9215141
[161,] -20.2898336 5.3198930
[162,] -4.8914547 -20.2898336
[163,] -20.6963180 -4.8914547
[164,] -35.4930758 -20.6963180
[165,] -8.9024059 -35.4930758
[166,] 39.2833997 -8.9024059
[167,] 34.4604104 39.2833997
[168,] 47.1009584 34.4604104
[169,] 44.3159596 47.1009584
[170,] 17.4080797 44.3159596
[171,] 46.3725172 17.4080797
[172,] 45.4917855 46.3725172
[173,] 44.8796289 45.4917855
[174,] 41.2245873 44.8796289
[175,] 57.0882141 41.2245873
[176,] 50.8852868 57.0882141
[177,] 33.4253214 50.8852868
[178,] -17.4125035 33.4253214
[179,] 5.3520561 -17.4125035
[180,] 20.9337551 5.3520561
[181,] 14.5091447 20.9337551
[182,] 20.7875657 14.5091447
[183,] -6.4552541 20.7875657
[184,] 16.2556169 -6.4552541
[185,] 12.1629628 16.2556169
[186,] -24.6080785 12.1629628
[187,] -17.4093094 -24.6080785
[188,] -11.9603593 -17.4093094
[189,] -5.8945115 -11.9603593
[190,] -18.0767784 -5.8945115
[191,] -33.9409522 -18.0767784
[192,] -25.8922728 -33.9409522
[193,] -20.8901130 -25.8922728
[194,] -20.3836886 -20.8901130
[195,] -12.3159489 -20.3836886
[196,] -10.9223112 -12.3159489
[197,] 27.7938371 -10.9223112
[198,] -19.9350235 27.7938371
[199,] -27.2261399 -19.9350235
[200,] -16.3282398 -27.2261399
[201,] -18.0461286 -16.3282398
[202,] -19.8428864 -18.0461286
[203,] 11.2799946 -19.8428864
[204,] -25.2294085 11.2799946
[205,] -24.9371832 -25.2294085
[206,] -28.9224018 -24.9371832
[207,] 28.0269862 -28.9224018
[208,] -13.8541669 28.0269862
[209,] -3.7439615 -13.8541669
[210,] 27.6898076 -3.7439615
[211,] 4.3767793 27.6898076
[212,] -31.1893552 4.3767793
[213,] 18.4225334 -31.1893552
[214,] 30.2548537 18.4225334
[215,] 34.7983132 30.2548537
[216,] 36.8407788 34.7983132
[217,] 21.0344430 36.8407788
[218,] 30.3427301 21.0344430
[219,] 3.8098062 30.3427301
[220,] -32.0924158 3.8098062
[221,] -28.4617139 -32.0924158
[222,] -33.1868295 -28.4617139
[223,] -21.3736862 -33.1868295
[224,] -45.7966518 -21.3736862
[225,] -13.2782412 -45.7966518
[226,] -50.0122145 -13.2782412
[227,] -35.7249985 -50.0122145
[228,] -16.9808781 -35.7249985
[229,] -37.1779067 -16.9808781
[230,] -72.1294954 -37.1779067
[231,] 49.0356959 -72.1294954
[232,] -84.6935087 49.0356959
[233,] 72.5014258 -84.6935087
[234,] -81.0318101 72.5014258
[235,] 45.7247179 -81.0318101
[236,] 5.8772932 45.7247179
[237,] -1.1845008 5.8772932
[238,] -39.7022598 -1.1845008
[239,] 16.9568356 -39.7022598
[240,] 0.5760466 16.9568356
[241,] 16.1991580 0.5760466
[242,] -18.3781973 16.1991580
[243,] 30.0619579 -18.3781973
[244,] 27.9670292 30.0619579
[245,] 29.7010354 27.9670292
[246,] 27.8941248 29.7010354
[247,] 25.7613183 27.8941248
[248,] 5.9096675 25.7613183
[249,] 61.9314258 5.9096675
[250,] 87.3603581 61.9314258
[251,] 56.8615242 87.3603581
[252,] 30.9722656 56.8615242
[253,] -2.4545046 30.9722656
[254,] -15.1725447 -2.4545046
[255,] 5.1113267 -15.1725447
[256,] 50.1555308 5.1113267
[257,] 39.7574199 50.1555308
[258,] 3.4215246 39.7574199
[259,] 6.5386330 3.4215246
[260,] 9.5214645 6.5386330
[261,] 19.1361857 9.5214645
[262,] -30.3995876 19.1361857
[263,] -25.1234410 -30.3995876
[264,] 4.2119183 -25.1234410
[265,] 0.5019836 4.2119183
[266,] -4.5485875 0.5019836
[267,] 2.9948719 -4.5485875
[268,] -0.4021302 2.9948719
[269,] 4.3049551 -0.4021302
[270,] 1.8850859 4.3049551
[271,] -0.8920734 1.8850859
[272,] -0.1592716 -0.8920734
[273,] 19.8143414 -0.1592716
[274,] 8.7368493 19.8143414
[275,] -4.3845790 8.7368493
[276,] -13.2199956 -4.3845790
[277,] -32.4334856 -13.2199956
[278,] 11.5855769 -32.4334856
[279,] 21.2086883 11.5855769
[280,] 15.4238428 21.2086883
[281,] 12.7121290 15.4238428
[282,] 5.7655495 12.7121290
[283,] 29.1840419 5.7655495
[284,] -0.7357510 29.1840419
[285,] 73.3407280 -0.7357510
[286,] 83.6711928 73.3407280
[287,] 60.0445228 83.6711928
[288,] 13.7744753 60.0445228
[289,] -6.6225268 13.7744753
[290,] 21.5877142 -6.6225268
[291,] -7.9264892 21.5877142
[292,] 21.7572139 -7.9264892
[293,] -1.4384378 21.7572139
[294,] -27.8713091 -1.4384378
[295,] -47.1308212 -27.8713091
[296,] 77.6332783 -47.1308212
[297,] 33.9373202 77.6332783
[298,] 11.7326055 33.9373202
[299,] 13.5355770 11.7326055
[300,] 10.8743217 13.5355770
[301,] 31.4815078 10.8743217
[302,] 37.0032088 31.4815078
[303,] 21.3173458 37.0032088
[304,] 17.6620362 21.3173458
[305,] 9.8535405 17.6620362
[306,] 3.4234995 9.8535405
[307,] -11.8284220 3.4234995
[308,] -1.1416226 -11.8284220
[309,] -20.7370301 -1.1416226
[310,] -28.8878093 -20.7370301
[311,] -27.9078271 -28.8878093
[312,] -26.2110093 -27.9078271
[313,] -91.8350236 -26.2110093
[314,] 8.0194358 -91.8350236
[315,] 2.6944952 8.0194358
[316,] -10.2316838 2.6944952
[317,] -61.6243640 -10.2316838
[318,] -62.2703426 -61.6243640
[319,] -58.6625078 -62.2703426
[320,] -39.9644855 -58.6625078
[321,] -53.1539236 -39.9644855
[322,] 7.3175890 -53.1539236
[323,] -67.8319646 7.3175890
[324,] -35.5094541 -67.8319646
[325,] -13.8329221 -35.5094541
[326,] -18.0757807 -13.8329221
[327,] -6.9965480 -18.0757807
[328,] -19.0242520 -6.9965480
[329,] -39.2227669 -19.0242520
[330,] -32.4787422 -39.2227669
[331,] 147.4678338 -32.4787422
[332,] -23.9763233 147.4678338
[333,] -52.5895091 -23.9763233
[334,] -11.6582063 -52.5895091
[335,] -9.7616830 -11.6582063
[336,] -5.0506396 -9.7616830
[337,] 8.4644264 -5.0506396
[338,] 21.6080344 8.4644264
[339,] 51.6833258 21.6080344
[340,] -37.3804163 51.6833258
[341,] -52.9906213 -37.3804163
[342,] 21.6612474 -52.9906213
[343,] 162.3328045 21.6612474
[344,] -15.7982038 162.3328045
[345,] -17.4656727 -15.7982038
[346,] -33.2290278 -17.4656727
[347,] 75.1905346 -33.2290278
[348,] 1.1786223 75.1905346
[349,] -5.2109616 1.1786223
[350,] 10.2797647 -5.2109616
[351,] 32.3826665 10.2797647
[352,] 13.1268059 32.3826665
[353,] -10.2116708 13.1268059
[354,] -65.4046068 -10.2116708
[355,] -35.7939227 -65.4046068
[356,] -19.3802672 -35.7939227
[357,] -24.4518137 -19.3802672
[358,] 176.7585202 -24.4518137
[359,] -148.2630243 176.7585202
[360,] -39.3515859 -148.2630243
[361,] -39.0398360 -39.3515859
[362,] -12.9027097 -39.0398360
[363,] -14.8540304 -12.9027097
[364,] -39.8525931 -14.8540304
[365,] -13.6213788 -39.8525931
[366,] -27.8336840 -13.6213788
[367,] -21.1800931 -27.8336840
[368,] 44.8416653 -21.1800931
[369,] 20.2566977 44.8416653
[370,] -17.0816256 20.2566977
[371,] -21.1277237 -17.0816256
[372,] 18.1901604 -21.1277237
[373,] -27.0171356 18.1901604
[374,] -35.8175258 -27.0171356
[375,] -6.1075912 -35.8175258
[376,] 10.6142558 -6.1075912
[377,] 4.0083128 10.6142558
[378,] 5.0116772 4.0083128
[379,] 24.8818733 5.0116772
[380,] -6.9453232 24.8818733
[381,] -46.1361380 -6.9453232
[382,] -53.4923193 -46.1361380
[383,] -56.6998833 -53.4923193
[384,] -31.2217567 -56.6998833
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 2.2072778 248.4408433
2 -172.2249231 2.2072778
3 -153.3507988 -172.2249231
4 -1.8611243 -153.3507988
5 -28.5953986 -1.8611243
6 -11.0582485 -28.5953986
7 -18.1909017 -11.0582485
8 15.9715218 -18.1909017
9 13.7964627 15.9715218
10 99.7359660 13.7964627
11 -79.7777035 99.7359660
12 -5.0155718 -79.7777035
13 -2.0155718 -5.0155718
14 15.6216118 -2.0155718
15 -5.9801239 15.6216118
16 -26.9956686 -5.9801239
17 86.7843752 -26.9956686
18 -75.0300872 86.7843752
19 -45.3267877 -75.0300872
20 -48.5302979 -45.3267877
21 -55.0945772 -48.5302979
22 -24.9891131 -55.0945772
23 -22.7350316 -24.9891131
24 -27.8335730 -22.7350316
25 10.9654981 -27.8335730
26 30.2436511 10.9654981
27 31.7531929 30.2436511
28 10.8709908 31.7531929
29 46.1109803 10.8709908
30 5.1229665 46.1109803
31 -18.4095768 5.1229665
32 -62.4820013 -18.4095768
33 -3.0031853 -62.4820013
34 109.3744863 -3.0031853
35 9.3074946 109.3744863
36 44.1450922 9.3074946
37 109.4428437 44.1450922
38 87.6594708 109.4428437
39 73.1600533 87.6594708
40 90.6839142 73.1600533
41 99.8933701 90.6839142
42 109.1969139 99.8933701
43 63.6062476 109.1969139
44 58.1496114 63.6062476
45 65.8107628 58.1496114
46 18.9700690 65.8107628
47 51.4960897 18.9700690
48 42.7668268 51.4960897
49 156.7117550 42.7668268
50 22.5201030 156.7117550
51 101.8350033 22.5201030
52 12.3289054 101.8350033
53 10.4045057 12.3289054
54 7.1256632 10.4045057
55 5.7272843 7.1256632
56 11.3289054 5.7272843
57 147.9110581 11.3289054
58 -66.1171215 147.9110581
59 -42.0812909 -66.1171215
60 -34.9617545 -42.0812909
61 -47.6829120 -34.9617545
62 -42.6829120 -47.6829120
63 -163.9219848 -42.6829120
64 -112.5080766 -163.9219848
65 -74.9064555 -112.5080766
66 32.3723870 -74.9064555
67 35.7707659 32.3723870
68 41.1691448 35.7707659
69 44.2886812 41.1691448
70 39.5675237 44.2886812
71 43.7707659 39.5675237
72 42.0098387 43.7707659
73 -12.7379548 42.0098387
74 -7.0565671 -12.7379548
75 -12.4630515 -7.0565671
76 -7.8614304 -12.4630515
77 -38.8935491 -7.8614304
78 90.9434329 -38.8935491
79 -4.2517038 90.9434329
80 -20.8496130 -4.2517038
81 27.2699234 -20.8496130
82 25.2699234 27.2699234
83 24.2699234 25.2699234
84 21.7520081 24.2699234
85 30.4731656 21.7520081
86 63.3496901 30.4731656
87 54.6964819 63.3496901
88 -44.9409698 54.6964819
89 -20.5425909 -44.9409698
90 -25.4149490 -20.5425909
91 -77.1802698 -25.4149490
92 26.6962546 -77.1802698
93 8.9353274 26.6962546
94 20.1385696 8.9353274
95 2.8994968 20.1385696
96 143.1785666 2.8994968
97 -104.6940188 143.1785666
98 -100.1761035 -104.6940188
99 -115.6940188 -100.1761035
100 4.1506143 -115.6940188
101 1.3457510 4.1506143
102 51.7443571 1.3457510
103 7.2264418 51.7443571
104 -6.2483204 7.2264418
105 -7.8499415 -6.2483204
106 -7.9694779 -7.8499415
107 -84.8103991 -7.9694779
108 -35.8501688 -84.8103991
109 -27.0453055 -35.8501688
110 -25.6469266 -27.0453055
111 -27.2485477 -25.6469266
112 -26.0453055 -27.2485477
113 49.2735340 -26.0453055
114 0.9946915 49.2735340
115 -47.0771970 0.9946915
116 1.3214092 -47.0771970
117 -22.5590544 1.3214092
118 -19.8818330 -22.5590544
119 -11.1606755 -19.8818330
120 -6.5192847 -11.1606755
121 -45.1648419 -6.5192847
122 -33.6869236 -45.1648419
123 33.1537703 -33.6869236
124 8.9944642 33.1537703
125 -3.0055358 8.9944642
126 17.7556187 -3.0055358
127 -13.9216027 17.7556187
128 -43.8778939 -13.9216027
129 -14.9974303 -43.8778939
130 -21.7941881 -14.9974303
131 -4.1925670 -21.7941881
132 10.7993275 -4.1925670
133 3.5960853 10.7993275
134 16.5725308 3.5960853
135 92.5532897 16.5725308
136 -15.1315828 92.5532897
137 -34.5618531 -15.1315828
138 -19.3667164 -34.5618531
139 -12.2471800 -19.3667164
140 -75.6735113 -12.2471800
141 -5.5979110 -75.6735113
142 -36.2434681 -5.5979110
143 -32.9565201 -36.2434681
144 -6.3109630 -32.9565201
145 -23.3190685 -6.3109630
146 -9.3109630 -23.3190685
147 -85.0900760 -9.3109630
148 -62.6998026 -85.0900760
149 -64.7800238 -62.6998026
150 -41.1386330 -64.7800238
151 -31.4174755 -41.1386330
152 -20.4093700 -31.4174755
153 -17.6045067 -20.4093700
154 -15.0028856 -17.6045067
155 -7.8833492 -15.0028856
156 -1.2817281 -7.8833492
157 5.7182719 -1.2817281
158 0.9215141 5.7182719
159 6.9215141 0.9215141
160 5.3198930 6.9215141
161 -20.2898336 5.3198930
162 -4.8914547 -20.2898336
163 -20.6963180 -4.8914547
164 -35.4930758 -20.6963180
165 -8.9024059 -35.4930758
166 39.2833997 -8.9024059
167 34.4604104 39.2833997
168 47.1009584 34.4604104
169 44.3159596 47.1009584
170 17.4080797 44.3159596
171 46.3725172 17.4080797
172 45.4917855 46.3725172
173 44.8796289 45.4917855
174 41.2245873 44.8796289
175 57.0882141 41.2245873
176 50.8852868 57.0882141
177 33.4253214 50.8852868
178 -17.4125035 33.4253214
179 5.3520561 -17.4125035
180 20.9337551 5.3520561
181 14.5091447 20.9337551
182 20.7875657 14.5091447
183 -6.4552541 20.7875657
184 16.2556169 -6.4552541
185 12.1629628 16.2556169
186 -24.6080785 12.1629628
187 -17.4093094 -24.6080785
188 -11.9603593 -17.4093094
189 -5.8945115 -11.9603593
190 -18.0767784 -5.8945115
191 -33.9409522 -18.0767784
192 -25.8922728 -33.9409522
193 -20.8901130 -25.8922728
194 -20.3836886 -20.8901130
195 -12.3159489 -20.3836886
196 -10.9223112 -12.3159489
197 27.7938371 -10.9223112
198 -19.9350235 27.7938371
199 -27.2261399 -19.9350235
200 -16.3282398 -27.2261399
201 -18.0461286 -16.3282398
202 -19.8428864 -18.0461286
203 11.2799946 -19.8428864
204 -25.2294085 11.2799946
205 -24.9371832 -25.2294085
206 -28.9224018 -24.9371832
207 28.0269862 -28.9224018
208 -13.8541669 28.0269862
209 -3.7439615 -13.8541669
210 27.6898076 -3.7439615
211 4.3767793 27.6898076
212 -31.1893552 4.3767793
213 18.4225334 -31.1893552
214 30.2548537 18.4225334
215 34.7983132 30.2548537
216 36.8407788 34.7983132
217 21.0344430 36.8407788
218 30.3427301 21.0344430
219 3.8098062 30.3427301
220 -32.0924158 3.8098062
221 -28.4617139 -32.0924158
222 -33.1868295 -28.4617139
223 -21.3736862 -33.1868295
224 -45.7966518 -21.3736862
225 -13.2782412 -45.7966518
226 -50.0122145 -13.2782412
227 -35.7249985 -50.0122145
228 -16.9808781 -35.7249985
229 -37.1779067 -16.9808781
230 -72.1294954 -37.1779067
231 49.0356959 -72.1294954
232 -84.6935087 49.0356959
233 72.5014258 -84.6935087
234 -81.0318101 72.5014258
235 45.7247179 -81.0318101
236 5.8772932 45.7247179
237 -1.1845008 5.8772932
238 -39.7022598 -1.1845008
239 16.9568356 -39.7022598
240 0.5760466 16.9568356
241 16.1991580 0.5760466
242 -18.3781973 16.1991580
243 30.0619579 -18.3781973
244 27.9670292 30.0619579
245 29.7010354 27.9670292
246 27.8941248 29.7010354
247 25.7613183 27.8941248
248 5.9096675 25.7613183
249 61.9314258 5.9096675
250 87.3603581 61.9314258
251 56.8615242 87.3603581
252 30.9722656 56.8615242
253 -2.4545046 30.9722656
254 -15.1725447 -2.4545046
255 5.1113267 -15.1725447
256 50.1555308 5.1113267
257 39.7574199 50.1555308
258 3.4215246 39.7574199
259 6.5386330 3.4215246
260 9.5214645 6.5386330
261 19.1361857 9.5214645
262 -30.3995876 19.1361857
263 -25.1234410 -30.3995876
264 4.2119183 -25.1234410
265 0.5019836 4.2119183
266 -4.5485875 0.5019836
267 2.9948719 -4.5485875
268 -0.4021302 2.9948719
269 4.3049551 -0.4021302
270 1.8850859 4.3049551
271 -0.8920734 1.8850859
272 -0.1592716 -0.8920734
273 19.8143414 -0.1592716
274 8.7368493 19.8143414
275 -4.3845790 8.7368493
276 -13.2199956 -4.3845790
277 -32.4334856 -13.2199956
278 11.5855769 -32.4334856
279 21.2086883 11.5855769
280 15.4238428 21.2086883
281 12.7121290 15.4238428
282 5.7655495 12.7121290
283 29.1840419 5.7655495
284 -0.7357510 29.1840419
285 73.3407280 -0.7357510
286 83.6711928 73.3407280
287 60.0445228 83.6711928
288 13.7744753 60.0445228
289 -6.6225268 13.7744753
290 21.5877142 -6.6225268
291 -7.9264892 21.5877142
292 21.7572139 -7.9264892
293 -1.4384378 21.7572139
294 -27.8713091 -1.4384378
295 -47.1308212 -27.8713091
296 77.6332783 -47.1308212
297 33.9373202 77.6332783
298 11.7326055 33.9373202
299 13.5355770 11.7326055
300 10.8743217 13.5355770
301 31.4815078 10.8743217
302 37.0032088 31.4815078
303 21.3173458 37.0032088
304 17.6620362 21.3173458
305 9.8535405 17.6620362
306 3.4234995 9.8535405
307 -11.8284220 3.4234995
308 -1.1416226 -11.8284220
309 -20.7370301 -1.1416226
310 -28.8878093 -20.7370301
311 -27.9078271 -28.8878093
312 -26.2110093 -27.9078271
313 -91.8350236 -26.2110093
314 8.0194358 -91.8350236
315 2.6944952 8.0194358
316 -10.2316838 2.6944952
317 -61.6243640 -10.2316838
318 -62.2703426 -61.6243640
319 -58.6625078 -62.2703426
320 -39.9644855 -58.6625078
321 -53.1539236 -39.9644855
322 7.3175890 -53.1539236
323 -67.8319646 7.3175890
324 -35.5094541 -67.8319646
325 -13.8329221 -35.5094541
326 -18.0757807 -13.8329221
327 -6.9965480 -18.0757807
328 -19.0242520 -6.9965480
329 -39.2227669 -19.0242520
330 -32.4787422 -39.2227669
331 147.4678338 -32.4787422
332 -23.9763233 147.4678338
333 -52.5895091 -23.9763233
334 -11.6582063 -52.5895091
335 -9.7616830 -11.6582063
336 -5.0506396 -9.7616830
337 8.4644264 -5.0506396
338 21.6080344 8.4644264
339 51.6833258 21.6080344
340 -37.3804163 51.6833258
341 -52.9906213 -37.3804163
342 21.6612474 -52.9906213
343 162.3328045 21.6612474
344 -15.7982038 162.3328045
345 -17.4656727 -15.7982038
346 -33.2290278 -17.4656727
347 75.1905346 -33.2290278
348 1.1786223 75.1905346
349 -5.2109616 1.1786223
350 10.2797647 -5.2109616
351 32.3826665 10.2797647
352 13.1268059 32.3826665
353 -10.2116708 13.1268059
354 -65.4046068 -10.2116708
355 -35.7939227 -65.4046068
356 -19.3802672 -35.7939227
357 -24.4518137 -19.3802672
358 176.7585202 -24.4518137
359 -148.2630243 176.7585202
360 -39.3515859 -148.2630243
361 -39.0398360 -39.3515859
362 -12.9027097 -39.0398360
363 -14.8540304 -12.9027097
364 -39.8525931 -14.8540304
365 -13.6213788 -39.8525931
366 -27.8336840 -13.6213788
367 -21.1800931 -27.8336840
368 44.8416653 -21.1800931
369 20.2566977 44.8416653
370 -17.0816256 20.2566977
371 -21.1277237 -17.0816256
372 18.1901604 -21.1277237
373 -27.0171356 18.1901604
374 -35.8175258 -27.0171356
375 -6.1075912 -35.8175258
376 10.6142558 -6.1075912
377 4.0083128 10.6142558
378 5.0116772 4.0083128
379 24.8818733 5.0116772
380 -6.9453232 24.8818733
381 -46.1361380 -6.9453232
382 -53.4923193 -46.1361380
383 -56.6998833 -53.4923193
384 -31.2217567 -56.6998833
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/7bs7h1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/8o8fg1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/html/rcomp/tmp/9br4w1199799857.ps",horizontal=F,pagecentre=F,paper="special",width=14.222222222222,height=10.666666666667)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> load(file='/var/www/html/rcomp/createtable')
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/10ensn1199799857.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/11uyrp1199799857.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/12z4pm1199799857.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/13mukx1199799857.tab")
>
> system("convert tmp/1c4su1199799857.ps tmp/1c4su1199799857.png")
> system("convert tmp/25nxz1199799857.ps tmp/25nxz1199799857.png")
> system("convert tmp/3wn5w1199799857.ps tmp/3wn5w1199799857.png")
> system("convert tmp/4wx1t1199799857.ps tmp/4wx1t1199799857.png")
> system("convert tmp/506m21199799857.ps tmp/506m21199799857.png")
> system("convert tmp/6kmav1199799857.ps tmp/6kmav1199799857.png")
> system("convert tmp/7bs7h1199799857.ps tmp/7bs7h1199799857.png")
> system("convert tmp/8o8fg1199799857.ps tmp/8o8fg1199799857.png")
> system("convert tmp/9br4w1199799857.ps tmp/9br4w1199799857.png")
>
>
> proc.time()
user system elapsed
6.140 2.549 7.056