R version 2.10.1 (2009-12-14)
Copyright (C) 2009 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(349
+ ,0
+ ,0
+ ,0
+ ,0
+ ,336
+ ,349
+ ,0
+ ,0
+ ,0
+ ,331
+ ,336
+ ,349
+ ,0
+ ,0
+ ,327
+ ,331
+ ,336
+ ,349
+ ,0
+ ,323
+ ,327
+ ,331
+ ,336
+ ,349
+ ,322
+ ,323
+ ,327
+ ,331
+ ,336
+ ,385
+ ,322
+ ,323
+ ,327
+ ,331
+ ,405
+ ,385
+ ,322
+ ,323
+ ,327
+ ,412
+ ,405
+ ,385
+ ,322
+ ,323
+ ,411
+ ,412
+ ,405
+ ,385
+ ,322
+ ,410
+ ,411
+ ,412
+ ,405
+ ,385
+ ,415
+ ,410
+ ,411
+ ,412
+ ,405
+ ,414
+ ,415
+ ,410
+ ,411
+ ,412
+ ,411
+ ,414
+ ,415
+ ,410
+ ,411
+ ,408
+ ,411
+ ,414
+ ,415
+ ,410
+ ,410
+ ,408
+ ,411
+ ,414
+ ,415
+ ,411
+ ,410
+ ,408
+ ,411
+ ,414
+ ,416
+ ,411
+ ,410
+ ,408
+ ,411
+ ,479
+ ,416
+ ,411
+ ,410
+ ,408
+ ,498
+ ,479
+ ,416
+ ,411
+ ,410
+ ,502
+ ,498
+ ,479
+ ,416
+ ,411
+ ,498
+ ,502
+ ,498
+ ,479
+ ,416
+ ,499
+ ,498
+ ,502
+ ,498
+ ,479
+ ,506
+ ,499
+ ,498
+ ,502
+ ,498
+ ,510
+ ,506
+ ,499
+ ,498
+ ,502
+ ,509
+ ,510
+ ,506
+ ,499
+ ,498
+ ,502
+ ,509
+ ,510
+ ,506
+ ,499
+ ,495
+ ,502
+ ,509
+ ,510
+ ,506
+ ,490
+ ,495
+ ,502
+ ,509
+ ,510
+ ,490
+ ,490
+ ,495
+ ,502
+ ,509
+ ,553
+ ,490
+ ,490
+ ,495
+ ,502
+ ,570
+ ,553
+ ,490
+ ,490
+ ,495
+ ,573
+ ,570
+ ,553
+ ,490
+ ,490
+ ,572
+ ,573
+ ,570
+ ,553
+ ,490
+ ,575
+ ,572
+ ,573
+ ,570
+ ,553
+ ,580
+ ,575
+ ,572
+ ,573
+ ,570
+ ,580
+ ,580
+ ,575
+ ,572
+ ,573
+ ,574
+ ,580
+ ,580
+ ,575
+ ,572
+ ,563
+ ,574
+ ,580
+ ,580
+ ,575
+ ,556
+ ,563
+ ,574
+ ,580
+ ,580
+ ,546
+ ,556
+ ,563
+ ,574
+ ,580
+ ,545
+ ,546
+ ,556
+ ,563
+ ,574
+ ,605
+ ,545
+ ,546
+ ,556
+ ,563
+ ,628
+ ,605
+ ,545
+ ,546
+ ,556
+ ,631
+ ,628
+ ,605
+ ,545
+ ,546
+ ,626
+ ,631
+ ,628
+ ,605
+ ,545
+ ,614
+ ,626
+ ,631
+ ,628
+ ,605
+ ,606
+ ,614
+ ,626
+ ,631
+ ,628
+ ,602
+ ,606
+ ,614
+ ,626
+ ,631
+ ,589
+ ,602
+ ,606
+ ,614
+ ,626
+ ,574
+ ,589
+ ,602
+ ,606
+ ,614
+ ,558
+ ,574
+ ,589
+ ,602
+ ,606
+ ,552
+ ,558
+ ,574
+ ,589
+ ,602
+ ,546
+ ,552
+ ,558
+ ,574
+ ,589
+ ,607
+ ,546
+ ,552
+ ,558
+ ,574
+ ,636
+ ,607
+ ,546
+ ,552
+ ,558
+ ,631
+ ,636
+ ,607
+ ,546
+ ,552
+ ,623
+ ,631
+ ,636
+ ,607
+ ,546
+ ,618
+ ,623
+ ,631
+ ,636
+ ,607
+ ,605
+ ,618
+ ,623
+ ,631
+ ,636
+ ,619
+ ,605
+ ,618
+ ,623
+ ,631
+ ,596
+ ,619
+ ,605
+ ,618
+ ,623
+ ,570
+ ,596
+ ,619
+ ,605
+ ,618
+ ,546
+ ,570
+ ,596
+ ,619
+ ,605
+ ,528
+ ,546
+ ,570
+ ,596
+ ,619
+ ,506
+ ,528
+ ,546
+ ,570
+ ,596
+ ,555
+ ,506
+ ,528
+ ,546
+ ,570
+ ,568
+ ,555
+ ,506
+ ,528
+ ,546
+ ,564
+ ,568
+ ,555
+ ,506
+ ,528
+ ,553
+ ,564
+ ,568
+ ,555
+ ,506
+ ,541
+ ,553
+ ,564
+ ,568
+ ,555
+ ,542
+ ,541
+ ,553
+ ,564
+ ,568
+ ,540
+ ,542
+ ,541
+ ,553
+ ,564
+ ,521
+ ,540
+ ,542
+ ,541
+ ,553
+ ,505
+ ,521
+ ,540
+ ,542
+ ,541
+ ,491
+ ,505
+ ,521
+ ,540
+ ,542
+ ,482
+ ,491
+ ,505
+ ,521
+ ,540
+ ,478
+ ,482
+ ,491
+ ,505
+ ,521
+ ,523
+ ,478
+ ,482
+ ,491
+ ,505
+ ,531
+ ,523
+ ,478
+ ,482
+ ,491
+ ,532
+ ,531
+ ,523
+ ,478
+ ,482
+ ,540
+ ,532
+ ,531
+ ,523
+ ,478
+ ,525
+ ,540
+ ,532
+ ,531
+ ,523
+ ,533
+ ,525
+ ,540
+ ,532
+ ,531
+ ,531
+ ,533
+ ,525
+ ,540
+ ,532
+ ,508
+ ,531
+ ,533
+ ,525
+ ,540
+ ,495
+ ,508
+ ,531
+ ,533
+ ,525
+ ,482
+ ,495
+ ,508
+ ,531
+ ,533
+ ,470
+ ,482
+ ,495
+ ,508
+ ,531
+ ,466
+ ,470
+ ,482
+ ,495
+ ,508
+ ,515
+ ,466
+ ,470
+ ,482
+ ,495
+ ,518
+ ,515
+ ,466
+ ,470
+ ,482
+ ,516
+ ,518
+ ,515
+ ,466
+ ,470
+ ,511
+ ,516
+ ,518
+ ,515
+ ,466
+ ,500
+ ,511
+ ,516
+ ,518
+ ,515
+ ,498
+ ,500
+ ,511
+ ,516
+ ,518
+ ,494
+ ,498
+ ,500
+ ,511
+ ,516
+ ,476
+ ,494
+ ,498
+ ,500
+ ,511
+ ,458
+ ,476
+ ,494
+ ,498
+ ,500
+ ,443
+ ,458
+ ,476
+ ,494
+ ,498
+ ,430
+ ,443
+ ,458
+ ,476
+ ,494
+ ,424
+ ,430
+ ,443
+ ,458
+ ,476
+ ,476
+ ,424
+ ,430
+ ,443
+ ,458
+ ,481
+ ,476
+ ,424
+ ,430
+ ,443
+ ,470
+ ,481
+ ,476
+ ,424
+ ,430
+ ,460
+ ,470
+ ,481
+ ,476
+ ,424
+ ,451
+ ,460
+ ,470
+ ,481
+ ,476
+ ,450
+ ,451
+ ,460
+ ,470
+ ,481
+ ,444
+ ,450
+ ,451
+ ,460
+ ,470
+ ,429
+ ,444
+ ,450
+ ,451
+ ,460
+ ,421
+ ,429
+ ,444
+ ,450
+ ,451
+ ,400
+ ,421
+ ,429
+ ,444
+ ,450
+ ,389
+ ,400
+ ,421
+ ,429
+ ,444
+ ,384
+ ,389
+ ,400
+ ,421
+ ,429
+ ,432
+ ,384
+ ,389
+ ,400
+ ,421
+ ,446
+ ,432
+ ,384
+ ,389
+ ,400
+ ,431
+ ,446
+ ,432
+ ,384
+ ,389
+ ,423
+ ,431
+ ,446
+ ,432
+ ,384
+ ,416
+ ,423
+ ,431
+ ,446
+ ,432
+ ,416
+ ,416
+ ,423
+ ,431
+ ,446
+ ,413
+ ,416
+ ,416
+ ,423
+ ,431
+ ,399
+ ,413
+ ,416
+ ,416
+ ,423
+ ,386
+ ,399
+ ,413
+ ,416
+ ,416
+ ,374
+ ,386
+ ,399
+ ,413
+ ,416
+ ,365
+ ,374
+ ,386
+ ,399
+ ,413
+ ,365
+ ,365
+ ,374
+ ,386
+ ,399
+ ,418
+ ,365
+ ,365
+ ,374
+ ,386
+ ,428
+ ,418
+ ,365
+ ,365
+ ,374
+ ,424
+ ,428
+ ,418
+ ,365
+ ,365
+ ,421
+ ,424
+ ,428
+ ,418
+ ,365
+ ,417
+ ,421
+ ,424
+ ,428
+ ,418
+ ,423
+ ,417
+ ,421
+ ,424
+ ,428
+ ,423
+ ,423
+ ,417
+ ,421
+ ,424
+ ,419
+ ,423
+ ,423
+ ,417
+ ,421
+ ,406
+ ,419
+ ,423
+ ,423
+ ,417
+ ,398
+ ,406
+ ,419
+ ,423
+ ,423
+ ,390
+ ,398
+ ,406
+ ,419
+ ,423
+ ,391
+ ,390
+ ,398
+ ,406
+ ,419
+ ,444
+ ,391
+ ,390
+ ,398
+ ,406
+ ,460
+ ,444
+ ,391
+ ,390
+ ,398
+ ,455
+ ,460
+ ,444
+ ,391
+ ,390
+ ,456
+ ,455
+ ,460
+ ,444
+ ,391
+ ,452
+ ,456
+ ,455
+ ,460
+ ,444
+ ,459
+ ,452
+ ,456
+ ,455
+ ,460
+ ,461
+ ,459
+ ,452
+ ,456
+ ,455
+ ,451
+ ,461
+ ,459
+ ,452
+ ,456
+ ,443
+ ,451
+ ,461
+ ,459
+ ,452
+ ,439
+ ,443
+ ,451
+ ,461
+ ,459
+ ,430
+ ,439
+ ,443
+ ,451
+ ,461
+ ,436
+ ,430
+ ,439
+ ,443
+ ,451
+ ,488
+ ,436
+ ,430
+ ,439
+ ,443
+ ,506
+ ,488
+ ,436
+ ,430
+ ,439
+ ,502
+ ,506
+ ,488
+ ,436
+ ,430
+ ,501
+ ,502
+ ,506
+ ,488
+ ,436
+ ,501
+ ,501
+ ,502
+ ,506
+ ,488
+ ,515
+ ,501
+ ,501
+ ,502
+ ,506
+ ,521
+ ,515
+ ,501
+ ,501
+ ,502
+ ,520
+ ,521
+ ,515
+ ,501
+ ,501
+ ,512
+ ,520
+ ,521
+ ,515
+ ,501
+ ,509
+ ,512
+ ,520
+ ,521
+ ,515
+ ,505
+ ,509
+ ,512
+ ,520
+ ,521
+ ,511
+ ,505
+ ,509
+ ,512
+ ,520
+ ,570
+ ,511
+ ,505
+ ,509
+ ,512
+ ,592
+ ,570
+ ,511
+ ,505
+ ,509
+ ,594
+ ,592
+ ,570
+ ,511
+ ,505
+ ,586
+ ,594
+ ,592
+ ,570
+ ,511
+ ,586
+ ,586
+ ,594
+ ,592
+ ,570
+ ,592
+ ,586
+ ,586
+ ,594
+ ,592
+ ,594
+ ,592
+ ,586
+ ,586
+ ,594
+ ,594
+ ,594
+ ,592
+ ,586
+ ,586
+ ,586
+ ,594
+ ,594
+ ,592
+ ,586
+ ,586
+ ,586
+ ,594
+ ,594
+ ,592
+ ,572
+ ,586
+ ,586
+ ,594
+ ,594
+ ,572
+ ,572
+ ,586
+ ,586
+ ,594
+ ,563
+ ,572
+ ,572
+ ,586
+ ,586
+ ,563
+ ,563
+ ,572
+ ,572
+ ,586
+ ,555
+ ,563
+ ,563
+ ,572
+ ,572
+ ,555
+ ,555
+ ,563
+ ,563
+ ,572
+ ,554
+ ,555
+ ,555
+ ,563
+ ,563
+ ,554
+ ,554
+ ,555
+ ,555
+ ,563
+ ,601
+ ,554
+ ,554
+ ,555
+ ,555
+ ,601
+ ,601
+ ,554
+ ,554
+ ,555
+ ,622
+ ,601
+ ,601
+ ,554
+ ,554
+ ,622
+ ,622
+ ,601
+ ,601
+ ,554
+ ,617
+ ,622
+ ,622
+ ,601
+ ,601
+ ,617
+ ,617
+ ,622
+ ,622
+ ,601
+ ,606
+ ,617
+ ,617
+ ,622
+ ,622
+ ,606
+ ,606
+ ,617
+ ,617
+ ,622
+ ,595
+ ,606
+ ,606
+ ,617
+ ,617
+ ,595
+ ,595
+ ,606
+ ,606
+ ,617
+ ,599
+ ,595
+ ,595
+ ,606
+ ,606
+ ,599
+ ,599
+ ,595
+ ,595
+ ,606
+ ,600
+ ,599
+ ,599
+ ,595
+ ,595
+ ,600
+ ,600
+ ,599
+ ,599
+ ,595
+ ,592
+ ,600
+ ,600
+ ,599
+ ,599
+ ,592
+ ,592
+ ,600
+ ,600
+ ,599
+ ,575
+ ,592
+ ,592
+ ,600
+ ,600
+ ,575
+ ,575
+ ,592
+ ,592
+ ,600
+ ,567
+ ,575
+ ,575
+ ,592
+ ,592
+ ,567
+ ,567
+ ,575
+ ,575
+ ,592
+ ,555
+ ,567
+ ,567
+ ,575
+ ,575
+ ,555
+ ,555
+ ,567
+ ,567
+ ,575
+ ,555
+ ,555
+ ,555
+ ,567
+ ,567
+ ,555
+ ,555
+ ,555
+ ,555
+ ,567
+ ,608
+ ,555
+ ,555
+ ,555
+ ,555
+ ,608
+ ,608
+ ,555
+ ,555
+ ,555
+ ,631
+ ,608
+ ,608
+ ,555
+ ,555
+ ,631
+ ,631
+ ,608
+ ,608
+ ,555
+ ,629
+ ,631
+ ,631
+ ,608
+ ,608
+ ,629
+ ,629
+ ,631
+ ,631
+ ,608
+ ,624
+ ,629
+ ,629
+ ,631
+ ,631
+ ,624
+ ,624
+ ,629
+ ,629
+ ,631
+ ,610
+ ,624
+ ,624
+ ,629
+ ,629
+ ,610
+ ,610
+ ,624
+ ,624
+ ,629
+ ,616
+ ,610
+ ,610
+ ,624
+ ,624
+ ,616
+ ,616
+ ,610
+ ,610
+ ,624
+ ,621
+ ,616
+ ,616
+ ,610
+ ,610
+ ,621
+ ,621
+ ,616
+ ,616
+ ,610
+ ,604
+ ,621
+ ,621
+ ,616
+ ,616
+ ,604
+ ,604
+ ,621
+ ,621
+ ,616
+ ,584
+ ,604
+ ,604
+ ,621
+ ,621
+ ,584
+ ,584
+ ,604
+ ,604
+ ,621
+ ,574
+ ,584
+ ,584
+ ,604
+ ,604
+ ,574
+ ,574
+ ,584
+ ,584
+ ,604
+ ,555
+ ,574
+ ,574
+ ,584
+ ,584
+ ,555
+ ,555
+ ,574
+ ,574
+ ,584
+ ,545
+ ,555
+ ,555
+ ,574
+ ,574
+ ,545
+ ,545
+ ,555
+ ,555
+ ,574
+ ,599
+ ,545
+ ,545
+ ,555
+ ,555
+ ,599
+ ,599
+ ,545
+ ,545
+ ,555
+ ,620
+ ,599
+ ,599
+ ,545
+ ,545
+ ,620
+ ,620
+ ,599
+ ,599
+ ,545
+ ,608
+ ,620
+ ,620
+ ,599
+ ,599
+ ,608
+ ,608
+ ,620
+ ,620
+ ,599
+ ,590
+ ,608
+ ,608
+ ,620
+ ,620
+ ,590
+ ,590
+ ,608
+ ,608
+ ,620
+ ,579
+ ,590
+ ,590
+ ,608
+ ,608
+ ,579
+ ,579
+ ,590
+ ,590
+ ,608
+ ,580
+ ,579
+ ,579
+ ,590
+ ,590
+ ,580
+ ,580
+ ,579
+ ,579
+ ,590
+ ,579
+ ,580
+ ,580
+ ,579
+ ,579
+ ,579
+ ,579
+ ,580
+ ,580
+ ,579
+ ,572
+ ,579
+ ,579
+ ,580
+ ,580
+ ,572
+ ,572
+ ,579
+ ,579
+ ,580
+ ,560
+ ,572
+ ,572
+ ,579
+ ,579
+ ,560
+ ,560
+ ,572
+ ,572
+ ,579
+ ,551
+ ,560
+ ,560
+ ,572
+ ,572
+ ,551
+ ,551
+ ,560
+ ,560
+ ,572
+ ,537
+ ,551
+ ,551
+ ,560
+ ,560
+ ,537
+ ,537
+ ,551
+ ,551
+ ,560
+ ,541
+ ,537
+ ,537
+ ,551
+ ,551
+ ,541
+ ,541
+ ,537
+ ,537
+ ,551
+ ,588
+ ,541
+ ,541
+ ,537
+ ,537
+ ,588
+ ,588
+ ,541
+ ,541
+ ,537
+ ,607
+ ,588
+ ,588
+ ,541
+ ,541
+ ,607
+ ,607
+ ,588
+ ,588
+ ,541
+ ,599
+ ,607
+ ,607
+ ,588
+ ,588
+ ,599
+ ,599
+ ,607
+ ,607
+ ,588
+ ,578
+ ,599
+ ,599
+ ,607
+ ,607
+ ,578
+ ,578
+ ,599
+ ,599
+ ,607
+ ,563
+ ,578
+ ,578
+ ,599
+ ,599
+ ,563
+ ,563
+ ,578
+ ,578
+ ,599
+ ,566
+ ,563
+ ,563
+ ,578
+ ,578
+ ,566
+ ,566
+ ,563
+ ,563
+ ,578
+ ,561
+ ,566
+ ,566
+ ,563
+ ,563
+ ,561
+ ,561
+ ,566
+ ,566
+ ,563
+ ,554
+ ,561
+ ,561
+ ,566
+ ,566
+ ,554
+ ,554
+ ,561
+ ,561
+ ,566
+ ,540
+ ,554
+ ,554
+ ,561
+ ,561
+ ,540
+ ,540
+ ,554
+ ,554
+ ,561
+ ,526
+ ,540
+ ,540
+ ,554
+ ,554
+ ,526
+ ,526
+ ,540
+ ,540
+ ,554
+ ,512
+ ,526
+ ,526
+ ,540
+ ,540
+ ,512
+ ,512
+ ,526
+ ,526
+ ,540
+ ,505
+ ,512
+ ,512
+ ,526
+ ,526
+ ,505
+ ,505
+ ,512
+ ,512
+ ,526
+ ,554
+ ,505
+ ,505
+ ,512
+ ,512
+ ,554
+ ,554
+ ,505
+ ,505
+ ,512
+ ,584
+ ,554
+ ,554
+ ,505
+ ,505
+ ,584
+ ,584
+ ,554
+ ,554
+ ,505
+ ,569
+ ,584
+ ,584
+ ,554
+ ,554
+ ,569
+ ,569
+ ,584
+ ,584
+ ,554
+ ,540
+ ,569
+ ,569
+ ,584
+ ,584
+ ,540
+ ,540
+ ,569
+ ,569
+ ,584
+ ,522
+ ,540
+ ,540
+ ,569
+ ,569
+ ,522
+ ,522
+ ,540
+ ,540
+ ,569
+ ,526
+ ,522
+ ,522
+ ,540
+ ,540
+ ,526
+ ,526
+ ,522
+ ,522
+ ,540
+ ,527
+ ,526
+ ,526
+ ,522
+ ,522
+ ,527
+ ,527
+ ,526
+ ,526
+ ,522
+ ,516
+ ,527
+ ,527
+ ,526
+ ,526
+ ,516
+ ,516
+ ,527
+ ,527
+ ,526
+ ,503
+ ,516
+ ,516
+ ,527
+ ,527
+ ,503
+ ,503
+ ,516
+ ,516
+ ,527
+ ,489
+ ,503
+ ,503
+ ,516
+ ,516
+ ,489
+ ,489
+ ,503
+ ,503
+ ,516
+ ,479
+ ,489
+ ,489
+ ,503
+ ,503
+ ,479
+ ,479
+ ,489
+ ,489
+ ,503
+ ,475
+ ,479
+ ,479
+ ,489
+ ,489
+ ,475
+ ,475
+ ,479
+ ,479
+ ,489
+ ,524
+ ,475
+ ,475
+ ,479
+ ,479
+ ,524
+ ,524
+ ,475
+ ,475
+ ,479
+ ,552
+ ,524
+ ,524
+ ,475
+ ,475
+ ,552
+ ,552
+ ,524
+ ,524
+ ,475
+ ,532
+ ,552
+ ,552
+ ,524
+ ,524
+ ,532
+ ,532
+ ,552
+ ,552
+ ,524
+ ,511
+ ,532
+ ,532
+ ,552
+ ,552
+ ,511
+ ,511
+ ,532
+ ,532
+ ,552
+ ,492
+ ,511
+ ,511
+ ,532
+ ,532
+ ,492
+ ,492
+ ,511
+ ,511
+ ,532
+ ,492
+ ,492
+ ,492
+ ,511
+ ,511
+ ,492
+ ,492
+ ,492
+ ,492
+ ,511
+ ,493
+ ,492
+ ,492
+ ,492
+ ,492
+ ,493
+ ,493
+ ,492
+ ,492
+ ,492
+ ,481
+ ,493
+ ,493
+ ,492
+ ,492
+ ,481
+ ,481
+ ,493
+ ,493
+ ,492
+ ,462
+ ,481
+ ,481
+ ,493
+ ,493
+ ,462
+ ,462
+ ,481
+ ,481
+ ,493
+ ,457
+ ,462
+ ,462
+ ,481
+ ,481
+ ,457
+ ,457
+ ,462
+ ,462
+ ,481
+ ,442
+ ,457
+ ,457
+ ,462
+ ,462
+ ,442
+ ,442
+ ,457
+ ,457
+ ,462
+ ,439
+ ,442
+ ,442
+ ,457
+ ,457
+ ,439
+ ,439
+ ,442
+ ,442
+ ,457
+ ,488
+ ,439
+ ,439
+ ,442
+ ,442
+ ,488
+ ,488
+ ,439
+ ,439
+ ,442
+ ,521
+ ,488
+ ,488
+ ,439
+ ,439
+ ,521
+ ,521
+ ,488
+ ,488
+ ,439
+ ,501
+ ,521
+ ,521
+ ,488
+ ,488
+ ,501
+ ,501
+ ,521
+ ,521
+ ,488
+ ,485
+ ,501
+ ,501
+ ,521
+ ,521
+ ,485
+ ,485
+ ,501
+ ,501
+ ,521
+ ,464
+ ,485
+ ,485
+ ,501
+ ,501
+ ,464
+ ,464
+ ,485
+ ,485
+ ,501
+ ,460
+ ,464
+ ,464
+ ,485
+ ,485
+ ,460
+ ,460
+ ,464
+ ,464
+ ,485
+ ,467
+ ,460
+ ,460
+ ,464
+ ,464
+ ,467
+ ,467
+ ,460
+ ,460
+ ,464
+ ,460
+ ,467
+ ,467
+ ,460
+ ,460
+ ,460
+ ,460
+ ,467
+ ,467
+ ,460
+ ,448
+ ,460
+ ,460
+ ,467
+ ,467
+ ,448
+ ,448
+ ,460
+ ,460
+ ,467
+ ,443
+ ,448
+ ,448
+ ,460
+ ,460
+ ,443
+ ,443
+ ,448
+ ,448
+ ,460
+ ,436
+ ,443
+ ,443
+ ,448
+ ,448
+ ,436
+ ,436
+ ,443
+ ,443
+ ,448
+ ,431
+ ,436
+ ,436
+ ,443
+ ,443
+ ,431
+ ,431
+ ,436
+ ,436
+ ,443
+ ,484
+ ,431
+ ,431
+ ,436
+ ,436
+ ,484
+ ,484
+ ,431
+ ,431
+ ,436
+ ,510
+ ,484
+ ,484
+ ,431
+ ,431
+ ,510
+ ,510
+ ,484
+ ,484
+ ,431
+ ,513
+ ,510
+ ,510
+ ,484
+ ,484
+ ,513
+ ,513
+ ,510
+ ,510
+ ,484
+ ,503
+ ,513
+ ,513
+ ,510
+ ,510
+ ,503
+ ,503
+ ,513
+ ,513
+ ,510
+ ,471
+ ,503
+ ,503
+ ,513
+ ,513
+ ,471
+ ,471
+ ,503
+ ,503
+ ,513
+ ,471
+ ,471
+ ,471
+ ,503
+ ,503
+ ,471
+ ,471
+ ,471
+ ,471
+ ,503
+ ,476
+ ,471
+ ,471
+ ,471
+ ,471
+ ,476
+ ,476
+ ,471
+ ,471
+ ,471
+ ,475
+ ,476
+ ,476
+ ,471
+ ,471
+ ,475
+ ,475
+ ,476
+ ,476
+ ,471
+ ,470
+ ,475
+ ,475
+ ,476
+ ,476
+ ,470
+ ,470
+ ,475
+ ,475
+ ,476
+ ,461
+ ,470
+ ,470
+ ,475
+ ,475
+ ,461
+ ,461
+ ,470
+ ,470
+ ,475
+ ,455
+ ,461
+ ,461
+ ,470
+ ,470
+ ,455
+ ,455
+ ,461
+ ,461
+ ,470
+ ,456
+ ,455
+ ,455
+ ,461
+ ,461
+ ,456
+ ,456
+ ,455
+ ,455
+ ,461
+ ,517
+ ,456
+ ,456
+ ,455
+ ,455
+ ,517
+ ,517
+ ,456
+ ,456
+ ,455
+ ,525
+ ,517
+ ,517
+ ,456
+ ,456
+ ,525
+ ,525
+ ,517
+ ,517
+ ,456
+ ,523
+ ,525
+ ,525
+ ,517
+ ,517
+ ,523
+ ,523
+ ,525
+ ,525
+ ,517
+ ,519
+ ,523
+ ,523
+ ,525
+ ,525
+ ,519
+ ,519
+ ,523
+ ,523
+ ,525
+ ,509
+ ,519
+ ,519
+ ,523
+ ,523
+ ,509
+ ,509
+ ,519
+ ,519
+ ,523
+ ,512
+ ,509
+ ,509
+ ,519
+ ,519
+ ,512
+ ,512
+ ,509
+ ,509
+ ,519
+ ,519
+ ,512
+ ,512
+ ,509
+ ,509
+ ,519
+ ,519
+ ,512
+ ,512
+ ,509
+ ,517
+ ,519
+ ,519
+ ,512
+ ,512
+ ,517
+ ,517
+ ,519
+ ,519
+ ,512
+ ,510
+ ,517
+ ,517
+ ,519
+ ,519
+ ,510
+ ,510
+ ,517
+ ,517
+ ,519
+ ,509
+ ,510
+ ,510
+ ,517
+ ,517
+ ,509
+ ,509
+ ,510
+ ,510
+ ,517
+ ,501
+ ,509
+ ,509
+ ,510
+ ,510
+ ,501
+ ,501
+ ,509
+ ,509
+ ,510
+ ,507
+ ,501
+ ,501
+ ,509
+ ,509
+ ,507
+ ,507
+ ,501
+ ,501
+ ,509
+ ,569
+ ,507
+ ,507
+ ,501
+ ,501
+ ,569
+ ,569
+ ,507
+ ,507
+ ,501
+ ,580
+ ,569
+ ,569
+ ,507
+ ,507
+ ,580
+ ,580
+ ,569
+ ,569
+ ,507
+ ,578
+ ,580
+ ,580
+ ,569
+ ,569
+ ,578
+ ,578
+ ,580
+ ,580
+ ,569
+ ,565
+ ,578
+ ,578
+ ,580
+ ,580
+ ,565
+ ,565
+ ,578
+ ,578
+ ,580
+ ,547
+ ,565
+ ,565
+ ,578
+ ,578
+ ,547
+ ,547
+ ,565
+ ,565
+ ,578
+ ,555
+ ,547
+ ,547
+ ,565
+ ,565
+ ,555
+ ,555
+ ,547
+ ,547
+ ,565
+ ,562
+ ,555
+ ,555
+ ,547
+ ,547
+ ,561
+ ,562
+ ,555
+ ,555
+ ,547
+ ,555
+ ,561
+ ,562
+ ,555
+ ,555
+ ,544
+ ,555
+ ,561
+ ,562
+ ,555
+ ,537
+ ,544
+ ,555
+ ,561
+ ,562
+ ,543
+ ,537
+ ,544
+ ,555
+ ,561
+ ,594
+ ,543
+ ,537
+ ,544
+ ,555
+ ,611
+ ,594
+ ,543
+ ,537
+ ,544
+ ,613
+ ,611
+ ,594
+ ,543
+ ,537
+ ,611
+ ,613
+ ,611
+ ,594
+ ,543
+ ,594
+ ,611
+ ,613
+ ,611
+ ,594
+ ,595
+ ,594
+ ,611
+ ,613
+ ,611
+ ,591
+ ,595
+ ,594
+ ,611
+ ,613
+ ,589
+ ,591
+ ,595
+ ,594
+ ,611
+ ,584
+ ,589
+ ,591
+ ,595
+ ,594
+ ,573
+ ,584
+ ,589
+ ,591
+ ,595
+ ,567
+ ,573
+ ,584
+ ,589
+ ,591
+ ,569
+ ,567
+ ,573
+ ,584
+ ,589
+ ,621
+ ,569
+ ,567
+ ,573
+ ,584
+ ,629
+ ,621
+ ,569
+ ,567
+ ,573
+ ,628
+ ,629
+ ,621
+ ,569
+ ,567
+ ,612
+ ,628
+ ,629
+ ,621
+ ,569
+ ,595
+ ,612
+ ,628
+ ,629
+ ,621
+ ,597
+ ,595
+ ,612
+ ,628
+ ,629
+ ,593
+ ,597
+ ,595
+ ,612
+ ,628
+ ,590
+ ,593
+ ,597
+ ,595
+ ,612
+ ,580
+ ,590
+ ,593
+ ,597
+ ,595
+ ,574
+ ,580
+ ,590
+ ,593
+ ,597
+ ,573
+ ,574
+ ,580
+ ,590
+ ,593
+ ,573
+ ,573
+ ,574
+ ,580
+ ,590
+ ,620
+ ,573
+ ,573
+ ,574
+ ,580
+ ,626
+ ,620
+ ,573
+ ,573
+ ,574
+ ,620
+ ,626
+ ,620
+ ,573
+ ,573
+ ,588
+ ,620
+ ,626
+ ,620
+ ,573
+ ,566
+ ,588
+ ,620
+ ,626
+ ,620
+ ,557
+ ,566
+ ,588
+ ,620
+ ,626
+ ,561
+ ,557
+ ,566
+ ,588
+ ,620
+ ,549
+ ,561
+ ,557
+ ,566
+ ,588
+ ,532
+ ,549
+ ,561
+ ,557
+ ,566
+ ,526
+ ,532
+ ,549
+ ,561
+ ,557
+ ,511
+ ,526
+ ,532
+ ,549
+ ,561
+ ,499
+ ,511
+ ,526
+ ,532
+ ,549
+ ,555
+ ,499
+ ,511
+ ,526
+ ,532
+ ,565
+ ,555
+ ,499
+ ,511
+ ,526
+ ,542
+ ,565
+ ,555
+ ,499
+ ,511
+ ,527
+ ,542
+ ,565
+ ,555
+ ,499
+ ,510
+ ,527
+ ,542
+ ,565
+ ,555
+ ,514
+ ,510
+ ,527
+ ,542
+ ,565
+ ,517
+ ,514
+ ,510
+ ,527
+ ,542
+ ,508
+ ,517
+ ,514
+ ,510
+ ,527
+ ,493
+ ,508
+ ,517
+ ,514
+ ,510
+ ,490
+ ,493
+ ,508
+ ,517
+ ,514
+ ,469
+ ,490
+ ,493
+ ,508
+ ,517
+ ,478
+ ,469
+ ,490
+ ,493
+ ,508
+ ,528
+ ,478
+ ,469
+ ,490
+ ,493
+ ,534
+ ,528
+ ,478
+ ,469
+ ,490
+ ,518
+ ,534
+ ,528
+ ,478
+ ,469
+ ,506
+ ,518
+ ,534
+ ,528
+ ,478
+ ,502
+ ,506
+ ,518
+ ,534
+ ,528
+ ,516
+ ,502
+ ,506
+ ,518
+ ,534
+ ,528
+ ,516
+ ,502
+ ,506
+ ,518
+ ,533
+ ,528
+ ,516
+ ,502
+ ,506
+ ,536
+ ,533
+ ,528
+ ,516
+ ,502
+ ,537
+ ,536
+ ,533
+ ,528
+ ,516
+ ,524
+ ,537
+ ,536
+ ,533
+ ,528
+ ,536
+ ,524
+ ,537
+ ,536
+ ,533
+ ,587
+ ,536
+ ,524
+ ,537
+ ,536
+ ,597
+ ,587
+ ,536
+ ,524
+ ,537
+ ,581
+ ,597
+ ,587
+ ,536
+ ,524
+ ,564
+ ,581
+ ,597
+ ,587
+ ,536
+ ,558
+ ,564
+ ,581
+ ,597
+ ,587
+ ,575
+ ,558
+ ,564
+ ,581
+ ,597
+ ,580
+ ,575
+ ,558
+ ,564
+ ,581
+ ,575
+ ,580
+ ,575
+ ,558
+ ,564
+ ,563
+ ,575
+ ,580
+ ,575
+ ,558
+ ,552
+ ,563
+ ,575
+ ,580
+ ,575
+ ,537
+ ,552
+ ,563
+ ,575
+ ,580
+ ,545
+ ,537
+ ,552
+ ,563
+ ,575
+ ,601
+ ,545
+ ,537
+ ,552
+ ,563
+ ,604
+ ,601
+ ,545
+ ,537
+ ,552
+ ,586
+ ,604
+ ,601
+ ,545
+ ,537
+ ,564
+ ,586
+ ,604
+ ,601
+ ,545
+ ,549
+ ,564
+ ,586
+ ,604
+ ,601
+ ,551
+ ,549
+ ,564
+ ,586
+ ,604)
+ ,dim=c(5
+ ,492)
+ ,dimnames=list(c('werkloosheid'
+ ,'y1'
+ ,'y2'
+ ,'y3'
+ ,'y4')
+ ,1:492))
> y <- array(NA,dim=c(5,492),dimnames=list(c('werkloosheid','y1','y2','y3','y4'),1:492))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '1'
> 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
werkloosheid y1 y2 y3 y4 t
1 349 0 0 0 0 1
2 336 349 0 0 0 2
3 331 336 349 0 0 3
4 327 331 336 349 0 4
5 323 327 331 336 349 5
6 322 323 327 331 336 6
7 385 322 323 327 331 7
8 405 385 322 323 327 8
9 412 405 385 322 323 9
10 411 412 405 385 322 10
11 410 411 412 405 385 11
12 415 410 411 412 405 12
13 414 415 410 411 412 13
14 411 414 415 410 411 14
15 408 411 414 415 410 15
16 410 408 411 414 415 16
17 411 410 408 411 414 17
18 416 411 410 408 411 18
19 479 416 411 410 408 19
20 498 479 416 411 410 20
21 502 498 479 416 411 21
22 498 502 498 479 416 22
23 499 498 502 498 479 23
24 506 499 498 502 498 24
25 510 506 499 498 502 25
26 509 510 506 499 498 26
27 502 509 510 506 499 27
28 495 502 509 510 506 28
29 490 495 502 509 510 29
30 490 490 495 502 509 30
31 553 490 490 495 502 31
32 570 553 490 490 495 32
33 573 570 553 490 490 33
34 572 573 570 553 490 34
35 575 572 573 570 553 35
36 580 575 572 573 570 36
37 580 580 575 572 573 37
38 574 580 580 575 572 38
39 563 574 580 580 575 39
40 556 563 574 580 580 40
41 546 556 563 574 580 41
42 545 546 556 563 574 42
43 605 545 546 556 563 43
44 628 605 545 546 556 44
45 631 628 605 545 546 45
46 626 631 628 605 545 46
47 614 626 631 628 605 47
48 606 614 626 631 628 48
49 602 606 614 626 631 49
50 589 602 606 614 626 50
51 574 589 602 606 614 51
52 558 574 589 602 606 52
53 552 558 574 589 602 53
54 546 552 558 574 589 54
55 607 546 552 558 574 55
56 636 607 546 552 558 56
57 631 636 607 546 552 57
58 623 631 636 607 546 58
59 618 623 631 636 607 59
60 605 618 623 631 636 60
61 619 605 618 623 631 61
62 596 619 605 618 623 62
63 570 596 619 605 618 63
64 546 570 596 619 605 64
65 528 546 570 596 619 65
66 506 528 546 570 596 66
67 555 506 528 546 570 67
68 568 555 506 528 546 68
69 564 568 555 506 528 69
70 553 564 568 555 506 70
71 541 553 564 568 555 71
72 542 541 553 564 568 72
73 540 542 541 553 564 73
74 521 540 542 541 553 74
75 505 521 540 542 541 75
76 491 505 521 540 542 76
77 482 491 505 521 540 77
78 478 482 491 505 521 78
79 523 478 482 491 505 79
80 531 523 478 482 491 80
81 532 531 523 478 482 81
82 540 532 531 523 478 82
83 525 540 532 531 523 83
84 533 525 540 532 531 84
85 531 533 525 540 532 85
86 508 531 533 525 540 86
87 495 508 531 533 525 87
88 482 495 508 531 533 88
89 470 482 495 508 531 89
90 466 470 482 495 508 90
91 515 466 470 482 495 91
92 518 515 466 470 482 92
93 516 518 515 466 470 93
94 511 516 518 515 466 94
95 500 511 516 518 515 95
96 498 500 511 516 518 96
97 494 498 500 511 516 97
98 476 494 498 500 511 98
99 458 476 494 498 500 99
100 443 458 476 494 498 100
101 430 443 458 476 494 101
102 424 430 443 458 476 102
103 476 424 430 443 458 103
104 481 476 424 430 443 104
105 470 481 476 424 430 105
106 460 470 481 476 424 106
107 451 460 470 481 476 107
108 450 451 460 470 481 108
109 444 450 451 460 470 109
110 429 444 450 451 460 110
111 421 429 444 450 451 111
112 400 421 429 444 450 112
113 389 400 421 429 444 113
114 384 389 400 421 429 114
115 432 384 389 400 421 115
116 446 432 384 389 400 116
117 431 446 432 384 389 117
118 423 431 446 432 384 118
119 416 423 431 446 432 119
120 416 416 423 431 446 120
121 413 416 416 423 431 121
122 399 413 416 416 423 122
123 386 399 413 416 416 123
124 374 386 399 413 416 124
125 365 374 386 399 413 125
126 365 365 374 386 399 126
127 418 365 365 374 386 127
128 428 418 365 365 374 128
129 424 428 418 365 365 129
130 421 424 428 418 365 130
131 417 421 424 428 418 131
132 423 417 421 424 428 132
133 423 423 417 421 424 133
134 419 423 423 417 421 134
135 406 419 423 423 417 135
136 398 406 419 423 423 136
137 390 398 406 419 423 137
138 391 390 398 406 419 138
139 444 391 390 398 406 139
140 460 444 391 390 398 140
141 455 460 444 391 390 141
142 456 455 460 444 391 142
143 452 456 455 460 444 143
144 459 452 456 455 460 144
145 461 459 452 456 455 145
146 451 461 459 452 456 146
147 443 451 461 459 452 147
148 439 443 451 461 459 148
149 430 439 443 451 461 149
150 436 430 439 443 451 150
151 488 436 430 439 443 151
152 506 488 436 430 439 152
153 502 506 488 436 430 153
154 501 502 506 488 436 154
155 501 501 502 506 488 155
156 515 501 501 502 506 156
157 521 515 501 501 502 157
158 520 521 515 501 501 158
159 512 520 521 515 501 159
160 509 512 520 521 515 160
161 505 509 512 520 521 161
162 511 505 509 512 520 162
163 570 511 505 509 512 163
164 592 570 511 505 509 164
165 594 592 570 511 505 165
166 586 594 592 570 511 166
167 586 586 594 592 570 167
168 592 586 586 594 592 168
169 594 592 586 586 594 169
170 594 594 592 586 586 170
171 586 594 594 592 586 171
172 586 586 594 594 592 172
173 572 586 586 594 594 173
174 572 572 586 586 594 174
175 563 572 572 586 586 175
176 563 563 572 572 586 176
177 555 563 563 572 572 177
178 555 555 563 563 572 178
179 554 555 555 563 563 179
180 554 554 555 555 563 180
181 601 554 554 555 555 181
182 601 601 554 554 555 182
183 622 601 601 554 554 183
184 622 622 601 601 554 184
185 617 622 622 601 601 185
186 617 617 622 622 601 186
187 606 617 617 622 622 187
188 606 606 617 617 622 188
189 595 606 606 617 617 189
190 595 595 606 606 617 190
191 599 595 595 606 606 191
192 599 599 595 595 606 192
193 600 599 599 595 595 193
194 600 600 599 599 595 194
195 592 600 600 599 599 195
196 592 592 600 600 599 196
197 575 592 592 600 600 197
198 575 575 592 592 600 198
199 567 575 575 592 592 199
200 567 567 575 575 592 200
201 555 567 567 575 575 201
202 555 555 567 567 575 202
203 555 555 555 567 567 203
204 555 555 555 555 567 204
205 608 555 555 555 555 205
206 608 608 555 555 555 206
207 631 608 608 555 555 207
208 631 631 608 608 555 208
209 629 631 631 608 608 209
210 629 629 631 631 608 210
211 624 629 629 631 631 211
212 624 624 629 629 631 212
213 610 624 624 629 629 213
214 610 610 624 624 629 214
215 616 610 610 624 624 215
216 616 616 610 610 624 216
217 621 616 616 610 610 217
218 621 621 616 616 610 218
219 604 621 621 616 616 219
220 604 604 621 621 616 220
221 584 604 604 621 621 221
222 584 584 604 604 621 222
223 574 584 584 604 604 223
224 574 574 584 584 604 224
225 555 574 574 584 584 225
226 555 555 574 574 584 226
227 545 555 555 574 574 227
228 545 545 555 555 574 228
229 599 545 545 555 555 229
230 599 599 545 545 555 230
231 620 599 599 545 545 231
232 620 620 599 599 545 232
233 608 620 620 599 599 233
234 608 608 620 620 599 234
235 590 608 608 620 620 235
236 590 590 608 608 620 236
237 579 590 590 608 608 237
238 579 579 590 590 608 238
239 580 579 579 590 590 239
240 580 580 579 579 590 240
241 579 580 580 579 579 241
242 579 579 580 580 579 242
243 572 579 579 580 580 243
244 572 572 579 579 580 244
245 560 572 572 579 579 245
246 560 560 572 572 579 246
247 551 560 560 572 572 247
248 551 551 560 560 572 248
249 537 551 551 560 560 249
250 537 537 551 551 560 250
251 541 537 537 551 551 251
252 541 541 537 537 551 252
253 588 541 541 537 537 253
254 588 588 541 541 537 254
255 607 588 588 541 541 255
256 607 607 588 588 541 256
257 599 607 607 588 588 257
258 599 599 607 607 588 258
259 578 599 599 607 607 259
260 578 578 599 599 607 260
261 563 578 578 599 599 261
262 563 563 578 578 599 262
263 566 563 563 578 578 263
264 566 566 563 563 578 264
265 561 566 566 563 563 265
266 561 561 566 566 563 266
267 554 561 561 566 566 267
268 554 554 561 561 566 268
269 540 554 554 561 561 269
270 540 540 554 554 561 270
271 526 540 540 554 554 271
272 526 526 540 540 554 272
273 512 526 526 540 540 273
274 512 512 526 526 540 274
275 505 512 512 526 526 275
276 505 505 512 512 526 276
277 554 505 505 512 512 277
278 554 554 505 505 512 278
279 584 554 554 505 505 279
280 584 584 554 554 505 280
281 569 584 584 554 554 281
282 569 569 584 584 554 282
283 540 569 569 584 584 283
284 540 540 569 569 584 284
285 522 540 540 569 569 285
286 522 522 540 540 569 286
287 526 522 522 540 540 287
288 526 526 522 522 540 288
289 527 526 526 522 522 289
290 527 527 526 526 522 290
291 516 527 527 526 526 291
292 516 516 527 527 526 292
293 503 516 516 527 527 293
294 503 503 516 516 527 294
295 489 503 503 516 516 295
296 489 489 503 503 516 296
297 479 489 489 503 503 297
298 479 479 489 489 503 298
299 475 479 479 489 489 299
300 475 475 479 479 489 300
301 524 475 475 479 479 301
302 524 524 475 475 479 302
303 552 524 524 475 475 303
304 552 552 524 524 475 304
305 532 552 552 524 524 305
306 532 532 552 552 524 306
307 511 532 532 552 552 307
308 511 511 532 532 552 308
309 492 511 511 532 532 309
310 492 492 511 511 532 310
311 492 492 492 511 511 311
312 492 492 492 492 511 312
313 493 492 492 492 492 313
314 493 493 492 492 492 314
315 481 493 493 492 492 315
316 481 481 493 493 492 316
317 462 481 481 493 493 317
318 462 462 481 481 493 318
319 457 462 462 481 481 319
320 457 457 462 462 481 320
321 442 457 457 462 462 321
322 442 442 457 457 462 322
323 439 442 442 457 457 323
324 439 439 442 442 457 324
325 488 439 439 442 442 325
326 488 488 439 439 442 326
327 521 488 488 439 439 327
328 521 521 488 488 439 328
329 501 521 521 488 488 329
330 501 501 521 521 488 330
331 485 501 501 521 521 331
332 485 485 501 501 521 332
333 464 485 485 501 501 333
334 464 464 485 485 501 334
335 460 464 464 485 485 335
336 460 460 464 464 485 336
337 467 460 460 464 464 337
338 467 467 460 460 464 338
339 460 467 467 460 460 339
340 460 460 467 467 460 340
341 448 460 460 467 467 341
342 448 448 460 460 467 342
343 443 448 448 460 460 343
344 443 443 448 448 460 344
345 436 443 443 448 448 345
346 436 436 443 443 448 346
347 431 436 436 443 443 347
348 431 431 436 436 443 348
349 484 431 431 436 436 349
350 484 484 431 431 436 350
351 510 484 484 431 431 351
352 510 510 484 484 431 352
353 513 510 510 484 484 353
354 513 513 510 510 484 354
355 503 513 513 510 510 355
356 503 503 513 513 510 356
357 471 503 503 513 513 357
358 471 471 503 503 513 358
359 471 471 471 503 503 359
360 471 471 471 471 503 360
361 476 471 471 471 471 361
362 476 476 471 471 471 362
363 475 476 476 471 471 363
364 475 475 476 476 471 364
365 470 475 475 476 476 365
366 470 470 475 475 476 366
367 461 470 470 475 475 367
368 461 461 470 470 475 368
369 455 461 461 470 470 369
370 455 455 461 461 470 370
371 456 455 455 461 461 371
372 456 456 455 455 461 372
373 517 456 456 455 455 373
374 517 517 456 456 455 374
375 525 517 517 456 456 375
376 525 525 517 517 456 376
377 523 525 525 517 517 377
378 523 523 525 525 517 378
379 519 523 523 525 525 379
380 519 519 523 523 525 380
381 509 519 519 523 523 381
382 509 509 519 519 523 382
383 512 509 509 519 519 383
384 512 512 509 509 519 384
385 519 512 512 509 509 385
386 519 519 512 512 509 386
387 517 519 519 512 512 387
388 517 517 519 519 512 388
389 510 517 517 519 519 389
390 510 510 517 517 519 390
391 509 510 510 517 517 391
392 509 509 510 510 517 392
393 501 509 509 510 510 393
394 501 501 509 509 510 394
395 507 501 501 509 509 395
396 507 507 501 501 509 396
397 569 507 507 501 501 397
398 569 569 507 507 501 398
399 580 569 569 507 507 399
400 580 580 569 569 507 400
401 578 580 580 569 569 401
402 578 578 580 580 569 402
403 565 578 578 580 580 403
404 565 565 578 578 580 404
405 547 565 565 578 578 405
406 547 547 565 565 578 406
407 555 547 547 565 565 407
408 555 555 547 547 565 408
409 562 555 555 547 547 409
410 561 562 555 555 547 410
411 555 561 562 555 555 411
412 544 555 561 562 555 412
413 537 544 555 561 562 413
414 543 537 544 555 561 414
415 594 543 537 544 555 415
416 611 594 543 537 544 416
417 613 611 594 543 537 417
418 611 613 611 594 543 418
419 594 611 613 611 594 419
420 595 594 611 613 611 420
421 591 595 594 611 613 421
422 589 591 595 594 611 422
423 584 589 591 595 594 423
424 573 584 589 591 595 424
425 567 573 584 589 591 425
426 569 567 573 584 589 426
427 621 569 567 573 584 427
428 629 621 569 567 573 428
429 628 629 621 569 567 429
430 612 628 629 621 569 430
431 595 612 628 629 621 431
432 597 595 612 628 629 432
433 593 597 595 612 628 433
434 590 593 597 595 612 434
435 580 590 593 597 595 435
436 574 580 590 593 597 436
437 573 574 580 590 593 437
438 573 573 574 580 590 438
439 620 573 573 574 580 439
440 626 620 573 573 574 440
441 620 626 620 573 573 441
442 588 620 626 620 573 442
443 566 588 620 626 620 443
444 557 566 588 620 626 444
445 561 557 566 588 620 445
446 549 561 557 566 588 446
447 532 549 561 557 566 447
448 526 532 549 561 557 448
449 511 526 532 549 561 449
450 499 511 526 532 549 450
451 555 499 511 526 532 451
452 565 555 499 511 526 452
453 542 565 555 499 511 453
454 527 542 565 555 499 454
455 510 527 542 565 555 455
456 514 510 527 542 565 456
457 517 514 510 527 542 457
458 508 517 514 510 527 458
459 493 508 517 514 510 459
460 490 493 508 517 514 460
461 469 490 493 508 517 461
462 478 469 490 493 508 462
463 528 478 469 490 493 463
464 534 528 478 469 490 464
465 518 534 528 478 469 465
466 506 518 534 528 478 466
467 502 506 518 534 528 467
468 516 502 506 518 534 468
469 528 516 502 506 518 469
470 533 528 516 502 506 470
471 536 533 528 516 502 471
472 537 536 533 528 516 472
473 524 537 536 533 528 473
474 536 524 537 536 533 474
475 587 536 524 537 536 475
476 597 587 536 524 537 476
477 581 597 587 536 524 477
478 564 581 597 587 536 478
479 558 564 581 597 587 479
480 575 558 564 581 597 480
481 580 575 558 564 581 481
482 575 580 575 558 564 482
483 563 575 580 575 558 483
484 552 563 575 580 575 484
485 537 552 563 575 580 485
486 545 537 552 563 575 486
487 601 545 537 552 563 487
488 604 601 545 537 552 488
489 586 604 601 545 537 489
490 564 586 604 601 545 490
491 549 564 586 604 601 491
492 551 549 564 586 604 492
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) y1 y2 y3 y4 t
60.0822007 0.9643136 -0.0623268 -0.0941034 0.0790286 0.0008752
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-60.629 -10.389 -2.535 4.978 288.917
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 60.0822007 7.3936888 8.126 3.7e-15 ***
y1 0.9643136 0.0448793 21.487 < 2e-16 ***
y2 -0.0623268 0.0625509 -0.996 0.3195
y3 -0.0941034 0.0625453 -1.505 0.1331
y4 0.0790286 0.0426348 1.854 0.0644 .
t 0.0008752 0.0070439 0.124 0.9012
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 21.38 on 486 degrees of freedom
Multiple R-squared: 0.8921, Adjusted R-squared: 0.8909
F-statistic: 803.3 on 5 and 486 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,] 0.6442869 7.114261e-01 3.557131e-01
[2,] 0.6731429 6.537142e-01 3.268571e-01
[3,] 0.7730450 4.539101e-01 2.269550e-01
[4,] 0.8594473 2.811054e-01 1.405527e-01
[5,] 0.9400684 1.198632e-01 5.993159e-02
[6,] 0.9802034 3.959321e-02 1.979660e-02
[7,] 0.9933503 1.329937e-02 6.649686e-03
[8,] 0.9963798 7.240411e-03 3.620206e-03
[9,] 0.9973763 5.247460e-03 2.623730e-03
[10,] 0.9971396 5.720793e-03 2.860397e-03
[11,] 0.9993031 1.393760e-03 6.968800e-04
[12,] 0.9997757 4.485913e-04 2.242956e-04
[13,] 0.9998289 3.421298e-04 1.710649e-04
[14,] 0.9997436 5.128634e-04 2.564317e-04
[15,] 0.9995589 8.822373e-04 4.411186e-04
[16,] 0.9992567 1.486629e-03 7.433143e-04
[17,] 0.9987962 2.407605e-03 1.203802e-03
[18,] 0.9983640 3.272061e-03 1.636031e-03
[19,] 0.9986478 2.704317e-03 1.352158e-03
[20,] 0.9994018 1.196437e-03 5.982185e-04
[21,] 0.9998408 3.184896e-04 1.592448e-04
[22,] 0.9999575 8.502786e-05 4.251393e-05
[23,] 0.9999740 5.205098e-05 2.602549e-05
[24,] 0.9999722 5.568406e-05 2.784203e-05
[25,] 0.9999627 7.468818e-05 3.734409e-05
[26,] 0.9999427 1.145629e-04 5.728143e-05
[27,] 0.9999086 1.828512e-04 9.142562e-05
[28,] 0.9998554 2.892813e-04 1.446407e-04
[29,] 0.9997855 4.290821e-04 2.145411e-04
[30,] 0.9997683 4.634193e-04 2.317096e-04
[31,] 0.9998830 2.339035e-04 1.169517e-04
[32,] 0.9999689 6.227338e-05 3.113669e-05
[33,] 0.9999971 5.893663e-06 2.946832e-06
[34,] 0.9999997 6.137501e-07 3.068751e-07
[35,] 0.9999998 3.783187e-07 1.891593e-07
[36,] 0.9999997 5.122081e-07 2.561040e-07
[37,] 0.9999996 8.205417e-07 4.102708e-07
[38,] 0.9999993 1.364524e-06 6.822618e-07
[39,] 0.9999992 1.667139e-06 8.335693e-07
[40,] 0.9999994 1.254976e-06 6.274880e-07
[41,] 0.9999997 6.474981e-07 3.237491e-07
[42,] 1.0000000 9.822434e-08 4.911217e-08
[43,] 1.0000000 3.962436e-09 1.981218e-09
[44,] 1.0000000 4.180720e-11 2.090360e-11
[45,] 1.0000000 6.659044e-13 3.329522e-13
[46,] 1.0000000 1.320397e-14 6.601985e-15
[47,] 1.0000000 4.038580e-15 2.019290e-15
[48,] 1.0000000 6.201508e-15 3.100754e-15
[49,] 1.0000000 1.054307e-14 5.271535e-15
[50,] 1.0000000 1.596455e-14 7.982276e-15
[51,] 1.0000000 1.835257e-14 9.176284e-15
[52,] 1.0000000 8.556829e-15 4.278415e-15
[53,] 1.0000000 9.306008e-15 4.653004e-15
[54,] 1.0000000 1.483415e-15 7.417075e-16
[55,] 1.0000000 4.259344e-17 2.129672e-17
[56,] 1.0000000 4.612755e-19 2.306378e-19
[57,] 1.0000000 2.994232e-21 1.497116e-21
[58,] 1.0000000 9.095889e-24 4.547945e-24
[59,] 1.0000000 2.601539e-24 1.300770e-24
[60,] 1.0000000 3.194582e-24 1.597291e-24
[61,] 1.0000000 2.607218e-24 1.303609e-24
[62,] 1.0000000 1.943840e-24 9.719201e-25
[63,] 1.0000000 8.167071e-25 4.083536e-25
[64,] 1.0000000 6.282972e-25 3.141486e-25
[65,] 1.0000000 4.760925e-25 2.380463e-25
[66,] 1.0000000 9.413622e-26 4.706811e-26
[67,] 1.0000000 2.279401e-26 1.139701e-26
[68,] 1.0000000 5.286594e-27 2.643297e-27
[69,] 1.0000000 1.672162e-27 8.360809e-28
[70,] 1.0000000 9.694797e-28 4.847399e-28
[71,] 1.0000000 4.124221e-28 2.062111e-28
[72,] 1.0000000 8.945184e-28 4.472592e-28
[73,] 1.0000000 1.882857e-27 9.414286e-28
[74,] 1.0000000 3.564088e-27 1.782044e-27
[75,] 1.0000000 4.043804e-27 2.021902e-27
[76,] 1.0000000 7.852253e-27 3.926126e-27
[77,] 1.0000000 1.483397e-26 7.416983e-27
[78,] 1.0000000 6.686301e-27 3.343151e-27
[79,] 1.0000000 6.552067e-27 3.276034e-27
[80,] 1.0000000 4.869514e-27 2.434757e-27
[81,] 1.0000000 3.016346e-27 1.508173e-27
[82,] 1.0000000 3.408058e-27 1.704029e-27
[83,] 1.0000000 6.458540e-28 3.229270e-28
[84,] 1.0000000 1.376899e-27 6.884497e-28
[85,] 1.0000000 2.911761e-27 1.455881e-27
[86,] 1.0000000 6.202383e-27 3.101192e-27
[87,] 1.0000000 9.180949e-27 4.590474e-27
[88,] 1.0000000 1.705764e-26 8.528819e-27
[89,] 1.0000000 2.951733e-26 1.475866e-26
[90,] 1.0000000 2.206533e-26 1.103266e-26
[91,] 1.0000000 1.573948e-26 7.869738e-27
[92,] 1.0000000 1.138027e-26 5.690136e-27
[93,] 1.0000000 7.507751e-27 3.753876e-27
[94,] 1.0000000 7.870357e-27 3.935178e-27
[95,] 1.0000000 1.036722e-27 5.183611e-28
[96,] 1.0000000 2.175810e-27 1.087905e-27
[97,] 1.0000000 3.116648e-27 1.558324e-27
[98,] 1.0000000 6.125806e-27 3.062903e-27
[99,] 1.0000000 9.297748e-27 4.648874e-27
[100,] 1.0000000 1.678807e-26 8.394037e-27
[101,] 1.0000000 2.569590e-26 1.284795e-26
[102,] 1.0000000 2.327395e-26 1.163698e-26
[103,] 1.0000000 3.217494e-26 1.608747e-26
[104,] 1.0000000 1.510728e-26 7.553642e-27
[105,] 1.0000000 1.335358e-26 6.676790e-27
[106,] 1.0000000 1.694794e-26 8.473970e-27
[107,] 1.0000000 4.264699e-27 2.132349e-27
[108,] 1.0000000 8.607229e-27 4.303615e-27
[109,] 1.0000000 8.492303e-27 4.246151e-27
[110,] 1.0000000 1.717474e-26 8.587372e-27
[111,] 1.0000000 2.993892e-26 1.496946e-26
[112,] 1.0000000 5.582797e-26 2.791398e-26
[113,] 1.0000000 9.773582e-26 4.886791e-26
[114,] 1.0000000 9.722873e-26 4.861437e-26
[115,] 1.0000000 1.055283e-25 5.276417e-26
[116,] 1.0000000 1.084804e-25 5.424021e-26
[117,] 1.0000000 1.137150e-25 5.685748e-26
[118,] 1.0000000 1.847574e-25 9.237871e-26
[119,] 1.0000000 2.068855e-26 1.034428e-26
[120,] 1.0000000 3.869782e-26 1.934891e-26
[121,] 1.0000000 6.013664e-26 3.006832e-26
[122,] 1.0000000 1.205327e-25 6.026634e-26
[123,] 1.0000000 2.285256e-25 1.142628e-25
[124,] 1.0000000 4.607886e-25 2.303943e-25
[125,] 1.0000000 8.818825e-25 4.409412e-25
[126,] 1.0000000 1.529579e-24 7.647894e-25
[127,] 1.0000000 1.950619e-24 9.753094e-25
[128,] 1.0000000 2.965507e-24 1.482753e-24
[129,] 1.0000000 4.172509e-24 2.086255e-24
[130,] 1.0000000 7.487718e-24 3.743859e-24
[131,] 1.0000000 5.141538e-25 2.570769e-25
[132,] 1.0000000 7.116056e-25 3.558028e-25
[133,] 1.0000000 9.224613e-25 4.612306e-25
[134,] 1.0000000 1.484920e-24 7.424599e-25
[135,] 1.0000000 2.736693e-24 1.368346e-24
[136,] 1.0000000 4.813755e-24 2.406877e-24
[137,] 1.0000000 8.665832e-24 4.332916e-24
[138,] 1.0000000 1.239248e-23 6.196238e-24
[139,] 1.0000000 2.063081e-23 1.031541e-23
[140,] 1.0000000 3.712042e-23 1.856021e-23
[141,] 1.0000000 5.263872e-23 2.631936e-23
[142,] 1.0000000 9.762208e-23 4.881104e-23
[143,] 1.0000000 3.868930e-24 1.934465e-24
[144,] 1.0000000 3.779760e-24 1.889880e-24
[145,] 1.0000000 4.491027e-24 2.245513e-24
[146,] 1.0000000 6.299844e-24 3.149922e-24
[147,] 1.0000000 1.026950e-23 5.134750e-24
[148,] 1.0000000 1.106473e-23 5.532363e-24
[149,] 1.0000000 1.550279e-23 7.751393e-24
[150,] 1.0000000 2.362996e-23 1.181498e-23
[151,] 1.0000000 3.832442e-23 1.916221e-23
[152,] 1.0000000 6.541163e-23 3.270581e-23
[153,] 1.0000000 1.137549e-22 5.687747e-23
[154,] 1.0000000 1.765326e-22 8.826628e-23
[155,] 1.0000000 5.555158e-25 2.777579e-25
[156,] 1.0000000 2.622436e-25 1.311218e-25
[157,] 1.0000000 2.603097e-25 1.301549e-25
[158,] 1.0000000 3.516323e-25 1.758162e-25
[159,] 1.0000000 4.559123e-25 2.279561e-25
[160,] 1.0000000 5.158197e-25 2.579099e-25
[161,] 1.0000000 7.277067e-25 3.638534e-25
[162,] 1.0000000 1.062816e-24 5.314078e-25
[163,] 1.0000000 1.792801e-24 8.964005e-25
[164,] 1.0000000 2.741894e-24 1.370947e-24
[165,] 1.0000000 4.595708e-24 2.297854e-24
[166,] 1.0000000 7.627418e-24 3.813709e-24
[167,] 1.0000000 1.381378e-23 6.906890e-24
[168,] 1.0000000 2.394589e-23 1.197295e-23
[169,] 1.0000000 4.310318e-23 2.155159e-23
[170,] 1.0000000 7.562643e-23 3.781322e-23
[171,] 1.0000000 1.341263e-22 6.706316e-23
[172,] 1.0000000 2.356970e-22 1.178485e-22
[173,] 1.0000000 5.883591e-24 2.941795e-24
[174,] 1.0000000 9.580491e-24 4.790245e-24
[175,] 1.0000000 5.171512e-24 2.585756e-24
[176,] 1.0000000 6.966835e-24 3.483417e-24
[177,] 1.0000000 1.173021e-23 5.865105e-24
[178,] 1.0000000 1.708672e-23 8.543362e-24
[179,] 1.0000000 3.079821e-23 1.539910e-23
[180,] 1.0000000 5.069639e-23 2.534819e-23
[181,] 1.0000000 9.225969e-23 4.612985e-23
[182,] 1.0000000 1.589532e-22 7.947658e-23
[183,] 1.0000000 2.408342e-22 1.204171e-22
[184,] 1.0000000 4.194358e-22 2.097179e-22
[185,] 1.0000000 6.990845e-22 3.495423e-22
[186,] 1.0000000 1.179374e-21 5.896871e-22
[187,] 1.0000000 2.130377e-21 1.065189e-21
[188,] 1.0000000 3.628286e-21 1.814143e-21
[189,] 1.0000000 5.539468e-21 2.769734e-21
[190,] 1.0000000 9.759176e-21 4.879588e-21
[191,] 1.0000000 1.774032e-20 8.870160e-21
[192,] 1.0000000 3.193585e-20 1.596792e-20
[193,] 1.0000000 5.264404e-20 2.632202e-20
[194,] 1.0000000 9.442659e-20 4.721330e-20
[195,] 1.0000000 1.684850e-19 8.424249e-20
[196,] 1.0000000 3.014770e-19 1.507385e-19
[197,] 1.0000000 4.878838e-21 2.439419e-21
[198,] 1.0000000 8.523613e-21 4.261806e-21
[199,] 1.0000000 5.046824e-21 2.523412e-21
[200,] 1.0000000 7.240772e-21 3.620386e-21
[201,] 1.0000000 1.240076e-20 6.200382e-21
[202,] 1.0000000 1.851544e-20 9.257719e-21
[203,] 1.0000000 3.299056e-20 1.649528e-20
[204,] 1.0000000 5.387277e-20 2.693639e-20
[205,] 1.0000000 9.272803e-20 4.636402e-20
[206,] 1.0000000 1.536329e-19 7.681644e-20
[207,] 1.0000000 2.086997e-19 1.043499e-19
[208,] 1.0000000 3.603466e-19 1.801733e-19
[209,] 1.0000000 5.206210e-19 2.603105e-19
[210,] 1.0000000 8.469443e-19 4.234721e-19
[211,] 1.0000000 1.317393e-18 6.586966e-19
[212,] 1.0000000 2.094402e-18 1.047201e-18
[213,] 1.0000000 2.962016e-18 1.481008e-18
[214,] 1.0000000 5.088466e-18 2.544233e-18
[215,] 1.0000000 8.792440e-18 4.396220e-18
[216,] 1.0000000 1.537094e-17 7.685472e-18
[217,] 1.0000000 1.979942e-17 9.899708e-18
[218,] 1.0000000 3.422159e-17 1.711079e-17
[219,] 1.0000000 5.592711e-17 2.796355e-17
[220,] 1.0000000 9.730656e-17 4.865328e-17
[221,] 1.0000000 1.665185e-18 8.325927e-19
[222,] 1.0000000 2.916531e-18 1.458265e-18
[223,] 1.0000000 2.227794e-18 1.113897e-18
[224,] 1.0000000 3.190209e-18 1.595104e-18
[225,] 1.0000000 5.423151e-18 2.711576e-18
[226,] 1.0000000 7.952319e-18 3.976159e-18
[227,] 1.0000000 1.205771e-17 6.028854e-18
[228,] 1.0000000 2.007647e-17 1.003823e-17
[229,] 1.0000000 3.425715e-17 1.712858e-17
[230,] 1.0000000 5.879677e-17 2.939839e-17
[231,] 1.0000000 9.641203e-17 4.820601e-17
[232,] 1.0000000 1.649694e-16 8.248470e-17
[233,] 1.0000000 2.795021e-16 1.397511e-16
[234,] 1.0000000 4.637208e-16 2.318604e-16
[235,] 1.0000000 7.904856e-16 3.952428e-16
[236,] 1.0000000 1.300131e-15 6.500657e-16
[237,] 1.0000000 2.051433e-15 1.025716e-15
[238,] 1.0000000 3.393740e-15 1.696870e-15
[239,] 1.0000000 5.513391e-15 2.756696e-15
[240,] 1.0000000 9.191511e-15 4.595755e-15
[241,] 1.0000000 1.301983e-14 6.509917e-15
[242,] 1.0000000 2.151255e-14 1.075628e-14
[243,] 1.0000000 3.373992e-14 1.686996e-14
[244,] 1.0000000 5.590824e-14 2.795412e-14
[245,] 1.0000000 3.651767e-15 1.825883e-15
[246,] 1.0000000 6.186531e-15 3.093266e-15
[247,] 1.0000000 5.065805e-15 2.532902e-15
[248,] 1.0000000 6.978253e-15 3.489126e-15
[249,] 1.0000000 1.177193e-14 5.885964e-15
[250,] 1.0000000 1.614680e-14 8.073401e-15
[251,] 1.0000000 2.183328e-14 1.091664e-14
[252,] 1.0000000 3.295316e-14 1.647658e-14
[253,] 1.0000000 5.018000e-14 2.509000e-14
[254,] 1.0000000 8.055599e-14 4.027799e-14
[255,] 1.0000000 1.181178e-13 5.905890e-14
[256,] 1.0000000 1.912022e-13 9.560110e-14
[257,] 1.0000000 3.112481e-13 1.556240e-13
[258,] 1.0000000 4.769987e-13 2.384994e-13
[259,] 1.0000000 7.656940e-13 3.828470e-13
[260,] 1.0000000 1.186588e-12 5.932942e-13
[261,] 1.0000000 1.708173e-12 8.540867e-13
[262,] 1.0000000 2.638607e-12 1.319304e-12
[263,] 1.0000000 3.676567e-12 1.838283e-12
[264,] 1.0000000 5.755703e-12 2.877851e-12
[265,] 1.0000000 7.678394e-12 3.839197e-12
[266,] 1.0000000 1.200167e-11 6.000837e-12
[267,] 1.0000000 1.793705e-11 8.968525e-12
[268,] 1.0000000 2.793485e-11 1.396743e-11
[269,] 1.0000000 1.998952e-12 9.994758e-13
[270,] 1.0000000 3.209516e-12 1.604758e-12
[271,] 1.0000000 1.216999e-12 6.084994e-13
[272,] 1.0000000 1.683004e-12 8.415022e-13
[273,] 1.0000000 2.541452e-12 1.270726e-12
[274,] 1.0000000 3.049810e-12 1.524905e-12
[275,] 1.0000000 2.782045e-12 1.391023e-12
[276,] 1.0000000 4.071915e-12 2.035958e-12
[277,] 1.0000000 5.278314e-12 2.639157e-12
[278,] 1.0000000 8.291346e-12 4.145673e-12
[279,] 1.0000000 1.207118e-11 6.035591e-12
[280,] 1.0000000 1.907479e-11 9.537397e-12
[281,] 1.0000000 2.912793e-11 1.456396e-11
[282,] 1.0000000 4.437590e-11 2.218795e-11
[283,] 1.0000000 6.440174e-11 3.220087e-11
[284,] 1.0000000 9.642838e-11 4.821419e-11
[285,] 1.0000000 1.310944e-10 6.554718e-11
[286,] 1.0000000 1.997699e-10 9.988494e-11
[287,] 1.0000000 2.534214e-10 1.267107e-10
[288,] 1.0000000 3.859910e-10 1.929955e-10
[289,] 1.0000000 5.216933e-10 2.608467e-10
[290,] 1.0000000 7.901429e-10 3.950714e-10
[291,] 1.0000000 1.159471e-09 5.797357e-10
[292,] 1.0000000 1.741687e-09 8.708434e-10
[293,] 1.0000000 1.816399e-10 9.081994e-11
[294,] 1.0000000 2.770585e-10 1.385293e-10
[295,] 1.0000000 1.329155e-10 6.645773e-11
[296,] 1.0000000 1.882594e-10 9.412968e-11
[297,] 1.0000000 2.346365e-10 1.173183e-10
[298,] 1.0000000 2.923301e-10 1.461651e-10
[299,] 1.0000000 3.347065e-10 1.673533e-10
[300,] 1.0000000 5.042312e-10 2.521156e-10
[301,] 1.0000000 5.640688e-10 2.820344e-10
[302,] 1.0000000 8.617174e-10 4.308587e-10
[303,] 1.0000000 1.318988e-09 6.594938e-10
[304,] 1.0000000 2.013262e-09 1.006631e-09
[305,] 1.0000000 3.053097e-09 1.526548e-09
[306,] 1.0000000 4.626967e-09 2.313484e-09
[307,] 1.0000000 6.005260e-09 3.002630e-09
[308,] 1.0000000 8.965272e-09 4.482636e-09
[309,] 1.0000000 8.617456e-09 4.308728e-09
[310,] 1.0000000 1.286459e-08 6.432294e-09
[311,] 1.0000000 1.798536e-08 8.992680e-09
[312,] 1.0000000 2.612676e-08 1.306338e-08
[313,] 1.0000000 2.634612e-08 1.317306e-08
[314,] 1.0000000 3.864378e-08 1.932189e-08
[315,] 1.0000000 5.313586e-08 2.656793e-08
[316,] 1.0000000 7.431179e-08 3.715590e-08
[317,] 1.0000000 1.361263e-08 6.806314e-09
[318,] 1.0000000 1.848697e-08 9.243484e-09
[319,] 1.0000000 7.534663e-09 3.767331e-09
[320,] 1.0000000 1.139946e-08 5.699731e-09
[321,] 1.0000000 1.266234e-08 6.331170e-09
[322,] 1.0000000 1.699700e-08 8.498501e-09
[323,] 1.0000000 1.939933e-08 9.699664e-09
[324,] 1.0000000 2.938157e-08 1.469079e-08
[325,] 1.0000000 2.457442e-08 1.228721e-08
[326,] 1.0000000 3.702543e-08 1.851271e-08
[327,] 1.0000000 5.198220e-08 2.599110e-08
[328,] 1.0000000 7.511186e-08 3.755593e-08
[329,] 0.9999999 1.116187e-07 5.580935e-08
[330,] 0.9999999 1.620981e-07 8.104903e-08
[331,] 0.9999999 2.153381e-07 1.076690e-07
[332,] 0.9999998 3.176933e-07 1.588467e-07
[333,] 0.9999998 3.483838e-07 1.741919e-07
[334,] 0.9999997 5.031939e-07 2.515969e-07
[335,] 0.9999997 6.451666e-07 3.225833e-07
[336,] 0.9999996 8.927978e-07 4.463989e-07
[337,] 0.9999995 1.058697e-06 5.293486e-07
[338,] 0.9999993 1.460311e-06 7.301557e-07
[339,] 0.9999991 1.753759e-06 8.768793e-07
[340,] 0.9999988 2.326384e-06 1.163192e-06
[341,] 0.9999998 4.339704e-07 2.169852e-07
[342,] 0.9999997 5.025791e-07 2.512896e-07
[343,] 0.9999998 3.467978e-07 1.733989e-07
[344,] 0.9999997 5.164641e-07 2.582321e-07
[345,] 0.9999996 7.491212e-07 3.745606e-07
[346,] 0.9999995 1.096021e-06 5.480103e-07
[347,] 0.9999993 1.455381e-06 7.276906e-07
[348,] 0.9999989 2.131971e-06 1.065986e-06
[349,] 0.9999996 8.672224e-07 4.336112e-07
[350,] 0.9999994 1.282705e-06 6.413524e-07
[351,] 0.9999991 1.767885e-06 8.839426e-07
[352,] 0.9999988 2.325589e-06 1.162795e-06
[353,] 0.9999983 3.410412e-06 1.705206e-06
[354,] 0.9999976 4.781144e-06 2.390572e-06
[355,] 0.9999966 6.709856e-06 3.354928e-06
[356,] 0.9999952 9.502120e-06 4.751060e-06
[357,] 0.9999939 1.212719e-05 6.063596e-06
[358,] 0.9999916 1.684687e-05 8.423437e-06
[359,] 0.9999908 1.835902e-05 9.179509e-06
[360,] 0.9999876 2.486403e-05 1.243202e-05
[361,] 0.9999862 2.754364e-05 1.377182e-05
[362,] 0.9999825 3.505377e-05 1.752689e-05
[363,] 0.9999779 4.423476e-05 2.211738e-05
[364,] 0.9999736 5.283767e-05 2.641884e-05
[365,] 0.9999972 5.540350e-06 2.770175e-06
[366,] 0.9999972 5.520091e-06 2.760045e-06
[367,] 0.9999965 6.969046e-06 3.484523e-06
[368,] 0.9999950 9.971272e-06 4.985636e-06
[369,] 0.9999928 1.438631e-05 7.193156e-06
[370,] 0.9999896 2.086138e-05 1.043069e-05
[371,] 0.9999857 2.859270e-05 1.429635e-05
[372,] 0.9999796 4.076361e-05 2.038181e-05
[373,] 0.9999762 4.757631e-05 2.378816e-05
[374,] 0.9999664 6.723322e-05 3.361661e-05
[375,] 0.9999531 9.377825e-05 4.688912e-05
[376,] 0.9999376 1.248160e-04 6.240802e-05
[377,] 0.9999130 1.739964e-04 8.699818e-05
[378,] 0.9998831 2.338961e-04 1.169481e-04
[379,] 0.9998430 3.139651e-04 1.569825e-04
[380,] 0.9997862 4.276639e-04 2.138319e-04
[381,] 0.9997473 5.054916e-04 2.527458e-04
[382,] 0.9996632 6.735660e-04 3.367830e-04
[383,] 0.9995762 8.476951e-04 4.238476e-04
[384,] 0.9994647 1.070523e-03 5.352614e-04
[385,] 0.9994350 1.129914e-03 5.649572e-04
[386,] 0.9992891 1.421864e-03 7.109318e-04
[387,] 0.9990881 1.823756e-03 9.118779e-04
[388,] 0.9989814 2.037178e-03 1.018589e-03
[389,] 0.9998842 2.315326e-04 1.157663e-04
[390,] 0.9998882 2.235911e-04 1.117956e-04
[391,] 0.9998978 2.043141e-04 1.021571e-04
[392,] 0.9998595 2.810849e-04 1.405424e-04
[393,] 0.9998016 3.967301e-04 1.983650e-04
[394,] 0.9997246 5.508045e-04 2.754023e-04
[395,] 0.9996670 6.659664e-04 3.329832e-04
[396,] 0.9995410 9.180526e-04 4.590263e-04
[397,] 0.9995894 8.211366e-04 4.105683e-04
[398,] 0.9994292 1.141547e-03 5.707733e-04
[399,] 0.9992219 1.556172e-03 7.780861e-04
[400,] 0.9989971 2.005731e-03 1.002865e-03
[401,] 0.9986737 2.652552e-03 1.326276e-03
[402,] 0.9982358 3.528424e-03 1.764212e-03
[403,] 0.9977031 4.593782e-03 2.296891e-03
[404,] 0.9972430 5.513970e-03 2.756985e-03
[405,] 0.9966063 6.787383e-03 3.393691e-03
[406,] 0.9955292 8.941627e-03 4.470813e-03
[407,] 0.9985889 2.822242e-03 1.411121e-03
[408,] 0.9981428 3.714315e-03 1.857157e-03
[409,] 0.9978240 4.352070e-03 2.176035e-03
[410,] 0.9975059 4.988140e-03 2.494070e-03
[411,] 0.9967694 6.461124e-03 3.230562e-03
[412,] 0.9960033 7.993332e-03 3.996666e-03
[413,] 0.9947424 1.051522e-02 5.257610e-03
[414,] 0.9930112 1.397751e-02 6.988753e-03
[415,] 0.9907949 1.841029e-02 9.205146e-03
[416,] 0.9884345 2.313109e-02 1.156555e-02
[417,] 0.9850184 2.996323e-02 1.498161e-02
[418,] 0.9807134 3.857328e-02 1.928664e-02
[419,] 0.9964874 7.025264e-03 3.512632e-03
[420,] 0.9952612 9.477534e-03 4.738767e-03
[421,] 0.9956094 8.781170e-03 4.390585e-03
[422,] 0.9941341 1.173184e-02 5.865920e-03
[423,] 0.9920729 1.585420e-02 7.927098e-03
[424,] 0.9904269 1.914613e-02 9.573065e-03
[425,] 0.9872870 2.542595e-02 1.271298e-02
[426,] 0.9835144 3.297112e-02 1.648556e-02
[427,] 0.9783962 4.320752e-02 2.160376e-02
[428,] 0.9720888 5.582231e-02 2.791115e-02
[429,] 0.9644750 7.104994e-02 3.552497e-02
[430,] 0.9550913 8.981731e-02 4.490865e-02
[431,] 0.9944588 1.108236e-02 5.541179e-03
[432,] 0.9928503 1.429945e-02 7.149724e-03
[433,] 0.9958998 8.200390e-03 4.100195e-03
[434,] 0.9943261 1.134774e-02 5.673870e-03
[435,] 0.9927661 1.446778e-02 7.233888e-03
[436,] 0.9901981 1.960372e-02 9.801861e-03
[437,] 0.9889125 2.217508e-02 1.108754e-02
[438,] 0.9847295 3.054109e-02 1.527055e-02
[439,] 0.9795068 4.098630e-02 2.049315e-02
[440,] 0.9736318 5.273639e-02 2.636819e-02
[441,] 0.9678286 6.434273e-02 3.217136e-02
[442,] 0.9591972 8.160558e-02 4.080279e-02
[443,] 0.9979292 4.141588e-03 2.070794e-03
[444,] 0.9968171 6.365828e-03 3.182914e-03
[445,] 0.9958446 8.310773e-03 4.155387e-03
[446,] 0.9958704 8.259185e-03 4.129593e-03
[447,] 0.9937740 1.245205e-02 6.226025e-03
[448,] 0.9933510 1.329799e-02 6.648997e-03
[449,] 0.9904966 1.900678e-02 9.503389e-03
[450,] 0.9859285 2.814302e-02 1.407151e-02
[451,] 0.9801282 3.974350e-02 1.987175e-02
[452,] 0.9713513 5.729746e-02 2.864873e-02
[453,] 0.9888588 2.228246e-02 1.114123e-02
[454,] 0.9833883 3.322344e-02 1.661172e-02
[455,] 0.9865105 2.697903e-02 1.348951e-02
[456,] 0.9872707 2.545856e-02 1.272928e-02
[457,] 0.9812193 3.756135e-02 1.878067e-02
[458,] 0.9723462 5.530768e-02 2.765384e-02
[459,] 0.9680344 6.393111e-02 3.196555e-02
[460,] 0.9538030 9.239392e-02 4.619696e-02
[461,] 0.9427481 1.145039e-01 5.725193e-02
[462,] 0.9279861 1.440277e-01 7.201386e-02
[463,] 0.9067703 1.864593e-01 9.322966e-02
[464,] 0.8926706 2.146589e-01 1.073294e-01
[465,] 0.9664504 6.709913e-02 3.354957e-02
[466,] 0.9786619 4.267622e-02 2.133811e-02
[467,] 0.9665418 6.691638e-02 3.345819e-02
[468,] 0.9537358 9.252834e-02 4.626417e-02
[469,] 0.9236197 1.527607e-01 7.638033e-02
[470,] 0.8742392 2.515217e-01 1.257608e-01
[471,] 0.8065094 3.869811e-01 1.934906e-01
[472,] 0.8357402 3.285196e-01 1.642598e-01
[473,] 0.7889518 4.220965e-01 2.110482e-01
[474,] 0.7411891 5.176219e-01 2.588109e-01
[475,] 0.6205193 7.589613e-01 3.794807e-01
> postscript(file="/var/yougetitorg/rcomp/tmp/1q0qf1296726865.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/yougetitorg/rcomp/tmp/2k3ct1296726865.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/yougetitorg/rcomp/tmp/3jcwx1296726865.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/yougetitorg/rcomp/tmp/4al2h1296726865.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/yougetitorg/rcomp/tmp/5lee41296726865.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 = 492
Frequency = 1
1 2 3 4 5 6
288.91692410 -60.62938294 -31.34212800 1.51039047 -27.74919833 -24.58527119
7 8 9 10 11 12
39.14758963 -1.72766549 -9.86621219 -10.36320565 -13.06021574 -8.08095312
13 14 15 16 17 18
-14.61302664 -16.35302900 -15.97374487 -11.75790628 -13.07767046 -8.96342982
19 20 21 22 23 24
49.70174660 8.19679736 -1.80795872 -2.94851005 -1.03366323 3.62671040
25 26 27 28 29 30
0.24543914 -3.76618480 -8.97374430 -9.46353830 -8.56072406 -4.75601400
31 32 33 34 35 36
57.82595363 14.15600781 5.08353378 8.17778552 8.94915813 9.93183893
37 38 39 40 41 42
4.96518712 -0.36271534 -5.34427823 -2.50680822 -7.00770351 0.63730401
43 44 45 46 47 48
61.18806568 25.87821692 11.13392088 10.39885197 0.82918488 2.55309010
49 50 51 52 53 54
4.81119905 -5.56513357 -8.08372313 -10.17432783 -2.58831736 -4.18471837
55 56 57 58 59 60
61.90610251 32.40797726 3.15349536 7.99614222 8.30639294 -3.13387582
61 62 63 64 65 66
22.73200745 -14.41779381 -18.19508254 -16.21250243 -15.96112721 -22.72923071
67 68 69 70 71 72
46.15917345 10.73857069 -3.39212624 -3.37580415 -7.66759629 2.81391082
73 74 75 76 77 78
-1.61822206 -18.88806884 -15.64919316 -15.67249601 -13.80011689 -9.99885525
79 80 81 82 83 84
38.24359347 2.85877142 -0.71706199 11.36712959 -14.08938846 8.33492850
85 86 87 88 89 90
-1.64155889 -24.25897202 -13.26703268 -15.98578378 -18.26715128 -10.91219756
91 92 93 94 95 96
41.00028824 -3.60312677 -4.87099930 -2.82908759 -12.72314112 -4.85349378
97 98 99 100 101 102
-7.92379625 -22.83206470 -23.04349490 -22.02696469 -23.06276500 -17.73381128
103 104 105 106 107 108
39.25191123 -6.30514417 -18.42384152 -12.13808684 -15.82039294 -10.19599426
109 110 111 112 113 114
-15.86521588 -25.19918055 -18.49215892 -33.19901924 -25.38530284 -20.65498927
115 116 117 118 119 120
30.13616685 -1.83892904 -26.94970950 -14.70120146 -17.39839710 -13.66564302
121 122 123 124 125 126
-16.67020336 -27.80463243 -26.93889784 -27.55858209 -26.87830424 -19.06522202
127 128 129 130 131 132
33.27109332 -7.73698724 -17.36641989 -10.89929450 -15.50401969 -7.00132076
133 134 135 136 137 138
-13.00358010 -16.76982206 -25.03270833 -21.22098622 -22.69401483 -15.38622521
139 140 141 142 143 144
36.42451683 1.25675194 -15.14348732 -3.41711618 -11.37680222 -2.19307120
145 146 147 148 149 150
-6.70420202 -18.65285879 -15.91110674 -13.18573513 -19.92706140 -5.46096244
151 152 153 154 155 156
40.44715534 8.14512015 -8.69652758 -0.29906288 -2.00055965 10.13730965
157 158 159 160 161 162
2.85805578 -2.97709691 -8.32225062 -4.21272471 -6.38754875 2.60805159
163 164 165 166 167 168
55.92190675 21.26116485 6.60340738 3.12302148 8.36889352 12.31898090
169 170 171 172 173 174
7.62134019 6.69802769 -0.61357370 6.81409455 -7.84345230 4.90323541
175 176 177 178 179 180
-4.33798601 3.02251372 -4.43290194 2.43380106 1.64556907 1.85618052
181 182 183 184 185 186
49.42520753 4.00749174 28.01500486 12.18640309 4.78004537 11.57690863
187 188 189 190 191 192
-1.39520169 8.74085544 -2.55047145 7.02096549 11.20381036 6.31054392
193 194 195 196 197 198
8.42829081 7.83951552 -0.41514736 7.39258928 -10.18592895 5.45369943
199 200 201 202 203 204
-2.97450239 3.13937369 -8.01662930 2.80143130 2.68486348 1.55474792
205 206 207 208 209 210
55.50221622 4.39272244 30.69516776 13.50255906 8.74668319 12.83881251
211 212 213 214 215 216
5.89562534 10.52811121 -3.62634074 9.40265707 14.92434978 7.82014613
217 218 219 220 221 222
14.29963250 10.04180971 -7.12160322 9.74136890 -11.71420504 5.97143373
223 224 225 226 227 228
-3.93249087 3.82770223 -14.21586848 3.16418029 -7.23061789 0.62367857
229 230 231 232 233 234
55.50107924 2.48623825 27.64129661 12.47141840 -2.48813970 11.05891847
235 236 237 238 239 240
-9.34947946 6.87804902 -4.29636511 4.61634828 6.35239352 4.35206776
241 242 243 244 245 246
4.28283424 5.34037598 -1.80185464 4.85336171 -7.50477246 3.40739149
247 248 249 250 251 252
-5.78820495 1.76050151 -11.85297140 0.79961295 4.63742016 -0.53815637
253 254 255 256 257 258
47.81667639 2.86947742 24.48184742 10.58187277 0.05086144 9.55245865
259 260 261 262 263 264
-13.44857483 6.04830779 -9.62920124 2.85845629 6.58228019 2.27691385
265 266 267 268 269 270
-1.35155156 3.75145114 -3.79814393 2.48065896 -11.56136072 1.27943036
271 272 273 274 275 276
-13.04081969 -0.85875217 -14.62580183 -2.44373431 -9.21078398 -3.77891137
277 278 279 280 281 282
45.89032658 -2.02063652 31.58570195 7.26648489 -9.73698886 7.54994026
283 284 285 286 287 288
-24.75669570 1.79597183 -16.82695124 -2.19917996 2.96989254 -2.58209745
289 290 291 292 293 294
0.08884981 -0.49992548 -11.75458836 -1.05391105 -14.81940968 -3.31934563
295 296 297 298 299 300
-17.26115437 -4.98498348 -14.83106177 -6.50624849 -10.02399095 -7.10864555
301 302 303 304 305 306
42.43145830 -5.19719471 26.17205790 3.78146795 -18.34665941 3.57363078
307 308 309 310 311 312
-20.88658194 -2.51893970 -21.24810523 -4.90319348 -4.42867678 -6.21751590
313 314 315 316 317 318
-3.71684722 -4.68203596 -16.62058435 -4.95559348 -24.78341891 -7.59157687
319 320 321 322 323 324
-12.82831780 -9.79558913 -23.60655445 -9.61324309 -13.15387719 -11.67336217
325 326 327 328 329 330
38.32421161 -9.21033804 27.07988594 -0.13227180 -21.94876515 0.44204187
331 332 333 334 335 336
-19.41331398 -5.86723953 -26.28477106 -7.54071536 -11.58599539 -9.70578700
337 338 339 340 341 342
-1.29636828 -8.42385184 -14.67232491 -7.26428164 -20.25464482 -9.34248086
343 344 345 346 347 348
-14.53807730 -10.84662507 -17.21079077 -10.93198788 -15.97400756 -11.81203851
349 350 351 352 353 354
41.42865267 -10.15135794 19.54623050 -0.53931888 -0.10821434 -0.55534272
355 356 357 358 359 360
-12.42398175 -2.49941126 -35.36064034 -5.44451532 -6.64956192 -9.66174478
361 362 363 364 365 366
-2.13370397 -6.95614695 -7.64538813 -6.21143294 -11.66977806 -6.94318882
367 368 369 370 371 372
-16.17666939 -7.96923938 -14.13591266 -9.19783679 -7.86141516 -9.39122410
373 374 375 376 377 378
52.14439926 -6.58549962 5.13653148 3.16145307 -3.16155382 -0.48097497
379 380 381 382 383 384
-5.23873277 -1.57056045 -11.66268560 -2.39683867 0.29513263 -3.53971689
385 386 387 388 389 390
4.43667458 -2.03208543 -3.83375888 -1.24728340 -8.92601256 -2.36489957
391 392 393 394 395 396
-3.64400513 -3.33929031 -10.84929193 -3.22976202 2.34977700 -4.18980646
397 398 399 400 401 402
58.81550817 -0.40818745 13.98102733 9.20711161 2.99205650 5.95494544
403 404 405 406 407 408
-8.03989823 4.30709611 -14.34597025 1.78745486 9.69206936 0.28282514
409 410 411 412 413 414
9.20307962 2.20483644 -3.02766658 -7.64626367 -5.06095427 6.51717907
415 416 417 418 419 420
50.73316965 19.13685512 9.03913690 10.49429008 -6.88400706 9.22851474
421 422 423 424 425 426
2.85750638 3.33451228 1.45054698 -5.30885610 -0.88600839 5.90094338
427 428 429 430 431 432
54.95748637 13.24165444 8.42964296 -1.37298654 -6.36383317 10.30506093
433 434 435 436 437 438
1.88937778 2.53511123 -3.29043714 -0.36962786 3.82591468 3.71144447
439 440 441 442 443 444
50.87390854 11.93036449 3.15199627 -18.26617862 -12.93270593 -3.75193244
445 446 447 448 449 450
5.01768883 -10.94273984 -15.23084566 -4.49864089 -16.21854524 -14.78009157
451 452 453 454 455 456
52.63476035 6.94702555 -22.15049532 -8.13075873 -15.58501633 0.91787330
457 458 459 460 461 462
-0.59370385 -12.65254034 -17.06771301 -6.19863045 -26.32548315 2.03695313
463 464 465 466 467 468
42.95151235 -0.44318432 -16.60706936 -8.81105619 -5.62420857 9.50442327
469 470 471 472 473 474
7.88906868 2.76093605 3.31997630 1.76063407 -12.49540094 11.98929390
475 476 477 478 479 480
50.46342507 10.72810767 -9.58062371 -6.67828584 -0.37248559 19.05702485
481 482 483 484 485 486
6.95355915 -1.03046176 -5.82420620 -6.43792250 -12.44493012 8.59920600
487 488 489 490 491 492
55.86212680 4.81607116 -10.64918751 -10.46787881 -9.51903106 3.64266105
> postscript(file="/var/yougetitorg/rcomp/tmp/6y5t41296726865.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 = 492
Frequency = 1
lag(myerror, k = 1) myerror
0 288.91692410 NA
1 -60.62938294 288.91692410
2 -31.34212800 -60.62938294
3 1.51039047 -31.34212800
4 -27.74919833 1.51039047
5 -24.58527119 -27.74919833
6 39.14758963 -24.58527119
7 -1.72766549 39.14758963
8 -9.86621219 -1.72766549
9 -10.36320565 -9.86621219
10 -13.06021574 -10.36320565
11 -8.08095312 -13.06021574
12 -14.61302664 -8.08095312
13 -16.35302900 -14.61302664
14 -15.97374487 -16.35302900
15 -11.75790628 -15.97374487
16 -13.07767046 -11.75790628
17 -8.96342982 -13.07767046
18 49.70174660 -8.96342982
19 8.19679736 49.70174660
20 -1.80795872 8.19679736
21 -2.94851005 -1.80795872
22 -1.03366323 -2.94851005
23 3.62671040 -1.03366323
24 0.24543914 3.62671040
25 -3.76618480 0.24543914
26 -8.97374430 -3.76618480
27 -9.46353830 -8.97374430
28 -8.56072406 -9.46353830
29 -4.75601400 -8.56072406
30 57.82595363 -4.75601400
31 14.15600781 57.82595363
32 5.08353378 14.15600781
33 8.17778552 5.08353378
34 8.94915813 8.17778552
35 9.93183893 8.94915813
36 4.96518712 9.93183893
37 -0.36271534 4.96518712
38 -5.34427823 -0.36271534
39 -2.50680822 -5.34427823
40 -7.00770351 -2.50680822
41 0.63730401 -7.00770351
42 61.18806568 0.63730401
43 25.87821692 61.18806568
44 11.13392088 25.87821692
45 10.39885197 11.13392088
46 0.82918488 10.39885197
47 2.55309010 0.82918488
48 4.81119905 2.55309010
49 -5.56513357 4.81119905
50 -8.08372313 -5.56513357
51 -10.17432783 -8.08372313
52 -2.58831736 -10.17432783
53 -4.18471837 -2.58831736
54 61.90610251 -4.18471837
55 32.40797726 61.90610251
56 3.15349536 32.40797726
57 7.99614222 3.15349536
58 8.30639294 7.99614222
59 -3.13387582 8.30639294
60 22.73200745 -3.13387582
61 -14.41779381 22.73200745
62 -18.19508254 -14.41779381
63 -16.21250243 -18.19508254
64 -15.96112721 -16.21250243
65 -22.72923071 -15.96112721
66 46.15917345 -22.72923071
67 10.73857069 46.15917345
68 -3.39212624 10.73857069
69 -3.37580415 -3.39212624
70 -7.66759629 -3.37580415
71 2.81391082 -7.66759629
72 -1.61822206 2.81391082
73 -18.88806884 -1.61822206
74 -15.64919316 -18.88806884
75 -15.67249601 -15.64919316
76 -13.80011689 -15.67249601
77 -9.99885525 -13.80011689
78 38.24359347 -9.99885525
79 2.85877142 38.24359347
80 -0.71706199 2.85877142
81 11.36712959 -0.71706199
82 -14.08938846 11.36712959
83 8.33492850 -14.08938846
84 -1.64155889 8.33492850
85 -24.25897202 -1.64155889
86 -13.26703268 -24.25897202
87 -15.98578378 -13.26703268
88 -18.26715128 -15.98578378
89 -10.91219756 -18.26715128
90 41.00028824 -10.91219756
91 -3.60312677 41.00028824
92 -4.87099930 -3.60312677
93 -2.82908759 -4.87099930
94 -12.72314112 -2.82908759
95 -4.85349378 -12.72314112
96 -7.92379625 -4.85349378
97 -22.83206470 -7.92379625
98 -23.04349490 -22.83206470
99 -22.02696469 -23.04349490
100 -23.06276500 -22.02696469
101 -17.73381128 -23.06276500
102 39.25191123 -17.73381128
103 -6.30514417 39.25191123
104 -18.42384152 -6.30514417
105 -12.13808684 -18.42384152
106 -15.82039294 -12.13808684
107 -10.19599426 -15.82039294
108 -15.86521588 -10.19599426
109 -25.19918055 -15.86521588
110 -18.49215892 -25.19918055
111 -33.19901924 -18.49215892
112 -25.38530284 -33.19901924
113 -20.65498927 -25.38530284
114 30.13616685 -20.65498927
115 -1.83892904 30.13616685
116 -26.94970950 -1.83892904
117 -14.70120146 -26.94970950
118 -17.39839710 -14.70120146
119 -13.66564302 -17.39839710
120 -16.67020336 -13.66564302
121 -27.80463243 -16.67020336
122 -26.93889784 -27.80463243
123 -27.55858209 -26.93889784
124 -26.87830424 -27.55858209
125 -19.06522202 -26.87830424
126 33.27109332 -19.06522202
127 -7.73698724 33.27109332
128 -17.36641989 -7.73698724
129 -10.89929450 -17.36641989
130 -15.50401969 -10.89929450
131 -7.00132076 -15.50401969
132 -13.00358010 -7.00132076
133 -16.76982206 -13.00358010
134 -25.03270833 -16.76982206
135 -21.22098622 -25.03270833
136 -22.69401483 -21.22098622
137 -15.38622521 -22.69401483
138 36.42451683 -15.38622521
139 1.25675194 36.42451683
140 -15.14348732 1.25675194
141 -3.41711618 -15.14348732
142 -11.37680222 -3.41711618
143 -2.19307120 -11.37680222
144 -6.70420202 -2.19307120
145 -18.65285879 -6.70420202
146 -15.91110674 -18.65285879
147 -13.18573513 -15.91110674
148 -19.92706140 -13.18573513
149 -5.46096244 -19.92706140
150 40.44715534 -5.46096244
151 8.14512015 40.44715534
152 -8.69652758 8.14512015
153 -0.29906288 -8.69652758
154 -2.00055965 -0.29906288
155 10.13730965 -2.00055965
156 2.85805578 10.13730965
157 -2.97709691 2.85805578
158 -8.32225062 -2.97709691
159 -4.21272471 -8.32225062
160 -6.38754875 -4.21272471
161 2.60805159 -6.38754875
162 55.92190675 2.60805159
163 21.26116485 55.92190675
164 6.60340738 21.26116485
165 3.12302148 6.60340738
166 8.36889352 3.12302148
167 12.31898090 8.36889352
168 7.62134019 12.31898090
169 6.69802769 7.62134019
170 -0.61357370 6.69802769
171 6.81409455 -0.61357370
172 -7.84345230 6.81409455
173 4.90323541 -7.84345230
174 -4.33798601 4.90323541
175 3.02251372 -4.33798601
176 -4.43290194 3.02251372
177 2.43380106 -4.43290194
178 1.64556907 2.43380106
179 1.85618052 1.64556907
180 49.42520753 1.85618052
181 4.00749174 49.42520753
182 28.01500486 4.00749174
183 12.18640309 28.01500486
184 4.78004537 12.18640309
185 11.57690863 4.78004537
186 -1.39520169 11.57690863
187 8.74085544 -1.39520169
188 -2.55047145 8.74085544
189 7.02096549 -2.55047145
190 11.20381036 7.02096549
191 6.31054392 11.20381036
192 8.42829081 6.31054392
193 7.83951552 8.42829081
194 -0.41514736 7.83951552
195 7.39258928 -0.41514736
196 -10.18592895 7.39258928
197 5.45369943 -10.18592895
198 -2.97450239 5.45369943
199 3.13937369 -2.97450239
200 -8.01662930 3.13937369
201 2.80143130 -8.01662930
202 2.68486348 2.80143130
203 1.55474792 2.68486348
204 55.50221622 1.55474792
205 4.39272244 55.50221622
206 30.69516776 4.39272244
207 13.50255906 30.69516776
208 8.74668319 13.50255906
209 12.83881251 8.74668319
210 5.89562534 12.83881251
211 10.52811121 5.89562534
212 -3.62634074 10.52811121
213 9.40265707 -3.62634074
214 14.92434978 9.40265707
215 7.82014613 14.92434978
216 14.29963250 7.82014613
217 10.04180971 14.29963250
218 -7.12160322 10.04180971
219 9.74136890 -7.12160322
220 -11.71420504 9.74136890
221 5.97143373 -11.71420504
222 -3.93249087 5.97143373
223 3.82770223 -3.93249087
224 -14.21586848 3.82770223
225 3.16418029 -14.21586848
226 -7.23061789 3.16418029
227 0.62367857 -7.23061789
228 55.50107924 0.62367857
229 2.48623825 55.50107924
230 27.64129661 2.48623825
231 12.47141840 27.64129661
232 -2.48813970 12.47141840
233 11.05891847 -2.48813970
234 -9.34947946 11.05891847
235 6.87804902 -9.34947946
236 -4.29636511 6.87804902
237 4.61634828 -4.29636511
238 6.35239352 4.61634828
239 4.35206776 6.35239352
240 4.28283424 4.35206776
241 5.34037598 4.28283424
242 -1.80185464 5.34037598
243 4.85336171 -1.80185464
244 -7.50477246 4.85336171
245 3.40739149 -7.50477246
246 -5.78820495 3.40739149
247 1.76050151 -5.78820495
248 -11.85297140 1.76050151
249 0.79961295 -11.85297140
250 4.63742016 0.79961295
251 -0.53815637 4.63742016
252 47.81667639 -0.53815637
253 2.86947742 47.81667639
254 24.48184742 2.86947742
255 10.58187277 24.48184742
256 0.05086144 10.58187277
257 9.55245865 0.05086144
258 -13.44857483 9.55245865
259 6.04830779 -13.44857483
260 -9.62920124 6.04830779
261 2.85845629 -9.62920124
262 6.58228019 2.85845629
263 2.27691385 6.58228019
264 -1.35155156 2.27691385
265 3.75145114 -1.35155156
266 -3.79814393 3.75145114
267 2.48065896 -3.79814393
268 -11.56136072 2.48065896
269 1.27943036 -11.56136072
270 -13.04081969 1.27943036
271 -0.85875217 -13.04081969
272 -14.62580183 -0.85875217
273 -2.44373431 -14.62580183
274 -9.21078398 -2.44373431
275 -3.77891137 -9.21078398
276 45.89032658 -3.77891137
277 -2.02063652 45.89032658
278 31.58570195 -2.02063652
279 7.26648489 31.58570195
280 -9.73698886 7.26648489
281 7.54994026 -9.73698886
282 -24.75669570 7.54994026
283 1.79597183 -24.75669570
284 -16.82695124 1.79597183
285 -2.19917996 -16.82695124
286 2.96989254 -2.19917996
287 -2.58209745 2.96989254
288 0.08884981 -2.58209745
289 -0.49992548 0.08884981
290 -11.75458836 -0.49992548
291 -1.05391105 -11.75458836
292 -14.81940968 -1.05391105
293 -3.31934563 -14.81940968
294 -17.26115437 -3.31934563
295 -4.98498348 -17.26115437
296 -14.83106177 -4.98498348
297 -6.50624849 -14.83106177
298 -10.02399095 -6.50624849
299 -7.10864555 -10.02399095
300 42.43145830 -7.10864555
301 -5.19719471 42.43145830
302 26.17205790 -5.19719471
303 3.78146795 26.17205790
304 -18.34665941 3.78146795
305 3.57363078 -18.34665941
306 -20.88658194 3.57363078
307 -2.51893970 -20.88658194
308 -21.24810523 -2.51893970
309 -4.90319348 -21.24810523
310 -4.42867678 -4.90319348
311 -6.21751590 -4.42867678
312 -3.71684722 -6.21751590
313 -4.68203596 -3.71684722
314 -16.62058435 -4.68203596
315 -4.95559348 -16.62058435
316 -24.78341891 -4.95559348
317 -7.59157687 -24.78341891
318 -12.82831780 -7.59157687
319 -9.79558913 -12.82831780
320 -23.60655445 -9.79558913
321 -9.61324309 -23.60655445
322 -13.15387719 -9.61324309
323 -11.67336217 -13.15387719
324 38.32421161 -11.67336217
325 -9.21033804 38.32421161
326 27.07988594 -9.21033804
327 -0.13227180 27.07988594
328 -21.94876515 -0.13227180
329 0.44204187 -21.94876515
330 -19.41331398 0.44204187
331 -5.86723953 -19.41331398
332 -26.28477106 -5.86723953
333 -7.54071536 -26.28477106
334 -11.58599539 -7.54071536
335 -9.70578700 -11.58599539
336 -1.29636828 -9.70578700
337 -8.42385184 -1.29636828
338 -14.67232491 -8.42385184
339 -7.26428164 -14.67232491
340 -20.25464482 -7.26428164
341 -9.34248086 -20.25464482
342 -14.53807730 -9.34248086
343 -10.84662507 -14.53807730
344 -17.21079077 -10.84662507
345 -10.93198788 -17.21079077
346 -15.97400756 -10.93198788
347 -11.81203851 -15.97400756
348 41.42865267 -11.81203851
349 -10.15135794 41.42865267
350 19.54623050 -10.15135794
351 -0.53931888 19.54623050
352 -0.10821434 -0.53931888
353 -0.55534272 -0.10821434
354 -12.42398175 -0.55534272
355 -2.49941126 -12.42398175
356 -35.36064034 -2.49941126
357 -5.44451532 -35.36064034
358 -6.64956192 -5.44451532
359 -9.66174478 -6.64956192
360 -2.13370397 -9.66174478
361 -6.95614695 -2.13370397
362 -7.64538813 -6.95614695
363 -6.21143294 -7.64538813
364 -11.66977806 -6.21143294
365 -6.94318882 -11.66977806
366 -16.17666939 -6.94318882
367 -7.96923938 -16.17666939
368 -14.13591266 -7.96923938
369 -9.19783679 -14.13591266
370 -7.86141516 -9.19783679
371 -9.39122410 -7.86141516
372 52.14439926 -9.39122410
373 -6.58549962 52.14439926
374 5.13653148 -6.58549962
375 3.16145307 5.13653148
376 -3.16155382 3.16145307
377 -0.48097497 -3.16155382
378 -5.23873277 -0.48097497
379 -1.57056045 -5.23873277
380 -11.66268560 -1.57056045
381 -2.39683867 -11.66268560
382 0.29513263 -2.39683867
383 -3.53971689 0.29513263
384 4.43667458 -3.53971689
385 -2.03208543 4.43667458
386 -3.83375888 -2.03208543
387 -1.24728340 -3.83375888
388 -8.92601256 -1.24728340
389 -2.36489957 -8.92601256
390 -3.64400513 -2.36489957
391 -3.33929031 -3.64400513
392 -10.84929193 -3.33929031
393 -3.22976202 -10.84929193
394 2.34977700 -3.22976202
395 -4.18980646 2.34977700
396 58.81550817 -4.18980646
397 -0.40818745 58.81550817
398 13.98102733 -0.40818745
399 9.20711161 13.98102733
400 2.99205650 9.20711161
401 5.95494544 2.99205650
402 -8.03989823 5.95494544
403 4.30709611 -8.03989823
404 -14.34597025 4.30709611
405 1.78745486 -14.34597025
406 9.69206936 1.78745486
407 0.28282514 9.69206936
408 9.20307962 0.28282514
409 2.20483644 9.20307962
410 -3.02766658 2.20483644
411 -7.64626367 -3.02766658
412 -5.06095427 -7.64626367
413 6.51717907 -5.06095427
414 50.73316965 6.51717907
415 19.13685512 50.73316965
416 9.03913690 19.13685512
417 10.49429008 9.03913690
418 -6.88400706 10.49429008
419 9.22851474 -6.88400706
420 2.85750638 9.22851474
421 3.33451228 2.85750638
422 1.45054698 3.33451228
423 -5.30885610 1.45054698
424 -0.88600839 -5.30885610
425 5.90094338 -0.88600839
426 54.95748637 5.90094338
427 13.24165444 54.95748637
428 8.42964296 13.24165444
429 -1.37298654 8.42964296
430 -6.36383317 -1.37298654
431 10.30506093 -6.36383317
432 1.88937778 10.30506093
433 2.53511123 1.88937778
434 -3.29043714 2.53511123
435 -0.36962786 -3.29043714
436 3.82591468 -0.36962786
437 3.71144447 3.82591468
438 50.87390854 3.71144447
439 11.93036449 50.87390854
440 3.15199627 11.93036449
441 -18.26617862 3.15199627
442 -12.93270593 -18.26617862
443 -3.75193244 -12.93270593
444 5.01768883 -3.75193244
445 -10.94273984 5.01768883
446 -15.23084566 -10.94273984
447 -4.49864089 -15.23084566
448 -16.21854524 -4.49864089
449 -14.78009157 -16.21854524
450 52.63476035 -14.78009157
451 6.94702555 52.63476035
452 -22.15049532 6.94702555
453 -8.13075873 -22.15049532
454 -15.58501633 -8.13075873
455 0.91787330 -15.58501633
456 -0.59370385 0.91787330
457 -12.65254034 -0.59370385
458 -17.06771301 -12.65254034
459 -6.19863045 -17.06771301
460 -26.32548315 -6.19863045
461 2.03695313 -26.32548315
462 42.95151235 2.03695313
463 -0.44318432 42.95151235
464 -16.60706936 -0.44318432
465 -8.81105619 -16.60706936
466 -5.62420857 -8.81105619
467 9.50442327 -5.62420857
468 7.88906868 9.50442327
469 2.76093605 7.88906868
470 3.31997630 2.76093605
471 1.76063407 3.31997630
472 -12.49540094 1.76063407
473 11.98929390 -12.49540094
474 50.46342507 11.98929390
475 10.72810767 50.46342507
476 -9.58062371 10.72810767
477 -6.67828584 -9.58062371
478 -0.37248559 -6.67828584
479 19.05702485 -0.37248559
480 6.95355915 19.05702485
481 -1.03046176 6.95355915
482 -5.82420620 -1.03046176
483 -6.43792250 -5.82420620
484 -12.44493012 -6.43792250
485 8.59920600 -12.44493012
486 55.86212680 8.59920600
487 4.81607116 55.86212680
488 -10.64918751 4.81607116
489 -10.46787881 -10.64918751
490 -9.51903106 -10.46787881
491 3.64266105 -9.51903106
492 NA 3.64266105
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -60.62938294 288.91692410
[2,] -31.34212800 -60.62938294
[3,] 1.51039047 -31.34212800
[4,] -27.74919833 1.51039047
[5,] -24.58527119 -27.74919833
[6,] 39.14758963 -24.58527119
[7,] -1.72766549 39.14758963
[8,] -9.86621219 -1.72766549
[9,] -10.36320565 -9.86621219
[10,] -13.06021574 -10.36320565
[11,] -8.08095312 -13.06021574
[12,] -14.61302664 -8.08095312
[13,] -16.35302900 -14.61302664
[14,] -15.97374487 -16.35302900
[15,] -11.75790628 -15.97374487
[16,] -13.07767046 -11.75790628
[17,] -8.96342982 -13.07767046
[18,] 49.70174660 -8.96342982
[19,] 8.19679736 49.70174660
[20,] -1.80795872 8.19679736
[21,] -2.94851005 -1.80795872
[22,] -1.03366323 -2.94851005
[23,] 3.62671040 -1.03366323
[24,] 0.24543914 3.62671040
[25,] -3.76618480 0.24543914
[26,] -8.97374430 -3.76618480
[27,] -9.46353830 -8.97374430
[28,] -8.56072406 -9.46353830
[29,] -4.75601400 -8.56072406
[30,] 57.82595363 -4.75601400
[31,] 14.15600781 57.82595363
[32,] 5.08353378 14.15600781
[33,] 8.17778552 5.08353378
[34,] 8.94915813 8.17778552
[35,] 9.93183893 8.94915813
[36,] 4.96518712 9.93183893
[37,] -0.36271534 4.96518712
[38,] -5.34427823 -0.36271534
[39,] -2.50680822 -5.34427823
[40,] -7.00770351 -2.50680822
[41,] 0.63730401 -7.00770351
[42,] 61.18806568 0.63730401
[43,] 25.87821692 61.18806568
[44,] 11.13392088 25.87821692
[45,] 10.39885197 11.13392088
[46,] 0.82918488 10.39885197
[47,] 2.55309010 0.82918488
[48,] 4.81119905 2.55309010
[49,] -5.56513357 4.81119905
[50,] -8.08372313 -5.56513357
[51,] -10.17432783 -8.08372313
[52,] -2.58831736 -10.17432783
[53,] -4.18471837 -2.58831736
[54,] 61.90610251 -4.18471837
[55,] 32.40797726 61.90610251
[56,] 3.15349536 32.40797726
[57,] 7.99614222 3.15349536
[58,] 8.30639294 7.99614222
[59,] -3.13387582 8.30639294
[60,] 22.73200745 -3.13387582
[61,] -14.41779381 22.73200745
[62,] -18.19508254 -14.41779381
[63,] -16.21250243 -18.19508254
[64,] -15.96112721 -16.21250243
[65,] -22.72923071 -15.96112721
[66,] 46.15917345 -22.72923071
[67,] 10.73857069 46.15917345
[68,] -3.39212624 10.73857069
[69,] -3.37580415 -3.39212624
[70,] -7.66759629 -3.37580415
[71,] 2.81391082 -7.66759629
[72,] -1.61822206 2.81391082
[73,] -18.88806884 -1.61822206
[74,] -15.64919316 -18.88806884
[75,] -15.67249601 -15.64919316
[76,] -13.80011689 -15.67249601
[77,] -9.99885525 -13.80011689
[78,] 38.24359347 -9.99885525
[79,] 2.85877142 38.24359347
[80,] -0.71706199 2.85877142
[81,] 11.36712959 -0.71706199
[82,] -14.08938846 11.36712959
[83,] 8.33492850 -14.08938846
[84,] -1.64155889 8.33492850
[85,] -24.25897202 -1.64155889
[86,] -13.26703268 -24.25897202
[87,] -15.98578378 -13.26703268
[88,] -18.26715128 -15.98578378
[89,] -10.91219756 -18.26715128
[90,] 41.00028824 -10.91219756
[91,] -3.60312677 41.00028824
[92,] -4.87099930 -3.60312677
[93,] -2.82908759 -4.87099930
[94,] -12.72314112 -2.82908759
[95,] -4.85349378 -12.72314112
[96,] -7.92379625 -4.85349378
[97,] -22.83206470 -7.92379625
[98,] -23.04349490 -22.83206470
[99,] -22.02696469 -23.04349490
[100,] -23.06276500 -22.02696469
[101,] -17.73381128 -23.06276500
[102,] 39.25191123 -17.73381128
[103,] -6.30514417 39.25191123
[104,] -18.42384152 -6.30514417
[105,] -12.13808684 -18.42384152
[106,] -15.82039294 -12.13808684
[107,] -10.19599426 -15.82039294
[108,] -15.86521588 -10.19599426
[109,] -25.19918055 -15.86521588
[110,] -18.49215892 -25.19918055
[111,] -33.19901924 -18.49215892
[112,] -25.38530284 -33.19901924
[113,] -20.65498927 -25.38530284
[114,] 30.13616685 -20.65498927
[115,] -1.83892904 30.13616685
[116,] -26.94970950 -1.83892904
[117,] -14.70120146 -26.94970950
[118,] -17.39839710 -14.70120146
[119,] -13.66564302 -17.39839710
[120,] -16.67020336 -13.66564302
[121,] -27.80463243 -16.67020336
[122,] -26.93889784 -27.80463243
[123,] -27.55858209 -26.93889784
[124,] -26.87830424 -27.55858209
[125,] -19.06522202 -26.87830424
[126,] 33.27109332 -19.06522202
[127,] -7.73698724 33.27109332
[128,] -17.36641989 -7.73698724
[129,] -10.89929450 -17.36641989
[130,] -15.50401969 -10.89929450
[131,] -7.00132076 -15.50401969
[132,] -13.00358010 -7.00132076
[133,] -16.76982206 -13.00358010
[134,] -25.03270833 -16.76982206
[135,] -21.22098622 -25.03270833
[136,] -22.69401483 -21.22098622
[137,] -15.38622521 -22.69401483
[138,] 36.42451683 -15.38622521
[139,] 1.25675194 36.42451683
[140,] -15.14348732 1.25675194
[141,] -3.41711618 -15.14348732
[142,] -11.37680222 -3.41711618
[143,] -2.19307120 -11.37680222
[144,] -6.70420202 -2.19307120
[145,] -18.65285879 -6.70420202
[146,] -15.91110674 -18.65285879
[147,] -13.18573513 -15.91110674
[148,] -19.92706140 -13.18573513
[149,] -5.46096244 -19.92706140
[150,] 40.44715534 -5.46096244
[151,] 8.14512015 40.44715534
[152,] -8.69652758 8.14512015
[153,] -0.29906288 -8.69652758
[154,] -2.00055965 -0.29906288
[155,] 10.13730965 -2.00055965
[156,] 2.85805578 10.13730965
[157,] -2.97709691 2.85805578
[158,] -8.32225062 -2.97709691
[159,] -4.21272471 -8.32225062
[160,] -6.38754875 -4.21272471
[161,] 2.60805159 -6.38754875
[162,] 55.92190675 2.60805159
[163,] 21.26116485 55.92190675
[164,] 6.60340738 21.26116485
[165,] 3.12302148 6.60340738
[166,] 8.36889352 3.12302148
[167,] 12.31898090 8.36889352
[168,] 7.62134019 12.31898090
[169,] 6.69802769 7.62134019
[170,] -0.61357370 6.69802769
[171,] 6.81409455 -0.61357370
[172,] -7.84345230 6.81409455
[173,] 4.90323541 -7.84345230
[174,] -4.33798601 4.90323541
[175,] 3.02251372 -4.33798601
[176,] -4.43290194 3.02251372
[177,] 2.43380106 -4.43290194
[178,] 1.64556907 2.43380106
[179,] 1.85618052 1.64556907
[180,] 49.42520753 1.85618052
[181,] 4.00749174 49.42520753
[182,] 28.01500486 4.00749174
[183,] 12.18640309 28.01500486
[184,] 4.78004537 12.18640309
[185,] 11.57690863 4.78004537
[186,] -1.39520169 11.57690863
[187,] 8.74085544 -1.39520169
[188,] -2.55047145 8.74085544
[189,] 7.02096549 -2.55047145
[190,] 11.20381036 7.02096549
[191,] 6.31054392 11.20381036
[192,] 8.42829081 6.31054392
[193,] 7.83951552 8.42829081
[194,] -0.41514736 7.83951552
[195,] 7.39258928 -0.41514736
[196,] -10.18592895 7.39258928
[197,] 5.45369943 -10.18592895
[198,] -2.97450239 5.45369943
[199,] 3.13937369 -2.97450239
[200,] -8.01662930 3.13937369
[201,] 2.80143130 -8.01662930
[202,] 2.68486348 2.80143130
[203,] 1.55474792 2.68486348
[204,] 55.50221622 1.55474792
[205,] 4.39272244 55.50221622
[206,] 30.69516776 4.39272244
[207,] 13.50255906 30.69516776
[208,] 8.74668319 13.50255906
[209,] 12.83881251 8.74668319
[210,] 5.89562534 12.83881251
[211,] 10.52811121 5.89562534
[212,] -3.62634074 10.52811121
[213,] 9.40265707 -3.62634074
[214,] 14.92434978 9.40265707
[215,] 7.82014613 14.92434978
[216,] 14.29963250 7.82014613
[217,] 10.04180971 14.29963250
[218,] -7.12160322 10.04180971
[219,] 9.74136890 -7.12160322
[220,] -11.71420504 9.74136890
[221,] 5.97143373 -11.71420504
[222,] -3.93249087 5.97143373
[223,] 3.82770223 -3.93249087
[224,] -14.21586848 3.82770223
[225,] 3.16418029 -14.21586848
[226,] -7.23061789 3.16418029
[227,] 0.62367857 -7.23061789
[228,] 55.50107924 0.62367857
[229,] 2.48623825 55.50107924
[230,] 27.64129661 2.48623825
[231,] 12.47141840 27.64129661
[232,] -2.48813970 12.47141840
[233,] 11.05891847 -2.48813970
[234,] -9.34947946 11.05891847
[235,] 6.87804902 -9.34947946
[236,] -4.29636511 6.87804902
[237,] 4.61634828 -4.29636511
[238,] 6.35239352 4.61634828
[239,] 4.35206776 6.35239352
[240,] 4.28283424 4.35206776
[241,] 5.34037598 4.28283424
[242,] -1.80185464 5.34037598
[243,] 4.85336171 -1.80185464
[244,] -7.50477246 4.85336171
[245,] 3.40739149 -7.50477246
[246,] -5.78820495 3.40739149
[247,] 1.76050151 -5.78820495
[248,] -11.85297140 1.76050151
[249,] 0.79961295 -11.85297140
[250,] 4.63742016 0.79961295
[251,] -0.53815637 4.63742016
[252,] 47.81667639 -0.53815637
[253,] 2.86947742 47.81667639
[254,] 24.48184742 2.86947742
[255,] 10.58187277 24.48184742
[256,] 0.05086144 10.58187277
[257,] 9.55245865 0.05086144
[258,] -13.44857483 9.55245865
[259,] 6.04830779 -13.44857483
[260,] -9.62920124 6.04830779
[261,] 2.85845629 -9.62920124
[262,] 6.58228019 2.85845629
[263,] 2.27691385 6.58228019
[264,] -1.35155156 2.27691385
[265,] 3.75145114 -1.35155156
[266,] -3.79814393 3.75145114
[267,] 2.48065896 -3.79814393
[268,] -11.56136072 2.48065896
[269,] 1.27943036 -11.56136072
[270,] -13.04081969 1.27943036
[271,] -0.85875217 -13.04081969
[272,] -14.62580183 -0.85875217
[273,] -2.44373431 -14.62580183
[274,] -9.21078398 -2.44373431
[275,] -3.77891137 -9.21078398
[276,] 45.89032658 -3.77891137
[277,] -2.02063652 45.89032658
[278,] 31.58570195 -2.02063652
[279,] 7.26648489 31.58570195
[280,] -9.73698886 7.26648489
[281,] 7.54994026 -9.73698886
[282,] -24.75669570 7.54994026
[283,] 1.79597183 -24.75669570
[284,] -16.82695124 1.79597183
[285,] -2.19917996 -16.82695124
[286,] 2.96989254 -2.19917996
[287,] -2.58209745 2.96989254
[288,] 0.08884981 -2.58209745
[289,] -0.49992548 0.08884981
[290,] -11.75458836 -0.49992548
[291,] -1.05391105 -11.75458836
[292,] -14.81940968 -1.05391105
[293,] -3.31934563 -14.81940968
[294,] -17.26115437 -3.31934563
[295,] -4.98498348 -17.26115437
[296,] -14.83106177 -4.98498348
[297,] -6.50624849 -14.83106177
[298,] -10.02399095 -6.50624849
[299,] -7.10864555 -10.02399095
[300,] 42.43145830 -7.10864555
[301,] -5.19719471 42.43145830
[302,] 26.17205790 -5.19719471
[303,] 3.78146795 26.17205790
[304,] -18.34665941 3.78146795
[305,] 3.57363078 -18.34665941
[306,] -20.88658194 3.57363078
[307,] -2.51893970 -20.88658194
[308,] -21.24810523 -2.51893970
[309,] -4.90319348 -21.24810523
[310,] -4.42867678 -4.90319348
[311,] -6.21751590 -4.42867678
[312,] -3.71684722 -6.21751590
[313,] -4.68203596 -3.71684722
[314,] -16.62058435 -4.68203596
[315,] -4.95559348 -16.62058435
[316,] -24.78341891 -4.95559348
[317,] -7.59157687 -24.78341891
[318,] -12.82831780 -7.59157687
[319,] -9.79558913 -12.82831780
[320,] -23.60655445 -9.79558913
[321,] -9.61324309 -23.60655445
[322,] -13.15387719 -9.61324309
[323,] -11.67336217 -13.15387719
[324,] 38.32421161 -11.67336217
[325,] -9.21033804 38.32421161
[326,] 27.07988594 -9.21033804
[327,] -0.13227180 27.07988594
[328,] -21.94876515 -0.13227180
[329,] 0.44204187 -21.94876515
[330,] -19.41331398 0.44204187
[331,] -5.86723953 -19.41331398
[332,] -26.28477106 -5.86723953
[333,] -7.54071536 -26.28477106
[334,] -11.58599539 -7.54071536
[335,] -9.70578700 -11.58599539
[336,] -1.29636828 -9.70578700
[337,] -8.42385184 -1.29636828
[338,] -14.67232491 -8.42385184
[339,] -7.26428164 -14.67232491
[340,] -20.25464482 -7.26428164
[341,] -9.34248086 -20.25464482
[342,] -14.53807730 -9.34248086
[343,] -10.84662507 -14.53807730
[344,] -17.21079077 -10.84662507
[345,] -10.93198788 -17.21079077
[346,] -15.97400756 -10.93198788
[347,] -11.81203851 -15.97400756
[348,] 41.42865267 -11.81203851
[349,] -10.15135794 41.42865267
[350,] 19.54623050 -10.15135794
[351,] -0.53931888 19.54623050
[352,] -0.10821434 -0.53931888
[353,] -0.55534272 -0.10821434
[354,] -12.42398175 -0.55534272
[355,] -2.49941126 -12.42398175
[356,] -35.36064034 -2.49941126
[357,] -5.44451532 -35.36064034
[358,] -6.64956192 -5.44451532
[359,] -9.66174478 -6.64956192
[360,] -2.13370397 -9.66174478
[361,] -6.95614695 -2.13370397
[362,] -7.64538813 -6.95614695
[363,] -6.21143294 -7.64538813
[364,] -11.66977806 -6.21143294
[365,] -6.94318882 -11.66977806
[366,] -16.17666939 -6.94318882
[367,] -7.96923938 -16.17666939
[368,] -14.13591266 -7.96923938
[369,] -9.19783679 -14.13591266
[370,] -7.86141516 -9.19783679
[371,] -9.39122410 -7.86141516
[372,] 52.14439926 -9.39122410
[373,] -6.58549962 52.14439926
[374,] 5.13653148 -6.58549962
[375,] 3.16145307 5.13653148
[376,] -3.16155382 3.16145307
[377,] -0.48097497 -3.16155382
[378,] -5.23873277 -0.48097497
[379,] -1.57056045 -5.23873277
[380,] -11.66268560 -1.57056045
[381,] -2.39683867 -11.66268560
[382,] 0.29513263 -2.39683867
[383,] -3.53971689 0.29513263
[384,] 4.43667458 -3.53971689
[385,] -2.03208543 4.43667458
[386,] -3.83375888 -2.03208543
[387,] -1.24728340 -3.83375888
[388,] -8.92601256 -1.24728340
[389,] -2.36489957 -8.92601256
[390,] -3.64400513 -2.36489957
[391,] -3.33929031 -3.64400513
[392,] -10.84929193 -3.33929031
[393,] -3.22976202 -10.84929193
[394,] 2.34977700 -3.22976202
[395,] -4.18980646 2.34977700
[396,] 58.81550817 -4.18980646
[397,] -0.40818745 58.81550817
[398,] 13.98102733 -0.40818745
[399,] 9.20711161 13.98102733
[400,] 2.99205650 9.20711161
[401,] 5.95494544 2.99205650
[402,] -8.03989823 5.95494544
[403,] 4.30709611 -8.03989823
[404,] -14.34597025 4.30709611
[405,] 1.78745486 -14.34597025
[406,] 9.69206936 1.78745486
[407,] 0.28282514 9.69206936
[408,] 9.20307962 0.28282514
[409,] 2.20483644 9.20307962
[410,] -3.02766658 2.20483644
[411,] -7.64626367 -3.02766658
[412,] -5.06095427 -7.64626367
[413,] 6.51717907 -5.06095427
[414,] 50.73316965 6.51717907
[415,] 19.13685512 50.73316965
[416,] 9.03913690 19.13685512
[417,] 10.49429008 9.03913690
[418,] -6.88400706 10.49429008
[419,] 9.22851474 -6.88400706
[420,] 2.85750638 9.22851474
[421,] 3.33451228 2.85750638
[422,] 1.45054698 3.33451228
[423,] -5.30885610 1.45054698
[424,] -0.88600839 -5.30885610
[425,] 5.90094338 -0.88600839
[426,] 54.95748637 5.90094338
[427,] 13.24165444 54.95748637
[428,] 8.42964296 13.24165444
[429,] -1.37298654 8.42964296
[430,] -6.36383317 -1.37298654
[431,] 10.30506093 -6.36383317
[432,] 1.88937778 10.30506093
[433,] 2.53511123 1.88937778
[434,] -3.29043714 2.53511123
[435,] -0.36962786 -3.29043714
[436,] 3.82591468 -0.36962786
[437,] 3.71144447 3.82591468
[438,] 50.87390854 3.71144447
[439,] 11.93036449 50.87390854
[440,] 3.15199627 11.93036449
[441,] -18.26617862 3.15199627
[442,] -12.93270593 -18.26617862
[443,] -3.75193244 -12.93270593
[444,] 5.01768883 -3.75193244
[445,] -10.94273984 5.01768883
[446,] -15.23084566 -10.94273984
[447,] -4.49864089 -15.23084566
[448,] -16.21854524 -4.49864089
[449,] -14.78009157 -16.21854524
[450,] 52.63476035 -14.78009157
[451,] 6.94702555 52.63476035
[452,] -22.15049532 6.94702555
[453,] -8.13075873 -22.15049532
[454,] -15.58501633 -8.13075873
[455,] 0.91787330 -15.58501633
[456,] -0.59370385 0.91787330
[457,] -12.65254034 -0.59370385
[458,] -17.06771301 -12.65254034
[459,] -6.19863045 -17.06771301
[460,] -26.32548315 -6.19863045
[461,] 2.03695313 -26.32548315
[462,] 42.95151235 2.03695313
[463,] -0.44318432 42.95151235
[464,] -16.60706936 -0.44318432
[465,] -8.81105619 -16.60706936
[466,] -5.62420857 -8.81105619
[467,] 9.50442327 -5.62420857
[468,] 7.88906868 9.50442327
[469,] 2.76093605 7.88906868
[470,] 3.31997630 2.76093605
[471,] 1.76063407 3.31997630
[472,] -12.49540094 1.76063407
[473,] 11.98929390 -12.49540094
[474,] 50.46342507 11.98929390
[475,] 10.72810767 50.46342507
[476,] -9.58062371 10.72810767
[477,] -6.67828584 -9.58062371
[478,] -0.37248559 -6.67828584
[479,] 19.05702485 -0.37248559
[480,] 6.95355915 19.05702485
[481,] -1.03046176 6.95355915
[482,] -5.82420620 -1.03046176
[483,] -6.43792250 -5.82420620
[484,] -12.44493012 -6.43792250
[485,] 8.59920600 -12.44493012
[486,] 55.86212680 8.59920600
[487,] 4.81607116 55.86212680
[488,] -10.64918751 4.81607116
[489,] -10.46787881 -10.64918751
[490,] -9.51903106 -10.46787881
[491,] 3.64266105 -9.51903106
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -60.62938294 288.91692410
2 -31.34212800 -60.62938294
3 1.51039047 -31.34212800
4 -27.74919833 1.51039047
5 -24.58527119 -27.74919833
6 39.14758963 -24.58527119
7 -1.72766549 39.14758963
8 -9.86621219 -1.72766549
9 -10.36320565 -9.86621219
10 -13.06021574 -10.36320565
11 -8.08095312 -13.06021574
12 -14.61302664 -8.08095312
13 -16.35302900 -14.61302664
14 -15.97374487 -16.35302900
15 -11.75790628 -15.97374487
16 -13.07767046 -11.75790628
17 -8.96342982 -13.07767046
18 49.70174660 -8.96342982
19 8.19679736 49.70174660
20 -1.80795872 8.19679736
21 -2.94851005 -1.80795872
22 -1.03366323 -2.94851005
23 3.62671040 -1.03366323
24 0.24543914 3.62671040
25 -3.76618480 0.24543914
26 -8.97374430 -3.76618480
27 -9.46353830 -8.97374430
28 -8.56072406 -9.46353830
29 -4.75601400 -8.56072406
30 57.82595363 -4.75601400
31 14.15600781 57.82595363
32 5.08353378 14.15600781
33 8.17778552 5.08353378
34 8.94915813 8.17778552
35 9.93183893 8.94915813
36 4.96518712 9.93183893
37 -0.36271534 4.96518712
38 -5.34427823 -0.36271534
39 -2.50680822 -5.34427823
40 -7.00770351 -2.50680822
41 0.63730401 -7.00770351
42 61.18806568 0.63730401
43 25.87821692 61.18806568
44 11.13392088 25.87821692
45 10.39885197 11.13392088
46 0.82918488 10.39885197
47 2.55309010 0.82918488
48 4.81119905 2.55309010
49 -5.56513357 4.81119905
50 -8.08372313 -5.56513357
51 -10.17432783 -8.08372313
52 -2.58831736 -10.17432783
53 -4.18471837 -2.58831736
54 61.90610251 -4.18471837
55 32.40797726 61.90610251
56 3.15349536 32.40797726
57 7.99614222 3.15349536
58 8.30639294 7.99614222
59 -3.13387582 8.30639294
60 22.73200745 -3.13387582
61 -14.41779381 22.73200745
62 -18.19508254 -14.41779381
63 -16.21250243 -18.19508254
64 -15.96112721 -16.21250243
65 -22.72923071 -15.96112721
66 46.15917345 -22.72923071
67 10.73857069 46.15917345
68 -3.39212624 10.73857069
69 -3.37580415 -3.39212624
70 -7.66759629 -3.37580415
71 2.81391082 -7.66759629
72 -1.61822206 2.81391082
73 -18.88806884 -1.61822206
74 -15.64919316 -18.88806884
75 -15.67249601 -15.64919316
76 -13.80011689 -15.67249601
77 -9.99885525 -13.80011689
78 38.24359347 -9.99885525
79 2.85877142 38.24359347
80 -0.71706199 2.85877142
81 11.36712959 -0.71706199
82 -14.08938846 11.36712959
83 8.33492850 -14.08938846
84 -1.64155889 8.33492850
85 -24.25897202 -1.64155889
86 -13.26703268 -24.25897202
87 -15.98578378 -13.26703268
88 -18.26715128 -15.98578378
89 -10.91219756 -18.26715128
90 41.00028824 -10.91219756
91 -3.60312677 41.00028824
92 -4.87099930 -3.60312677
93 -2.82908759 -4.87099930
94 -12.72314112 -2.82908759
95 -4.85349378 -12.72314112
96 -7.92379625 -4.85349378
97 -22.83206470 -7.92379625
98 -23.04349490 -22.83206470
99 -22.02696469 -23.04349490
100 -23.06276500 -22.02696469
101 -17.73381128 -23.06276500
102 39.25191123 -17.73381128
103 -6.30514417 39.25191123
104 -18.42384152 -6.30514417
105 -12.13808684 -18.42384152
106 -15.82039294 -12.13808684
107 -10.19599426 -15.82039294
108 -15.86521588 -10.19599426
109 -25.19918055 -15.86521588
110 -18.49215892 -25.19918055
111 -33.19901924 -18.49215892
112 -25.38530284 -33.19901924
113 -20.65498927 -25.38530284
114 30.13616685 -20.65498927
115 -1.83892904 30.13616685
116 -26.94970950 -1.83892904
117 -14.70120146 -26.94970950
118 -17.39839710 -14.70120146
119 -13.66564302 -17.39839710
120 -16.67020336 -13.66564302
121 -27.80463243 -16.67020336
122 -26.93889784 -27.80463243
123 -27.55858209 -26.93889784
124 -26.87830424 -27.55858209
125 -19.06522202 -26.87830424
126 33.27109332 -19.06522202
127 -7.73698724 33.27109332
128 -17.36641989 -7.73698724
129 -10.89929450 -17.36641989
130 -15.50401969 -10.89929450
131 -7.00132076 -15.50401969
132 -13.00358010 -7.00132076
133 -16.76982206 -13.00358010
134 -25.03270833 -16.76982206
135 -21.22098622 -25.03270833
136 -22.69401483 -21.22098622
137 -15.38622521 -22.69401483
138 36.42451683 -15.38622521
139 1.25675194 36.42451683
140 -15.14348732 1.25675194
141 -3.41711618 -15.14348732
142 -11.37680222 -3.41711618
143 -2.19307120 -11.37680222
144 -6.70420202 -2.19307120
145 -18.65285879 -6.70420202
146 -15.91110674 -18.65285879
147 -13.18573513 -15.91110674
148 -19.92706140 -13.18573513
149 -5.46096244 -19.92706140
150 40.44715534 -5.46096244
151 8.14512015 40.44715534
152 -8.69652758 8.14512015
153 -0.29906288 -8.69652758
154 -2.00055965 -0.29906288
155 10.13730965 -2.00055965
156 2.85805578 10.13730965
157 -2.97709691 2.85805578
158 -8.32225062 -2.97709691
159 -4.21272471 -8.32225062
160 -6.38754875 -4.21272471
161 2.60805159 -6.38754875
162 55.92190675 2.60805159
163 21.26116485 55.92190675
164 6.60340738 21.26116485
165 3.12302148 6.60340738
166 8.36889352 3.12302148
167 12.31898090 8.36889352
168 7.62134019 12.31898090
169 6.69802769 7.62134019
170 -0.61357370 6.69802769
171 6.81409455 -0.61357370
172 -7.84345230 6.81409455
173 4.90323541 -7.84345230
174 -4.33798601 4.90323541
175 3.02251372 -4.33798601
176 -4.43290194 3.02251372
177 2.43380106 -4.43290194
178 1.64556907 2.43380106
179 1.85618052 1.64556907
180 49.42520753 1.85618052
181 4.00749174 49.42520753
182 28.01500486 4.00749174
183 12.18640309 28.01500486
184 4.78004537 12.18640309
185 11.57690863 4.78004537
186 -1.39520169 11.57690863
187 8.74085544 -1.39520169
188 -2.55047145 8.74085544
189 7.02096549 -2.55047145
190 11.20381036 7.02096549
191 6.31054392 11.20381036
192 8.42829081 6.31054392
193 7.83951552 8.42829081
194 -0.41514736 7.83951552
195 7.39258928 -0.41514736
196 -10.18592895 7.39258928
197 5.45369943 -10.18592895
198 -2.97450239 5.45369943
199 3.13937369 -2.97450239
200 -8.01662930 3.13937369
201 2.80143130 -8.01662930
202 2.68486348 2.80143130
203 1.55474792 2.68486348
204 55.50221622 1.55474792
205 4.39272244 55.50221622
206 30.69516776 4.39272244
207 13.50255906 30.69516776
208 8.74668319 13.50255906
209 12.83881251 8.74668319
210 5.89562534 12.83881251
211 10.52811121 5.89562534
212 -3.62634074 10.52811121
213 9.40265707 -3.62634074
214 14.92434978 9.40265707
215 7.82014613 14.92434978
216 14.29963250 7.82014613
217 10.04180971 14.29963250
218 -7.12160322 10.04180971
219 9.74136890 -7.12160322
220 -11.71420504 9.74136890
221 5.97143373 -11.71420504
222 -3.93249087 5.97143373
223 3.82770223 -3.93249087
224 -14.21586848 3.82770223
225 3.16418029 -14.21586848
226 -7.23061789 3.16418029
227 0.62367857 -7.23061789
228 55.50107924 0.62367857
229 2.48623825 55.50107924
230 27.64129661 2.48623825
231 12.47141840 27.64129661
232 -2.48813970 12.47141840
233 11.05891847 -2.48813970
234 -9.34947946 11.05891847
235 6.87804902 -9.34947946
236 -4.29636511 6.87804902
237 4.61634828 -4.29636511
238 6.35239352 4.61634828
239 4.35206776 6.35239352
240 4.28283424 4.35206776
241 5.34037598 4.28283424
242 -1.80185464 5.34037598
243 4.85336171 -1.80185464
244 -7.50477246 4.85336171
245 3.40739149 -7.50477246
246 -5.78820495 3.40739149
247 1.76050151 -5.78820495
248 -11.85297140 1.76050151
249 0.79961295 -11.85297140
250 4.63742016 0.79961295
251 -0.53815637 4.63742016
252 47.81667639 -0.53815637
253 2.86947742 47.81667639
254 24.48184742 2.86947742
255 10.58187277 24.48184742
256 0.05086144 10.58187277
257 9.55245865 0.05086144
258 -13.44857483 9.55245865
259 6.04830779 -13.44857483
260 -9.62920124 6.04830779
261 2.85845629 -9.62920124
262 6.58228019 2.85845629
263 2.27691385 6.58228019
264 -1.35155156 2.27691385
265 3.75145114 -1.35155156
266 -3.79814393 3.75145114
267 2.48065896 -3.79814393
268 -11.56136072 2.48065896
269 1.27943036 -11.56136072
270 -13.04081969 1.27943036
271 -0.85875217 -13.04081969
272 -14.62580183 -0.85875217
273 -2.44373431 -14.62580183
274 -9.21078398 -2.44373431
275 -3.77891137 -9.21078398
276 45.89032658 -3.77891137
277 -2.02063652 45.89032658
278 31.58570195 -2.02063652
279 7.26648489 31.58570195
280 -9.73698886 7.26648489
281 7.54994026 -9.73698886
282 -24.75669570 7.54994026
283 1.79597183 -24.75669570
284 -16.82695124 1.79597183
285 -2.19917996 -16.82695124
286 2.96989254 -2.19917996
287 -2.58209745 2.96989254
288 0.08884981 -2.58209745
289 -0.49992548 0.08884981
290 -11.75458836 -0.49992548
291 -1.05391105 -11.75458836
292 -14.81940968 -1.05391105
293 -3.31934563 -14.81940968
294 -17.26115437 -3.31934563
295 -4.98498348 -17.26115437
296 -14.83106177 -4.98498348
297 -6.50624849 -14.83106177
298 -10.02399095 -6.50624849
299 -7.10864555 -10.02399095
300 42.43145830 -7.10864555
301 -5.19719471 42.43145830
302 26.17205790 -5.19719471
303 3.78146795 26.17205790
304 -18.34665941 3.78146795
305 3.57363078 -18.34665941
306 -20.88658194 3.57363078
307 -2.51893970 -20.88658194
308 -21.24810523 -2.51893970
309 -4.90319348 -21.24810523
310 -4.42867678 -4.90319348
311 -6.21751590 -4.42867678
312 -3.71684722 -6.21751590
313 -4.68203596 -3.71684722
314 -16.62058435 -4.68203596
315 -4.95559348 -16.62058435
316 -24.78341891 -4.95559348
317 -7.59157687 -24.78341891
318 -12.82831780 -7.59157687
319 -9.79558913 -12.82831780
320 -23.60655445 -9.79558913
321 -9.61324309 -23.60655445
322 -13.15387719 -9.61324309
323 -11.67336217 -13.15387719
324 38.32421161 -11.67336217
325 -9.21033804 38.32421161
326 27.07988594 -9.21033804
327 -0.13227180 27.07988594
328 -21.94876515 -0.13227180
329 0.44204187 -21.94876515
330 -19.41331398 0.44204187
331 -5.86723953 -19.41331398
332 -26.28477106 -5.86723953
333 -7.54071536 -26.28477106
334 -11.58599539 -7.54071536
335 -9.70578700 -11.58599539
336 -1.29636828 -9.70578700
337 -8.42385184 -1.29636828
338 -14.67232491 -8.42385184
339 -7.26428164 -14.67232491
340 -20.25464482 -7.26428164
341 -9.34248086 -20.25464482
342 -14.53807730 -9.34248086
343 -10.84662507 -14.53807730
344 -17.21079077 -10.84662507
345 -10.93198788 -17.21079077
346 -15.97400756 -10.93198788
347 -11.81203851 -15.97400756
348 41.42865267 -11.81203851
349 -10.15135794 41.42865267
350 19.54623050 -10.15135794
351 -0.53931888 19.54623050
352 -0.10821434 -0.53931888
353 -0.55534272 -0.10821434
354 -12.42398175 -0.55534272
355 -2.49941126 -12.42398175
356 -35.36064034 -2.49941126
357 -5.44451532 -35.36064034
358 -6.64956192 -5.44451532
359 -9.66174478 -6.64956192
360 -2.13370397 -9.66174478
361 -6.95614695 -2.13370397
362 -7.64538813 -6.95614695
363 -6.21143294 -7.64538813
364 -11.66977806 -6.21143294
365 -6.94318882 -11.66977806
366 -16.17666939 -6.94318882
367 -7.96923938 -16.17666939
368 -14.13591266 -7.96923938
369 -9.19783679 -14.13591266
370 -7.86141516 -9.19783679
371 -9.39122410 -7.86141516
372 52.14439926 -9.39122410
373 -6.58549962 52.14439926
374 5.13653148 -6.58549962
375 3.16145307 5.13653148
376 -3.16155382 3.16145307
377 -0.48097497 -3.16155382
378 -5.23873277 -0.48097497
379 -1.57056045 -5.23873277
380 -11.66268560 -1.57056045
381 -2.39683867 -11.66268560
382 0.29513263 -2.39683867
383 -3.53971689 0.29513263
384 4.43667458 -3.53971689
385 -2.03208543 4.43667458
386 -3.83375888 -2.03208543
387 -1.24728340 -3.83375888
388 -8.92601256 -1.24728340
389 -2.36489957 -8.92601256
390 -3.64400513 -2.36489957
391 -3.33929031 -3.64400513
392 -10.84929193 -3.33929031
393 -3.22976202 -10.84929193
394 2.34977700 -3.22976202
395 -4.18980646 2.34977700
396 58.81550817 -4.18980646
397 -0.40818745 58.81550817
398 13.98102733 -0.40818745
399 9.20711161 13.98102733
400 2.99205650 9.20711161
401 5.95494544 2.99205650
402 -8.03989823 5.95494544
403 4.30709611 -8.03989823
404 -14.34597025 4.30709611
405 1.78745486 -14.34597025
406 9.69206936 1.78745486
407 0.28282514 9.69206936
408 9.20307962 0.28282514
409 2.20483644 9.20307962
410 -3.02766658 2.20483644
411 -7.64626367 -3.02766658
412 -5.06095427 -7.64626367
413 6.51717907 -5.06095427
414 50.73316965 6.51717907
415 19.13685512 50.73316965
416 9.03913690 19.13685512
417 10.49429008 9.03913690
418 -6.88400706 10.49429008
419 9.22851474 -6.88400706
420 2.85750638 9.22851474
421 3.33451228 2.85750638
422 1.45054698 3.33451228
423 -5.30885610 1.45054698
424 -0.88600839 -5.30885610
425 5.90094338 -0.88600839
426 54.95748637 5.90094338
427 13.24165444 54.95748637
428 8.42964296 13.24165444
429 -1.37298654 8.42964296
430 -6.36383317 -1.37298654
431 10.30506093 -6.36383317
432 1.88937778 10.30506093
433 2.53511123 1.88937778
434 -3.29043714 2.53511123
435 -0.36962786 -3.29043714
436 3.82591468 -0.36962786
437 3.71144447 3.82591468
438 50.87390854 3.71144447
439 11.93036449 50.87390854
440 3.15199627 11.93036449
441 -18.26617862 3.15199627
442 -12.93270593 -18.26617862
443 -3.75193244 -12.93270593
444 5.01768883 -3.75193244
445 -10.94273984 5.01768883
446 -15.23084566 -10.94273984
447 -4.49864089 -15.23084566
448 -16.21854524 -4.49864089
449 -14.78009157 -16.21854524
450 52.63476035 -14.78009157
451 6.94702555 52.63476035
452 -22.15049532 6.94702555
453 -8.13075873 -22.15049532
454 -15.58501633 -8.13075873
455 0.91787330 -15.58501633
456 -0.59370385 0.91787330
457 -12.65254034 -0.59370385
458 -17.06771301 -12.65254034
459 -6.19863045 -17.06771301
460 -26.32548315 -6.19863045
461 2.03695313 -26.32548315
462 42.95151235 2.03695313
463 -0.44318432 42.95151235
464 -16.60706936 -0.44318432
465 -8.81105619 -16.60706936
466 -5.62420857 -8.81105619
467 9.50442327 -5.62420857
468 7.88906868 9.50442327
469 2.76093605 7.88906868
470 3.31997630 2.76093605
471 1.76063407 3.31997630
472 -12.49540094 1.76063407
473 11.98929390 -12.49540094
474 50.46342507 11.98929390
475 10.72810767 50.46342507
476 -9.58062371 10.72810767
477 -6.67828584 -9.58062371
478 -0.37248559 -6.67828584
479 19.05702485 -0.37248559
480 6.95355915 19.05702485
481 -1.03046176 6.95355915
482 -5.82420620 -1.03046176
483 -6.43792250 -5.82420620
484 -12.44493012 -6.43792250
485 8.59920600 -12.44493012
486 55.86212680 8.59920600
487 4.81607116 55.86212680
488 -10.64918751 4.81607116
489 -10.46787881 -10.64918751
490 -9.51903106 -10.46787881
491 3.64266105 -9.51903106
> 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/yougetitorg/rcomp/tmp/7h0ep1296726865.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/yougetitorg/rcomp/tmp/8ug2e1296726865.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/yougetitorg/rcomp/tmp/9iqh61296726865.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/yougetitorg/rcomp/tmp/10dhe61296726865.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/yougetitorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/yougetitorg/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/yougetitorg/rcomp/tmp/1118tn1296726865.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/yougetitorg/rcomp/tmp/12f4fb1296726865.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/yougetitorg/rcomp/tmp/136opj1296726865.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/yougetitorg/rcomp/tmp/14gfg41296726865.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/yougetitorg/rcomp/tmp/15vw731296726865.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/yougetitorg/rcomp/tmp/16zlgk1296726865.tab")
+ }
>
> try(system("convert tmp/1q0qf1296726865.ps tmp/1q0qf1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/2k3ct1296726865.ps tmp/2k3ct1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jcwx1296726865.ps tmp/3jcwx1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/4al2h1296726865.ps tmp/4al2h1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lee41296726865.ps tmp/5lee41296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/6y5t41296726865.ps tmp/6y5t41296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/7h0ep1296726865.ps tmp/7h0ep1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/8ug2e1296726865.ps tmp/8ug2e1296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/9iqh61296726865.ps tmp/9iqh61296726865.png",intern=TRUE))
character(0)
> try(system("convert tmp/10dhe61296726865.ps tmp/10dhe61296726865.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
20.020 2.490 21.097