R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
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(0
+ ,1281151
+ ,13
+ ,0
+ ,1281814
+ ,1298756
+ ,1707752
+ ,1867943
+ ,0
+ ,1164976
+ ,14
+ ,0
+ ,1281151
+ ,1281814
+ ,1298756
+ ,1707752
+ ,0
+ ,1454329
+ ,15
+ ,0
+ ,1164976
+ ,1281151
+ ,1281814
+ ,1298756
+ ,0
+ ,1645288
+ ,16
+ ,0
+ ,1454329
+ ,1164976
+ ,1281151
+ ,1281814
+ ,0
+ ,1817743
+ ,17
+ ,0
+ ,1645288
+ ,1454329
+ ,1164976
+ ,1281151
+ ,0
+ ,1895785
+ ,18
+ ,0
+ ,1817743
+ ,1645288
+ ,1454329
+ ,1164976
+ ,0
+ ,2236311
+ ,19
+ ,0
+ ,1895785
+ ,1817743
+ ,1645288
+ ,1454329
+ ,0
+ ,2295951
+ ,20
+ ,0
+ ,2236311
+ ,1895785
+ ,1817743
+ ,1645288
+ ,0
+ ,2087315
+ ,21
+ ,0
+ ,2295951
+ ,2236311
+ ,1895785
+ ,1817743
+ ,0
+ ,1980891
+ ,22
+ ,0
+ ,2087315
+ ,2295951
+ ,2236311
+ ,1895785
+ ,0
+ ,1465446
+ ,23
+ ,0
+ ,1980891
+ ,2087315
+ ,2295951
+ ,2236311
+ ,0
+ ,1445026
+ ,24
+ ,0
+ ,1465446
+ ,1980891
+ ,2087315
+ ,2295951
+ ,0
+ ,1488120
+ ,25
+ ,0
+ ,1445026
+ ,1465446
+ ,1980891
+ ,2087315
+ ,0
+ ,1338333
+ ,26
+ ,0
+ ,1488120
+ ,1445026
+ ,1465446
+ ,1980891
+ ,0
+ ,1715789
+ ,27
+ ,0
+ ,1338333
+ ,1488120
+ ,1445026
+ ,1465446
+ ,0
+ ,1806090
+ ,28
+ ,0
+ ,1715789
+ ,1338333
+ ,1488120
+ ,1445026
+ ,0
+ ,2083316
+ ,29
+ ,0
+ ,1806090
+ ,1715789
+ ,1338333
+ ,1488120
+ ,0
+ ,2092278
+ ,30
+ ,0
+ ,2083316
+ ,1806090
+ ,1715789
+ ,1338333
+ ,0
+ ,2430800
+ ,31
+ ,0
+ ,2092278
+ ,2083316
+ ,1806090
+ ,1715789
+ ,0
+ ,2424894
+ ,32
+ ,0
+ ,2430800
+ ,2092278
+ ,2083316
+ ,1806090
+ ,0
+ ,2299016
+ ,33
+ ,0
+ ,2424894
+ ,2430800
+ ,2092278
+ ,2083316
+ ,0
+ ,2130688
+ ,34
+ ,0
+ ,2299016
+ ,2424894
+ ,2430800
+ ,2092278
+ ,0
+ ,1652221
+ ,35
+ ,0
+ ,2130688
+ ,2299016
+ ,2424894
+ ,2430800
+ ,0
+ ,1608162
+ ,36
+ ,0
+ ,1652221
+ ,2130688
+ ,2299016
+ ,2424894
+ ,0
+ ,1647074
+ ,37
+ ,0
+ ,1608162
+ ,1652221
+ ,2130688
+ ,2299016
+ ,0
+ ,1479691
+ ,38
+ ,0
+ ,1647074
+ ,1608162
+ ,1652221
+ ,2130688
+ ,0
+ ,1884978
+ ,39
+ ,0
+ ,1479691
+ ,1647074
+ ,1608162
+ ,1652221
+ ,0
+ ,2007898
+ ,40
+ ,0
+ ,1884978
+ ,1479691
+ ,1647074
+ ,1608162
+ ,0
+ ,2208954
+ ,41
+ ,0
+ ,2007898
+ ,1884978
+ ,1479691
+ ,1647074
+ ,0
+ ,2217164
+ ,42
+ ,0
+ ,2208954
+ ,2007898
+ ,1884978
+ ,1479691
+ ,0
+ ,2534291
+ ,43
+ ,0
+ ,2217164
+ ,2208954
+ ,2007898
+ ,1884978
+ ,0
+ ,2560312
+ ,44
+ ,0
+ ,2534291
+ ,2217164
+ ,2208954
+ ,2007898
+ ,0
+ ,2429069
+ ,45
+ ,0
+ ,2560312
+ ,2534291
+ ,2217164
+ ,2208954
+ ,0
+ ,2315077
+ ,46
+ ,0
+ ,2429069
+ ,2560312
+ ,2534291
+ ,2217164
+ ,0
+ ,1799608
+ ,47
+ ,0
+ ,2315077
+ ,2429069
+ ,2560312
+ ,2534291
+ ,0
+ ,1772590
+ ,48
+ ,0
+ ,1799608
+ ,2315077
+ ,2429069
+ ,2560312
+ ,0
+ ,1744799
+ ,49
+ ,0
+ ,1772590
+ ,1799608
+ ,2315077
+ ,2429069
+ ,0
+ ,1659093
+ ,50
+ ,0
+ ,1744799
+ ,1772590
+ ,1799608
+ ,2315077
+ ,0
+ ,2099821
+ ,51
+ ,0
+ ,1659093
+ ,1744799
+ ,1772590
+ ,1799608
+ ,0
+ ,2135736
+ ,52
+ ,0
+ ,2099821
+ ,1659093
+ ,1744799
+ ,1772590
+ ,0
+ ,2427894
+ ,53
+ ,0
+ ,2135736
+ ,2099821
+ ,1659093
+ ,1744799
+ ,0
+ ,2468882
+ ,54
+ ,0
+ ,2427894
+ ,2135736
+ ,2099821
+ ,1659093
+ ,0
+ ,2703217
+ ,55
+ ,0
+ ,2468882
+ ,2427894
+ ,2135736
+ ,2099821
+ ,0
+ ,2766841
+ ,56
+ ,0
+ ,2703217
+ ,2468882
+ ,2427894
+ ,2135736
+ ,0
+ ,2655236
+ ,57
+ ,0
+ ,2766841
+ ,2703217
+ ,2468882
+ ,2427894
+ ,0
+ ,2550373
+ ,58
+ ,0
+ ,2655236
+ ,2766841
+ ,2703217
+ ,2468882
+ ,0
+ ,2052097
+ ,59
+ ,0
+ ,2550373
+ ,2655236
+ ,2766841
+ ,2703217
+ ,0
+ ,1998055
+ ,60
+ ,0
+ ,2052097
+ ,2550373
+ ,2655236
+ ,2766841
+ ,0
+ ,1920748
+ ,61
+ ,0
+ ,1998055
+ ,2052097
+ ,2550373
+ ,2655236
+ ,0
+ ,1876694
+ ,62
+ ,0
+ ,1920748
+ ,1998055
+ ,2052097
+ ,2550373
+ ,0
+ ,2380930
+ ,63
+ ,0
+ ,1876694
+ ,1920748
+ ,1998055
+ ,2052097
+ ,0
+ ,2467402
+ ,64
+ ,0
+ ,2380930
+ ,1876694
+ ,1920748
+ ,1998055
+ ,0
+ ,2770771
+ ,65
+ ,0
+ ,2467402
+ ,2380930
+ ,1876694
+ ,1920748
+ ,0
+ ,2781340
+ ,66
+ ,0
+ ,2770771
+ ,2467402
+ ,2380930
+ ,1876694
+ ,0
+ ,3143926
+ ,67
+ ,0
+ ,2781340
+ ,2770771
+ ,2467402
+ ,2380930
+ ,0
+ ,3172235
+ ,68
+ ,0
+ ,3143926
+ ,2781340
+ ,2770771
+ ,2467402
+ ,0
+ ,2952540
+ ,69
+ ,0
+ ,3172235
+ ,3143926
+ ,2781340
+ ,2770771
+ ,0
+ ,2920877
+ ,70
+ ,0
+ ,2952540
+ ,3172235
+ ,3143926
+ ,2781340
+ ,0
+ ,2384552
+ ,71
+ ,0
+ ,2920877
+ ,2952540
+ ,3172235
+ ,3143926
+ ,0
+ ,2248987
+ ,72
+ ,0
+ ,2384552
+ ,2920877
+ ,2952540
+ ,3172235
+ ,0
+ ,2208616
+ ,73
+ ,0
+ ,2248987
+ ,2384552
+ ,2920877
+ ,2952540
+ ,0
+ ,2178756
+ ,74
+ ,0
+ ,2208616
+ ,2248987
+ ,2384552
+ ,2920877
+ ,0
+ ,2632870
+ ,75
+ ,0
+ ,2178756
+ ,2208616
+ ,2248987
+ ,2384552
+ ,0
+ ,2706905
+ ,76
+ ,0
+ ,2632870
+ ,2178756
+ ,2208616
+ ,2248987
+ ,0
+ ,3029745
+ ,77
+ ,0
+ ,2706905
+ ,2632870
+ ,2178756
+ ,2208616
+ ,0
+ ,3015402
+ ,78
+ ,0
+ ,3029745
+ ,2706905
+ ,2632870
+ ,2178756
+ ,0
+ ,3391414
+ ,79
+ ,0
+ ,3015402
+ ,3029745
+ ,2706905
+ ,2632870
+ ,0
+ ,3507805
+ ,80
+ ,0
+ ,3391414
+ ,3015402
+ ,3029745
+ ,2706905
+ ,0
+ ,3177852
+ ,81
+ ,0
+ ,3507805
+ ,3391414
+ ,3015402
+ ,3029745
+ ,0
+ ,3142961
+ ,82
+ ,0
+ ,3177852
+ ,3507805
+ ,3391414
+ ,3015402
+ ,0
+ ,2545815
+ ,83
+ ,0
+ ,3142961
+ ,3177852
+ ,3507805
+ ,3391414
+ ,0
+ ,2414007
+ ,84
+ ,0
+ ,2545815
+ ,3142961
+ ,3177852
+ ,3507805
+ ,0
+ ,2372578
+ ,85
+ ,0
+ ,2414007
+ ,2545815
+ ,3142961
+ ,3177852
+ ,0
+ ,2332664
+ ,86
+ ,0
+ ,2372578
+ ,2414007
+ ,2545815
+ ,3142961
+ ,0
+ ,2825328
+ ,87
+ ,0
+ ,2332664
+ ,2372578
+ ,2414007
+ ,2545815
+ ,0
+ ,2901478
+ ,88
+ ,0
+ ,2825328
+ ,2332664
+ ,2372578
+ ,2414007
+ ,0
+ ,3263955
+ ,89
+ ,0
+ ,2901478
+ ,2825328
+ ,2332664
+ ,2372578
+ ,0
+ ,3226738
+ ,90
+ ,0
+ ,3263955
+ ,2901478
+ ,2825328
+ ,2332664
+ ,0
+ ,3610786
+ ,91
+ ,0
+ ,3226738
+ ,3263955
+ ,2901478
+ ,2825328
+ ,0
+ ,3709274
+ ,92
+ ,0
+ ,3610786
+ ,3226738
+ ,3263955
+ ,2901478
+ ,0
+ ,3467185
+ ,93
+ ,0
+ ,3709274
+ ,3610786
+ ,3226738
+ ,3263955
+ ,0
+ ,3449646
+ ,94
+ ,0
+ ,3467185
+ ,3709274
+ ,3610786
+ ,3226738
+ ,0
+ ,2802951
+ ,95
+ ,0
+ ,3449646
+ ,3467185
+ ,3709274
+ ,3610786
+ ,0
+ ,2462530
+ ,96
+ ,0
+ ,2802951
+ ,3449646
+ ,3467185
+ ,3709274
+ ,0
+ ,2490645
+ ,97
+ ,0
+ ,2462530
+ ,2802951
+ ,3449646
+ ,3467185
+ ,0
+ ,2561520
+ ,98
+ ,0
+ ,2490645
+ ,2462530
+ ,2802951
+ ,3449646
+ ,0
+ ,3067554
+ ,99
+ ,0
+ ,2561520
+ ,2490645
+ ,2462530
+ ,2802951
+ ,0
+ ,3226951
+ ,100
+ ,0
+ ,3067554
+ ,2561520
+ ,2490645
+ ,2462530
+ ,0
+ ,3546493
+ ,101
+ ,0
+ ,3226951
+ ,3067554
+ ,2561520
+ ,2490645
+ ,0
+ ,3492787
+ ,102
+ ,0
+ ,3546493
+ ,3226951
+ ,3067554
+ ,2561520
+ ,0
+ ,3952263
+ ,103
+ ,0
+ ,3492787
+ ,3546493
+ ,3226951
+ ,3067554
+ ,0
+ ,3932072
+ ,104
+ ,0
+ ,3952263
+ ,3492787
+ ,3546493
+ ,3226951
+ ,0
+ ,3720284
+ ,105
+ ,0
+ ,3932072
+ ,3952263
+ ,3492787
+ ,3546493
+ ,0
+ ,3651555
+ ,106
+ ,0
+ ,3720284
+ ,3932072
+ ,3952263
+ ,3492787
+ ,0
+ ,2914972
+ ,107
+ ,0
+ ,3651555
+ ,3720284
+ ,3932072
+ ,3952263
+ ,0
+ ,2713514
+ ,108
+ ,0
+ ,2914972
+ ,3651555
+ ,3720284
+ ,3932072
+ ,0
+ ,2703997
+ ,109
+ ,0
+ ,2713514
+ ,2914972
+ ,3651555
+ ,3720284
+ ,0
+ ,2591373
+ ,110
+ ,0
+ ,2703997
+ ,2713514
+ ,2914972
+ ,3651555
+ ,0
+ ,3163748
+ ,111
+ ,0
+ ,2591373
+ ,2703997
+ ,2713514
+ ,2914972
+ ,0
+ ,3355137
+ ,112
+ ,0
+ ,3163748
+ ,2591373
+ ,2703997
+ ,2713514
+ ,0
+ ,3613702
+ ,113
+ ,0
+ ,3355137
+ ,3163748
+ ,2591373
+ ,2703997
+ ,0
+ ,3686773
+ ,114
+ ,0
+ ,3613702
+ ,3355137
+ ,3163748
+ ,2591373
+ ,0
+ ,4098716
+ ,115
+ ,0
+ ,3686773
+ ,3613702
+ ,3355137
+ ,3163748
+ ,0
+ ,4063517
+ ,116
+ ,0
+ ,4098716
+ ,3686773
+ ,3613702
+ ,3355137
+ ,1
+ ,3551489
+ ,117
+ ,117
+ ,4063517
+ ,4098716
+ ,3686773
+ ,3613702
+ ,1
+ ,3226663
+ ,118
+ ,118
+ ,3551489
+ ,4063517
+ ,4098716
+ ,3686773
+ ,1
+ ,2656842
+ ,119
+ ,119
+ ,3226663
+ ,3551489
+ ,4063517
+ ,4098716
+ ,1
+ ,2597484
+ ,120
+ ,120
+ ,2656842
+ ,3226663
+ ,3551489
+ ,4063517
+ ,1
+ ,2572399
+ ,121
+ ,121
+ ,2597484
+ ,2656842
+ ,3226663
+ ,3551489
+ ,1
+ ,2596631
+ ,122
+ ,122
+ ,2572399
+ ,2597484
+ ,2656842
+ ,3226663
+ ,1
+ ,3165225
+ ,123
+ ,123
+ ,2596631
+ ,2572399
+ ,2597484
+ ,2656842
+ ,1
+ ,3303145
+ ,124
+ ,124
+ ,3165225
+ ,2596631
+ ,2572399
+ ,2597484
+ ,1
+ ,3698247
+ ,125
+ ,125
+ ,3303145
+ ,3165225
+ ,2596631
+ ,2572399
+ ,1
+ ,3668631
+ ,126
+ ,126
+ ,3698247
+ ,3303145
+ ,3165225
+ ,2596631
+ ,1
+ ,4130433
+ ,127
+ ,127
+ ,3668631
+ ,3698247
+ ,3303145
+ ,3165225
+ ,1
+ ,4131400
+ ,128
+ ,128
+ ,4130433
+ ,3668631
+ ,3698247
+ ,3303145
+ ,1
+ ,3864358
+ ,129
+ ,129
+ ,4131400
+ ,4130433
+ ,3668631
+ ,3698247
+ ,1
+ ,3721110
+ ,130
+ ,130
+ ,3864358
+ ,4131400
+ ,4130433
+ ,3668631
+ ,1
+ ,2892532
+ ,131
+ ,131
+ ,3721110
+ ,3864358
+ ,4131400
+ ,4130433
+ ,1
+ ,2843451
+ ,132
+ ,132
+ ,2892532
+ ,3721110
+ ,3864358
+ ,4131400
+ ,1
+ ,2747502
+ ,133
+ ,133
+ ,2843451
+ ,2892532
+ ,3721110
+ ,3864358
+ ,1
+ ,2668775
+ ,134
+ ,134
+ ,2747502
+ ,2843451
+ ,2892532
+ ,3721110
+ ,1
+ ,3018602
+ ,135
+ ,135
+ ,2668775
+ ,2747502
+ ,2843451
+ ,2892532
+ ,1
+ ,3013392
+ ,136
+ ,136
+ ,3018602
+ ,2668775
+ ,2747502
+ ,2843451
+ ,1
+ ,3393657
+ ,137
+ ,137
+ ,3013392
+ ,3018602
+ ,2668775
+ ,2747502
+ ,1
+ ,3544233
+ ,138
+ ,138
+ ,3393657
+ ,3013392
+ ,3018602
+ ,2668775
+ ,1
+ ,4075832
+ ,139
+ ,139
+ ,3544233
+ ,3393657
+ ,3013392
+ ,3018602
+ ,1
+ ,4032923
+ ,140
+ ,140
+ ,4075832
+ ,3544233
+ ,3393657
+ ,3013392
+ ,1
+ ,3734509
+ ,141
+ ,141
+ ,4032923
+ ,4075832
+ ,3544233
+ ,3393657
+ ,1
+ ,3761285
+ ,142
+ ,142
+ ,3734509
+ ,4032923
+ ,4075832
+ ,3544233
+ ,1
+ ,2970090
+ ,143
+ ,143
+ ,3761285
+ ,3734509
+ ,4032923
+ ,4075832
+ ,1
+ ,2847849
+ ,144
+ ,144
+ ,2970090
+ ,3761285
+ ,3734509
+ ,4032923
+ ,1
+ ,2741680
+ ,145
+ ,145
+ ,2847849
+ ,2970090
+ ,3761285
+ ,3734509
+ ,1
+ ,2830639
+ ,146
+ ,146
+ ,2741680
+ ,2847849
+ ,2970090
+ ,3761285
+ ,1
+ ,3257673
+ ,147
+ ,147
+ ,2830639
+ ,2741680
+ ,2847849
+ ,2970090
+ ,1
+ ,3480085
+ ,148
+ ,148
+ ,3257673
+ ,2830639
+ ,2741680
+ ,2847849
+ ,1
+ ,3843271
+ ,149
+ ,149
+ ,3480085
+ ,3257673
+ ,2830639
+ ,2741680
+ ,1
+ ,3796961
+ ,150
+ ,150
+ ,3843271
+ ,3480085
+ ,3257673
+ ,2830639
+ ,1
+ ,4337767
+ ,151
+ ,151
+ ,3796961
+ ,3843271
+ ,3480085
+ ,3257673
+ ,1
+ ,4243630
+ ,152
+ ,152
+ ,4337767
+ ,3796961
+ ,3843271
+ ,3480085
+ ,1
+ ,3927202
+ ,153
+ ,153
+ ,4243630
+ ,4337767
+ ,3796961
+ ,3843271
+ ,1
+ ,3915296
+ ,154
+ ,154
+ ,3927202
+ ,4243630
+ ,4337767
+ ,3796961
+ ,1
+ ,3087396
+ ,155
+ ,155
+ ,3915296
+ ,3927202
+ ,4243630
+ ,4337767
+ ,1
+ ,2963792
+ ,156
+ ,156
+ ,3087396
+ ,3915296
+ ,3927202
+ ,4243630
+ ,1
+ ,2955792
+ ,157
+ ,157
+ ,2963792
+ ,3087396
+ ,3915296
+ ,3927202
+ ,1
+ ,2829925
+ ,158
+ ,158
+ ,2955792
+ ,2963792
+ ,3087396
+ ,3915296
+ ,1
+ ,3281195
+ ,159
+ ,159
+ ,2829925
+ ,2955792
+ ,2963792
+ ,3087396
+ ,1
+ ,3548011
+ ,160
+ ,160
+ ,3281195
+ ,2829925
+ ,2955792
+ ,2963792
+ ,1
+ ,4059648
+ ,161
+ ,161
+ ,3548011
+ ,3281195
+ ,2829925
+ ,2955792
+ ,1
+ ,3941175
+ ,162
+ ,162
+ ,4059648
+ ,3548011
+ ,3281195
+ ,2829925
+ ,1
+ ,4528594
+ ,163
+ ,163
+ ,3941175
+ ,4059648
+ ,3548011
+ ,3281195
+ ,1
+ ,4433151
+ ,164
+ ,164
+ ,4528594
+ ,3941175
+ ,4059648
+ ,3548011
+ ,1
+ ,4145737
+ ,165
+ ,165
+ ,4433151
+ ,4528594
+ ,3941175
+ ,4059648
+ ,1
+ ,4077132
+ ,166
+ ,166
+ ,4145737
+ ,4433151
+ ,4528594
+ ,3941175
+ ,1
+ ,3198519
+ ,167
+ ,167
+ ,4077132
+ ,4145737
+ ,4433151
+ ,4528594
+ ,1
+ ,3078660
+ ,168
+ ,168
+ ,3198519
+ ,4077132
+ ,4145737
+ ,4433151
+ ,1
+ ,3028202
+ ,169
+ ,169
+ ,3078660
+ ,3198519
+ ,4077132
+ ,4145737
+ ,1
+ ,2858642
+ ,170
+ ,170
+ ,3028202
+ ,3078660
+ ,3198519
+ ,4077132
+ ,1
+ ,3398954
+ ,171
+ ,171
+ ,2858642
+ ,3028202
+ ,3078660
+ ,3198519
+ ,1
+ ,3808883
+ ,172
+ ,172
+ ,3398954
+ ,2858642
+ ,3028202
+ ,3078660
+ ,1
+ ,4175961
+ ,173
+ ,173
+ ,3808883
+ ,3398954
+ ,2858642
+ ,3028202
+ ,1
+ ,4227542
+ ,174
+ ,174
+ ,4175961
+ ,3808883
+ ,3398954
+ ,2858642
+ ,1
+ ,4744616
+ ,175
+ ,175
+ ,4227542
+ ,4175961
+ ,3808883
+ ,3398954
+ ,1
+ ,4608012
+ ,176
+ ,176
+ ,4744616
+ ,4227542
+ ,4175961
+ ,3808883
+ ,1
+ ,4295049
+ ,177
+ ,177
+ ,4608012
+ ,4744616
+ ,4227542
+ ,4175961
+ ,1
+ ,4201144
+ ,178
+ ,178
+ ,4295049
+ ,4608012
+ ,4744616
+ ,4227542
+ ,1
+ ,3353276
+ ,179
+ ,179
+ ,4201144
+ ,4295049
+ ,4608012
+ ,4744616
+ ,1
+ ,3286851
+ ,180
+ ,180
+ ,3353276
+ ,4201144
+ ,4295049
+ ,4608012
+ ,1
+ ,3169889
+ ,181
+ ,181
+ ,3286851
+ ,3353276
+ ,4201144
+ ,4295049
+ ,1
+ ,3051720
+ ,182
+ ,182
+ ,3169889
+ ,3286851
+ ,3353276
+ ,4201144
+ ,1
+ ,3695426
+ ,183
+ ,183
+ ,3051720
+ ,3169889
+ ,3286851
+ ,3353276
+ ,1
+ ,3905501
+ ,184
+ ,184
+ ,3695426
+ ,3051720
+ ,3169889
+ ,3286851
+ ,1
+ ,4296458
+ ,185
+ ,185
+ ,3905501
+ ,3695426
+ ,3051720
+ ,3169889
+ ,1
+ ,4246247
+ ,186
+ ,186
+ ,4296458
+ ,3905501
+ ,3695426
+ ,3051720
+ ,1
+ ,4921849
+ ,187
+ ,187
+ ,4246247
+ ,4296458
+ ,3905501
+ ,3695426
+ ,1
+ ,4821446
+ ,188
+ ,188
+ ,4921849
+ ,4246247
+ ,4296458
+ ,3905501
+ ,1
+ ,4425064
+ ,189
+ ,189
+ ,4821446
+ ,4921849
+ ,4246247
+ ,4296458
+ ,1
+ ,4379099
+ ,190
+ ,190
+ ,4425064
+ ,4821446
+ ,4921849
+ ,4246247
+ ,1
+ ,3472889
+ ,191
+ ,191
+ ,4379099
+ ,4425064
+ ,4821446
+ ,4921849
+ ,1
+ ,3359160
+ ,192
+ ,192
+ ,3472889
+ ,4379099
+ ,4425064
+ ,4821446
+ ,1
+ ,3200944
+ ,193
+ ,193
+ ,3359160
+ ,3472889
+ ,4379099
+ ,4425064
+ ,1
+ ,3153170
+ ,194
+ ,194
+ ,3200944
+ ,3359160
+ ,3472889
+ ,4379099
+ ,1
+ ,3741498
+ ,195
+ ,195
+ ,3153170
+ ,3200944
+ ,3359160
+ ,3472889
+ ,1
+ ,3918719
+ ,196
+ ,196
+ ,3741498
+ ,3153170
+ ,3200944
+ ,3359160
+ ,1
+ ,4403449
+ ,197
+ ,197
+ ,3918719
+ ,3741498
+ ,3153170
+ ,3200944
+ ,1
+ ,4400407
+ ,198
+ ,198
+ ,4403449
+ ,3918719
+ ,3741498
+ ,3153170
+ ,1
+ ,4847473
+ ,199
+ ,199
+ ,4400407
+ ,4403449
+ ,3918719
+ ,3741498
+ ,1
+ ,4716136
+ ,200
+ ,200
+ ,4847473
+ ,4400407
+ ,4403449
+ ,3918719
+ ,1
+ ,4297440
+ ,201
+ ,201
+ ,4716136
+ ,4847473
+ ,4400407
+ ,4403449
+ ,1
+ ,4272253
+ ,202
+ ,202
+ ,4297440
+ ,4716136
+ ,4847473
+ ,4400407
+ ,1
+ ,3271834
+ ,203
+ ,203
+ ,4272253
+ ,4297440
+ ,4716136
+ ,4847473
+ ,1
+ ,3168388
+ ,204
+ ,204
+ ,3271834
+ ,4272253
+ ,4297440
+ ,4716136
+ ,1
+ ,2911748
+ ,205
+ ,205
+ ,3168388
+ ,3271834
+ ,4272253
+ ,4297440
+ ,1
+ ,2720999
+ ,206
+ ,206
+ ,2911748
+ ,3168388
+ ,3271834
+ ,4272253
+ ,1
+ ,3199918
+ ,207
+ ,207
+ ,2720999
+ ,2911748
+ ,3168388
+ ,3271834
+ ,1
+ ,3672623
+ ,208
+ ,208
+ ,3199918
+ ,2720999
+ ,2911748
+ ,3168388
+ ,1
+ ,3892013
+ ,209
+ ,209
+ ,3672623
+ ,3199918
+ ,2720999
+ ,2911748
+ ,1
+ ,3850845
+ ,210
+ ,210
+ ,3892013
+ ,3672623
+ ,3199918
+ ,2720999
+ ,1
+ ,4532467
+ ,211
+ ,211
+ ,3850845
+ ,3892013
+ ,3672623
+ ,3199918
+ ,1
+ ,4484739
+ ,212
+ ,212
+ ,4532467
+ ,3850845
+ ,3892013
+ ,3672623
+ ,1
+ ,4014972
+ ,213
+ ,213
+ ,4484739
+ ,4532467
+ ,3850845
+ ,3892013
+ ,1
+ ,3983758
+ ,214
+ ,214
+ ,4014972
+ ,4484739
+ ,4532467
+ ,3850845
+ ,1
+ ,3158459
+ ,215
+ ,215
+ ,3983758
+ ,4014972
+ ,4484739
+ ,4532467
+ ,1
+ ,3100569
+ ,216
+ ,216
+ ,3158459
+ ,3983758
+ ,4014972
+ ,4484739
+ ,1
+ ,2935404
+ ,217
+ ,217
+ ,3100569
+ ,3158459
+ ,3983758
+ ,4014972
+ ,1
+ ,2855719
+ ,218
+ ,218
+ ,2935404
+ ,3100569
+ ,3158459
+ ,3983758
+ ,1
+ ,3465611
+ ,219
+ ,219
+ ,2855719
+ ,2935404
+ ,3100569
+ ,3158459
+ ,1
+ ,3006985
+ ,220
+ ,220
+ ,3465611
+ ,2855719
+ ,2935404
+ ,3100569
+ ,1
+ ,4095110
+ ,221
+ ,221
+ ,3006985
+ ,3465611
+ ,2855719
+ ,2935404
+ ,1
+ ,4104793
+ ,222
+ ,222
+ ,4095110
+ ,3006985
+ ,3465611
+ ,2855719
+ ,1
+ ,4730788
+ ,223
+ ,223
+ ,4104793
+ ,4095110
+ ,3006985
+ ,3465611
+ ,1
+ ,4642726
+ ,224
+ ,224
+ ,4730788
+ ,4104793
+ ,4095110
+ ,3006985
+ ,1
+ ,4246919
+ ,225
+ ,225
+ ,4642726
+ ,4730788
+ ,4104793
+ ,4095110
+ ,1
+ ,4308117
+ ,226
+ ,226
+ ,4246919
+ ,4642726
+ ,4730788
+ ,4104793)
+ ,dim=c(8
+ ,214)
+ ,dimnames=list(c('9/11'
+ ,'Passagiers'
+ ,'t'
+ ,'T_9/11'
+ ,'Passagiers-1'
+ ,'Passagiers-2'
+ ,'Passagiers-3'
+ ,'Passagiers-4')
+ ,1:214))
> y <- array(NA,dim=c(8,214),dimnames=list(c('9/11','Passagiers','t','T_9/11','Passagiers-1','Passagiers-2','Passagiers-3','Passagiers-4'),1:214))
> 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 = 'Include Monthly Dummies'
> par1 = '2'
> #'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)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> 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
Passagiers 9/11 t T_9/11 Passagiers-1 Passagiers-2 Passagiers-3
1 1281151 0 13 0 1281814 1298756 1707752
2 1164976 0 14 0 1281151 1281814 1298756
3 1454329 0 15 0 1164976 1281151 1281814
4 1645288 0 16 0 1454329 1164976 1281151
5 1817743 0 17 0 1645288 1454329 1164976
6 1895785 0 18 0 1817743 1645288 1454329
7 2236311 0 19 0 1895785 1817743 1645288
8 2295951 0 20 0 2236311 1895785 1817743
9 2087315 0 21 0 2295951 2236311 1895785
10 1980891 0 22 0 2087315 2295951 2236311
11 1465446 0 23 0 1980891 2087315 2295951
12 1445026 0 24 0 1465446 1980891 2087315
13 1488120 0 25 0 1445026 1465446 1980891
14 1338333 0 26 0 1488120 1445026 1465446
15 1715789 0 27 0 1338333 1488120 1445026
16 1806090 0 28 0 1715789 1338333 1488120
17 2083316 0 29 0 1806090 1715789 1338333
18 2092278 0 30 0 2083316 1806090 1715789
19 2430800 0 31 0 2092278 2083316 1806090
20 2424894 0 32 0 2430800 2092278 2083316
21 2299016 0 33 0 2424894 2430800 2092278
22 2130688 0 34 0 2299016 2424894 2430800
23 1652221 0 35 0 2130688 2299016 2424894
24 1608162 0 36 0 1652221 2130688 2299016
25 1647074 0 37 0 1608162 1652221 2130688
26 1479691 0 38 0 1647074 1608162 1652221
27 1884978 0 39 0 1479691 1647074 1608162
28 2007898 0 40 0 1884978 1479691 1647074
29 2208954 0 41 0 2007898 1884978 1479691
30 2217164 0 42 0 2208954 2007898 1884978
31 2534291 0 43 0 2217164 2208954 2007898
32 2560312 0 44 0 2534291 2217164 2208954
33 2429069 0 45 0 2560312 2534291 2217164
34 2315077 0 46 0 2429069 2560312 2534291
35 1799608 0 47 0 2315077 2429069 2560312
36 1772590 0 48 0 1799608 2315077 2429069
37 1744799 0 49 0 1772590 1799608 2315077
38 1659093 0 50 0 1744799 1772590 1799608
39 2099821 0 51 0 1659093 1744799 1772590
40 2135736 0 52 0 2099821 1659093 1744799
41 2427894 0 53 0 2135736 2099821 1659093
42 2468882 0 54 0 2427894 2135736 2099821
43 2703217 0 55 0 2468882 2427894 2135736
44 2766841 0 56 0 2703217 2468882 2427894
45 2655236 0 57 0 2766841 2703217 2468882
46 2550373 0 58 0 2655236 2766841 2703217
47 2052097 0 59 0 2550373 2655236 2766841
48 1998055 0 60 0 2052097 2550373 2655236
49 1920748 0 61 0 1998055 2052097 2550373
50 1876694 0 62 0 1920748 1998055 2052097
51 2380930 0 63 0 1876694 1920748 1998055
52 2467402 0 64 0 2380930 1876694 1920748
53 2770771 0 65 0 2467402 2380930 1876694
54 2781340 0 66 0 2770771 2467402 2380930
55 3143926 0 67 0 2781340 2770771 2467402
56 3172235 0 68 0 3143926 2781340 2770771
57 2952540 0 69 0 3172235 3143926 2781340
58 2920877 0 70 0 2952540 3172235 3143926
59 2384552 0 71 0 2920877 2952540 3172235
60 2248987 0 72 0 2384552 2920877 2952540
61 2208616 0 73 0 2248987 2384552 2920877
62 2178756 0 74 0 2208616 2248987 2384552
63 2632870 0 75 0 2178756 2208616 2248987
64 2706905 0 76 0 2632870 2178756 2208616
65 3029745 0 77 0 2706905 2632870 2178756
66 3015402 0 78 0 3029745 2706905 2632870
67 3391414 0 79 0 3015402 3029745 2706905
68 3507805 0 80 0 3391414 3015402 3029745
69 3177852 0 81 0 3507805 3391414 3015402
70 3142961 0 82 0 3177852 3507805 3391414
71 2545815 0 83 0 3142961 3177852 3507805
72 2414007 0 84 0 2545815 3142961 3177852
73 2372578 0 85 0 2414007 2545815 3142961
74 2332664 0 86 0 2372578 2414007 2545815
75 2825328 0 87 0 2332664 2372578 2414007
76 2901478 0 88 0 2825328 2332664 2372578
77 3263955 0 89 0 2901478 2825328 2332664
78 3226738 0 90 0 3263955 2901478 2825328
79 3610786 0 91 0 3226738 3263955 2901478
80 3709274 0 92 0 3610786 3226738 3263955
81 3467185 0 93 0 3709274 3610786 3226738
82 3449646 0 94 0 3467185 3709274 3610786
83 2802951 0 95 0 3449646 3467185 3709274
84 2462530 0 96 0 2802951 3449646 3467185
85 2490645 0 97 0 2462530 2802951 3449646
86 2561520 0 98 0 2490645 2462530 2802951
87 3067554 0 99 0 2561520 2490645 2462530
88 3226951 0 100 0 3067554 2561520 2490645
89 3546493 0 101 0 3226951 3067554 2561520
90 3492787 0 102 0 3546493 3226951 3067554
91 3952263 0 103 0 3492787 3546493 3226951
92 3932072 0 104 0 3952263 3492787 3546493
93 3720284 0 105 0 3932072 3952263 3492787
94 3651555 0 106 0 3720284 3932072 3952263
95 2914972 0 107 0 3651555 3720284 3932072
96 2713514 0 108 0 2914972 3651555 3720284
97 2703997 0 109 0 2713514 2914972 3651555
98 2591373 0 110 0 2703997 2713514 2914972
99 3163748 0 111 0 2591373 2703997 2713514
100 3355137 0 112 0 3163748 2591373 2703997
101 3613702 0 113 0 3355137 3163748 2591373
102 3686773 0 114 0 3613702 3355137 3163748
103 4098716 0 115 0 3686773 3613702 3355137
104 4063517 0 116 0 4098716 3686773 3613702
105 3551489 1 117 117 4063517 4098716 3686773
106 3226663 1 118 118 3551489 4063517 4098716
107 2656842 1 119 119 3226663 3551489 4063517
108 2597484 1 120 120 2656842 3226663 3551489
109 2572399 1 121 121 2597484 2656842 3226663
110 2596631 1 122 122 2572399 2597484 2656842
111 3165225 1 123 123 2596631 2572399 2597484
112 3303145 1 124 124 3165225 2596631 2572399
113 3698247 1 125 125 3303145 3165225 2596631
114 3668631 1 126 126 3698247 3303145 3165225
115 4130433 1 127 127 3668631 3698247 3303145
116 4131400 1 128 128 4130433 3668631 3698247
117 3864358 1 129 129 4131400 4130433 3668631
118 3721110 1 130 130 3864358 4131400 4130433
119 2892532 1 131 131 3721110 3864358 4131400
120 2843451 1 132 132 2892532 3721110 3864358
121 2747502 1 133 133 2843451 2892532 3721110
122 2668775 1 134 134 2747502 2843451 2892532
123 3018602 1 135 135 2668775 2747502 2843451
124 3013392 1 136 136 3018602 2668775 2747502
125 3393657 1 137 137 3013392 3018602 2668775
126 3544233 1 138 138 3393657 3013392 3018602
127 4075832 1 139 139 3544233 3393657 3013392
128 4032923 1 140 140 4075832 3544233 3393657
129 3734509 1 141 141 4032923 4075832 3544233
130 3761285 1 142 142 3734509 4032923 4075832
131 2970090 1 143 143 3761285 3734509 4032923
132 2847849 1 144 144 2970090 3761285 3734509
133 2741680 1 145 145 2847849 2970090 3761285
134 2830639 1 146 146 2741680 2847849 2970090
135 3257673 1 147 147 2830639 2741680 2847849
136 3480085 1 148 148 3257673 2830639 2741680
137 3843271 1 149 149 3480085 3257673 2830639
138 3796961 1 150 150 3843271 3480085 3257673
139 4337767 1 151 151 3796961 3843271 3480085
140 4243630 1 152 152 4337767 3796961 3843271
141 3927202 1 153 153 4243630 4337767 3796961
142 3915296 1 154 154 3927202 4243630 4337767
143 3087396 1 155 155 3915296 3927202 4243630
144 2963792 1 156 156 3087396 3915296 3927202
145 2955792 1 157 157 2963792 3087396 3915296
146 2829925 1 158 158 2955792 2963792 3087396
147 3281195 1 159 159 2829925 2955792 2963792
148 3548011 1 160 160 3281195 2829925 2955792
149 4059648 1 161 161 3548011 3281195 2829925
150 3941175 1 162 162 4059648 3548011 3281195
151 4528594 1 163 163 3941175 4059648 3548011
152 4433151 1 164 164 4528594 3941175 4059648
153 4145737 1 165 165 4433151 4528594 3941175
154 4077132 1 166 166 4145737 4433151 4528594
155 3198519 1 167 167 4077132 4145737 4433151
156 3078660 1 168 168 3198519 4077132 4145737
157 3028202 1 169 169 3078660 3198519 4077132
158 2858642 1 170 170 3028202 3078660 3198519
159 3398954 1 171 171 2858642 3028202 3078660
160 3808883 1 172 172 3398954 2858642 3028202
161 4175961 1 173 173 3808883 3398954 2858642
162 4227542 1 174 174 4175961 3808883 3398954
163 4744616 1 175 175 4227542 4175961 3808883
164 4608012 1 176 176 4744616 4227542 4175961
165 4295049 1 177 177 4608012 4744616 4227542
166 4201144 1 178 178 4295049 4608012 4744616
167 3353276 1 179 179 4201144 4295049 4608012
168 3286851 1 180 180 3353276 4201144 4295049
169 3169889 1 181 181 3286851 3353276 4201144
170 3051720 1 182 182 3169889 3286851 3353276
171 3695426 1 183 183 3051720 3169889 3286851
172 3905501 1 184 184 3695426 3051720 3169889
173 4296458 1 185 185 3905501 3695426 3051720
174 4246247 1 186 186 4296458 3905501 3695426
175 4921849 1 187 187 4246247 4296458 3905501
176 4821446 1 188 188 4921849 4246247 4296458
177 4425064 1 189 189 4821446 4921849 4246247
178 4379099 1 190 190 4425064 4821446 4921849
179 3472889 1 191 191 4379099 4425064 4821446
180 3359160 1 192 192 3472889 4379099 4425064
181 3200944 1 193 193 3359160 3472889 4379099
182 3153170 1 194 194 3200944 3359160 3472889
183 3741498 1 195 195 3153170 3200944 3359160
184 3918719 1 196 196 3741498 3153170 3200944
185 4403449 1 197 197 3918719 3741498 3153170
186 4400407 1 198 198 4403449 3918719 3741498
187 4847473 1 199 199 4400407 4403449 3918719
188 4716136 1 200 200 4847473 4400407 4403449
189 4297440 1 201 201 4716136 4847473 4400407
190 4272253 1 202 202 4297440 4716136 4847473
191 3271834 1 203 203 4272253 4297440 4716136
192 3168388 1 204 204 3271834 4272253 4297440
193 2911748 1 205 205 3168388 3271834 4272253
194 2720999 1 206 206 2911748 3168388 3271834
195 3199918 1 207 207 2720999 2911748 3168388
196 3672623 1 208 208 3199918 2720999 2911748
197 3892013 1 209 209 3672623 3199918 2720999
198 3850845 1 210 210 3892013 3672623 3199918
199 4532467 1 211 211 3850845 3892013 3672623
200 4484739 1 212 212 4532467 3850845 3892013
201 4014972 1 213 213 4484739 4532467 3850845
202 3983758 1 214 214 4014972 4484739 4532467
203 3158459 1 215 215 3983758 4014972 4484739
204 3100569 1 216 216 3158459 3983758 4014972
205 2935404 1 217 217 3100569 3158459 3983758
206 2855719 1 218 218 2935404 3100569 3158459
207 3465611 1 219 219 2855719 2935404 3100569
208 3006985 1 220 220 3465611 2855719 2935404
209 4095110 1 221 221 3006985 3465611 2855719
210 4104793 1 222 222 4095110 3006985 3465611
211 4730788 1 223 223 4104793 4095110 3006985
212 4642726 1 224 224 4730788 4104793 4095110
213 4246919 1 225 225 4642726 4730788 4104793
214 4308117 1 226 226 4246919 4642726 4730788
Passagiers-4 M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 M11
1 1867943 1 0 0 0 0 0 0 0 0 0 0
2 1707752 0 1 0 0 0 0 0 0 0 0 0
3 1298756 0 0 1 0 0 0 0 0 0 0 0
4 1281814 0 0 0 1 0 0 0 0 0 0 0
5 1281151 0 0 0 0 1 0 0 0 0 0 0
6 1164976 0 0 0 0 0 1 0 0 0 0 0
7 1454329 0 0 0 0 0 0 1 0 0 0 0
8 1645288 0 0 0 0 0 0 0 1 0 0 0
9 1817743 0 0 0 0 0 0 0 0 1 0 0
10 1895785 0 0 0 0 0 0 0 0 0 1 0
11 2236311 0 0 0 0 0 0 0 0 0 0 1
12 2295951 0 0 0 0 0 0 0 0 0 0 0
13 2087315 1 0 0 0 0 0 0 0 0 0 0
14 1980891 0 1 0 0 0 0 0 0 0 0 0
15 1465446 0 0 1 0 0 0 0 0 0 0 0
16 1445026 0 0 0 1 0 0 0 0 0 0 0
17 1488120 0 0 0 0 1 0 0 0 0 0 0
18 1338333 0 0 0 0 0 1 0 0 0 0 0
19 1715789 0 0 0 0 0 0 1 0 0 0 0
20 1806090 0 0 0 0 0 0 0 1 0 0 0
21 2083316 0 0 0 0 0 0 0 0 1 0 0
22 2092278 0 0 0 0 0 0 0 0 0 1 0
23 2430800 0 0 0 0 0 0 0 0 0 0 1
24 2424894 0 0 0 0 0 0 0 0 0 0 0
25 2299016 1 0 0 0 0 0 0 0 0 0 0
26 2130688 0 1 0 0 0 0 0 0 0 0 0
27 1652221 0 0 1 0 0 0 0 0 0 0 0
28 1608162 0 0 0 1 0 0 0 0 0 0 0
29 1647074 0 0 0 0 1 0 0 0 0 0 0
30 1479691 0 0 0 0 0 1 0 0 0 0 0
31 1884978 0 0 0 0 0 0 1 0 0 0 0
32 2007898 0 0 0 0 0 0 0 1 0 0 0
33 2208954 0 0 0 0 0 0 0 0 1 0 0
34 2217164 0 0 0 0 0 0 0 0 0 1 0
35 2534291 0 0 0 0 0 0 0 0 0 0 1
36 2560312 0 0 0 0 0 0 0 0 0 0 0
37 2429069 1 0 0 0 0 0 0 0 0 0 0
38 2315077 0 1 0 0 0 0 0 0 0 0 0
39 1799608 0 0 1 0 0 0 0 0 0 0 0
40 1772590 0 0 0 1 0 0 0 0 0 0 0
41 1744799 0 0 0 0 1 0 0 0 0 0 0
42 1659093 0 0 0 0 0 1 0 0 0 0 0
43 2099821 0 0 0 0 0 0 1 0 0 0 0
44 2135736 0 0 0 0 0 0 0 1 0 0 0
45 2427894 0 0 0 0 0 0 0 0 1 0 0
46 2468882 0 0 0 0 0 0 0 0 0 1 0
47 2703217 0 0 0 0 0 0 0 0 0 0 1
48 2766841 0 0 0 0 0 0 0 0 0 0 0
49 2655236 1 0 0 0 0 0 0 0 0 0 0
50 2550373 0 1 0 0 0 0 0 0 0 0 0
51 2052097 0 0 1 0 0 0 0 0 0 0 0
52 1998055 0 0 0 1 0 0 0 0 0 0 0
53 1920748 0 0 0 0 1 0 0 0 0 0 0
54 1876694 0 0 0 0 0 1 0 0 0 0 0
55 2380930 0 0 0 0 0 0 1 0 0 0 0
56 2467402 0 0 0 0 0 0 0 1 0 0 0
57 2770771 0 0 0 0 0 0 0 0 1 0 0
58 2781340 0 0 0 0 0 0 0 0 0 1 0
59 3143926 0 0 0 0 0 0 0 0 0 0 1
60 3172235 0 0 0 0 0 0 0 0 0 0 0
61 2952540 1 0 0 0 0 0 0 0 0 0 0
62 2920877 0 1 0 0 0 0 0 0 0 0 0
63 2384552 0 0 1 0 0 0 0 0 0 0 0
64 2248987 0 0 0 1 0 0 0 0 0 0 0
65 2208616 0 0 0 0 1 0 0 0 0 0 0
66 2178756 0 0 0 0 0 1 0 0 0 0 0
67 2632870 0 0 0 0 0 0 1 0 0 0 0
68 2706905 0 0 0 0 0 0 0 1 0 0 0
69 3029745 0 0 0 0 0 0 0 0 1 0 0
70 3015402 0 0 0 0 0 0 0 0 0 1 0
71 3391414 0 0 0 0 0 0 0 0 0 0 1
72 3507805 0 0 0 0 0 0 0 0 0 0 0
73 3177852 1 0 0 0 0 0 0 0 0 0 0
74 3142961 0 1 0 0 0 0 0 0 0 0 0
75 2545815 0 0 1 0 0 0 0 0 0 0 0
76 2414007 0 0 0 1 0 0 0 0 0 0 0
77 2372578 0 0 0 0 1 0 0 0 0 0 0
78 2332664 0 0 0 0 0 1 0 0 0 0 0
79 2825328 0 0 0 0 0 0 1 0 0 0 0
80 2901478 0 0 0 0 0 0 0 1 0 0 0
81 3263955 0 0 0 0 0 0 0 0 1 0 0
82 3226738 0 0 0 0 0 0 0 0 0 1 0
83 3610786 0 0 0 0 0 0 0 0 0 0 1
84 3709274 0 0 0 0 0 0 0 0 0 0 0
85 3467185 1 0 0 0 0 0 0 0 0 0 0
86 3449646 0 1 0 0 0 0 0 0 0 0 0
87 2802951 0 0 1 0 0 0 0 0 0 0 0
88 2462530 0 0 0 1 0 0 0 0 0 0 0
89 2490645 0 0 0 0 1 0 0 0 0 0 0
90 2561520 0 0 0 0 0 1 0 0 0 0 0
91 3067554 0 0 0 0 0 0 1 0 0 0 0
92 3226951 0 0 0 0 0 0 0 1 0 0 0
93 3546493 0 0 0 0 0 0 0 0 1 0 0
94 3492787 0 0 0 0 0 0 0 0 0 1 0
95 3952263 0 0 0 0 0 0 0 0 0 0 1
96 3932072 0 0 0 0 0 0 0 0 0 0 0
97 3720284 1 0 0 0 0 0 0 0 0 0 0
98 3651555 0 1 0 0 0 0 0 0 0 0 0
99 2914972 0 0 1 0 0 0 0 0 0 0 0
100 2713514 0 0 0 1 0 0 0 0 0 0 0
101 2703997 0 0 0 0 1 0 0 0 0 0 0
102 2591373 0 0 0 0 0 1 0 0 0 0 0
103 3163748 0 0 0 0 0 0 1 0 0 0 0
104 3355137 0 0 0 0 0 0 0 1 0 0 0
105 3613702 0 0 0 0 0 0 0 0 1 0 0
106 3686773 0 0 0 0 0 0 0 0 0 1 0
107 4098716 0 0 0 0 0 0 0 0 0 0 1
108 4063517 0 0 0 0 0 0 0 0 0 0 0
109 3551489 1 0 0 0 0 0 0 0 0 0 0
110 3226663 0 1 0 0 0 0 0 0 0 0 0
111 2656842 0 0 1 0 0 0 0 0 0 0 0
112 2597484 0 0 0 1 0 0 0 0 0 0 0
113 2572399 0 0 0 0 1 0 0 0 0 0 0
114 2596631 0 0 0 0 0 1 0 0 0 0 0
115 3165225 0 0 0 0 0 0 1 0 0 0 0
116 3303145 0 0 0 0 0 0 0 1 0 0 0
117 3698247 0 0 0 0 0 0 0 0 1 0 0
118 3668631 0 0 0 0 0 0 0 0 0 1 0
119 4130433 0 0 0 0 0 0 0 0 0 0 1
120 4131400 0 0 0 0 0 0 0 0 0 0 0
121 3864358 1 0 0 0 0 0 0 0 0 0 0
122 3721110 0 1 0 0 0 0 0 0 0 0 0
123 2892532 0 0 1 0 0 0 0 0 0 0 0
124 2843451 0 0 0 1 0 0 0 0 0 0 0
125 2747502 0 0 0 0 1 0 0 0 0 0 0
126 2668775 0 0 0 0 0 1 0 0 0 0 0
127 3018602 0 0 0 0 0 0 1 0 0 0 0
128 3013392 0 0 0 0 0 0 0 1 0 0 0
129 3393657 0 0 0 0 0 0 0 0 1 0 0
130 3544233 0 0 0 0 0 0 0 0 0 1 0
131 4075832 0 0 0 0 0 0 0 0 0 0 1
132 4032923 0 0 0 0 0 0 0 0 0 0 0
133 3734509 1 0 0 0 0 0 0 0 0 0 0
134 3761285 0 1 0 0 0 0 0 0 0 0 0
135 2970090 0 0 1 0 0 0 0 0 0 0 0
136 2847849 0 0 0 1 0 0 0 0 0 0 0
137 2741680 0 0 0 0 1 0 0 0 0 0 0
138 2830639 0 0 0 0 0 1 0 0 0 0 0
139 3257673 0 0 0 0 0 0 1 0 0 0 0
140 3480085 0 0 0 0 0 0 0 1 0 0 0
141 3843271 0 0 0 0 0 0 0 0 1 0 0
142 3796961 0 0 0 0 0 0 0 0 0 1 0
143 4337767 0 0 0 0 0 0 0 0 0 0 1
144 4243630 0 0 0 0 0 0 0 0 0 0 0
145 3927202 1 0 0 0 0 0 0 0 0 0 0
146 3915296 0 1 0 0 0 0 0 0 0 0 0
147 3087396 0 0 1 0 0 0 0 0 0 0 0
148 2963792 0 0 0 1 0 0 0 0 0 0 0
149 2955792 0 0 0 0 1 0 0 0 0 0 0
150 2829925 0 0 0 0 0 1 0 0 0 0 0
151 3281195 0 0 0 0 0 0 1 0 0 0 0
152 3548011 0 0 0 0 0 0 0 1 0 0 0
153 4059648 0 0 0 0 0 0 0 0 1 0 0
154 3941175 0 0 0 0 0 0 0 0 0 1 0
155 4528594 0 0 0 0 0 0 0 0 0 0 1
156 4433151 0 0 0 0 0 0 0 0 0 0 0
157 4145737 1 0 0 0 0 0 0 0 0 0 0
158 4077132 0 1 0 0 0 0 0 0 0 0 0
159 3198519 0 0 1 0 0 0 0 0 0 0 0
160 3078660 0 0 0 1 0 0 0 0 0 0 0
161 3028202 0 0 0 0 1 0 0 0 0 0 0
162 2858642 0 0 0 0 0 1 0 0 0 0 0
163 3398954 0 0 0 0 0 0 1 0 0 0 0
164 3808883 0 0 0 0 0 0 0 1 0 0 0
165 4175961 0 0 0 0 0 0 0 0 1 0 0
166 4227542 0 0 0 0 0 0 0 0 0 1 0
167 4744616 0 0 0 0 0 0 0 0 0 0 1
168 4608012 0 0 0 0 0 0 0 0 0 0 0
169 4295049 1 0 0 0 0 0 0 0 0 0 0
170 4201144 0 1 0 0 0 0 0 0 0 0 0
171 3353276 0 0 1 0 0 0 0 0 0 0 0
172 3286851 0 0 0 1 0 0 0 0 0 0 0
173 3169889 0 0 0 0 1 0 0 0 0 0 0
174 3051720 0 0 0 0 0 1 0 0 0 0 0
175 3695426 0 0 0 0 0 0 1 0 0 0 0
176 3905501 0 0 0 0 0 0 0 1 0 0 0
177 4296458 0 0 0 0 0 0 0 0 1 0 0
178 4246247 0 0 0 0 0 0 0 0 0 1 0
179 4921849 0 0 0 0 0 0 0 0 0 0 1
180 4821446 0 0 0 0 0 0 0 0 0 0 0
181 4425064 1 0 0 0 0 0 0 0 0 0 0
182 4379099 0 1 0 0 0 0 0 0 0 0 0
183 3472889 0 0 1 0 0 0 0 0 0 0 0
184 3359160 0 0 0 1 0 0 0 0 0 0 0
185 3200944 0 0 0 0 1 0 0 0 0 0 0
186 3153170 0 0 0 0 0 1 0 0 0 0 0
187 3741498 0 0 0 0 0 0 1 0 0 0 0
188 3918719 0 0 0 0 0 0 0 1 0 0 0
189 4403449 0 0 0 0 0 0 0 0 1 0 0
190 4400407 0 0 0 0 0 0 0 0 0 1 0
191 4847473 0 0 0 0 0 0 0 0 0 0 1
192 4716136 0 0 0 0 0 0 0 0 0 0 0
193 4297440 1 0 0 0 0 0 0 0 0 0 0
194 4272253 0 1 0 0 0 0 0 0 0 0 0
195 3271834 0 0 1 0 0 0 0 0 0 0 0
196 3168388 0 0 0 1 0 0 0 0 0 0 0
197 2911748 0 0 0 0 1 0 0 0 0 0 0
198 2720999 0 0 0 0 0 1 0 0 0 0 0
199 3199918 0 0 0 0 0 0 1 0 0 0 0
200 3672623 0 0 0 0 0 0 0 1 0 0 0
201 3892013 0 0 0 0 0 0 0 0 1 0 0
202 3850845 0 0 0 0 0 0 0 0 0 1 0
203 4532467 0 0 0 0 0 0 0 0 0 0 1
204 4484739 0 0 0 0 0 0 0 0 0 0 0
205 4014972 1 0 0 0 0 0 0 0 0 0 0
206 3983758 0 1 0 0 0 0 0 0 0 0 0
207 3158459 0 0 1 0 0 0 0 0 0 0 0
208 3100569 0 0 0 1 0 0 0 0 0 0 0
209 2935404 0 0 0 0 1 0 0 0 0 0 0
210 2855719 0 0 0 0 0 1 0 0 0 0 0
211 3465611 0 0 0 0 0 0 1 0 0 0 0
212 3006985 0 0 0 0 0 0 0 1 0 0 0
213 4095110 0 0 0 0 0 0 0 0 1 0 0
214 4104793 0 0 0 0 0 0 0 0 0 1 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) `9/11` t `T_9/11` `Passagiers-1`
1.940e+05 3.801e+05 5.405e+03 -3.722e+03 6.587e-01
`Passagiers-2` `Passagiers-3` `Passagiers-4` M1 M2
2.647e-01 -2.761e-03 -2.097e-01 1.343e+05 1.003e+05
M3 M4 M5 M6 M7
5.012e+05 3.151e+05 4.537e+05 1.558e+05 6.148e+05
M8 M9 M10 M11
3.144e+05 -7.045e+03 1.172e+05 -3.740e+05
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-632906 -50188 -2490 62624 420778
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.940e+05 9.745e+04 1.991 0.047893 *
`9/11` 3.801e+05 9.458e+04 4.019 8.36e-05 ***
t 5.405e+03 1.010e+03 5.352 2.43e-07 ***
`T_9/11` -3.722e+03 8.183e+02 -4.548 9.50e-06 ***
`Passagiers-1` 6.587e-01 6.991e-02 9.422 < 2e-16 ***
`Passagiers-2` 2.647e-01 8.454e-02 3.131 0.002013 **
`Passagiers-3` -2.761e-03 8.369e-02 -0.033 0.973716
`Passagiers-4` -2.097e-01 6.920e-02 -3.031 0.002768 **
M1 1.343e+05 6.410e+04 2.095 0.037451 *
M2 1.003e+05 6.629e+04 1.513 0.131935
M3 5.012e+05 6.780e+04 7.392 4.15e-12 ***
M4 3.151e+05 9.234e+04 3.413 0.000782 ***
M5 4.537e+05 8.790e+04 5.162 6.00e-07 ***
M6 1.558e+05 9.935e+04 1.569 0.118379
M7 6.148e+05 7.784e+04 7.898 2.02e-13 ***
M8 3.144e+05 9.604e+04 3.274 0.001257 **
M9 -7.045e+03 7.884e+04 -0.089 0.928892
M10 1.172e+05 6.842e+04 1.713 0.088297 .
M11 -3.740e+05 6.017e+04 -6.215 3.04e-09 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 111900 on 195 degrees of freedom
Multiple R-squared: 0.9848, Adjusted R-squared: 0.9833
F-statistic: 699.7 on 18 and 195 DF, p-value: < 2.2e-16
> if (n > n25) {
+ kp3 <- k + 3
+ nmkm3 <- n - k - 3
+ gqarr <- array(NA, dim=c(nmkm3-kp3+1,3))
+ numgqtests <- 0
+ numsignificant1 <- 0
+ numsignificant5 <- 0
+ numsignificant10 <- 0
+ for (mypoint in kp3:nmkm3) {
+ j <- 0
+ numgqtests <- numgqtests + 1
+ for (myalt in c('greater', 'two.sided', 'less')) {
+ j <- j + 1
+ gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value
+ }
+ if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1
+ if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1
+ }
+ gqarr
+ }
[,1] [,2] [,3]
[1,] 5.503968e-02 1.100794e-01 0.9449603
[2,] 2.268071e-02 4.536142e-02 0.9773193
[3,] 6.910215e-03 1.382043e-02 0.9930898
[4,] 2.244675e-03 4.489349e-03 0.9977553
[5,] 8.240686e-04 1.648137e-03 0.9991759
[6,] 2.592132e-04 5.184264e-04 0.9997408
[7,] 8.320187e-05 1.664037e-04 0.9999168
[8,] 3.276403e-05 6.552806e-05 0.9999672
[9,] 2.100580e-05 4.201159e-05 0.9999790
[10,] 1.997461e-05 3.994922e-05 0.9999800
[11,] 7.475639e-06 1.495128e-05 0.9999925
[12,] 2.297465e-06 4.594930e-06 0.9999977
[13,] 7.538132e-07 1.507626e-06 0.9999992
[14,] 2.227786e-07 4.455572e-07 0.9999998
[15,] 5.942688e-08 1.188538e-07 0.9999999
[16,] 4.056417e-08 8.112834e-08 1.0000000
[17,] 1.061834e-08 2.123668e-08 1.0000000
[18,] 7.000681e-08 1.400136e-07 0.9999999
[19,] 3.927747e-08 7.855494e-08 1.0000000
[20,] 1.721160e-08 3.442320e-08 1.0000000
[21,] 1.618468e-08 3.236937e-08 1.0000000
[22,] 1.635081e-08 3.270162e-08 1.0000000
[23,] 5.746189e-09 1.149238e-08 1.0000000
[24,] 3.444419e-09 6.888838e-09 1.0000000
[25,] 2.783449e-09 5.566898e-09 1.0000000
[26,] 1.854487e-09 3.708973e-09 1.0000000
[27,] 6.421296e-10 1.284259e-09 1.0000000
[28,] 7.131577e-10 1.426315e-09 1.0000000
[29,] 3.194761e-10 6.389521e-10 1.0000000
[30,] 9.305914e-09 1.861183e-08 1.0000000
[31,] 7.979599e-09 1.595920e-08 1.0000000
[32,] 3.219159e-08 6.438319e-08 1.0000000
[33,] 1.535458e-08 3.070916e-08 1.0000000
[34,] 1.264126e-08 2.528253e-08 1.0000000
[35,] 5.135807e-09 1.027161e-08 1.0000000
[36,] 4.393649e-09 8.787298e-09 1.0000000
[37,] 3.208311e-09 6.416621e-09 1.0000000
[38,] 1.809115e-09 3.618231e-09 1.0000000
[39,] 2.702841e-09 5.405682e-09 1.0000000
[40,] 2.452507e-09 4.905014e-09 1.0000000
[41,] 1.528883e-09 3.057766e-09 1.0000000
[42,] 1.089229e-09 2.178459e-09 1.0000000
[43,] 4.700153e-10 9.400306e-10 1.0000000
[44,] 5.599262e-10 1.119852e-09 1.0000000
[45,] 2.424227e-10 4.848454e-10 1.0000000
[46,] 2.206127e-10 4.412254e-10 1.0000000
[47,] 3.785197e-10 7.570395e-10 1.0000000
[48,] 1.062501e-09 2.125001e-09 1.0000000
[49,] 4.803793e-10 9.607586e-10 1.0000000
[50,] 5.200693e-10 1.040139e-09 1.0000000
[51,] 6.021107e-10 1.204221e-09 1.0000000
[52,] 5.524302e-10 1.104860e-09 1.0000000
[53,] 2.812170e-10 5.624341e-10 1.0000000
[54,] 2.667208e-10 5.334417e-10 1.0000000
[55,] 1.264864e-10 2.529728e-10 1.0000000
[56,] 1.606567e-10 3.213133e-10 1.0000000
[57,] 8.500874e-11 1.700175e-10 1.0000000
[58,] 1.017013e-10 2.034027e-10 1.0000000
[59,] 8.019629e-11 1.603926e-10 1.0000000
[60,] 3.762807e-11 7.525614e-11 1.0000000
[61,] 2.849166e-11 5.698331e-11 1.0000000
[62,] 3.560369e-11 7.120738e-11 1.0000000
[63,] 2.926880e-08 5.853760e-08 1.0000000
[64,] 1.715960e-08 3.431920e-08 1.0000000
[65,] 3.547913e-08 7.095826e-08 1.0000000
[66,] 2.926445e-08 5.852890e-08 1.0000000
[67,] 1.972926e-08 3.945853e-08 1.0000000
[68,] 1.480418e-08 2.960836e-08 1.0000000
[69,] 9.973404e-09 1.994681e-08 1.0000000
[70,] 3.350715e-08 6.701430e-08 1.0000000
[71,] 2.081733e-08 4.163467e-08 1.0000000
[72,] 1.154438e-08 2.308876e-08 1.0000000
[73,] 6.461862e-09 1.292372e-08 1.0000000
[74,] 1.644439e-08 3.288878e-08 1.0000000
[75,] 2.264896e-08 4.529793e-08 1.0000000
[76,] 1.381039e-08 2.762079e-08 1.0000000
[77,] 1.007923e-08 2.015847e-08 1.0000000
[78,] 7.304676e-09 1.460935e-08 1.0000000
[79,] 5.413549e-09 1.082710e-08 1.0000000
[80,] 3.725478e-09 7.450956e-09 1.0000000
[81,] 2.432209e-09 4.864419e-09 1.0000000
[82,] 1.446065e-09 2.892131e-09 1.0000000
[83,] 9.020298e-10 1.804060e-09 1.0000000
[84,] 7.721606e-10 1.544321e-09 1.0000000
[85,] 4.881690e-09 9.763379e-09 1.0000000
[86,] 4.504987e-09 9.009974e-09 1.0000000
[87,] 2.498152e-09 4.996305e-09 1.0000000
[88,] 3.898742e-09 7.797485e-09 1.0000000
[89,] 2.861289e-09 5.722578e-09 1.0000000
[90,] 1.872157e-09 3.744314e-09 1.0000000
[91,] 1.292028e-09 2.584057e-09 1.0000000
[92,] 7.587223e-10 1.517445e-09 1.0000000
[93,] 8.640420e-10 1.728084e-09 1.0000000
[94,] 8.080094e-10 1.616019e-09 1.0000000
[95,] 6.342273e-10 1.268455e-09 1.0000000
[96,] 5.907779e-10 1.181556e-09 1.0000000
[97,] 5.575577e-10 1.115115e-09 1.0000000
[98,] 1.218475e-08 2.436950e-08 1.0000000
[99,] 6.830575e-09 1.366115e-08 1.0000000
[100,] 8.468299e-09 1.693660e-08 1.0000000
[101,] 5.355428e-09 1.071086e-08 1.0000000
[102,] 3.806004e-08 7.612008e-08 1.0000000
[103,] 4.703977e-07 9.407955e-07 0.9999995
[104,] 1.688519e-06 3.377038e-06 0.9999983
[105,] 1.214998e-06 2.429995e-06 0.9999988
[106,] 1.073698e-06 2.147396e-06 0.9999989
[107,] 1.177217e-06 2.354433e-06 0.9999988
[108,] 9.138129e-07 1.827626e-06 0.9999991
[109,] 9.617738e-07 1.923548e-06 0.9999990
[110,] 1.452408e-06 2.904815e-06 0.9999985
[111,] 8.506275e-07 1.701255e-06 0.9999991
[112,] 6.795914e-07 1.359183e-06 0.9999993
[113,] 1.806171e-06 3.612342e-06 0.9999982
[114,] 1.115100e-06 2.230199e-06 0.9999989
[115,] 8.652097e-07 1.730419e-06 0.9999991
[116,] 1.031221e-06 2.062443e-06 0.9999990
[117,] 9.103058e-07 1.820612e-06 0.9999991
[118,] 2.387738e-06 4.775477e-06 0.9999976
[119,] 1.895158e-06 3.790316e-06 0.9999981
[120,] 1.147474e-06 2.294948e-06 0.9999989
[121,] 9.688313e-07 1.937663e-06 0.9999990
[122,] 1.002739e-06 2.005477e-06 0.9999990
[123,] 5.902345e-07 1.180469e-06 0.9999994
[124,] 4.166478e-07 8.332955e-07 0.9999996
[125,] 2.653059e-07 5.306119e-07 0.9999997
[126,] 2.167876e-07 4.335751e-07 0.9999998
[127,] 1.625023e-07 3.250047e-07 0.9999998
[128,] 2.830914e-07 5.661829e-07 0.9999997
[129,] 3.353270e-07 6.706540e-07 0.9999997
[130,] 5.903468e-07 1.180694e-06 0.9999994
[131,] 4.385780e-07 8.771560e-07 0.9999996
[132,] 2.534046e-07 5.068092e-07 0.9999997
[133,] 1.602647e-07 3.205293e-07 0.9999998
[134,] 1.576951e-07 3.153902e-07 0.9999998
[135,] 1.028323e-07 2.056646e-07 0.9999999
[136,] 5.671141e-08 1.134228e-07 0.9999999
[137,] 4.142025e-08 8.284050e-08 1.0000000
[138,] 2.926407e-08 5.852813e-08 1.0000000
[139,] 1.828126e-07 3.656252e-07 0.9999998
[140,] 1.096685e-07 2.193370e-07 0.9999999
[141,] 5.846766e-08 1.169353e-07 0.9999999
[142,] 5.294914e-08 1.058983e-07 0.9999999
[143,] 4.315926e-08 8.631852e-08 1.0000000
[144,] 2.153324e-08 4.306647e-08 1.0000000
[145,] 1.673689e-08 3.347378e-08 1.0000000
[146,] 9.525291e-09 1.905058e-08 1.0000000
[147,] 4.921374e-09 9.842748e-09 1.0000000
[148,] 2.928608e-09 5.857215e-09 1.0000000
[149,] 1.433057e-09 2.866113e-09 1.0000000
[150,] 1.109205e-09 2.218410e-09 1.0000000
[151,] 2.290698e-09 4.581396e-09 1.0000000
[152,] 1.124213e-09 2.248427e-09 1.0000000
[153,] 6.853806e-10 1.370761e-09 1.0000000
[154,] 1.809688e-09 3.619375e-09 1.0000000
[155,] 7.887702e-10 1.577540e-09 1.0000000
[156,] 4.459495e-10 8.918991e-10 1.0000000
[157,] 1.927033e-10 3.854065e-10 1.0000000
[158,] 1.137755e-10 2.275510e-10 1.0000000
[159,] 4.399068e-11 8.798135e-11 1.0000000
[160,] 3.241201e-11 6.482401e-11 1.0000000
[161,] 3.058442e-11 6.116883e-11 1.0000000
[162,] 4.442718e-11 8.885435e-11 1.0000000
[163,] 2.284143e-09 4.568286e-09 1.0000000
[164,] 4.847850e-09 9.695700e-09 1.0000000
[165,] 2.040623e-08 4.081245e-08 1.0000000
[166,] 8.609384e-09 1.721877e-08 1.0000000
[167,] 3.584265e-09 7.168530e-09 1.0000000
[168,] 2.127695e-09 4.255390e-09 1.0000000
[169,] 1.149321e-09 2.298643e-09 1.0000000
[170,] 2.347826e-09 4.695651e-09 1.0000000
[171,] 1.289481e-09 2.578961e-09 1.0000000
> postscript(file="/var/www/rcomp/tmp/191aa1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> 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/rcomp/tmp/2karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> 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/rcomp/tmp/3karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/4karv1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/5vj8x1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> qqnorm(mysum$resid, main='Residual Normal Q-Q Plot')
> qqline(mysum$resid)
> grid()
> dev.off()
null device
1
> (myerror <- as.ts(mysum$resid))
Time Series:
Start = 1
End = 214
Frequency = 1
1 2 3 4 5 6
90945.6890 -26431.6816 -152473.9555 55710.4893 -118660.0363 64123.7765
7 8 9 10 11 12
-95510.1935 54659.5163 69004.7594 -28123.8249 139124.4708 118966.5259
13 14 15 16 17 18
128188.4557 -39721.7390 -89436.6614 -31658.9483 -49185.9588 15337.1862
19 20 21 22 23 24
-110325.6568 -26920.0787 135674.2027 -75019.3021 147479.8459 82190.0877
25 26 27 28 29 30
110200.8293 -79175.1075 -80670.8124 -8919.5067 -132395.6533 -30695.8873
31 32 33 34 35 36
-151170.7823 -14914.5713 110968.5848 -50521.9376 96203.1530 64628.7423
37 38 39 40 41 42
23529.3743 -33446.4849 -43363.3649 -100191.3284 -98397.1124 16325.1713
43 44 45 46 47 48
-225453.3180 -23728.1457 138133.4213 -30470.8909 104975.4532 40576.0430
49 50 51 52 53 54
-32644.3365 -6235.5921 36555.0189 -28379.6699 -75757.2988 -3309.1513
55 56 57 58 59 60
-86306.5507 14306.0401 59663.0798 38782.5595 143370.5497 -4567.2274
61 62 63 64 65 66
453.0063 53549.4485 18886.6483 -46220.9533 -44882.8762 -4042.5790
67 68 69 70 71 72
-72903.7388 110985.6216 -11476.5180 8544.5494 86679.7326 1587.4990
73 74 75 76 77 78
-3960.9673 37939.0716 35983.3656 -48957.3089 -19828.3947 -30534.3341
79 80 81 82 83 84
-78691.5335 88608.4850 71926.7697 51384.1134 46920.3948 -222250.1779
85 86 87 88 89 90
10756.1608 176351.7003 85410.7913 2025.4305 -55265.3548 -52937.9954
91 92 93 94 95 96
-407.4535 20232.0154 83014.6718 19483.2485 -33675.7533 -115929.6101
97 98 99 100 101 102
17910.1839 -22968.5277 64799.9742 47314.7810 -117985.4785 4516.1914
103 104 105 106 107 108
-43857.0717 -33927.4872 -206026.0388 -293724.3211 61755.6975 79277.8501
109 110 111 112 113 114
-150.4262 18941.2861 55987.1307 -15217.2609 -6923.6792 -30484.1118
115 116 117 118 119 120
5287.7493 38602.6563 51213.7185 -47258.6252 -124432.2239 34021.0692
121 122 123 124 125 126
-2670.9085 -5209.2059 -154593.4596 -195608.0717 -65111.2965 116979.0506
127 128 129 130 131 132
61498.3824 -72793.0986 -83739.3648 58076.0874 -70893.9582 -64516.3751
133 134 135 136 137 138
-79240.1214 147766.5629 -24527.6230 51466.0230 -7173.8251 -35582.5990
139 140 141 142 143 144
69180.9646 -22599.5075 -24376.7115 62913.8086 -70722.6475 -42085.0054
145 146 147 148 149 150
48088.1086 -12251.9365 -52482.0974 108795.0843 182937.3208 -72166.1508
151 152 153 154 155 156
92660.8092 -2308.4495 44387.2132 41205.5935 -103700.9420 -23098.8464
157 158 159 160 161 162
41500.5171 -47586.9472 30618.6918 288586.6540 91303.7015 54683.3950
163 164 165 166 167 168
94476.7178 -10707.7616 26322.3494 61036.8960 -44563.8256 67212.4307
169 170 171 172 173 174
16543.1990 3291.0507 175235.5645 162661.7182 79736.1003 -10445.5213
175 176 177 178 179 180
269742.7892 81423.3615 -26043.8208 81068.5625 -59026.1680 38551.5690
181 182 183 184 185 186
-24132.5768 82599.4785 151346.5140 113730.9565 152419.6456 70942.5853
187 188 189 190 191 192
54998.2251 -32827.0687 -61945.2561 98093.5329 -192008.4879 -34138.0828
193 194 195 196 197 198
-181707.6964 -151739.3894 -91898.8012 277768.0944 -135608.0500 -188919.2089
199 200 201 202 203 204
102905.7018 15510.8803 -237591.8634 -79415.9118 -127485.2911 -20426.4917
205 206 207 208 209 210
-163608.4908 -95671.9868 34623.0762 -632906.1831 420778.2462 116210.1826
211 212 213 214
113874.9593 -183602.4075 -139109.1972 83945.8618
> postscript(file="/var/www/rcomp/tmp/6vj8x1292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> dum <- cbind(lag(myerror,k=1),myerror)
> dum
Time Series:
Start = 0
End = 214
Frequency = 1
lag(myerror, k = 1) myerror
0 90945.6890 NA
1 -26431.6816 90945.6890
2 -152473.9555 -26431.6816
3 55710.4893 -152473.9555
4 -118660.0363 55710.4893
5 64123.7765 -118660.0363
6 -95510.1935 64123.7765
7 54659.5163 -95510.1935
8 69004.7594 54659.5163
9 -28123.8249 69004.7594
10 139124.4708 -28123.8249
11 118966.5259 139124.4708
12 128188.4557 118966.5259
13 -39721.7390 128188.4557
14 -89436.6614 -39721.7390
15 -31658.9483 -89436.6614
16 -49185.9588 -31658.9483
17 15337.1862 -49185.9588
18 -110325.6568 15337.1862
19 -26920.0787 -110325.6568
20 135674.2027 -26920.0787
21 -75019.3021 135674.2027
22 147479.8459 -75019.3021
23 82190.0877 147479.8459
24 110200.8293 82190.0877
25 -79175.1075 110200.8293
26 -80670.8124 -79175.1075
27 -8919.5067 -80670.8124
28 -132395.6533 -8919.5067
29 -30695.8873 -132395.6533
30 -151170.7823 -30695.8873
31 -14914.5713 -151170.7823
32 110968.5848 -14914.5713
33 -50521.9376 110968.5848
34 96203.1530 -50521.9376
35 64628.7423 96203.1530
36 23529.3743 64628.7423
37 -33446.4849 23529.3743
38 -43363.3649 -33446.4849
39 -100191.3284 -43363.3649
40 -98397.1124 -100191.3284
41 16325.1713 -98397.1124
42 -225453.3180 16325.1713
43 -23728.1457 -225453.3180
44 138133.4213 -23728.1457
45 -30470.8909 138133.4213
46 104975.4532 -30470.8909
47 40576.0430 104975.4532
48 -32644.3365 40576.0430
49 -6235.5921 -32644.3365
50 36555.0189 -6235.5921
51 -28379.6699 36555.0189
52 -75757.2988 -28379.6699
53 -3309.1513 -75757.2988
54 -86306.5507 -3309.1513
55 14306.0401 -86306.5507
56 59663.0798 14306.0401
57 38782.5595 59663.0798
58 143370.5497 38782.5595
59 -4567.2274 143370.5497
60 453.0063 -4567.2274
61 53549.4485 453.0063
62 18886.6483 53549.4485
63 -46220.9533 18886.6483
64 -44882.8762 -46220.9533
65 -4042.5790 -44882.8762
66 -72903.7388 -4042.5790
67 110985.6216 -72903.7388
68 -11476.5180 110985.6216
69 8544.5494 -11476.5180
70 86679.7326 8544.5494
71 1587.4990 86679.7326
72 -3960.9673 1587.4990
73 37939.0716 -3960.9673
74 35983.3656 37939.0716
75 -48957.3089 35983.3656
76 -19828.3947 -48957.3089
77 -30534.3341 -19828.3947
78 -78691.5335 -30534.3341
79 88608.4850 -78691.5335
80 71926.7697 88608.4850
81 51384.1134 71926.7697
82 46920.3948 51384.1134
83 -222250.1779 46920.3948
84 10756.1608 -222250.1779
85 176351.7003 10756.1608
86 85410.7913 176351.7003
87 2025.4305 85410.7913
88 -55265.3548 2025.4305
89 -52937.9954 -55265.3548
90 -407.4535 -52937.9954
91 20232.0154 -407.4535
92 83014.6718 20232.0154
93 19483.2485 83014.6718
94 -33675.7533 19483.2485
95 -115929.6101 -33675.7533
96 17910.1839 -115929.6101
97 -22968.5277 17910.1839
98 64799.9742 -22968.5277
99 47314.7810 64799.9742
100 -117985.4785 47314.7810
101 4516.1914 -117985.4785
102 -43857.0717 4516.1914
103 -33927.4872 -43857.0717
104 -206026.0388 -33927.4872
105 -293724.3211 -206026.0388
106 61755.6975 -293724.3211
107 79277.8501 61755.6975
108 -150.4262 79277.8501
109 18941.2861 -150.4262
110 55987.1307 18941.2861
111 -15217.2609 55987.1307
112 -6923.6792 -15217.2609
113 -30484.1118 -6923.6792
114 5287.7493 -30484.1118
115 38602.6563 5287.7493
116 51213.7185 38602.6563
117 -47258.6252 51213.7185
118 -124432.2239 -47258.6252
119 34021.0692 -124432.2239
120 -2670.9085 34021.0692
121 -5209.2059 -2670.9085
122 -154593.4596 -5209.2059
123 -195608.0717 -154593.4596
124 -65111.2965 -195608.0717
125 116979.0506 -65111.2965
126 61498.3824 116979.0506
127 -72793.0986 61498.3824
128 -83739.3648 -72793.0986
129 58076.0874 -83739.3648
130 -70893.9582 58076.0874
131 -64516.3751 -70893.9582
132 -79240.1214 -64516.3751
133 147766.5629 -79240.1214
134 -24527.6230 147766.5629
135 51466.0230 -24527.6230
136 -7173.8251 51466.0230
137 -35582.5990 -7173.8251
138 69180.9646 -35582.5990
139 -22599.5075 69180.9646
140 -24376.7115 -22599.5075
141 62913.8086 -24376.7115
142 -70722.6475 62913.8086
143 -42085.0054 -70722.6475
144 48088.1086 -42085.0054
145 -12251.9365 48088.1086
146 -52482.0974 -12251.9365
147 108795.0843 -52482.0974
148 182937.3208 108795.0843
149 -72166.1508 182937.3208
150 92660.8092 -72166.1508
151 -2308.4495 92660.8092
152 44387.2132 -2308.4495
153 41205.5935 44387.2132
154 -103700.9420 41205.5935
155 -23098.8464 -103700.9420
156 41500.5171 -23098.8464
157 -47586.9472 41500.5171
158 30618.6918 -47586.9472
159 288586.6540 30618.6918
160 91303.7015 288586.6540
161 54683.3950 91303.7015
162 94476.7178 54683.3950
163 -10707.7616 94476.7178
164 26322.3494 -10707.7616
165 61036.8960 26322.3494
166 -44563.8256 61036.8960
167 67212.4307 -44563.8256
168 16543.1990 67212.4307
169 3291.0507 16543.1990
170 175235.5645 3291.0507
171 162661.7182 175235.5645
172 79736.1003 162661.7182
173 -10445.5213 79736.1003
174 269742.7892 -10445.5213
175 81423.3615 269742.7892
176 -26043.8208 81423.3615
177 81068.5625 -26043.8208
178 -59026.1680 81068.5625
179 38551.5690 -59026.1680
180 -24132.5768 38551.5690
181 82599.4785 -24132.5768
182 151346.5140 82599.4785
183 113730.9565 151346.5140
184 152419.6456 113730.9565
185 70942.5853 152419.6456
186 54998.2251 70942.5853
187 -32827.0687 54998.2251
188 -61945.2561 -32827.0687
189 98093.5329 -61945.2561
190 -192008.4879 98093.5329
191 -34138.0828 -192008.4879
192 -181707.6964 -34138.0828
193 -151739.3894 -181707.6964
194 -91898.8012 -151739.3894
195 277768.0944 -91898.8012
196 -135608.0500 277768.0944
197 -188919.2089 -135608.0500
198 102905.7018 -188919.2089
199 15510.8803 102905.7018
200 -237591.8634 15510.8803
201 -79415.9118 -237591.8634
202 -127485.2911 -79415.9118
203 -20426.4917 -127485.2911
204 -163608.4908 -20426.4917
205 -95671.9868 -163608.4908
206 34623.0762 -95671.9868
207 -632906.1831 34623.0762
208 420778.2462 -632906.1831
209 116210.1826 420778.2462
210 113874.9593 116210.1826
211 -183602.4075 113874.9593
212 -139109.1972 -183602.4075
213 83945.8618 -139109.1972
214 NA 83945.8618
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -26431.6816 90945.6890
[2,] -152473.9555 -26431.6816
[3,] 55710.4893 -152473.9555
[4,] -118660.0363 55710.4893
[5,] 64123.7765 -118660.0363
[6,] -95510.1935 64123.7765
[7,] 54659.5163 -95510.1935
[8,] 69004.7594 54659.5163
[9,] -28123.8249 69004.7594
[10,] 139124.4708 -28123.8249
[11,] 118966.5259 139124.4708
[12,] 128188.4557 118966.5259
[13,] -39721.7390 128188.4557
[14,] -89436.6614 -39721.7390
[15,] -31658.9483 -89436.6614
[16,] -49185.9588 -31658.9483
[17,] 15337.1862 -49185.9588
[18,] -110325.6568 15337.1862
[19,] -26920.0787 -110325.6568
[20,] 135674.2027 -26920.0787
[21,] -75019.3021 135674.2027
[22,] 147479.8459 -75019.3021
[23,] 82190.0877 147479.8459
[24,] 110200.8293 82190.0877
[25,] -79175.1075 110200.8293
[26,] -80670.8124 -79175.1075
[27,] -8919.5067 -80670.8124
[28,] -132395.6533 -8919.5067
[29,] -30695.8873 -132395.6533
[30,] -151170.7823 -30695.8873
[31,] -14914.5713 -151170.7823
[32,] 110968.5848 -14914.5713
[33,] -50521.9376 110968.5848
[34,] 96203.1530 -50521.9376
[35,] 64628.7423 96203.1530
[36,] 23529.3743 64628.7423
[37,] -33446.4849 23529.3743
[38,] -43363.3649 -33446.4849
[39,] -100191.3284 -43363.3649
[40,] -98397.1124 -100191.3284
[41,] 16325.1713 -98397.1124
[42,] -225453.3180 16325.1713
[43,] -23728.1457 -225453.3180
[44,] 138133.4213 -23728.1457
[45,] -30470.8909 138133.4213
[46,] 104975.4532 -30470.8909
[47,] 40576.0430 104975.4532
[48,] -32644.3365 40576.0430
[49,] -6235.5921 -32644.3365
[50,] 36555.0189 -6235.5921
[51,] -28379.6699 36555.0189
[52,] -75757.2988 -28379.6699
[53,] -3309.1513 -75757.2988
[54,] -86306.5507 -3309.1513
[55,] 14306.0401 -86306.5507
[56,] 59663.0798 14306.0401
[57,] 38782.5595 59663.0798
[58,] 143370.5497 38782.5595
[59,] -4567.2274 143370.5497
[60,] 453.0063 -4567.2274
[61,] 53549.4485 453.0063
[62,] 18886.6483 53549.4485
[63,] -46220.9533 18886.6483
[64,] -44882.8762 -46220.9533
[65,] -4042.5790 -44882.8762
[66,] -72903.7388 -4042.5790
[67,] 110985.6216 -72903.7388
[68,] -11476.5180 110985.6216
[69,] 8544.5494 -11476.5180
[70,] 86679.7326 8544.5494
[71,] 1587.4990 86679.7326
[72,] -3960.9673 1587.4990
[73,] 37939.0716 -3960.9673
[74,] 35983.3656 37939.0716
[75,] -48957.3089 35983.3656
[76,] -19828.3947 -48957.3089
[77,] -30534.3341 -19828.3947
[78,] -78691.5335 -30534.3341
[79,] 88608.4850 -78691.5335
[80,] 71926.7697 88608.4850
[81,] 51384.1134 71926.7697
[82,] 46920.3948 51384.1134
[83,] -222250.1779 46920.3948
[84,] 10756.1608 -222250.1779
[85,] 176351.7003 10756.1608
[86,] 85410.7913 176351.7003
[87,] 2025.4305 85410.7913
[88,] -55265.3548 2025.4305
[89,] -52937.9954 -55265.3548
[90,] -407.4535 -52937.9954
[91,] 20232.0154 -407.4535
[92,] 83014.6718 20232.0154
[93,] 19483.2485 83014.6718
[94,] -33675.7533 19483.2485
[95,] -115929.6101 -33675.7533
[96,] 17910.1839 -115929.6101
[97,] -22968.5277 17910.1839
[98,] 64799.9742 -22968.5277
[99,] 47314.7810 64799.9742
[100,] -117985.4785 47314.7810
[101,] 4516.1914 -117985.4785
[102,] -43857.0717 4516.1914
[103,] -33927.4872 -43857.0717
[104,] -206026.0388 -33927.4872
[105,] -293724.3211 -206026.0388
[106,] 61755.6975 -293724.3211
[107,] 79277.8501 61755.6975
[108,] -150.4262 79277.8501
[109,] 18941.2861 -150.4262
[110,] 55987.1307 18941.2861
[111,] -15217.2609 55987.1307
[112,] -6923.6792 -15217.2609
[113,] -30484.1118 -6923.6792
[114,] 5287.7493 -30484.1118
[115,] 38602.6563 5287.7493
[116,] 51213.7185 38602.6563
[117,] -47258.6252 51213.7185
[118,] -124432.2239 -47258.6252
[119,] 34021.0692 -124432.2239
[120,] -2670.9085 34021.0692
[121,] -5209.2059 -2670.9085
[122,] -154593.4596 -5209.2059
[123,] -195608.0717 -154593.4596
[124,] -65111.2965 -195608.0717
[125,] 116979.0506 -65111.2965
[126,] 61498.3824 116979.0506
[127,] -72793.0986 61498.3824
[128,] -83739.3648 -72793.0986
[129,] 58076.0874 -83739.3648
[130,] -70893.9582 58076.0874
[131,] -64516.3751 -70893.9582
[132,] -79240.1214 -64516.3751
[133,] 147766.5629 -79240.1214
[134,] -24527.6230 147766.5629
[135,] 51466.0230 -24527.6230
[136,] -7173.8251 51466.0230
[137,] -35582.5990 -7173.8251
[138,] 69180.9646 -35582.5990
[139,] -22599.5075 69180.9646
[140,] -24376.7115 -22599.5075
[141,] 62913.8086 -24376.7115
[142,] -70722.6475 62913.8086
[143,] -42085.0054 -70722.6475
[144,] 48088.1086 -42085.0054
[145,] -12251.9365 48088.1086
[146,] -52482.0974 -12251.9365
[147,] 108795.0843 -52482.0974
[148,] 182937.3208 108795.0843
[149,] -72166.1508 182937.3208
[150,] 92660.8092 -72166.1508
[151,] -2308.4495 92660.8092
[152,] 44387.2132 -2308.4495
[153,] 41205.5935 44387.2132
[154,] -103700.9420 41205.5935
[155,] -23098.8464 -103700.9420
[156,] 41500.5171 -23098.8464
[157,] -47586.9472 41500.5171
[158,] 30618.6918 -47586.9472
[159,] 288586.6540 30618.6918
[160,] 91303.7015 288586.6540
[161,] 54683.3950 91303.7015
[162,] 94476.7178 54683.3950
[163,] -10707.7616 94476.7178
[164,] 26322.3494 -10707.7616
[165,] 61036.8960 26322.3494
[166,] -44563.8256 61036.8960
[167,] 67212.4307 -44563.8256
[168,] 16543.1990 67212.4307
[169,] 3291.0507 16543.1990
[170,] 175235.5645 3291.0507
[171,] 162661.7182 175235.5645
[172,] 79736.1003 162661.7182
[173,] -10445.5213 79736.1003
[174,] 269742.7892 -10445.5213
[175,] 81423.3615 269742.7892
[176,] -26043.8208 81423.3615
[177,] 81068.5625 -26043.8208
[178,] -59026.1680 81068.5625
[179,] 38551.5690 -59026.1680
[180,] -24132.5768 38551.5690
[181,] 82599.4785 -24132.5768
[182,] 151346.5140 82599.4785
[183,] 113730.9565 151346.5140
[184,] 152419.6456 113730.9565
[185,] 70942.5853 152419.6456
[186,] 54998.2251 70942.5853
[187,] -32827.0687 54998.2251
[188,] -61945.2561 -32827.0687
[189,] 98093.5329 -61945.2561
[190,] -192008.4879 98093.5329
[191,] -34138.0828 -192008.4879
[192,] -181707.6964 -34138.0828
[193,] -151739.3894 -181707.6964
[194,] -91898.8012 -151739.3894
[195,] 277768.0944 -91898.8012
[196,] -135608.0500 277768.0944
[197,] -188919.2089 -135608.0500
[198,] 102905.7018 -188919.2089
[199,] 15510.8803 102905.7018
[200,] -237591.8634 15510.8803
[201,] -79415.9118 -237591.8634
[202,] -127485.2911 -79415.9118
[203,] -20426.4917 -127485.2911
[204,] -163608.4908 -20426.4917
[205,] -95671.9868 -163608.4908
[206,] 34623.0762 -95671.9868
[207,] -632906.1831 34623.0762
[208,] 420778.2462 -632906.1831
[209,] 116210.1826 420778.2462
[210,] 113874.9593 116210.1826
[211,] -183602.4075 113874.9593
[212,] -139109.1972 -183602.4075
[213,] 83945.8618 -139109.1972
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -26431.6816 90945.6890
2 -152473.9555 -26431.6816
3 55710.4893 -152473.9555
4 -118660.0363 55710.4893
5 64123.7765 -118660.0363
6 -95510.1935 64123.7765
7 54659.5163 -95510.1935
8 69004.7594 54659.5163
9 -28123.8249 69004.7594
10 139124.4708 -28123.8249
11 118966.5259 139124.4708
12 128188.4557 118966.5259
13 -39721.7390 128188.4557
14 -89436.6614 -39721.7390
15 -31658.9483 -89436.6614
16 -49185.9588 -31658.9483
17 15337.1862 -49185.9588
18 -110325.6568 15337.1862
19 -26920.0787 -110325.6568
20 135674.2027 -26920.0787
21 -75019.3021 135674.2027
22 147479.8459 -75019.3021
23 82190.0877 147479.8459
24 110200.8293 82190.0877
25 -79175.1075 110200.8293
26 -80670.8124 -79175.1075
27 -8919.5067 -80670.8124
28 -132395.6533 -8919.5067
29 -30695.8873 -132395.6533
30 -151170.7823 -30695.8873
31 -14914.5713 -151170.7823
32 110968.5848 -14914.5713
33 -50521.9376 110968.5848
34 96203.1530 -50521.9376
35 64628.7423 96203.1530
36 23529.3743 64628.7423
37 -33446.4849 23529.3743
38 -43363.3649 -33446.4849
39 -100191.3284 -43363.3649
40 -98397.1124 -100191.3284
41 16325.1713 -98397.1124
42 -225453.3180 16325.1713
43 -23728.1457 -225453.3180
44 138133.4213 -23728.1457
45 -30470.8909 138133.4213
46 104975.4532 -30470.8909
47 40576.0430 104975.4532
48 -32644.3365 40576.0430
49 -6235.5921 -32644.3365
50 36555.0189 -6235.5921
51 -28379.6699 36555.0189
52 -75757.2988 -28379.6699
53 -3309.1513 -75757.2988
54 -86306.5507 -3309.1513
55 14306.0401 -86306.5507
56 59663.0798 14306.0401
57 38782.5595 59663.0798
58 143370.5497 38782.5595
59 -4567.2274 143370.5497
60 453.0063 -4567.2274
61 53549.4485 453.0063
62 18886.6483 53549.4485
63 -46220.9533 18886.6483
64 -44882.8762 -46220.9533
65 -4042.5790 -44882.8762
66 -72903.7388 -4042.5790
67 110985.6216 -72903.7388
68 -11476.5180 110985.6216
69 8544.5494 -11476.5180
70 86679.7326 8544.5494
71 1587.4990 86679.7326
72 -3960.9673 1587.4990
73 37939.0716 -3960.9673
74 35983.3656 37939.0716
75 -48957.3089 35983.3656
76 -19828.3947 -48957.3089
77 -30534.3341 -19828.3947
78 -78691.5335 -30534.3341
79 88608.4850 -78691.5335
80 71926.7697 88608.4850
81 51384.1134 71926.7697
82 46920.3948 51384.1134
83 -222250.1779 46920.3948
84 10756.1608 -222250.1779
85 176351.7003 10756.1608
86 85410.7913 176351.7003
87 2025.4305 85410.7913
88 -55265.3548 2025.4305
89 -52937.9954 -55265.3548
90 -407.4535 -52937.9954
91 20232.0154 -407.4535
92 83014.6718 20232.0154
93 19483.2485 83014.6718
94 -33675.7533 19483.2485
95 -115929.6101 -33675.7533
96 17910.1839 -115929.6101
97 -22968.5277 17910.1839
98 64799.9742 -22968.5277
99 47314.7810 64799.9742
100 -117985.4785 47314.7810
101 4516.1914 -117985.4785
102 -43857.0717 4516.1914
103 -33927.4872 -43857.0717
104 -206026.0388 -33927.4872
105 -293724.3211 -206026.0388
106 61755.6975 -293724.3211
107 79277.8501 61755.6975
108 -150.4262 79277.8501
109 18941.2861 -150.4262
110 55987.1307 18941.2861
111 -15217.2609 55987.1307
112 -6923.6792 -15217.2609
113 -30484.1118 -6923.6792
114 5287.7493 -30484.1118
115 38602.6563 5287.7493
116 51213.7185 38602.6563
117 -47258.6252 51213.7185
118 -124432.2239 -47258.6252
119 34021.0692 -124432.2239
120 -2670.9085 34021.0692
121 -5209.2059 -2670.9085
122 -154593.4596 -5209.2059
123 -195608.0717 -154593.4596
124 -65111.2965 -195608.0717
125 116979.0506 -65111.2965
126 61498.3824 116979.0506
127 -72793.0986 61498.3824
128 -83739.3648 -72793.0986
129 58076.0874 -83739.3648
130 -70893.9582 58076.0874
131 -64516.3751 -70893.9582
132 -79240.1214 -64516.3751
133 147766.5629 -79240.1214
134 -24527.6230 147766.5629
135 51466.0230 -24527.6230
136 -7173.8251 51466.0230
137 -35582.5990 -7173.8251
138 69180.9646 -35582.5990
139 -22599.5075 69180.9646
140 -24376.7115 -22599.5075
141 62913.8086 -24376.7115
142 -70722.6475 62913.8086
143 -42085.0054 -70722.6475
144 48088.1086 -42085.0054
145 -12251.9365 48088.1086
146 -52482.0974 -12251.9365
147 108795.0843 -52482.0974
148 182937.3208 108795.0843
149 -72166.1508 182937.3208
150 92660.8092 -72166.1508
151 -2308.4495 92660.8092
152 44387.2132 -2308.4495
153 41205.5935 44387.2132
154 -103700.9420 41205.5935
155 -23098.8464 -103700.9420
156 41500.5171 -23098.8464
157 -47586.9472 41500.5171
158 30618.6918 -47586.9472
159 288586.6540 30618.6918
160 91303.7015 288586.6540
161 54683.3950 91303.7015
162 94476.7178 54683.3950
163 -10707.7616 94476.7178
164 26322.3494 -10707.7616
165 61036.8960 26322.3494
166 -44563.8256 61036.8960
167 67212.4307 -44563.8256
168 16543.1990 67212.4307
169 3291.0507 16543.1990
170 175235.5645 3291.0507
171 162661.7182 175235.5645
172 79736.1003 162661.7182
173 -10445.5213 79736.1003
174 269742.7892 -10445.5213
175 81423.3615 269742.7892
176 -26043.8208 81423.3615
177 81068.5625 -26043.8208
178 -59026.1680 81068.5625
179 38551.5690 -59026.1680
180 -24132.5768 38551.5690
181 82599.4785 -24132.5768
182 151346.5140 82599.4785
183 113730.9565 151346.5140
184 152419.6456 113730.9565
185 70942.5853 152419.6456
186 54998.2251 70942.5853
187 -32827.0687 54998.2251
188 -61945.2561 -32827.0687
189 98093.5329 -61945.2561
190 -192008.4879 98093.5329
191 -34138.0828 -192008.4879
192 -181707.6964 -34138.0828
193 -151739.3894 -181707.6964
194 -91898.8012 -151739.3894
195 277768.0944 -91898.8012
196 -135608.0500 277768.0944
197 -188919.2089 -135608.0500
198 102905.7018 -188919.2089
199 15510.8803 102905.7018
200 -237591.8634 15510.8803
201 -79415.9118 -237591.8634
202 -127485.2911 -79415.9118
203 -20426.4917 -127485.2911
204 -163608.4908 -20426.4917
205 -95671.9868 -163608.4908
206 34623.0762 -95671.9868
207 -632906.1831 34623.0762
208 420778.2462 -632906.1831
209 116210.1826 420778.2462
210 113874.9593 116210.1826
211 -183602.4075 113874.9593
212 -139109.1972 -183602.4075
213 83945.8618 -139109.1972
> 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/rcomp/tmp/7ns801292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/8ns801292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/9g2p41292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> 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
> if (n > n25) {
+ postscript(file="/var/www/rcomp/tmp/10g2p41292011607.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/11uc841292011608.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/rcomp/tmp/12xvps1292011608.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/rcomp/tmp/134em41292011608.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/rcomp/tmp/14f5l71292011608.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/150n1u1292011608.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/16exz31292011608.tab")
+ }
> try(system("convert tmp/191aa1292011607.ps tmp/191aa1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/2karv1292011607.ps tmp/2karv1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/3karv1292011607.ps tmp/3karv1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/4karv1292011607.ps tmp/4karv1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/5vj8x1292011607.ps tmp/5vj8x1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/6vj8x1292011607.ps tmp/6vj8x1292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/7ns801292011607.ps tmp/7ns801292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ns801292011607.ps tmp/8ns801292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/9g2p41292011607.ps tmp/9g2p41292011607.png",intern=TRUE))
character(0)
> try(system("convert tmp/10g2p41292011607.ps tmp/10g2p41292011607.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
6.60 0.89 7.49