R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(56
+ ,3
+ ,30
+ ,112285
+ ,24188
+ ,146283
+ ,144
+ ,145
+ ,1
+ ,56
+ ,4
+ ,28
+ ,84786
+ ,18273
+ ,98364
+ ,103
+ ,101
+ ,1
+ ,54
+ ,12
+ ,38
+ ,83123
+ ,14130
+ ,86146
+ ,98
+ ,98
+ ,1
+ ,89
+ ,2
+ ,30
+ ,101193
+ ,32287
+ ,96933
+ ,135
+ ,132
+ ,1
+ ,40
+ ,1
+ ,22
+ ,38361
+ ,8654
+ ,79234
+ ,61
+ ,60
+ ,1
+ ,25
+ ,3
+ ,26
+ ,68504
+ ,9245
+ ,42551
+ ,39
+ ,38
+ ,1
+ ,92
+ ,0
+ ,25
+ ,119182
+ ,33251
+ ,195663
+ ,150
+ ,144
+ ,1
+ ,18
+ ,0
+ ,18
+ ,22807
+ ,1271
+ ,6853
+ ,5
+ ,5
+ ,1
+ ,63
+ ,0
+ ,11
+ ,17140
+ ,5279
+ ,21529
+ ,28
+ ,28
+ ,0
+ ,44
+ ,5
+ ,26
+ ,116174
+ ,27101
+ ,95757
+ ,84
+ ,84
+ ,1
+ ,33
+ ,0
+ ,25
+ ,57635
+ ,16373
+ ,85584
+ ,80
+ ,79
+ ,1
+ ,84
+ ,0
+ ,38
+ ,66198
+ ,19716
+ ,143983
+ ,130
+ ,127
+ ,1
+ ,88
+ ,7
+ ,44
+ ,71701
+ ,17753
+ ,75851
+ ,82
+ ,78
+ ,1
+ ,55
+ ,7
+ ,30
+ ,57793
+ ,9028
+ ,59238
+ ,60
+ ,60
+ ,1
+ ,60
+ ,3
+ ,40
+ ,80444
+ ,18653
+ ,93163
+ ,131
+ ,131
+ ,1
+ ,66
+ ,9
+ ,34
+ ,53855
+ ,8828
+ ,96037
+ ,84
+ ,84
+ ,1
+ ,154
+ ,0
+ ,47
+ ,97668
+ ,29498
+ ,151511
+ ,140
+ ,133
+ ,1
+ ,53
+ ,4
+ ,30
+ ,133824
+ ,27563
+ ,136368
+ ,151
+ ,150
+ ,1
+ ,119
+ ,3
+ ,31
+ ,101481
+ ,18293
+ ,112642
+ ,91
+ ,91
+ ,1
+ ,41
+ ,0
+ ,23
+ ,99645
+ ,22530
+ ,94728
+ ,138
+ ,132
+ ,1
+ ,61
+ ,7
+ ,36
+ ,114789
+ ,15977
+ ,105499
+ ,150
+ ,136
+ ,1
+ ,58
+ ,0
+ ,36
+ ,99052
+ ,35082
+ ,121527
+ ,124
+ ,124
+ ,1
+ ,75
+ ,1
+ ,30
+ ,67654
+ ,16116
+ ,127766
+ ,119
+ ,118
+ ,1
+ ,33
+ ,5
+ ,25
+ ,65553
+ ,15849
+ ,98958
+ ,73
+ ,70
+ ,1
+ ,40
+ ,7
+ ,39
+ ,97500
+ ,16026
+ ,77900
+ ,110
+ ,107
+ ,1
+ ,92
+ ,0
+ ,34
+ ,69112
+ ,26569
+ ,85646
+ ,123
+ ,119
+ ,1
+ ,100
+ ,0
+ ,31
+ ,82753
+ ,24785
+ ,98579
+ ,90
+ ,89
+ ,1
+ ,112
+ ,5
+ ,31
+ ,85323
+ ,17569
+ ,130767
+ ,116
+ ,112
+ ,1
+ ,73
+ ,0
+ ,33
+ ,72654
+ ,23825
+ ,131741
+ ,113
+ ,108
+ ,1
+ ,40
+ ,0
+ ,25
+ ,30727
+ ,7869
+ ,53907
+ ,56
+ ,52
+ ,1
+ ,45
+ ,0
+ ,33
+ ,77873
+ ,14975
+ ,178812
+ ,115
+ ,112
+ ,1
+ ,60
+ ,3
+ ,35
+ ,117478
+ ,37791
+ ,146761
+ ,119
+ ,116
+ ,1
+ ,62
+ ,4
+ ,42
+ ,74007
+ ,9605
+ ,82036
+ ,129
+ ,123
+ ,1
+ ,75
+ ,1
+ ,43
+ ,90183
+ ,27295
+ ,163253
+ ,127
+ ,125
+ ,1
+ ,31
+ ,4
+ ,30
+ ,61542
+ ,2746
+ ,27032
+ ,27
+ ,27
+ ,1
+ ,77
+ ,2
+ ,33
+ ,101494
+ ,34461
+ ,171975
+ ,175
+ ,162
+ ,1
+ ,34
+ ,0
+ ,13
+ ,27570
+ ,8098
+ ,65990
+ ,35
+ ,32
+ ,0
+ ,46
+ ,0
+ ,32
+ ,55813
+ ,4787
+ ,86572
+ ,64
+ ,64
+ ,1
+ ,99
+ ,0
+ ,36
+ ,79215
+ ,24919
+ ,159676
+ ,96
+ ,92
+ ,1
+ ,17
+ ,0
+ ,0
+ ,1423
+ ,603
+ ,1929
+ ,0
+ ,0
+ ,0
+ ,66
+ ,2
+ ,28
+ ,55461
+ ,16329
+ ,85371
+ ,84
+ ,83
+ ,1
+ ,30
+ ,1
+ ,14
+ ,31081
+ ,12558
+ ,58391
+ ,41
+ ,41
+ ,1
+ ,76
+ ,0
+ ,17
+ ,22996
+ ,7784
+ ,31580
+ ,47
+ ,47
+ ,0
+ ,146
+ ,2
+ ,32
+ ,83122
+ ,28522
+ ,136815
+ ,126
+ ,120
+ ,1
+ ,67
+ ,10
+ ,30
+ ,70106
+ ,22265
+ ,120642
+ ,105
+ ,105
+ ,1
+ ,56
+ ,6
+ ,35
+ ,60578
+ ,14459
+ ,69107
+ ,80
+ ,79
+ ,1
+ ,107
+ ,0
+ ,20
+ ,39992
+ ,14526
+ ,50495
+ ,70
+ ,65
+ ,0
+ ,58
+ ,5
+ ,28
+ ,79892
+ ,22240
+ ,108016
+ ,73
+ ,70
+ ,1
+ ,34
+ ,4
+ ,28
+ ,49810
+ ,11802
+ ,46341
+ ,57
+ ,55
+ ,1
+ ,61
+ ,1
+ ,39
+ ,71570
+ ,7623
+ ,78348
+ ,40
+ ,39
+ ,1
+ ,119
+ ,2
+ ,34
+ ,100708
+ ,11912
+ ,79336
+ ,68
+ ,67
+ ,1
+ ,42
+ ,2
+ ,26
+ ,33032
+ ,7935
+ ,56968
+ ,21
+ ,21
+ ,1
+ ,66
+ ,0
+ ,39
+ ,82875
+ ,18220
+ ,93176
+ ,127
+ ,127
+ ,1
+ ,89
+ ,8
+ ,39
+ ,139077
+ ,19199
+ ,161632
+ ,154
+ ,152
+ ,1
+ ,44
+ ,3
+ ,33
+ ,71595
+ ,19918
+ ,87850
+ ,116
+ ,113
+ ,1
+ ,66
+ ,0
+ ,28
+ ,72260
+ ,21884
+ ,127969
+ ,102
+ ,99
+ ,1
+ ,24
+ ,0
+ ,4
+ ,5950
+ ,2694
+ ,15049
+ ,7
+ ,7
+ ,1
+ ,259
+ ,8
+ ,39
+ ,115762
+ ,15808
+ ,155135
+ ,148
+ ,141
+ ,1
+ ,17
+ ,5
+ ,18
+ ,32551
+ ,3597
+ ,25109
+ ,21
+ ,21
+ ,1
+ ,64
+ ,3
+ ,14
+ ,31701
+ ,5296
+ ,45824
+ ,35
+ ,35
+ ,1
+ ,41
+ ,1
+ ,29
+ ,80670
+ ,25239
+ ,102996
+ ,112
+ ,109
+ ,1
+ ,68
+ ,5
+ ,44
+ ,143558
+ ,29801
+ ,160604
+ ,137
+ ,133
+ ,1
+ ,168
+ ,1
+ ,21
+ ,117105
+ ,18450
+ ,158051
+ ,135
+ ,123
+ ,0
+ ,43
+ ,1
+ ,16
+ ,23789
+ ,7132
+ ,44547
+ ,26
+ ,26
+ ,0
+ ,132
+ ,5
+ ,28
+ ,120733
+ ,34861
+ ,162647
+ ,230
+ ,230
+ ,1
+ ,105
+ ,0
+ ,35
+ ,105195
+ ,35940
+ ,174141
+ ,181
+ ,166
+ ,1
+ ,71
+ ,12
+ ,28
+ ,73107
+ ,16688
+ ,60622
+ ,71
+ ,68
+ ,1
+ ,112
+ ,8
+ ,38
+ ,132068
+ ,24683
+ ,179566
+ ,147
+ ,147
+ ,1
+ ,94
+ ,8
+ ,23
+ ,149193
+ ,46230
+ ,184301
+ ,190
+ ,179
+ ,1
+ ,82
+ ,8
+ ,36
+ ,46821
+ ,10387
+ ,75661
+ ,64
+ ,61
+ ,1
+ ,70
+ ,8
+ ,32
+ ,87011
+ ,21436
+ ,96144
+ ,105
+ ,101
+ ,1
+ ,57
+ ,2
+ ,29
+ ,95260
+ ,30546
+ ,129847
+ ,107
+ ,108
+ ,1
+ ,53
+ ,0
+ ,25
+ ,55183
+ ,19746
+ ,117286
+ ,94
+ ,90
+ ,1
+ ,103
+ ,5
+ ,27
+ ,106671
+ ,15977
+ ,71180
+ ,116
+ ,114
+ ,1
+ ,121
+ ,8
+ ,36
+ ,73511
+ ,22583
+ ,109377
+ ,106
+ ,103
+ ,1
+ ,62
+ ,2
+ ,28
+ ,92945
+ ,17274
+ ,85298
+ ,143
+ ,142
+ ,1
+ ,52
+ ,5
+ ,23
+ ,78664
+ ,16469
+ ,73631
+ ,81
+ ,79
+ ,1
+ ,52
+ ,12
+ ,40
+ ,70054
+ ,14251
+ ,86767
+ ,89
+ ,88
+ ,1
+ ,32
+ ,6
+ ,23
+ ,22618
+ ,3007
+ ,23824
+ ,26
+ ,25
+ ,1
+ ,62
+ ,7
+ ,40
+ ,74011
+ ,16851
+ ,93487
+ ,84
+ ,83
+ ,1
+ ,45
+ ,2
+ ,28
+ ,83737
+ ,21113
+ ,82981
+ ,113
+ ,113
+ ,1
+ ,46
+ ,0
+ ,34
+ ,69094
+ ,17401
+ ,73815
+ ,120
+ ,118
+ ,1
+ ,63
+ ,4
+ ,33
+ ,93133
+ ,23958
+ ,94552
+ ,110
+ ,110
+ ,1
+ ,75
+ ,3
+ ,28
+ ,95536
+ ,23567
+ ,132190
+ ,134
+ ,129
+ ,1
+ ,88
+ ,6
+ ,34
+ ,225920
+ ,13065
+ ,128754
+ ,54
+ ,51
+ ,1
+ ,46
+ ,2
+ ,30
+ ,62133
+ ,15358
+ ,66363
+ ,96
+ ,93
+ ,1
+ ,53
+ ,0
+ ,33
+ ,61370
+ ,14587
+ ,67808
+ ,78
+ ,76
+ ,1
+ ,37
+ ,1
+ ,22
+ ,43836
+ ,12770
+ ,61724
+ ,51
+ ,49
+ ,1
+ ,90
+ ,0
+ ,38
+ ,106117
+ ,24021
+ ,131722
+ ,121
+ ,118
+ ,1
+ ,63
+ ,5
+ ,26
+ ,38692
+ ,9648
+ ,68580
+ ,38
+ ,38
+ ,1
+ ,78
+ ,2
+ ,35
+ ,84651
+ ,20537
+ ,106175
+ ,145
+ ,141
+ ,1
+ ,25
+ ,0
+ ,8
+ ,56622
+ ,7905
+ ,55792
+ ,59
+ ,58
+ ,1
+ ,45
+ ,0
+ ,24
+ ,15986
+ ,4527
+ ,25157
+ ,27
+ ,27
+ ,1
+ ,46
+ ,5
+ ,29
+ ,95364
+ ,30495
+ ,76669
+ ,91
+ ,91
+ ,1
+ ,41
+ ,0
+ ,20
+ ,26706
+ ,7117
+ ,57283
+ ,48
+ ,48
+ ,0
+ ,144
+ ,1
+ ,29
+ ,89691
+ ,17719
+ ,105805
+ ,68
+ ,63
+ ,1
+ ,82
+ ,0
+ ,45
+ ,67267
+ ,27056
+ ,129484
+ ,58
+ ,56
+ ,1
+ ,91
+ ,1
+ ,37
+ ,126846
+ ,33473
+ ,72413
+ ,150
+ ,144
+ ,1
+ ,71
+ ,1
+ ,33
+ ,41140
+ ,9758
+ ,87831
+ ,74
+ ,73
+ ,1
+ ,63
+ ,2
+ ,33
+ ,102860
+ ,21115
+ ,96971
+ ,181
+ ,168
+ ,1
+ ,53
+ ,6
+ ,25
+ ,51715
+ ,7236
+ ,71299
+ ,65
+ ,64
+ ,1
+ ,62
+ ,1
+ ,32
+ ,55801
+ ,13790
+ ,77494
+ ,97
+ ,97
+ ,1
+ ,63
+ ,4
+ ,29
+ ,111813
+ ,32902
+ ,120336
+ ,121
+ ,117
+ ,1
+ ,32
+ ,2
+ ,28
+ ,120293
+ ,25131
+ ,93913
+ ,99
+ ,100
+ ,1
+ ,39
+ ,3
+ ,28
+ ,138599
+ ,30910
+ ,136048
+ ,152
+ ,149
+ ,1
+ ,62
+ ,0
+ ,31
+ ,161647
+ ,35947
+ ,181248
+ ,188
+ ,187
+ ,1
+ ,117
+ ,10
+ ,52
+ ,115929
+ ,29848
+ ,146123
+ ,138
+ ,127
+ ,1
+ ,34
+ ,0
+ ,21
+ ,24266
+ ,6943
+ ,32036
+ ,40
+ ,37
+ ,0
+ ,92
+ ,9
+ ,24
+ ,162901
+ ,42705
+ ,186646
+ ,254
+ ,245
+ ,1
+ ,93
+ ,7
+ ,41
+ ,109825
+ ,31808
+ ,102255
+ ,87
+ ,87
+ ,1
+ ,54
+ ,0
+ ,33
+ ,129838
+ ,26675
+ ,168237
+ ,178
+ ,177
+ ,1
+ ,144
+ ,0
+ ,32
+ ,37510
+ ,8435
+ ,64219
+ ,51
+ ,49
+ ,1
+ ,14
+ ,4
+ ,19
+ ,43750
+ ,7409
+ ,19630
+ ,49
+ ,49
+ ,1
+ ,61
+ ,4
+ ,20
+ ,40652
+ ,14993
+ ,76825
+ ,73
+ ,73
+ ,1
+ ,109
+ ,0
+ ,31
+ ,87771
+ ,36867
+ ,115338
+ ,176
+ ,177
+ ,1
+ ,38
+ ,0
+ ,31
+ ,85872
+ ,33835
+ ,109427
+ ,94
+ ,94
+ ,1
+ ,73
+ ,0
+ ,32
+ ,89275
+ ,24164
+ ,118168
+ ,120
+ ,117
+ ,1
+ ,75
+ ,1
+ ,18
+ ,44418
+ ,12607
+ ,84845
+ ,66
+ ,60
+ ,0
+ ,50
+ ,0
+ ,23
+ ,192565
+ ,22609
+ ,153197
+ ,56
+ ,55
+ ,1
+ ,61
+ ,1
+ ,17
+ ,35232
+ ,5892
+ ,29877
+ ,39
+ ,39
+ ,0
+ ,55
+ ,0
+ ,20
+ ,40909
+ ,17014
+ ,63506
+ ,66
+ ,64
+ ,0
+ ,77
+ ,0
+ ,12
+ ,13294
+ ,5394
+ ,22445
+ ,27
+ ,26
+ ,0
+ ,75
+ ,4
+ ,17
+ ,32387
+ ,9178
+ ,47695
+ ,65
+ ,64
+ ,0
+ ,72
+ ,0
+ ,30
+ ,140867
+ ,6440
+ ,68370
+ ,58
+ ,58
+ ,1
+ ,50
+ ,4
+ ,31
+ ,120662
+ ,21916
+ ,146304
+ ,98
+ ,95
+ ,1
+ ,32
+ ,4
+ ,10
+ ,21233
+ ,4011
+ ,38233
+ ,25
+ ,25
+ ,0
+ ,53
+ ,3
+ ,13
+ ,44332
+ ,5818
+ ,42071
+ ,26
+ ,26
+ ,0
+ ,42
+ ,0
+ ,22
+ ,61056
+ ,18647
+ ,50517
+ ,77
+ ,76
+ ,0
+ ,71
+ ,0
+ ,42
+ ,101338
+ ,20556
+ ,103950
+ ,130
+ ,129
+ ,1
+ ,10
+ ,0
+ ,1
+ ,1168
+ ,238
+ ,5841
+ ,11
+ ,11
+ ,1
+ ,35
+ ,5
+ ,9
+ ,13497
+ ,70
+ ,2341
+ ,2
+ ,2
+ ,0
+ ,65
+ ,0
+ ,32
+ ,65567
+ ,22392
+ ,84396
+ ,101
+ ,101
+ ,1
+ ,25
+ ,4
+ ,11
+ ,25162
+ ,3913
+ ,24610
+ ,31
+ ,28
+ ,1
+ ,66
+ ,0
+ ,25
+ ,32334
+ ,12237
+ ,35753
+ ,36
+ ,36
+ ,0
+ ,41
+ ,0
+ ,36
+ ,40735
+ ,8388
+ ,55515
+ ,120
+ ,89
+ ,1
+ ,86
+ ,1
+ ,31
+ ,91413
+ ,22120
+ ,209056
+ ,195
+ ,193
+ ,1
+ ,16
+ ,0
+ ,0
+ ,855
+ ,338
+ ,6622
+ ,4
+ ,4
+ ,1
+ ,42
+ ,5
+ ,24
+ ,97068
+ ,11727
+ ,115814
+ ,89
+ ,84
+ ,1
+ ,19
+ ,0
+ ,13
+ ,44339
+ ,3704
+ ,11609
+ ,24
+ ,23
+ ,0
+ ,19
+ ,0
+ ,8
+ ,14116
+ ,3988
+ ,13155
+ ,39
+ ,39
+ ,1
+ ,45
+ ,0
+ ,13
+ ,10288
+ ,3030
+ ,18274
+ ,14
+ ,14
+ ,0
+ ,65
+ ,0
+ ,19
+ ,65622
+ ,13520
+ ,72875
+ ,78
+ ,78
+ ,0
+ ,35
+ ,0
+ ,18
+ ,16563
+ ,1421
+ ,10112
+ ,15
+ ,14
+ ,0
+ ,95
+ ,2
+ ,33
+ ,76643
+ ,20923
+ ,142775
+ ,106
+ ,101
+ ,1
+ ,49
+ ,7
+ ,40
+ ,110681
+ ,20237
+ ,68847
+ ,83
+ ,82
+ ,1
+ ,37
+ ,1
+ ,22
+ ,29011
+ ,3219
+ ,17659
+ ,24
+ ,24
+ ,0
+ ,64
+ ,8
+ ,38
+ ,92696
+ ,3769
+ ,20112
+ ,37
+ ,36
+ ,1
+ ,38
+ ,2
+ ,24
+ ,94785
+ ,12252
+ ,61023
+ ,77
+ ,75
+ ,1
+ ,34
+ ,0
+ ,8
+ ,8773
+ ,1888
+ ,13983
+ ,16
+ ,16
+ ,1
+ ,32
+ ,2
+ ,35
+ ,83209
+ ,14497
+ ,65176
+ ,56
+ ,55
+ ,1
+ ,65
+ ,0
+ ,43
+ ,93815
+ ,28864
+ ,132432
+ ,132
+ ,131
+ ,1
+ ,52
+ ,0
+ ,43
+ ,86687
+ ,21721
+ ,112494
+ ,144
+ ,131
+ ,1
+ ,62
+ ,1
+ ,14
+ ,34553
+ ,4821
+ ,45109
+ ,40
+ ,39
+ ,0
+ ,65
+ ,3
+ ,41
+ ,105547
+ ,33644
+ ,170875
+ ,153
+ ,144
+ ,1
+ ,83
+ ,0
+ ,38
+ ,103487
+ ,15923
+ ,180759
+ ,143
+ ,139
+ ,1
+ ,95
+ ,3
+ ,45
+ ,213688
+ ,42935
+ ,214921
+ ,220
+ ,211
+ ,1
+ ,29
+ ,0
+ ,31
+ ,71220
+ ,18864
+ ,100226
+ ,79
+ ,78
+ ,1
+ ,18
+ ,0
+ ,13
+ ,23517
+ ,4977
+ ,32043
+ ,50
+ ,50
+ ,0
+ ,33
+ ,0
+ ,28
+ ,56926
+ ,7785
+ ,54454
+ ,39
+ ,39
+ ,1
+ ,247
+ ,4
+ ,31
+ ,91721
+ ,17939
+ ,78876
+ ,95
+ ,90
+ ,1
+ ,139
+ ,4
+ ,40
+ ,115168
+ ,23436
+ ,170745
+ ,169
+ ,166
+ ,1
+ ,29
+ ,11
+ ,30
+ ,111194
+ ,325
+ ,6940
+ ,12
+ ,12
+ ,1
+ ,118
+ ,0
+ ,16
+ ,51009
+ ,13539
+ ,49025
+ ,63
+ ,57
+ ,0
+ ,110
+ ,0
+ ,37
+ ,135777
+ ,34538
+ ,122037
+ ,134
+ ,133
+ ,1
+ ,67
+ ,4
+ ,30
+ ,51513
+ ,12198
+ ,53782
+ ,69
+ ,69
+ ,1
+ ,42
+ ,0
+ ,35
+ ,74163
+ ,26924
+ ,127748
+ ,119
+ ,119
+ ,1
+ ,65
+ ,1
+ ,32
+ ,51633
+ ,12716
+ ,86839
+ ,119
+ ,119
+ ,1
+ ,94
+ ,0
+ ,27
+ ,75345
+ ,8172
+ ,44830
+ ,75
+ ,65
+ ,1
+ ,64
+ ,0
+ ,20
+ ,33416
+ ,10855
+ ,77395
+ ,63
+ ,61
+ ,0
+ ,81
+ ,0
+ ,18
+ ,83305
+ ,11932
+ ,89324
+ ,55
+ ,49
+ ,0
+ ,95
+ ,9
+ ,31
+ ,98952
+ ,14300
+ ,103300
+ ,103
+ ,101
+ ,1
+ ,67
+ ,1
+ ,31
+ ,102372
+ ,25515
+ ,112283
+ ,197
+ ,196
+ ,1
+ ,63
+ ,3
+ ,21
+ ,37238
+ ,2805
+ ,10901
+ ,16
+ ,15
+ ,1
+ ,83
+ ,10
+ ,39
+ ,103772
+ ,29402
+ ,120691
+ ,140
+ ,136
+ ,1
+ ,45
+ ,5
+ ,41
+ ,123969
+ ,16440
+ ,58106
+ ,89
+ ,89
+ ,1
+ ,30
+ ,0
+ ,13
+ ,27142
+ ,11221
+ ,57140
+ ,40
+ ,40
+ ,0
+ ,70
+ ,2
+ ,32
+ ,135400
+ ,28732
+ ,122422
+ ,125
+ ,123
+ ,1
+ ,32
+ ,0
+ ,18
+ ,21399
+ ,5250
+ ,25899
+ ,21
+ ,21
+ ,0
+ ,83
+ ,1
+ ,39
+ ,130115
+ ,28608
+ ,139296
+ ,167
+ ,163
+ ,1
+ ,31
+ ,2
+ ,14
+ ,24874
+ ,8092
+ ,52678
+ ,32
+ ,29
+ ,0
+ ,67
+ ,4
+ ,7
+ ,34988
+ ,4473
+ ,23853
+ ,36
+ ,35
+ ,0
+ ,66
+ ,0
+ ,17
+ ,45549
+ ,1572
+ ,17306
+ ,13
+ ,13
+ ,0
+ ,10
+ ,0
+ ,0
+ ,6023
+ ,2065
+ ,7953
+ ,5
+ ,5
+ ,1
+ ,70
+ ,2
+ ,30
+ ,64466
+ ,14817
+ ,89455
+ ,96
+ ,96
+ ,1
+ ,103
+ ,1
+ ,37
+ ,54990
+ ,16714
+ ,147866
+ ,151
+ ,151
+ ,1
+ ,5
+ ,0
+ ,0
+ ,1644
+ ,556
+ ,4245
+ ,6
+ ,6
+ ,1
+ ,20
+ ,0
+ ,5
+ ,6179
+ ,2089
+ ,21509
+ ,13
+ ,13
+ ,1
+ ,5
+ ,0
+ ,1
+ ,3926
+ ,2658
+ ,7670
+ ,3
+ ,3
+ ,1
+ ,36
+ ,1
+ ,16
+ ,32755
+ ,10695
+ ,66675
+ ,57
+ ,56
+ ,0
+ ,34
+ ,0
+ ,32
+ ,34777
+ ,1669
+ ,14336
+ ,23
+ ,23
+ ,1
+ ,48
+ ,2
+ ,24
+ ,73224
+ ,16267
+ ,53608
+ ,61
+ ,57
+ ,1
+ ,40
+ ,0
+ ,17
+ ,27114
+ ,7768
+ ,30059
+ ,21
+ ,14
+ ,0
+ ,43
+ ,3
+ ,11
+ ,20760
+ ,7252
+ ,29668
+ ,43
+ ,43
+ ,0
+ ,31
+ ,6
+ ,24
+ ,37636
+ ,6387
+ ,22097
+ ,20
+ ,20
+ ,0
+ ,42
+ ,0
+ ,22
+ ,65461
+ ,18715
+ ,96841
+ ,82
+ ,72
+ ,0
+ ,46
+ ,2
+ ,12
+ ,30080
+ ,7936
+ ,41907
+ ,90
+ ,87
+ ,0
+ ,33
+ ,0
+ ,19
+ ,24094
+ ,8643
+ ,27080
+ ,25
+ ,21
+ ,0
+ ,18
+ ,2
+ ,13
+ ,69008
+ ,7294
+ ,35885
+ ,60
+ ,56
+ ,0
+ ,55
+ ,1
+ ,17
+ ,54968
+ ,4570
+ ,41247
+ ,61
+ ,59
+ ,0
+ ,35
+ ,1
+ ,15
+ ,46090
+ ,7185
+ ,28313
+ ,85
+ ,82
+ ,0
+ ,59
+ ,2
+ ,16
+ ,27507
+ ,10058
+ ,36845
+ ,43
+ ,43
+ ,0
+ ,19
+ ,1
+ ,24
+ ,10672
+ ,2342
+ ,16548
+ ,25
+ ,25
+ ,0
+ ,66
+ ,0
+ ,15
+ ,34029
+ ,8509
+ ,36134
+ ,41
+ ,38
+ ,0
+ ,60
+ ,1
+ ,17
+ ,46300
+ ,13275
+ ,55764
+ ,26
+ ,25
+ ,0
+ ,36
+ ,3
+ ,18
+ ,24760
+ ,6816
+ ,28910
+ ,38
+ ,38
+ ,0
+ ,25
+ ,0
+ ,20
+ ,18779
+ ,1930
+ ,13339
+ ,12
+ ,12
+ ,0
+ ,47
+ ,0
+ ,16
+ ,21280
+ ,8086
+ ,25319
+ ,29
+ ,29
+ ,0
+ ,54
+ ,0
+ ,16
+ ,40662
+ ,10737
+ ,66956
+ ,49
+ ,47
+ ,0
+ ,53
+ ,0
+ ,18
+ ,28987
+ ,8033
+ ,47487
+ ,46
+ ,45
+ ,0
+ ,40
+ ,1
+ ,22
+ ,22827
+ ,7058
+ ,52785
+ ,41
+ ,40
+ ,0
+ ,40
+ ,4
+ ,8
+ ,18513
+ ,6782
+ ,44683
+ ,31
+ ,30
+ ,0
+ ,39
+ ,0
+ ,17
+ ,30594
+ ,5401
+ ,35619
+ ,41
+ ,41
+ ,0
+ ,14
+ ,0
+ ,18
+ ,24006
+ ,6521
+ ,21920
+ ,26
+ ,25
+ ,0
+ ,45
+ ,0
+ ,16
+ ,27913
+ ,10856
+ ,45608
+ ,23
+ ,23
+ ,0
+ ,36
+ ,7
+ ,23
+ ,42744
+ ,2154
+ ,7721
+ ,14
+ ,14
+ ,0
+ ,28
+ ,2
+ ,22
+ ,12934
+ ,6117
+ ,20634
+ ,16
+ ,16
+ ,0
+ ,44
+ ,0
+ ,13
+ ,22574
+ ,5238
+ ,29788
+ ,25
+ ,26
+ ,0
+ ,30
+ ,7
+ ,13
+ ,41385
+ ,4820
+ ,31931
+ ,21
+ ,21
+ ,0
+ ,22
+ ,3
+ ,16
+ ,18653
+ ,5615
+ ,37754
+ ,32
+ ,27
+ ,0
+ ,17
+ ,0
+ ,16
+ ,18472
+ ,4272
+ ,32505
+ ,9
+ ,9
+ ,0
+ ,31
+ ,0
+ ,20
+ ,30976
+ ,8702
+ ,40557
+ ,35
+ ,33
+ ,0
+ ,55
+ ,6
+ ,22
+ ,63339
+ ,15340
+ ,94238
+ ,42
+ ,42
+ ,0
+ ,54
+ ,2
+ ,17
+ ,25568
+ ,8030
+ ,44197
+ ,68
+ ,68
+ ,0
+ ,21
+ ,0
+ ,18
+ ,33747
+ ,9526
+ ,43228
+ ,32
+ ,32
+ ,0
+ ,14
+ ,0
+ ,17
+ ,4154
+ ,1278
+ ,4103
+ ,6
+ ,6
+ ,0
+ ,81
+ ,3
+ ,12
+ ,19474
+ ,4236
+ ,44144
+ ,68
+ ,67
+ ,0
+ ,35
+ ,0
+ ,7
+ ,35130
+ ,3023
+ ,32868
+ ,33
+ ,33
+ ,0
+ ,43
+ ,1
+ ,17
+ ,39067
+ ,7196
+ ,27640
+ ,84
+ ,77
+ ,0
+ ,46
+ ,1
+ ,14
+ ,13310
+ ,3394
+ ,14063
+ ,46
+ ,46
+ ,0
+ ,30
+ ,0
+ ,23
+ ,65892
+ ,6371
+ ,28990
+ ,30
+ ,30
+ ,0
+ ,23
+ ,1
+ ,17
+ ,4143
+ ,1574
+ ,4694
+ ,0
+ ,0
+ ,0
+ ,38
+ ,0
+ ,14
+ ,28579
+ ,9620
+ ,42648
+ ,36
+ ,36
+ ,0
+ ,54
+ ,0
+ ,15
+ ,51776
+ ,6978
+ ,64329
+ ,47
+ ,46
+ ,0
+ ,20
+ ,0
+ ,17
+ ,21152
+ ,4911
+ ,21928
+ ,20
+ ,18
+ ,0
+ ,53
+ ,0
+ ,21
+ ,38084
+ ,8645
+ ,25836
+ ,50
+ ,48
+ ,0
+ ,45
+ ,0
+ ,18
+ ,27717
+ ,8987
+ ,22779
+ ,30
+ ,29
+ ,0
+ ,39
+ ,0
+ ,18
+ ,32928
+ ,5544
+ ,40820
+ ,30
+ ,28
+ ,0
+ ,20
+ ,0
+ ,17
+ ,11342
+ ,3083
+ ,27530
+ ,34
+ ,34
+ ,0
+ ,24
+ ,0
+ ,17
+ ,19499
+ ,6909
+ ,32378
+ ,33
+ ,33
+ ,0
+ ,31
+ ,0
+ ,16
+ ,16380
+ ,3189
+ ,10824
+ ,34
+ ,34
+ ,0
+ ,35
+ ,0
+ ,15
+ ,36874
+ ,6745
+ ,39613
+ ,37
+ ,33
+ ,0
+ ,151
+ ,0
+ ,21
+ ,48259
+ ,16724
+ ,60865
+ ,83
+ ,80
+ ,0
+ ,52
+ ,0
+ ,16
+ ,16734
+ ,4850
+ ,19787
+ ,32
+ ,32
+ ,0
+ ,30
+ ,2
+ ,14
+ ,28207
+ ,7025
+ ,20107
+ ,30
+ ,30
+ ,0
+ ,31
+ ,0
+ ,15
+ ,30143
+ ,6047
+ ,36605
+ ,43
+ ,41
+ ,0
+ ,29
+ ,1
+ ,17
+ ,41369
+ ,7377
+ ,40961
+ ,41
+ ,41
+ ,0
+ ,57
+ ,1
+ ,15
+ ,45833
+ ,9078
+ ,48231
+ ,51
+ ,51
+ ,0
+ ,40
+ ,0
+ ,15
+ ,29156
+ ,4605
+ ,39725
+ ,19
+ ,18
+ ,0
+ ,44
+ ,0
+ ,10
+ ,35944
+ ,3238
+ ,21455
+ ,37
+ ,34
+ ,0
+ ,25
+ ,0
+ ,6
+ ,36278
+ ,8100
+ ,23430
+ ,33
+ ,31
+ ,0
+ ,77
+ ,0
+ ,22
+ ,45588
+ ,9653
+ ,62991
+ ,41
+ ,39
+ ,0
+ ,35
+ ,0
+ ,21
+ ,45097
+ ,8914
+ ,49363
+ ,54
+ ,54
+ ,0
+ ,11
+ ,0
+ ,1
+ ,3895
+ ,786
+ ,9604
+ ,14
+ ,14
+ ,0
+ ,63
+ ,1
+ ,18
+ ,28394
+ ,6700
+ ,24552
+ ,25
+ ,24
+ ,0
+ ,44
+ ,0
+ ,17
+ ,18632
+ ,5788
+ ,31493
+ ,25
+ ,24
+ ,0
+ ,19
+ ,0
+ ,4
+ ,2325
+ ,593
+ ,3439
+ ,8
+ ,8
+ ,0
+ ,13
+ ,0
+ ,10
+ ,25139
+ ,4506
+ ,19555
+ ,26
+ ,26
+ ,0
+ ,42
+ ,0
+ ,16
+ ,27975
+ ,6382
+ ,21228
+ ,20
+ ,19
+ ,0
+ ,38
+ ,1
+ ,16
+ ,14483
+ ,5621
+ ,23177
+ ,11
+ ,11
+ ,0
+ ,29
+ ,0
+ ,9
+ ,13127
+ ,3997
+ ,22094
+ ,14
+ ,14
+ ,0
+ ,20
+ ,0
+ ,16
+ ,5839
+ ,520
+ ,2342
+ ,3
+ ,1
+ ,0
+ ,27
+ ,0
+ ,17
+ ,24069
+ ,8891
+ ,38798
+ ,40
+ ,39
+ ,0
+ ,20
+ ,0
+ ,7
+ ,3738
+ ,999
+ ,3255
+ ,5
+ ,5
+ ,0
+ ,19
+ ,0
+ ,15
+ ,18625
+ ,7067
+ ,24261
+ ,38
+ ,37
+ ,0
+ ,37
+ ,0
+ ,14
+ ,36341
+ ,4639
+ ,18511
+ ,32
+ ,32
+ ,0
+ ,26
+ ,0
+ ,14
+ ,24548
+ ,5654
+ ,40798
+ ,41
+ ,38
+ ,0
+ ,42
+ ,0
+ ,18
+ ,21792
+ ,6928
+ ,28893
+ ,46
+ ,47
+ ,0
+ ,49
+ ,0
+ ,12
+ ,26263
+ ,1514
+ ,21425
+ ,47
+ ,47
+ ,0
+ ,30
+ ,0
+ ,16
+ ,23686
+ ,9238
+ ,50276
+ ,37
+ ,37
+ ,0
+ ,49
+ ,0
+ ,21
+ ,49303
+ ,8204
+ ,37643
+ ,51
+ ,51
+ ,0
+ ,67
+ ,1
+ ,19
+ ,25659
+ ,5926
+ ,30377
+ ,49
+ ,45
+ ,0
+ ,28
+ ,0
+ ,16
+ ,28904
+ ,5785
+ ,27126
+ ,21
+ ,21
+ ,0
+ ,19
+ ,0
+ ,1
+ ,2781
+ ,4
+ ,13
+ ,1
+ ,1
+ ,0
+ ,49
+ ,1
+ ,16
+ ,29236
+ ,5930
+ ,42097
+ ,44
+ ,42
+ ,0
+ ,27
+ ,0
+ ,10
+ ,19546
+ ,3710
+ ,24451
+ ,26
+ ,26
+ ,0
+ ,30
+ ,6
+ ,19
+ ,22818
+ ,705
+ ,14335
+ ,21
+ ,21
+ ,0
+ ,22
+ ,3
+ ,12
+ ,32689
+ ,443
+ ,5084
+ ,4
+ ,4
+ ,0
+ ,12
+ ,1
+ ,2
+ ,5752
+ ,2416
+ ,9927
+ ,10
+ ,10
+ ,0
+ ,31
+ ,2
+ ,14
+ ,22197
+ ,7747
+ ,43527
+ ,43
+ ,43
+ ,0
+ ,20
+ ,0
+ ,17
+ ,20055
+ ,5432
+ ,27184
+ ,34
+ ,34
+ ,0
+ ,20
+ ,0
+ ,19
+ ,25272
+ ,4913
+ ,21610
+ ,32
+ ,31
+ ,0
+ ,39
+ ,0
+ ,14
+ ,82206
+ ,2650
+ ,20484
+ ,20
+ ,19
+ ,0
+ ,29
+ ,3
+ ,11
+ ,32073
+ ,2370
+ ,20156
+ ,34
+ ,34
+ ,0
+ ,16
+ ,1
+ ,4
+ ,5444
+ ,775
+ ,6012
+ ,6
+ ,6
+ ,0
+ ,27
+ ,0
+ ,16
+ ,20154
+ ,5576
+ ,18475
+ ,12
+ ,11
+ ,0
+ ,21
+ ,0
+ ,20
+ ,36944
+ ,1352
+ ,12645
+ ,24
+ ,24
+ ,0
+ ,19
+ ,1
+ ,12
+ ,8019
+ ,3080
+ ,11017
+ ,16
+ ,16
+ ,0
+ ,35
+ ,0
+ ,15
+ ,30884
+ ,10205
+ ,37623
+ ,72
+ ,72
+ ,0
+ ,14
+ ,0
+ ,16
+ ,19540
+ ,6095
+ ,35873
+ ,27
+ ,21
+ ,0)
+ ,dim=c(9
+ ,289)
+ ,dimnames=list(c('logins'
+ ,'shared_compendiums'
+ ,'compendiums_reviewed'
+ ,'totsize'
+ ,'totrevisions'
+ ,'totseconds'
+ ,'tothyperlinks'
+ ,'totblogs'
+ ,'course_cid')
+ ,1:289))
> y <- array(NA,dim=c(9,289),dimnames=list(c('logins','shared_compendiums','compendiums_reviewed','totsize','totrevisions','totseconds','tothyperlinks','totblogs','course_cid'),1:289))
> for (i in 1:dim(x)[1])
+ {
+ for (j in 1:dim(x)[2])
+ {
+ y[i,j] <- as.numeric(x[i,j])
+ }
+ }
> par3 = 'No Linear Trend'
> par2 = 'Do not include Seasonal Dummies'
> par1 = '4'
> 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
totsize logins shared_compendiums compendiums_reviewed totrevisions
1 112285 56 3 30 24188
2 84786 56 4 28 18273
3 83123 54 12 38 14130
4 101193 89 2 30 32287
5 38361 40 1 22 8654
6 68504 25 3 26 9245
7 119182 92 0 25 33251
8 22807 18 0 18 1271
9 17140 63 0 11 5279
10 116174 44 5 26 27101
11 57635 33 0 25 16373
12 66198 84 0 38 19716
13 71701 88 7 44 17753
14 57793 55 7 30 9028
15 80444 60 3 40 18653
16 53855 66 9 34 8828
17 97668 154 0 47 29498
18 133824 53 4 30 27563
19 101481 119 3 31 18293
20 99645 41 0 23 22530
21 114789 61 7 36 15977
22 99052 58 0 36 35082
23 67654 75 1 30 16116
24 65553 33 5 25 15849
25 97500 40 7 39 16026
26 69112 92 0 34 26569
27 82753 100 0 31 24785
28 85323 112 5 31 17569
29 72654 73 0 33 23825
30 30727 40 0 25 7869
31 77873 45 0 33 14975
32 117478 60 3 35 37791
33 74007 62 4 42 9605
34 90183 75 1 43 27295
35 61542 31 4 30 2746
36 101494 77 2 33 34461
37 27570 34 0 13 8098
38 55813 46 0 32 4787
39 79215 99 0 36 24919
40 1423 17 0 0 603
41 55461 66 2 28 16329
42 31081 30 1 14 12558
43 22996 76 0 17 7784
44 83122 146 2 32 28522
45 70106 67 10 30 22265
46 60578 56 6 35 14459
47 39992 107 0 20 14526
48 79892 58 5 28 22240
49 49810 34 4 28 11802
50 71570 61 1 39 7623
51 100708 119 2 34 11912
52 33032 42 2 26 7935
53 82875 66 0 39 18220
54 139077 89 8 39 19199
55 71595 44 3 33 19918
56 72260 66 0 28 21884
57 5950 24 0 4 2694
58 115762 259 8 39 15808
59 32551 17 5 18 3597
60 31701 64 3 14 5296
61 80670 41 1 29 25239
62 143558 68 5 44 29801
63 117105 168 1 21 18450
64 23789 43 1 16 7132
65 120733 132 5 28 34861
66 105195 105 0 35 35940
67 73107 71 12 28 16688
68 132068 112 8 38 24683
69 149193 94 8 23 46230
70 46821 82 8 36 10387
71 87011 70 8 32 21436
72 95260 57 2 29 30546
73 55183 53 0 25 19746
74 106671 103 5 27 15977
75 73511 121 8 36 22583
76 92945 62 2 28 17274
77 78664 52 5 23 16469
78 70054 52 12 40 14251
79 22618 32 6 23 3007
80 74011 62 7 40 16851
81 83737 45 2 28 21113
82 69094 46 0 34 17401
83 93133 63 4 33 23958
84 95536 75 3 28 23567
85 225920 88 6 34 13065
86 62133 46 2 30 15358
87 61370 53 0 33 14587
88 43836 37 1 22 12770
89 106117 90 0 38 24021
90 38692 63 5 26 9648
91 84651 78 2 35 20537
92 56622 25 0 8 7905
93 15986 45 0 24 4527
94 95364 46 5 29 30495
95 26706 41 0 20 7117
96 89691 144 1 29 17719
97 67267 82 0 45 27056
98 126846 91 1 37 33473
99 41140 71 1 33 9758
100 102860 63 2 33 21115
101 51715 53 6 25 7236
102 55801 62 1 32 13790
103 111813 63 4 29 32902
104 120293 32 2 28 25131
105 138599 39 3 28 30910
106 161647 62 0 31 35947
107 115929 117 10 52 29848
108 24266 34 0 21 6943
109 162901 92 9 24 42705
110 109825 93 7 41 31808
111 129838 54 0 33 26675
112 37510 144 0 32 8435
113 43750 14 4 19 7409
114 40652 61 4 20 14993
115 87771 109 0 31 36867
116 85872 38 0 31 33835
117 89275 73 0 32 24164
118 44418 75 1 18 12607
119 192565 50 0 23 22609
120 35232 61 1 17 5892
121 40909 55 0 20 17014
122 13294 77 0 12 5394
123 32387 75 4 17 9178
124 140867 72 0 30 6440
125 120662 50 4 31 21916
126 21233 32 4 10 4011
127 44332 53 3 13 5818
128 61056 42 0 22 18647
129 101338 71 0 42 20556
130 1168 10 0 1 238
131 13497 35 5 9 70
132 65567 65 0 32 22392
133 25162 25 4 11 3913
134 32334 66 0 25 12237
135 40735 41 0 36 8388
136 91413 86 1 31 22120
137 855 16 0 0 338
138 97068 42 5 24 11727
139 44339 19 0 13 3704
140 14116 19 0 8 3988
141 10288 45 0 13 3030
142 65622 65 0 19 13520
143 16563 35 0 18 1421
144 76643 95 2 33 20923
145 110681 49 7 40 20237
146 29011 37 1 22 3219
147 92696 64 8 38 3769
148 94785 38 2 24 12252
149 8773 34 0 8 1888
150 83209 32 2 35 14497
151 93815 65 0 43 28864
152 86687 52 0 43 21721
153 34553 62 1 14 4821
154 105547 65 3 41 33644
155 103487 83 0 38 15923
156 213688 95 3 45 42935
157 71220 29 0 31 18864
158 23517 18 0 13 4977
159 56926 33 0 28 7785
160 91721 247 4 31 17939
161 115168 139 4 40 23436
162 111194 29 11 30 325
163 51009 118 0 16 13539
164 135777 110 0 37 34538
165 51513 67 4 30 12198
166 74163 42 0 35 26924
167 51633 65 1 32 12716
168 75345 94 0 27 8172
169 33416 64 0 20 10855
170 83305 81 0 18 11932
171 98952 95 9 31 14300
172 102372 67 1 31 25515
173 37238 63 3 21 2805
174 103772 83 10 39 29402
175 123969 45 5 41 16440
176 27142 30 0 13 11221
177 135400 70 2 32 28732
178 21399 32 0 18 5250
179 130115 83 1 39 28608
180 24874 31 2 14 8092
181 34988 67 4 7 4473
182 45549 66 0 17 1572
183 6023 10 0 0 2065
184 64466 70 2 30 14817
185 54990 103 1 37 16714
186 1644 5 0 0 556
187 6179 20 0 5 2089
188 3926 5 0 1 2658
189 32755 36 1 16 10695
190 34777 34 0 32 1669
191 73224 48 2 24 16267
192 27114 40 0 17 7768
193 20760 43 3 11 7252
194 37636 31 6 24 6387
195 65461 42 0 22 18715
196 30080 46 2 12 7936
197 24094 33 0 19 8643
198 69008 18 2 13 7294
199 54968 55 1 17 4570
200 46090 35 1 15 7185
201 27507 59 2 16 10058
202 10672 19 1 24 2342
203 34029 66 0 15 8509
204 46300 60 1 17 13275
205 24760 36 3 18 6816
206 18779 25 0 20 1930
207 21280 47 0 16 8086
208 40662 54 0 16 10737
209 28987 53 0 18 8033
210 22827 40 1 22 7058
211 18513 40 4 8 6782
212 30594 39 0 17 5401
213 24006 14 0 18 6521
214 27913 45 0 16 10856
215 42744 36 7 23 2154
216 12934 28 2 22 6117
217 22574 44 0 13 5238
218 41385 30 7 13 4820
219 18653 22 3 16 5615
220 18472 17 0 16 4272
221 30976 31 0 20 8702
222 63339 55 6 22 15340
223 25568 54 2 17 8030
224 33747 21 0 18 9526
225 4154 14 0 17 1278
226 19474 81 3 12 4236
227 35130 35 0 7 3023
228 39067 43 1 17 7196
229 13310 46 1 14 3394
230 65892 30 0 23 6371
231 4143 23 1 17 1574
232 28579 38 0 14 9620
233 51776 54 0 15 6978
234 21152 20 0 17 4911
235 38084 53 0 21 8645
236 27717 45 0 18 8987
237 32928 39 0 18 5544
238 11342 20 0 17 3083
239 19499 24 0 17 6909
240 16380 31 0 16 3189
241 36874 35 0 15 6745
242 48259 151 0 21 16724
243 16734 52 0 16 4850
244 28207 30 2 14 7025
245 30143 31 0 15 6047
246 41369 29 1 17 7377
247 45833 57 1 15 9078
248 29156 40 0 15 4605
249 35944 44 0 10 3238
250 36278 25 0 6 8100
251 45588 77 0 22 9653
252 45097 35 0 21 8914
253 3895 11 0 1 786
254 28394 63 1 18 6700
255 18632 44 0 17 5788
256 2325 19 0 4 593
257 25139 13 0 10 4506
258 27975 42 0 16 6382
259 14483 38 1 16 5621
260 13127 29 0 9 3997
261 5839 20 0 16 520
262 24069 27 0 17 8891
263 3738 20 0 7 999
264 18625 19 0 15 7067
265 36341 37 0 14 4639
266 24548 26 0 14 5654
267 21792 42 0 18 6928
268 26263 49 0 12 1514
269 23686 30 0 16 9238
270 49303 49 0 21 8204
271 25659 67 1 19 5926
272 28904 28 0 16 5785
273 2781 19 0 1 4
274 29236 49 1 16 5930
275 19546 27 0 10 3710
276 22818 30 6 19 705
277 32689 22 3 12 443
278 5752 12 1 2 2416
279 22197 31 2 14 7747
280 20055 20 0 17 5432
281 25272 20 0 19 4913
282 82206 39 0 14 2650
283 32073 29 3 11 2370
284 5444 16 1 4 775
285 20154 27 0 16 5576
286 36944 21 0 20 1352
287 8019 19 1 12 3080
288 30884 35 0 15 10205
289 19540 14 0 16 6095
totseconds tothyperlinks totblogs course_cid
1 146283 144 145 1
2 98364 103 101 1
3 86146 98 98 1
4 96933 135 132 1
5 79234 61 60 1
6 42551 39 38 1
7 195663 150 144 1
8 6853 5 5 1
9 21529 28 28 0
10 95757 84 84 1
11 85584 80 79 1
12 143983 130 127 1
13 75851 82 78 1
14 59238 60 60 1
15 93163 131 131 1
16 96037 84 84 1
17 151511 140 133 1
18 136368 151 150 1
19 112642 91 91 1
20 94728 138 132 1
21 105499 150 136 1
22 121527 124 124 1
23 127766 119 118 1
24 98958 73 70 1
25 77900 110 107 1
26 85646 123 119 1
27 98579 90 89 1
28 130767 116 112 1
29 131741 113 108 1
30 53907 56 52 1
31 178812 115 112 1
32 146761 119 116 1
33 82036 129 123 1
34 163253 127 125 1
35 27032 27 27 1
36 171975 175 162 1
37 65990 35 32 0
38 86572 64 64 1
39 159676 96 92 1
40 1929 0 0 0
41 85371 84 83 1
42 58391 41 41 1
43 31580 47 47 0
44 136815 126 120 1
45 120642 105 105 1
46 69107 80 79 1
47 50495 70 65 0
48 108016 73 70 1
49 46341 57 55 1
50 78348 40 39 1
51 79336 68 67 1
52 56968 21 21 1
53 93176 127 127 1
54 161632 154 152 1
55 87850 116 113 1
56 127969 102 99 1
57 15049 7 7 1
58 155135 148 141 1
59 25109 21 21 1
60 45824 35 35 1
61 102996 112 109 1
62 160604 137 133 1
63 158051 135 123 0
64 44547 26 26 0
65 162647 230 230 1
66 174141 181 166 1
67 60622 71 68 1
68 179566 147 147 1
69 184301 190 179 1
70 75661 64 61 1
71 96144 105 101 1
72 129847 107 108 1
73 117286 94 90 1
74 71180 116 114 1
75 109377 106 103 1
76 85298 143 142 1
77 73631 81 79 1
78 86767 89 88 1
79 23824 26 25 1
80 93487 84 83 1
81 82981 113 113 1
82 73815 120 118 1
83 94552 110 110 1
84 132190 134 129 1
85 128754 54 51 1
86 66363 96 93 1
87 67808 78 76 1
88 61724 51 49 1
89 131722 121 118 1
90 68580 38 38 1
91 106175 145 141 1
92 55792 59 58 1
93 25157 27 27 1
94 76669 91 91 1
95 57283 48 48 0
96 105805 68 63 1
97 129484 58 56 1
98 72413 150 144 1
99 87831 74 73 1
100 96971 181 168 1
101 71299 65 64 1
102 77494 97 97 1
103 120336 121 117 1
104 93913 99 100 1
105 136048 152 149 1
106 181248 188 187 1
107 146123 138 127 1
108 32036 40 37 0
109 186646 254 245 1
110 102255 87 87 1
111 168237 178 177 1
112 64219 51 49 1
113 19630 49 49 1
114 76825 73 73 1
115 115338 176 177 1
116 109427 94 94 1
117 118168 120 117 1
118 84845 66 60 0
119 153197 56 55 1
120 29877 39 39 0
121 63506 66 64 0
122 22445 27 26 0
123 47695 65 64 0
124 68370 58 58 1
125 146304 98 95 1
126 38233 25 25 0
127 42071 26 26 0
128 50517 77 76 0
129 103950 130 129 1
130 5841 11 11 1
131 2341 2 2 0
132 84396 101 101 1
133 24610 31 28 1
134 35753 36 36 0
135 55515 120 89 1
136 209056 195 193 1
137 6622 4 4 1
138 115814 89 84 1
139 11609 24 23 0
140 13155 39 39 1
141 18274 14 14 0
142 72875 78 78 0
143 10112 15 14 0
144 142775 106 101 1
145 68847 83 82 1
146 17659 24 24 0
147 20112 37 36 1
148 61023 77 75 1
149 13983 16 16 1
150 65176 56 55 1
151 132432 132 131 1
152 112494 144 131 1
153 45109 40 39 0
154 170875 153 144 1
155 180759 143 139 1
156 214921 220 211 1
157 100226 79 78 1
158 32043 50 50 0
159 54454 39 39 1
160 78876 95 90 1
161 170745 169 166 1
162 6940 12 12 1
163 49025 63 57 0
164 122037 134 133 1
165 53782 69 69 1
166 127748 119 119 1
167 86839 119 119 1
168 44830 75 65 1
169 77395 63 61 0
170 89324 55 49 0
171 103300 103 101 1
172 112283 197 196 1
173 10901 16 15 1
174 120691 140 136 1
175 58106 89 89 1
176 57140 40 40 0
177 122422 125 123 1
178 25899 21 21 0
179 139296 167 163 1
180 52678 32 29 0
181 23853 36 35 0
182 17306 13 13 0
183 7953 5 5 1
184 89455 96 96 1
185 147866 151 151 1
186 4245 6 6 1
187 21509 13 13 1
188 7670 3 3 1
189 66675 57 56 0
190 14336 23 23 1
191 53608 61 57 1
192 30059 21 14 0
193 29668 43 43 0
194 22097 20 20 0
195 96841 82 72 0
196 41907 90 87 0
197 27080 25 21 0
198 35885 60 56 0
199 41247 61 59 0
200 28313 85 82 0
201 36845 43 43 0
202 16548 25 25 0
203 36134 41 38 0
204 55764 26 25 0
205 28910 38 38 0
206 13339 12 12 0
207 25319 29 29 0
208 66956 49 47 0
209 47487 46 45 0
210 52785 41 40 0
211 44683 31 30 0
212 35619 41 41 0
213 21920 26 25 0
214 45608 23 23 0
215 7721 14 14 0
216 20634 16 16 0
217 29788 25 26 0
218 31931 21 21 0
219 37754 32 27 0
220 32505 9 9 0
221 40557 35 33 0
222 94238 42 42 0
223 44197 68 68 0
224 43228 32 32 0
225 4103 6 6 0
226 44144 68 67 0
227 32868 33 33 0
228 27640 84 77 0
229 14063 46 46 0
230 28990 30 30 0
231 4694 0 0 0
232 42648 36 36 0
233 64329 47 46 0
234 21928 20 18 0
235 25836 50 48 0
236 22779 30 29 0
237 40820 30 28 0
238 27530 34 34 0
239 32378 33 33 0
240 10824 34 34 0
241 39613 37 33 0
242 60865 83 80 0
243 19787 32 32 0
244 20107 30 30 0
245 36605 43 41 0
246 40961 41 41 0
247 48231 51 51 0
248 39725 19 18 0
249 21455 37 34 0
250 23430 33 31 0
251 62991 41 39 0
252 49363 54 54 0
253 9604 14 14 0
254 24552 25 24 0
255 31493 25 24 0
256 3439 8 8 0
257 19555 26 26 0
258 21228 20 19 0
259 23177 11 11 0
260 22094 14 14 0
261 2342 3 1 0
262 38798 40 39 0
263 3255 5 5 0
264 24261 38 37 0
265 18511 32 32 0
266 40798 41 38 0
267 28893 46 47 0
268 21425 47 47 0
269 50276 37 37 0
270 37643 51 51 0
271 30377 49 45 0
272 27126 21 21 0
273 13 1 1 0
274 42097 44 42 0
275 24451 26 26 0
276 14335 21 21 0
277 5084 4 4 0
278 9927 10 10 0
279 43527 43 43 0
280 27184 34 34 0
281 21610 32 31 0
282 20484 20 19 0
283 20156 34 34 0
284 6012 6 6 0
285 18475 12 11 0
286 12645 24 24 0
287 11017 16 16 0
288 37623 72 72 0
289 35873 27 21 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) logins shared_compendiums
528.2800 9.5830 2014.6175
compendiums_reviewed totrevisions totseconds
539.6194 1.1471 0.2528
tothyperlinks totblogs course_cid
72.9181 4.3617 5358.7324
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-42081 -10786 -3571 7324 137064
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.283e+02 3.045e+03 0.173 0.86241
logins 9.583e+00 4.593e+01 0.209 0.83487
shared_compendiums 2.015e+03 4.824e+02 4.177 3.96e-05 ***
compendiums_reviewed 5.396e+02 1.919e+02 2.813 0.00526 **
totrevisions 1.147e+00 2.782e-01 4.123 4.92e-05 ***
totseconds 2.528e-01 6.111e-02 4.136 4.68e-05 ***
tothyperlinks 7.292e+01 4.051e+02 0.180 0.85730
totblogs 4.362e+00 4.171e+02 0.010 0.99166
course_cid 5.359e+03 3.454e+03 1.552 0.12190
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 20090 on 280 degrees of freedom
Multiple R-squared: 0.7598, Adjusted R-squared: 0.7529
F-statistic: 110.7 on 8 and 280 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.307162738 6.143255e-01 6.928373e-01
[2,] 0.168059014 3.361180e-01 8.319410e-01
[3,] 0.083835179 1.676704e-01 9.161648e-01
[4,] 0.043414971 8.682994e-02 9.565850e-01
[5,] 0.022473571 4.494714e-02 9.775264e-01
[6,] 0.023393471 4.678694e-02 9.766065e-01
[7,] 0.037387459 7.477492e-02 9.626125e-01
[8,] 0.050988930 1.019779e-01 9.490111e-01
[9,] 0.067598301 1.351966e-01 9.324017e-01
[10,] 0.101840233 2.036805e-01 8.981598e-01
[11,] 0.071925249 1.438505e-01 9.280748e-01
[12,] 0.049293410 9.858682e-02 9.507066e-01
[13,] 0.041486756 8.297351e-02 9.585132e-01
[14,] 0.037423428 7.484686e-02 9.625766e-01
[15,] 0.056035286 1.120706e-01 9.439647e-01
[16,] 0.038238467 7.647693e-02 9.617615e-01
[17,] 0.025780286 5.156057e-02 9.742197e-01
[18,] 0.020041284 4.008257e-02 9.799587e-01
[19,] 0.014719414 2.943883e-02 9.852806e-01
[20,] 0.012013663 2.402733e-02 9.879863e-01
[21,] 0.007785767 1.557153e-02 9.922142e-01
[22,] 0.005370678 1.074136e-02 9.946293e-01
[23,] 0.003445684 6.891368e-03 9.965543e-01
[24,] 0.009151172 1.830234e-02 9.908488e-01
[25,] 0.013286346 2.657269e-02 9.867137e-01
[26,] 0.008937157 1.787431e-02 9.910628e-01
[27,] 0.008307689 1.661538e-02 9.916923e-01
[28,] 0.005828222 1.165644e-02 9.941718e-01
[29,] 0.004057219 8.114439e-03 9.959428e-01
[30,] 0.003798831 7.597662e-03 9.962012e-01
[31,] 0.005377918 1.075584e-02 9.946221e-01
[32,] 0.003538664 7.077328e-03 9.964613e-01
[33,] 0.002736755 5.473511e-03 9.972632e-01
[34,] 0.008298208 1.659642e-02 9.917018e-01
[35,] 0.006915767 1.383153e-02 9.930842e-01
[36,] 0.004756638 9.513277e-03 9.952434e-01
[37,] 0.003255439 6.510878e-03 9.967446e-01
[38,] 0.002411606 4.823212e-03 9.975884e-01
[39,] 0.004319437 8.638875e-03 9.956806e-01
[40,] 0.021544202 4.308840e-02 9.784558e-01
[41,] 0.018071360 3.614272e-02 9.819286e-01
[42,] 0.013181950 2.636390e-02 9.868180e-01
[43,] 0.022806789 4.561358e-02 9.771932e-01
[44,] 0.019582198 3.916440e-02 9.804178e-01
[45,] 0.015556283 3.111257e-02 9.844437e-01
[46,] 0.013203490 2.640698e-02 9.867965e-01
[47,] 0.010250464 2.050093e-02 9.897495e-01
[48,] 0.007550920 1.510184e-02 9.924491e-01
[49,] 0.005802150 1.160430e-02 9.941979e-01
[50,] 0.004195665 8.391331e-03 9.958043e-01
[51,] 0.006780204 1.356041e-02 9.932198e-01
[52,] 0.014691704 2.938341e-02 9.853083e-01
[53,] 0.011344779 2.268956e-02 9.886552e-01
[54,] 0.010333444 2.066689e-02 9.896666e-01
[55,] 0.009761510 1.952302e-02 9.902385e-01
[56,] 0.007870389 1.574078e-02 9.921296e-01
[57,] 0.006238688 1.247738e-02 9.937613e-01
[58,] 0.004707051 9.414102e-03 9.952929e-01
[59,] 0.006766971 1.353394e-02 9.932330e-01
[60,] 0.005193656 1.038731e-02 9.948063e-01
[61,] 0.003889173 7.778345e-03 9.961108e-01
[62,] 0.004024209 8.048417e-03 9.959758e-01
[63,] 0.006293945 1.258789e-02 9.937061e-01
[64,] 0.008771261 1.754252e-02 9.912287e-01
[65,] 0.007218922 1.443784e-02 9.927811e-01
[66,] 0.005826155 1.165231e-02 9.941738e-01
[67,] 0.006802754 1.360551e-02 9.931972e-01
[68,] 0.006349927 1.269985e-02 9.936501e-01
[69,] 0.005527758 1.105552e-02 9.944722e-01
[70,] 0.004173140 8.346281e-03 9.958269e-01
[71,] 0.003145073 6.290147e-03 9.968549e-01
[72,] 0.002369561 4.739122e-03 9.976304e-01
[73,] 0.001748532 3.497063e-03 9.982515e-01
[74,] 0.978200843 4.359831e-02 2.179916e-02
[75,] 0.973424623 5.315075e-02 2.657538e-02
[76,] 0.967597075 6.480585e-02 3.240292e-02
[77,] 0.962095257 7.580949e-02 3.790474e-02
[78,] 0.956352198 8.729560e-02 4.364780e-02
[79,] 0.958957792 8.208442e-02 4.104221e-02
[80,] 0.951166037 9.766793e-02 4.883396e-02
[81,] 0.948990328 1.020193e-01 5.100967e-02
[82,] 0.945218764 1.095625e-01 5.478124e-02
[83,] 0.937889090 1.242218e-01 6.211091e-02
[84,] 0.928587651 1.428247e-01 7.141235e-02
[85,] 0.920234463 1.595311e-01 7.976554e-02
[86,] 0.936153370 1.276933e-01 6.384663e-02
[87,] 0.952870858 9.425828e-02 4.712914e-02
[88,] 0.957187258 8.562548e-02 4.281274e-02
[89,] 0.951312852 9.737430e-02 4.868715e-02
[90,] 0.946227406 1.075452e-01 5.377259e-02
[91,] 0.940748635 1.185027e-01 5.925137e-02
[92,] 0.931064576 1.378708e-01 6.893542e-02
[93,] 0.952931418 9.413716e-02 4.706858e-02
[94,] 0.961263931 7.747214e-02 3.873607e-02
[95,] 0.973853941 5.229212e-02 2.614606e-02
[96,] 0.976115960 4.776808e-02 2.388404e-02
[97,] 0.971118131 5.776374e-02 2.888187e-02
[98,] 0.967491990 6.501602e-02 3.250801e-02
[99,] 0.962943469 7.411306e-02 3.705653e-02
[100,] 0.961433403 7.713319e-02 3.856660e-02
[101,] 0.963393057 7.321389e-02 3.660694e-02
[102,] 0.956060944 8.787811e-02 4.393906e-02
[103,] 0.964534338 7.093132e-02 3.546566e-02
[104,] 0.965324811 6.935038e-02 3.467519e-02
[105,] 0.961986192 7.602762e-02 3.801381e-02
[106,] 0.954783016 9.043397e-02 4.521698e-02
[107,] 0.947522131 1.049557e-01 5.247787e-02
[108,] 0.999928435 1.431310e-04 7.156549e-05
[109,] 0.999904582 1.908363e-04 9.541813e-05
[110,] 0.999875948 2.481044e-04 1.240522e-04
[111,] 0.999835319 3.293629e-04 1.646814e-04
[112,] 0.999804111 3.917781e-04 1.958890e-04
[113,] 0.999999911 1.777208e-07 8.886040e-08
[114,] 0.999999917 1.667079e-07 8.335397e-08
[115,] 0.999999877 2.460083e-07 1.230042e-07
[116,] 0.999999839 3.213822e-07 1.606911e-07
[117,] 0.999999775 4.504507e-07 2.252253e-07
[118,] 0.999999696 6.073353e-07 3.036676e-07
[119,] 0.999999561 8.787621e-07 4.393811e-07
[120,] 0.999999357 1.286556e-06 6.432782e-07
[121,] 0.999999272 1.455034e-06 7.275168e-07
[122,] 0.999999019 1.961487e-06 9.807434e-07
[123,] 0.999998713 2.574002e-06 1.287001e-06
[124,] 0.999999449 1.101180e-06 5.505899e-07
[125,] 0.999999590 8.205944e-07 4.102972e-07
[126,] 0.999999403 1.193574e-06 5.967870e-07
[127,] 0.999999374 1.252054e-06 6.260271e-07
[128,] 0.999999623 7.533485e-07 3.766743e-07
[129,] 0.999999467 1.066331e-06 5.331655e-07
[130,] 0.999999205 1.589921e-06 7.949605e-07
[131,] 0.999999167 1.665042e-06 8.325210e-07
[132,] 0.999998735 2.530984e-06 1.265492e-06
[133,] 0.999998812 2.376830e-06 1.188415e-06
[134,] 0.999998705 2.589800e-06 1.294900e-06
[135,] 0.999998062 3.876975e-06 1.938488e-06
[136,] 0.999998948 2.103321e-06 1.051660e-06
[137,] 0.999999583 8.330324e-07 4.165162e-07
[138,] 0.999999416 1.167809e-06 5.839044e-07
[139,] 0.999999281 1.437485e-06 7.187424e-07
[140,] 0.999999102 1.795980e-06 8.979899e-07
[141,] 0.999999276 1.447216e-06 7.236082e-07
[142,] 0.999998929 2.141327e-06 1.070663e-06
[143,] 0.999999573 8.548993e-07 4.274497e-07
[144,] 0.999999379 1.242287e-06 6.211437e-07
[145,] 0.999999974 5.213414e-08 2.606707e-08
[146,] 0.999999958 8.419119e-08 4.209560e-08
[147,] 0.999999932 1.365238e-07 6.826188e-08
[148,] 0.999999901 1.975431e-07 9.877155e-08
[149,] 0.999999875 2.490769e-07 1.245384e-07
[150,] 0.999999803 3.948862e-07 1.974431e-07
[151,] 0.999999998 3.981023e-09 1.990512e-09
[152,] 0.999999997 6.097076e-09 3.048538e-09
[153,] 0.999999998 3.970667e-09 1.985334e-09
[154,] 0.999999998 4.918144e-09 2.459072e-09
[155,] 0.999999998 4.681161e-09 2.340581e-09
[156,] 0.999999998 4.010240e-09 2.005120e-09
[157,] 0.999999998 4.594541e-09 2.297271e-09
[158,] 0.999999997 5.876421e-09 2.938210e-09
[159,] 0.999999999 1.244701e-09 6.223507e-10
[160,] 0.999999999 1.687825e-09 8.439124e-10
[161,] 0.999999999 2.861579e-09 1.430790e-09
[162,] 0.999999998 4.929332e-09 2.464666e-09
[163,] 0.999999999 2.749608e-09 1.374804e-09
[164,] 1.000000000 3.130401e-10 1.565201e-10
[165,] 1.000000000 4.964079e-10 2.482040e-10
[166,] 1.000000000 4.429616e-11 2.214808e-11
[167,] 1.000000000 8.235144e-11 4.117572e-11
[168,] 1.000000000 1.582494e-11 7.912472e-12
[169,] 1.000000000 2.408434e-11 1.204217e-11
[170,] 1.000000000 4.147902e-11 2.073951e-11
[171,] 1.000000000 1.774128e-11 8.870640e-12
[172,] 1.000000000 3.426284e-11 1.713142e-11
[173,] 1.000000000 5.840963e-11 2.920481e-11
[174,] 1.000000000 1.101317e-11 5.506584e-12
[175,] 1.000000000 2.025930e-11 1.012965e-11
[176,] 1.000000000 2.628185e-11 1.314092e-11
[177,] 1.000000000 3.360839e-11 1.680420e-11
[178,] 1.000000000 4.958218e-11 2.479109e-11
[179,] 1.000000000 6.742550e-11 3.371275e-11
[180,] 1.000000000 1.251071e-10 6.255356e-11
[181,] 1.000000000 2.372460e-10 1.186230e-10
[182,] 1.000000000 3.806129e-10 1.903064e-10
[183,] 1.000000000 6.937333e-10 3.468667e-10
[184,] 0.999999999 1.313558e-09 6.567789e-10
[185,] 0.999999999 1.754780e-09 8.773898e-10
[186,] 0.999999998 3.133475e-09 1.566738e-09
[187,] 1.000000000 3.122058e-10 1.561029e-10
[188,] 1.000000000 2.473805e-10 1.236903e-10
[189,] 1.000000000 2.872610e-10 1.436305e-10
[190,] 1.000000000 4.977831e-10 2.488915e-10
[191,] 1.000000000 5.251394e-10 2.625697e-10
[192,] 0.999999999 1.014451e-09 5.072257e-10
[193,] 0.999999999 1.629296e-09 8.146479e-10
[194,] 0.999999999 2.806081e-09 1.403041e-09
[195,] 0.999999997 5.217100e-09 2.608550e-09
[196,] 0.999999995 9.326431e-09 4.663215e-09
[197,] 0.999999991 1.760038e-08 8.800191e-09
[198,] 0.999999985 2.931322e-08 1.465661e-08
[199,] 0.999999987 2.576705e-08 1.288353e-08
[200,] 0.999999983 3.459434e-08 1.729717e-08
[201,] 0.999999968 6.447059e-08 3.223530e-08
[202,] 0.999999940 1.193767e-07 5.968833e-08
[203,] 0.999999895 2.098091e-07 1.049045e-07
[204,] 0.999999880 2.396521e-07 1.198260e-07
[205,] 0.999999866 2.683968e-07 1.341984e-07
[206,] 0.999999762 4.761564e-07 2.380782e-07
[207,] 0.999999726 5.488897e-07 2.744448e-07
[208,] 0.999999630 7.405597e-07 3.702799e-07
[209,] 0.999999422 1.155673e-06 5.778363e-07
[210,] 0.999999000 1.999631e-06 9.998155e-07
[211,] 0.999998467 3.065058e-06 1.532529e-06
[212,] 0.999997980 4.039400e-06 2.019700e-06
[213,] 0.999996456 7.087053e-06 3.543526e-06
[214,] 0.999995767 8.465410e-06 4.232705e-06
[215,] 0.999996677 6.645797e-06 3.322898e-06
[216,] 0.999995389 9.222882e-06 4.611441e-06
[217,] 0.999992260 1.547981e-05 7.739903e-06
[218,] 0.999990521 1.895881e-05 9.479404e-06
[219,] 0.999998999 2.001142e-06 1.000571e-06
[220,] 0.999998684 2.632614e-06 1.316307e-06
[221,] 0.999997553 4.894663e-06 2.447331e-06
[222,] 0.999996517 6.966669e-06 3.483334e-06
[223,] 0.999993652 1.269537e-05 6.347683e-06
[224,] 0.999989180 2.163951e-05 1.081975e-05
[225,] 0.999980912 3.817602e-05 1.908801e-05
[226,] 0.999966311 6.737810e-05 3.368905e-05
[227,] 0.999968837 6.232667e-05 3.116333e-05
[228,] 0.999952054 9.589244e-05 4.794622e-05
[229,] 0.999928695 1.426106e-04 7.130529e-05
[230,] 0.999890713 2.185732e-04 1.092866e-04
[231,] 0.999830922 3.381565e-04 1.690782e-04
[232,] 0.999808224 3.835522e-04 1.917761e-04
[233,] 0.999692817 6.143654e-04 3.071827e-04
[234,] 0.999482175 1.035651e-03 5.178255e-04
[235,] 0.999304606 1.390787e-03 6.953936e-04
[236,] 0.999013760 1.972481e-03 9.862404e-04
[237,] 0.998387188 3.225624e-03 1.612812e-03
[238,] 0.997811582 4.376837e-03 2.188418e-03
[239,] 0.999144049 1.711902e-03 8.559511e-04
[240,] 0.998568058 2.863885e-03 1.431942e-03
[241,] 0.997739714 4.520572e-03 2.260286e-03
[242,] 0.996358273 7.283454e-03 3.641727e-03
[243,] 0.994129414 1.174117e-02 5.870586e-03
[244,] 0.992856107 1.428779e-02 7.143893e-03
[245,] 0.989780084 2.043983e-02 1.021992e-02
[246,] 0.986473726 2.705255e-02 1.352627e-02
[247,] 0.979766522 4.046696e-02 2.023348e-02
[248,] 0.974534305 5.093139e-02 2.546569e-02
[249,] 0.965133332 6.973334e-02 3.486667e-02
[250,] 0.967235260 6.552948e-02 3.276474e-02
[251,] 0.950804342 9.839132e-02 4.919566e-02
[252,] 0.941933261 1.161335e-01 5.806674e-02
[253,] 0.915618749 1.687625e-01 8.438125e-02
[254,] 0.888690903 2.226182e-01 1.113091e-01
[255,] 0.844006963 3.119861e-01 1.559930e-01
[256,] 0.816470330 3.670593e-01 1.835297e-01
[257,] 0.810857298 3.782854e-01 1.891427e-01
[258,] 0.754453280 4.910934e-01 2.455467e-01
[259,] 0.682422743 6.351545e-01 3.175773e-01
[260,] 0.765388242 4.692235e-01 2.346118e-01
[261,] 0.676564394 6.468712e-01 3.234356e-01
[262,] 0.655566810 6.888664e-01 3.444332e-01
[263,] 0.914071325 1.718573e-01 8.592867e-02
[264,] 0.944295850 1.114083e-01 5.570415e-02
[265,] 0.890838547 2.183229e-01 1.091615e-01
[266,] 0.905851460 1.882971e-01 9.414854e-02
> postscript(file="/var/wessaorg/rcomp/tmp/1euqg1324145294.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/wessaorg/rcomp/tmp/2njje1324145294.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/wessaorg/rcomp/tmp/365ss1324145294.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/wessaorg/rcomp/tmp/4av0q1324145294.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/wessaorg/rcomp/tmp/5y8n21324145294.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 = 289
Frequency = 1
1 2 3 4 5 6
7775.7772 1420.0847 -13518.6087 2278.8402 -16459.5432 17933.8491
7 8 9 10 11 12
-240.8192 3457.8234 -3588.8140 23979.8658 -8650.3435 -30042.2078
13 14 15 16 17 18
-18730.7895 -8877.7754 -8714.8850 -30035.6479 -17978.4392 25431.4449
19 20 21 22 23 24
15193.8497 20527.9996 18264.7705 -7359.1211 -17127.7298 -13035.1420
25 26 27 28 29 30
9521.1825 -17616.5121 -1118.7795 -10591.7092 -21079.2198 -15996.1021
31 32 33 34 35 36
-17501.3966 -3542.7005 -4892.9376 -24020.6381 19041.7336 -23433.2243
37 38 39 40 41 42
-8959.7482 -101.8572 -23392.6932 -447.4581 -16993.4841 -16995.3554
43 44 45 46 47 48
-7977.3759 -22471.0041 -36905.7478 -17051.3628 -7167.5393 -10175.0672
49 50 51 52 53 54
-9217.9742 10404.1422 32336.4386 -16441.1940 1044.6059 20405.2553
55 56 57 58 59 60
-12569.0572 -14686.7899 -9760.5008 1478.5744 -5380.7072 -8760.1904
61 62 63 64 65 66
-6900.4229 17854.2787 30126.6432 -9249.9531 -10475.1836 -19749.5174
67 68 69 70 71 72
-12684.3672 3424.0765 -372.0612 -31366.8962 -9918.9255 -6983.9013
73 74 75 76 77 78
-24245.0321 29879.8610 -30808.2930 14903.8073 6040.9775 -27243.5644
79 80 81 82 83 84
-19550.6279 -17603.8757 4354.7794 -3463.7447 894.6354 -3002.5614
85 86 87 88 89 90
137064.0312 -6209.3563 -2723.2346 -10474.1071 8675.7866 -23240.0762
91 92 93 94 95 96
-6482.1840 18453.4946 -16921.2513 1922.4606 -11359.7933 12458.5582
97 98 99 100 101 102
-31926.3117 29841.1763 -24357.4463 11870.4100 -11598.9544 -12864.5674
103 104 105 106 107 108
4123.7824 34740.7519 29607.5292 36866.5858 -21074.7208 -7059.8949
109 110 111 112 113 114
9296.9503 -2235.8083 18752.1998 -16865.1354 2170.4463 -26928.5223
115 116 117 118 119 120
-20936.8408 -10842.1287 -1426.3803 -9537.9737 104807.4306 5606.7126
121 122 123 124 125 126
-11598.8216 -8390.3988 -13694.2456 88950.5901 19829.2014 -9253.4243
127 128 129 130 131 132
10919.9404 8149.0512 12210.5116 -7953.9351 -3122.9209 -13033.5359
133 134 135 136 137 138
-8050.4768 -8173.1341 -17763.4138 -27316.8568 -7555.9711 18173.1399
139 140 141 142 143 144
27580.1378 -7183.5975 -6863.1242 14261.6272 645.3913 -20249.5024
145 146 147 148 149 150
21612.0802 4231.2019 37311.2998 36133.5522 -8693.3373 16672.8377
151 152 153 154 155 156
-12678.2859 -7323.5627 3842.4926 -22698.8930 1311.3232 56027.6916
157 158 159 160 161 162
-4745.7783 -1871.1133 9905.5222 10846.1812 -4782.4715 63625.3697
163 164 165 166 167 168
7951.5718 28054.6084 -12181.5878 -23383.3754 -19891.5878 27529.8044
169 170 171 172 173 174
-15392.0247 31798.3711 6830.0927 4231.6693 6166.3621 -19136.0587
175 176 177 178 179 180
45030.3670 -11094.2591 33992.2758 -3340.4312 19460.0563 -12592.3428
181 182 183 184 185 186
8044.1058 28032.5465 -4725.1808 -9335.6821 -42081.2029 -6465.3613
187 188 189 190 191 192
-11435.3391 -7768.0095 -12288.3663 3980.8218 12990.6559 -1071.7764
193 194 195 196 197 198
-11300.6262 -2685.2076 419.8854 -8031.6142 -5676.8469 35206.4092
199 200 201 202 203 204
22351.2908 13163.4787 -10423.2119 -13805.0176 2724.7412 2681.1442
205 206 207 208 209 210
-9932.7629 705.9302 -6248.6783 -2035.8463 -6530.2779 -16573.1385
211 212 213 214 215 216
-16239.0376 2151.4416 -1395.1613 -7438.5538 9852.8583 -17232.1273
217 218 219 220 221 222
-865.0432 4229.0820 -15198.6789 -4665.0070 -3570.9939 -6337.5012
223 224 225 226 227 228
-14317.9593 -1022.1312 -8648.7023 -15617.4573 16163.3042 5236.7946
229 230 231 232 233 234
-8231.0444 35710.9558 -10785.8056 -4464.9053 14743.8378 -1454.2300
235 236 237 238 239 240
5413.6810 -1336.1648 3325.9676 -11673.9977 -9092.1587 -2100.7122
241 242 243 244 245 246
7324.3709 -6017.6006 -5964.2079 348.4075 1720.2725 7390.7695
247 248 249 250 251 252
8104.1205 3362.8387 17614.3349 14517.3573 2295.9042 6025.9996
253 254 255 256 257 258
-1689.3718 -284.6652 -8018.6510 -2711.5422 6969.1586 4182.7887
259 260 261 262 263 264
-10214.0525 -3787.0978 -4926.4181 -8983.7503 -3114.3484 -7350.5978
265 266 267 268 269 270
15430.3350 -3737.3091 -7661.1836 8005.4396 -11927.6581 14106.4902
271 272 273 274 275 276
-6023.7573 4358.3133 1445.8700 -3244.7077 917.5384 -6393.1730
277 278 279 280 281 282
17328.2852 -4038.4484 -13423.7245 -5568.0219 732.8970 63990.7185
283 284 285 286 287 288
8846.3634 -2282.9797 -1255.8275 18820.3326 -8735.5727 -4853.7493
289
-7875.5501
> postscript(file="/var/wessaorg/rcomp/tmp/6sxlr1324145294.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 = 289
Frequency = 1
lag(myerror, k = 1) myerror
0 7775.7772 NA
1 1420.0847 7775.7772
2 -13518.6087 1420.0847
3 2278.8402 -13518.6087
4 -16459.5432 2278.8402
5 17933.8491 -16459.5432
6 -240.8192 17933.8491
7 3457.8234 -240.8192
8 -3588.8140 3457.8234
9 23979.8658 -3588.8140
10 -8650.3435 23979.8658
11 -30042.2078 -8650.3435
12 -18730.7895 -30042.2078
13 -8877.7754 -18730.7895
14 -8714.8850 -8877.7754
15 -30035.6479 -8714.8850
16 -17978.4392 -30035.6479
17 25431.4449 -17978.4392
18 15193.8497 25431.4449
19 20527.9996 15193.8497
20 18264.7705 20527.9996
21 -7359.1211 18264.7705
22 -17127.7298 -7359.1211
23 -13035.1420 -17127.7298
24 9521.1825 -13035.1420
25 -17616.5121 9521.1825
26 -1118.7795 -17616.5121
27 -10591.7092 -1118.7795
28 -21079.2198 -10591.7092
29 -15996.1021 -21079.2198
30 -17501.3966 -15996.1021
31 -3542.7005 -17501.3966
32 -4892.9376 -3542.7005
33 -24020.6381 -4892.9376
34 19041.7336 -24020.6381
35 -23433.2243 19041.7336
36 -8959.7482 -23433.2243
37 -101.8572 -8959.7482
38 -23392.6932 -101.8572
39 -447.4581 -23392.6932
40 -16993.4841 -447.4581
41 -16995.3554 -16993.4841
42 -7977.3759 -16995.3554
43 -22471.0041 -7977.3759
44 -36905.7478 -22471.0041
45 -17051.3628 -36905.7478
46 -7167.5393 -17051.3628
47 -10175.0672 -7167.5393
48 -9217.9742 -10175.0672
49 10404.1422 -9217.9742
50 32336.4386 10404.1422
51 -16441.1940 32336.4386
52 1044.6059 -16441.1940
53 20405.2553 1044.6059
54 -12569.0572 20405.2553
55 -14686.7899 -12569.0572
56 -9760.5008 -14686.7899
57 1478.5744 -9760.5008
58 -5380.7072 1478.5744
59 -8760.1904 -5380.7072
60 -6900.4229 -8760.1904
61 17854.2787 -6900.4229
62 30126.6432 17854.2787
63 -9249.9531 30126.6432
64 -10475.1836 -9249.9531
65 -19749.5174 -10475.1836
66 -12684.3672 -19749.5174
67 3424.0765 -12684.3672
68 -372.0612 3424.0765
69 -31366.8962 -372.0612
70 -9918.9255 -31366.8962
71 -6983.9013 -9918.9255
72 -24245.0321 -6983.9013
73 29879.8610 -24245.0321
74 -30808.2930 29879.8610
75 14903.8073 -30808.2930
76 6040.9775 14903.8073
77 -27243.5644 6040.9775
78 -19550.6279 -27243.5644
79 -17603.8757 -19550.6279
80 4354.7794 -17603.8757
81 -3463.7447 4354.7794
82 894.6354 -3463.7447
83 -3002.5614 894.6354
84 137064.0312 -3002.5614
85 -6209.3563 137064.0312
86 -2723.2346 -6209.3563
87 -10474.1071 -2723.2346
88 8675.7866 -10474.1071
89 -23240.0762 8675.7866
90 -6482.1840 -23240.0762
91 18453.4946 -6482.1840
92 -16921.2513 18453.4946
93 1922.4606 -16921.2513
94 -11359.7933 1922.4606
95 12458.5582 -11359.7933
96 -31926.3117 12458.5582
97 29841.1763 -31926.3117
98 -24357.4463 29841.1763
99 11870.4100 -24357.4463
100 -11598.9544 11870.4100
101 -12864.5674 -11598.9544
102 4123.7824 -12864.5674
103 34740.7519 4123.7824
104 29607.5292 34740.7519
105 36866.5858 29607.5292
106 -21074.7208 36866.5858
107 -7059.8949 -21074.7208
108 9296.9503 -7059.8949
109 -2235.8083 9296.9503
110 18752.1998 -2235.8083
111 -16865.1354 18752.1998
112 2170.4463 -16865.1354
113 -26928.5223 2170.4463
114 -20936.8408 -26928.5223
115 -10842.1287 -20936.8408
116 -1426.3803 -10842.1287
117 -9537.9737 -1426.3803
118 104807.4306 -9537.9737
119 5606.7126 104807.4306
120 -11598.8216 5606.7126
121 -8390.3988 -11598.8216
122 -13694.2456 -8390.3988
123 88950.5901 -13694.2456
124 19829.2014 88950.5901
125 -9253.4243 19829.2014
126 10919.9404 -9253.4243
127 8149.0512 10919.9404
128 12210.5116 8149.0512
129 -7953.9351 12210.5116
130 -3122.9209 -7953.9351
131 -13033.5359 -3122.9209
132 -8050.4768 -13033.5359
133 -8173.1341 -8050.4768
134 -17763.4138 -8173.1341
135 -27316.8568 -17763.4138
136 -7555.9711 -27316.8568
137 18173.1399 -7555.9711
138 27580.1378 18173.1399
139 -7183.5975 27580.1378
140 -6863.1242 -7183.5975
141 14261.6272 -6863.1242
142 645.3913 14261.6272
143 -20249.5024 645.3913
144 21612.0802 -20249.5024
145 4231.2019 21612.0802
146 37311.2998 4231.2019
147 36133.5522 37311.2998
148 -8693.3373 36133.5522
149 16672.8377 -8693.3373
150 -12678.2859 16672.8377
151 -7323.5627 -12678.2859
152 3842.4926 -7323.5627
153 -22698.8930 3842.4926
154 1311.3232 -22698.8930
155 56027.6916 1311.3232
156 -4745.7783 56027.6916
157 -1871.1133 -4745.7783
158 9905.5222 -1871.1133
159 10846.1812 9905.5222
160 -4782.4715 10846.1812
161 63625.3697 -4782.4715
162 7951.5718 63625.3697
163 28054.6084 7951.5718
164 -12181.5878 28054.6084
165 -23383.3754 -12181.5878
166 -19891.5878 -23383.3754
167 27529.8044 -19891.5878
168 -15392.0247 27529.8044
169 31798.3711 -15392.0247
170 6830.0927 31798.3711
171 4231.6693 6830.0927
172 6166.3621 4231.6693
173 -19136.0587 6166.3621
174 45030.3670 -19136.0587
175 -11094.2591 45030.3670
176 33992.2758 -11094.2591
177 -3340.4312 33992.2758
178 19460.0563 -3340.4312
179 -12592.3428 19460.0563
180 8044.1058 -12592.3428
181 28032.5465 8044.1058
182 -4725.1808 28032.5465
183 -9335.6821 -4725.1808
184 -42081.2029 -9335.6821
185 -6465.3613 -42081.2029
186 -11435.3391 -6465.3613
187 -7768.0095 -11435.3391
188 -12288.3663 -7768.0095
189 3980.8218 -12288.3663
190 12990.6559 3980.8218
191 -1071.7764 12990.6559
192 -11300.6262 -1071.7764
193 -2685.2076 -11300.6262
194 419.8854 -2685.2076
195 -8031.6142 419.8854
196 -5676.8469 -8031.6142
197 35206.4092 -5676.8469
198 22351.2908 35206.4092
199 13163.4787 22351.2908
200 -10423.2119 13163.4787
201 -13805.0176 -10423.2119
202 2724.7412 -13805.0176
203 2681.1442 2724.7412
204 -9932.7629 2681.1442
205 705.9302 -9932.7629
206 -6248.6783 705.9302
207 -2035.8463 -6248.6783
208 -6530.2779 -2035.8463
209 -16573.1385 -6530.2779
210 -16239.0376 -16573.1385
211 2151.4416 -16239.0376
212 -1395.1613 2151.4416
213 -7438.5538 -1395.1613
214 9852.8583 -7438.5538
215 -17232.1273 9852.8583
216 -865.0432 -17232.1273
217 4229.0820 -865.0432
218 -15198.6789 4229.0820
219 -4665.0070 -15198.6789
220 -3570.9939 -4665.0070
221 -6337.5012 -3570.9939
222 -14317.9593 -6337.5012
223 -1022.1312 -14317.9593
224 -8648.7023 -1022.1312
225 -15617.4573 -8648.7023
226 16163.3042 -15617.4573
227 5236.7946 16163.3042
228 -8231.0444 5236.7946
229 35710.9558 -8231.0444
230 -10785.8056 35710.9558
231 -4464.9053 -10785.8056
232 14743.8378 -4464.9053
233 -1454.2300 14743.8378
234 5413.6810 -1454.2300
235 -1336.1648 5413.6810
236 3325.9676 -1336.1648
237 -11673.9977 3325.9676
238 -9092.1587 -11673.9977
239 -2100.7122 -9092.1587
240 7324.3709 -2100.7122
241 -6017.6006 7324.3709
242 -5964.2079 -6017.6006
243 348.4075 -5964.2079
244 1720.2725 348.4075
245 7390.7695 1720.2725
246 8104.1205 7390.7695
247 3362.8387 8104.1205
248 17614.3349 3362.8387
249 14517.3573 17614.3349
250 2295.9042 14517.3573
251 6025.9996 2295.9042
252 -1689.3718 6025.9996
253 -284.6652 -1689.3718
254 -8018.6510 -284.6652
255 -2711.5422 -8018.6510
256 6969.1586 -2711.5422
257 4182.7887 6969.1586
258 -10214.0525 4182.7887
259 -3787.0978 -10214.0525
260 -4926.4181 -3787.0978
261 -8983.7503 -4926.4181
262 -3114.3484 -8983.7503
263 -7350.5978 -3114.3484
264 15430.3350 -7350.5978
265 -3737.3091 15430.3350
266 -7661.1836 -3737.3091
267 8005.4396 -7661.1836
268 -11927.6581 8005.4396
269 14106.4902 -11927.6581
270 -6023.7573 14106.4902
271 4358.3133 -6023.7573
272 1445.8700 4358.3133
273 -3244.7077 1445.8700
274 917.5384 -3244.7077
275 -6393.1730 917.5384
276 17328.2852 -6393.1730
277 -4038.4484 17328.2852
278 -13423.7245 -4038.4484
279 -5568.0219 -13423.7245
280 732.8970 -5568.0219
281 63990.7185 732.8970
282 8846.3634 63990.7185
283 -2282.9797 8846.3634
284 -1255.8275 -2282.9797
285 18820.3326 -1255.8275
286 -8735.5727 18820.3326
287 -4853.7493 -8735.5727
288 -7875.5501 -4853.7493
289 NA -7875.5501
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 1420.0847 7775.7772
[2,] -13518.6087 1420.0847
[3,] 2278.8402 -13518.6087
[4,] -16459.5432 2278.8402
[5,] 17933.8491 -16459.5432
[6,] -240.8192 17933.8491
[7,] 3457.8234 -240.8192
[8,] -3588.8140 3457.8234
[9,] 23979.8658 -3588.8140
[10,] -8650.3435 23979.8658
[11,] -30042.2078 -8650.3435
[12,] -18730.7895 -30042.2078
[13,] -8877.7754 -18730.7895
[14,] -8714.8850 -8877.7754
[15,] -30035.6479 -8714.8850
[16,] -17978.4392 -30035.6479
[17,] 25431.4449 -17978.4392
[18,] 15193.8497 25431.4449
[19,] 20527.9996 15193.8497
[20,] 18264.7705 20527.9996
[21,] -7359.1211 18264.7705
[22,] -17127.7298 -7359.1211
[23,] -13035.1420 -17127.7298
[24,] 9521.1825 -13035.1420
[25,] -17616.5121 9521.1825
[26,] -1118.7795 -17616.5121
[27,] -10591.7092 -1118.7795
[28,] -21079.2198 -10591.7092
[29,] -15996.1021 -21079.2198
[30,] -17501.3966 -15996.1021
[31,] -3542.7005 -17501.3966
[32,] -4892.9376 -3542.7005
[33,] -24020.6381 -4892.9376
[34,] 19041.7336 -24020.6381
[35,] -23433.2243 19041.7336
[36,] -8959.7482 -23433.2243
[37,] -101.8572 -8959.7482
[38,] -23392.6932 -101.8572
[39,] -447.4581 -23392.6932
[40,] -16993.4841 -447.4581
[41,] -16995.3554 -16993.4841
[42,] -7977.3759 -16995.3554
[43,] -22471.0041 -7977.3759
[44,] -36905.7478 -22471.0041
[45,] -17051.3628 -36905.7478
[46,] -7167.5393 -17051.3628
[47,] -10175.0672 -7167.5393
[48,] -9217.9742 -10175.0672
[49,] 10404.1422 -9217.9742
[50,] 32336.4386 10404.1422
[51,] -16441.1940 32336.4386
[52,] 1044.6059 -16441.1940
[53,] 20405.2553 1044.6059
[54,] -12569.0572 20405.2553
[55,] -14686.7899 -12569.0572
[56,] -9760.5008 -14686.7899
[57,] 1478.5744 -9760.5008
[58,] -5380.7072 1478.5744
[59,] -8760.1904 -5380.7072
[60,] -6900.4229 -8760.1904
[61,] 17854.2787 -6900.4229
[62,] 30126.6432 17854.2787
[63,] -9249.9531 30126.6432
[64,] -10475.1836 -9249.9531
[65,] -19749.5174 -10475.1836
[66,] -12684.3672 -19749.5174
[67,] 3424.0765 -12684.3672
[68,] -372.0612 3424.0765
[69,] -31366.8962 -372.0612
[70,] -9918.9255 -31366.8962
[71,] -6983.9013 -9918.9255
[72,] -24245.0321 -6983.9013
[73,] 29879.8610 -24245.0321
[74,] -30808.2930 29879.8610
[75,] 14903.8073 -30808.2930
[76,] 6040.9775 14903.8073
[77,] -27243.5644 6040.9775
[78,] -19550.6279 -27243.5644
[79,] -17603.8757 -19550.6279
[80,] 4354.7794 -17603.8757
[81,] -3463.7447 4354.7794
[82,] 894.6354 -3463.7447
[83,] -3002.5614 894.6354
[84,] 137064.0312 -3002.5614
[85,] -6209.3563 137064.0312
[86,] -2723.2346 -6209.3563
[87,] -10474.1071 -2723.2346
[88,] 8675.7866 -10474.1071
[89,] -23240.0762 8675.7866
[90,] -6482.1840 -23240.0762
[91,] 18453.4946 -6482.1840
[92,] -16921.2513 18453.4946
[93,] 1922.4606 -16921.2513
[94,] -11359.7933 1922.4606
[95,] 12458.5582 -11359.7933
[96,] -31926.3117 12458.5582
[97,] 29841.1763 -31926.3117
[98,] -24357.4463 29841.1763
[99,] 11870.4100 -24357.4463
[100,] -11598.9544 11870.4100
[101,] -12864.5674 -11598.9544
[102,] 4123.7824 -12864.5674
[103,] 34740.7519 4123.7824
[104,] 29607.5292 34740.7519
[105,] 36866.5858 29607.5292
[106,] -21074.7208 36866.5858
[107,] -7059.8949 -21074.7208
[108,] 9296.9503 -7059.8949
[109,] -2235.8083 9296.9503
[110,] 18752.1998 -2235.8083
[111,] -16865.1354 18752.1998
[112,] 2170.4463 -16865.1354
[113,] -26928.5223 2170.4463
[114,] -20936.8408 -26928.5223
[115,] -10842.1287 -20936.8408
[116,] -1426.3803 -10842.1287
[117,] -9537.9737 -1426.3803
[118,] 104807.4306 -9537.9737
[119,] 5606.7126 104807.4306
[120,] -11598.8216 5606.7126
[121,] -8390.3988 -11598.8216
[122,] -13694.2456 -8390.3988
[123,] 88950.5901 -13694.2456
[124,] 19829.2014 88950.5901
[125,] -9253.4243 19829.2014
[126,] 10919.9404 -9253.4243
[127,] 8149.0512 10919.9404
[128,] 12210.5116 8149.0512
[129,] -7953.9351 12210.5116
[130,] -3122.9209 -7953.9351
[131,] -13033.5359 -3122.9209
[132,] -8050.4768 -13033.5359
[133,] -8173.1341 -8050.4768
[134,] -17763.4138 -8173.1341
[135,] -27316.8568 -17763.4138
[136,] -7555.9711 -27316.8568
[137,] 18173.1399 -7555.9711
[138,] 27580.1378 18173.1399
[139,] -7183.5975 27580.1378
[140,] -6863.1242 -7183.5975
[141,] 14261.6272 -6863.1242
[142,] 645.3913 14261.6272
[143,] -20249.5024 645.3913
[144,] 21612.0802 -20249.5024
[145,] 4231.2019 21612.0802
[146,] 37311.2998 4231.2019
[147,] 36133.5522 37311.2998
[148,] -8693.3373 36133.5522
[149,] 16672.8377 -8693.3373
[150,] -12678.2859 16672.8377
[151,] -7323.5627 -12678.2859
[152,] 3842.4926 -7323.5627
[153,] -22698.8930 3842.4926
[154,] 1311.3232 -22698.8930
[155,] 56027.6916 1311.3232
[156,] -4745.7783 56027.6916
[157,] -1871.1133 -4745.7783
[158,] 9905.5222 -1871.1133
[159,] 10846.1812 9905.5222
[160,] -4782.4715 10846.1812
[161,] 63625.3697 -4782.4715
[162,] 7951.5718 63625.3697
[163,] 28054.6084 7951.5718
[164,] -12181.5878 28054.6084
[165,] -23383.3754 -12181.5878
[166,] -19891.5878 -23383.3754
[167,] 27529.8044 -19891.5878
[168,] -15392.0247 27529.8044
[169,] 31798.3711 -15392.0247
[170,] 6830.0927 31798.3711
[171,] 4231.6693 6830.0927
[172,] 6166.3621 4231.6693
[173,] -19136.0587 6166.3621
[174,] 45030.3670 -19136.0587
[175,] -11094.2591 45030.3670
[176,] 33992.2758 -11094.2591
[177,] -3340.4312 33992.2758
[178,] 19460.0563 -3340.4312
[179,] -12592.3428 19460.0563
[180,] 8044.1058 -12592.3428
[181,] 28032.5465 8044.1058
[182,] -4725.1808 28032.5465
[183,] -9335.6821 -4725.1808
[184,] -42081.2029 -9335.6821
[185,] -6465.3613 -42081.2029
[186,] -11435.3391 -6465.3613
[187,] -7768.0095 -11435.3391
[188,] -12288.3663 -7768.0095
[189,] 3980.8218 -12288.3663
[190,] 12990.6559 3980.8218
[191,] -1071.7764 12990.6559
[192,] -11300.6262 -1071.7764
[193,] -2685.2076 -11300.6262
[194,] 419.8854 -2685.2076
[195,] -8031.6142 419.8854
[196,] -5676.8469 -8031.6142
[197,] 35206.4092 -5676.8469
[198,] 22351.2908 35206.4092
[199,] 13163.4787 22351.2908
[200,] -10423.2119 13163.4787
[201,] -13805.0176 -10423.2119
[202,] 2724.7412 -13805.0176
[203,] 2681.1442 2724.7412
[204,] -9932.7629 2681.1442
[205,] 705.9302 -9932.7629
[206,] -6248.6783 705.9302
[207,] -2035.8463 -6248.6783
[208,] -6530.2779 -2035.8463
[209,] -16573.1385 -6530.2779
[210,] -16239.0376 -16573.1385
[211,] 2151.4416 -16239.0376
[212,] -1395.1613 2151.4416
[213,] -7438.5538 -1395.1613
[214,] 9852.8583 -7438.5538
[215,] -17232.1273 9852.8583
[216,] -865.0432 -17232.1273
[217,] 4229.0820 -865.0432
[218,] -15198.6789 4229.0820
[219,] -4665.0070 -15198.6789
[220,] -3570.9939 -4665.0070
[221,] -6337.5012 -3570.9939
[222,] -14317.9593 -6337.5012
[223,] -1022.1312 -14317.9593
[224,] -8648.7023 -1022.1312
[225,] -15617.4573 -8648.7023
[226,] 16163.3042 -15617.4573
[227,] 5236.7946 16163.3042
[228,] -8231.0444 5236.7946
[229,] 35710.9558 -8231.0444
[230,] -10785.8056 35710.9558
[231,] -4464.9053 -10785.8056
[232,] 14743.8378 -4464.9053
[233,] -1454.2300 14743.8378
[234,] 5413.6810 -1454.2300
[235,] -1336.1648 5413.6810
[236,] 3325.9676 -1336.1648
[237,] -11673.9977 3325.9676
[238,] -9092.1587 -11673.9977
[239,] -2100.7122 -9092.1587
[240,] 7324.3709 -2100.7122
[241,] -6017.6006 7324.3709
[242,] -5964.2079 -6017.6006
[243,] 348.4075 -5964.2079
[244,] 1720.2725 348.4075
[245,] 7390.7695 1720.2725
[246,] 8104.1205 7390.7695
[247,] 3362.8387 8104.1205
[248,] 17614.3349 3362.8387
[249,] 14517.3573 17614.3349
[250,] 2295.9042 14517.3573
[251,] 6025.9996 2295.9042
[252,] -1689.3718 6025.9996
[253,] -284.6652 -1689.3718
[254,] -8018.6510 -284.6652
[255,] -2711.5422 -8018.6510
[256,] 6969.1586 -2711.5422
[257,] 4182.7887 6969.1586
[258,] -10214.0525 4182.7887
[259,] -3787.0978 -10214.0525
[260,] -4926.4181 -3787.0978
[261,] -8983.7503 -4926.4181
[262,] -3114.3484 -8983.7503
[263,] -7350.5978 -3114.3484
[264,] 15430.3350 -7350.5978
[265,] -3737.3091 15430.3350
[266,] -7661.1836 -3737.3091
[267,] 8005.4396 -7661.1836
[268,] -11927.6581 8005.4396
[269,] 14106.4902 -11927.6581
[270,] -6023.7573 14106.4902
[271,] 4358.3133 -6023.7573
[272,] 1445.8700 4358.3133
[273,] -3244.7077 1445.8700
[274,] 917.5384 -3244.7077
[275,] -6393.1730 917.5384
[276,] 17328.2852 -6393.1730
[277,] -4038.4484 17328.2852
[278,] -13423.7245 -4038.4484
[279,] -5568.0219 -13423.7245
[280,] 732.8970 -5568.0219
[281,] 63990.7185 732.8970
[282,] 8846.3634 63990.7185
[283,] -2282.9797 8846.3634
[284,] -1255.8275 -2282.9797
[285,] 18820.3326 -1255.8275
[286,] -8735.5727 18820.3326
[287,] -4853.7493 -8735.5727
[288,] -7875.5501 -4853.7493
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 1420.0847 7775.7772
2 -13518.6087 1420.0847
3 2278.8402 -13518.6087
4 -16459.5432 2278.8402
5 17933.8491 -16459.5432
6 -240.8192 17933.8491
7 3457.8234 -240.8192
8 -3588.8140 3457.8234
9 23979.8658 -3588.8140
10 -8650.3435 23979.8658
11 -30042.2078 -8650.3435
12 -18730.7895 -30042.2078
13 -8877.7754 -18730.7895
14 -8714.8850 -8877.7754
15 -30035.6479 -8714.8850
16 -17978.4392 -30035.6479
17 25431.4449 -17978.4392
18 15193.8497 25431.4449
19 20527.9996 15193.8497
20 18264.7705 20527.9996
21 -7359.1211 18264.7705
22 -17127.7298 -7359.1211
23 -13035.1420 -17127.7298
24 9521.1825 -13035.1420
25 -17616.5121 9521.1825
26 -1118.7795 -17616.5121
27 -10591.7092 -1118.7795
28 -21079.2198 -10591.7092
29 -15996.1021 -21079.2198
30 -17501.3966 -15996.1021
31 -3542.7005 -17501.3966
32 -4892.9376 -3542.7005
33 -24020.6381 -4892.9376
34 19041.7336 -24020.6381
35 -23433.2243 19041.7336
36 -8959.7482 -23433.2243
37 -101.8572 -8959.7482
38 -23392.6932 -101.8572
39 -447.4581 -23392.6932
40 -16993.4841 -447.4581
41 -16995.3554 -16993.4841
42 -7977.3759 -16995.3554
43 -22471.0041 -7977.3759
44 -36905.7478 -22471.0041
45 -17051.3628 -36905.7478
46 -7167.5393 -17051.3628
47 -10175.0672 -7167.5393
48 -9217.9742 -10175.0672
49 10404.1422 -9217.9742
50 32336.4386 10404.1422
51 -16441.1940 32336.4386
52 1044.6059 -16441.1940
53 20405.2553 1044.6059
54 -12569.0572 20405.2553
55 -14686.7899 -12569.0572
56 -9760.5008 -14686.7899
57 1478.5744 -9760.5008
58 -5380.7072 1478.5744
59 -8760.1904 -5380.7072
60 -6900.4229 -8760.1904
61 17854.2787 -6900.4229
62 30126.6432 17854.2787
63 -9249.9531 30126.6432
64 -10475.1836 -9249.9531
65 -19749.5174 -10475.1836
66 -12684.3672 -19749.5174
67 3424.0765 -12684.3672
68 -372.0612 3424.0765
69 -31366.8962 -372.0612
70 -9918.9255 -31366.8962
71 -6983.9013 -9918.9255
72 -24245.0321 -6983.9013
73 29879.8610 -24245.0321
74 -30808.2930 29879.8610
75 14903.8073 -30808.2930
76 6040.9775 14903.8073
77 -27243.5644 6040.9775
78 -19550.6279 -27243.5644
79 -17603.8757 -19550.6279
80 4354.7794 -17603.8757
81 -3463.7447 4354.7794
82 894.6354 -3463.7447
83 -3002.5614 894.6354
84 137064.0312 -3002.5614
85 -6209.3563 137064.0312
86 -2723.2346 -6209.3563
87 -10474.1071 -2723.2346
88 8675.7866 -10474.1071
89 -23240.0762 8675.7866
90 -6482.1840 -23240.0762
91 18453.4946 -6482.1840
92 -16921.2513 18453.4946
93 1922.4606 -16921.2513
94 -11359.7933 1922.4606
95 12458.5582 -11359.7933
96 -31926.3117 12458.5582
97 29841.1763 -31926.3117
98 -24357.4463 29841.1763
99 11870.4100 -24357.4463
100 -11598.9544 11870.4100
101 -12864.5674 -11598.9544
102 4123.7824 -12864.5674
103 34740.7519 4123.7824
104 29607.5292 34740.7519
105 36866.5858 29607.5292
106 -21074.7208 36866.5858
107 -7059.8949 -21074.7208
108 9296.9503 -7059.8949
109 -2235.8083 9296.9503
110 18752.1998 -2235.8083
111 -16865.1354 18752.1998
112 2170.4463 -16865.1354
113 -26928.5223 2170.4463
114 -20936.8408 -26928.5223
115 -10842.1287 -20936.8408
116 -1426.3803 -10842.1287
117 -9537.9737 -1426.3803
118 104807.4306 -9537.9737
119 5606.7126 104807.4306
120 -11598.8216 5606.7126
121 -8390.3988 -11598.8216
122 -13694.2456 -8390.3988
123 88950.5901 -13694.2456
124 19829.2014 88950.5901
125 -9253.4243 19829.2014
126 10919.9404 -9253.4243
127 8149.0512 10919.9404
128 12210.5116 8149.0512
129 -7953.9351 12210.5116
130 -3122.9209 -7953.9351
131 -13033.5359 -3122.9209
132 -8050.4768 -13033.5359
133 -8173.1341 -8050.4768
134 -17763.4138 -8173.1341
135 -27316.8568 -17763.4138
136 -7555.9711 -27316.8568
137 18173.1399 -7555.9711
138 27580.1378 18173.1399
139 -7183.5975 27580.1378
140 -6863.1242 -7183.5975
141 14261.6272 -6863.1242
142 645.3913 14261.6272
143 -20249.5024 645.3913
144 21612.0802 -20249.5024
145 4231.2019 21612.0802
146 37311.2998 4231.2019
147 36133.5522 37311.2998
148 -8693.3373 36133.5522
149 16672.8377 -8693.3373
150 -12678.2859 16672.8377
151 -7323.5627 -12678.2859
152 3842.4926 -7323.5627
153 -22698.8930 3842.4926
154 1311.3232 -22698.8930
155 56027.6916 1311.3232
156 -4745.7783 56027.6916
157 -1871.1133 -4745.7783
158 9905.5222 -1871.1133
159 10846.1812 9905.5222
160 -4782.4715 10846.1812
161 63625.3697 -4782.4715
162 7951.5718 63625.3697
163 28054.6084 7951.5718
164 -12181.5878 28054.6084
165 -23383.3754 -12181.5878
166 -19891.5878 -23383.3754
167 27529.8044 -19891.5878
168 -15392.0247 27529.8044
169 31798.3711 -15392.0247
170 6830.0927 31798.3711
171 4231.6693 6830.0927
172 6166.3621 4231.6693
173 -19136.0587 6166.3621
174 45030.3670 -19136.0587
175 -11094.2591 45030.3670
176 33992.2758 -11094.2591
177 -3340.4312 33992.2758
178 19460.0563 -3340.4312
179 -12592.3428 19460.0563
180 8044.1058 -12592.3428
181 28032.5465 8044.1058
182 -4725.1808 28032.5465
183 -9335.6821 -4725.1808
184 -42081.2029 -9335.6821
185 -6465.3613 -42081.2029
186 -11435.3391 -6465.3613
187 -7768.0095 -11435.3391
188 -12288.3663 -7768.0095
189 3980.8218 -12288.3663
190 12990.6559 3980.8218
191 -1071.7764 12990.6559
192 -11300.6262 -1071.7764
193 -2685.2076 -11300.6262
194 419.8854 -2685.2076
195 -8031.6142 419.8854
196 -5676.8469 -8031.6142
197 35206.4092 -5676.8469
198 22351.2908 35206.4092
199 13163.4787 22351.2908
200 -10423.2119 13163.4787
201 -13805.0176 -10423.2119
202 2724.7412 -13805.0176
203 2681.1442 2724.7412
204 -9932.7629 2681.1442
205 705.9302 -9932.7629
206 -6248.6783 705.9302
207 -2035.8463 -6248.6783
208 -6530.2779 -2035.8463
209 -16573.1385 -6530.2779
210 -16239.0376 -16573.1385
211 2151.4416 -16239.0376
212 -1395.1613 2151.4416
213 -7438.5538 -1395.1613
214 9852.8583 -7438.5538
215 -17232.1273 9852.8583
216 -865.0432 -17232.1273
217 4229.0820 -865.0432
218 -15198.6789 4229.0820
219 -4665.0070 -15198.6789
220 -3570.9939 -4665.0070
221 -6337.5012 -3570.9939
222 -14317.9593 -6337.5012
223 -1022.1312 -14317.9593
224 -8648.7023 -1022.1312
225 -15617.4573 -8648.7023
226 16163.3042 -15617.4573
227 5236.7946 16163.3042
228 -8231.0444 5236.7946
229 35710.9558 -8231.0444
230 -10785.8056 35710.9558
231 -4464.9053 -10785.8056
232 14743.8378 -4464.9053
233 -1454.2300 14743.8378
234 5413.6810 -1454.2300
235 -1336.1648 5413.6810
236 3325.9676 -1336.1648
237 -11673.9977 3325.9676
238 -9092.1587 -11673.9977
239 -2100.7122 -9092.1587
240 7324.3709 -2100.7122
241 -6017.6006 7324.3709
242 -5964.2079 -6017.6006
243 348.4075 -5964.2079
244 1720.2725 348.4075
245 7390.7695 1720.2725
246 8104.1205 7390.7695
247 3362.8387 8104.1205
248 17614.3349 3362.8387
249 14517.3573 17614.3349
250 2295.9042 14517.3573
251 6025.9996 2295.9042
252 -1689.3718 6025.9996
253 -284.6652 -1689.3718
254 -8018.6510 -284.6652
255 -2711.5422 -8018.6510
256 6969.1586 -2711.5422
257 4182.7887 6969.1586
258 -10214.0525 4182.7887
259 -3787.0978 -10214.0525
260 -4926.4181 -3787.0978
261 -8983.7503 -4926.4181
262 -3114.3484 -8983.7503
263 -7350.5978 -3114.3484
264 15430.3350 -7350.5978
265 -3737.3091 15430.3350
266 -7661.1836 -3737.3091
267 8005.4396 -7661.1836
268 -11927.6581 8005.4396
269 14106.4902 -11927.6581
270 -6023.7573 14106.4902
271 4358.3133 -6023.7573
272 1445.8700 4358.3133
273 -3244.7077 1445.8700
274 917.5384 -3244.7077
275 -6393.1730 917.5384
276 17328.2852 -6393.1730
277 -4038.4484 17328.2852
278 -13423.7245 -4038.4484
279 -5568.0219 -13423.7245
280 732.8970 -5568.0219
281 63990.7185 732.8970
282 8846.3634 63990.7185
283 -2282.9797 8846.3634
284 -1255.8275 -2282.9797
285 18820.3326 -1255.8275
286 -8735.5727 18820.3326
287 -4853.7493 -8735.5727
288 -7875.5501 -4853.7493
> 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/wessaorg/rcomp/tmp/72zji1324145294.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/wessaorg/rcomp/tmp/8nhn11324145294.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/wessaorg/rcomp/tmp/9beiz1324145294.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/wessaorg/rcomp/tmp/10x1d31324145294.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/1166wh1324145294.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/wessaorg/rcomp/tmp/121x801324145294.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/wessaorg/rcomp/tmp/135v2y1324145294.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/wessaorg/rcomp/tmp/14tnpv1324145294.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/wessaorg/rcomp/tmp/15n9ea1324145294.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/wessaorg/rcomp/tmp/16q2c21324145294.tab")
+ }
>
> try(system("convert tmp/1euqg1324145294.ps tmp/1euqg1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/2njje1324145294.ps tmp/2njje1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/365ss1324145294.ps tmp/365ss1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/4av0q1324145294.ps tmp/4av0q1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/5y8n21324145294.ps tmp/5y8n21324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sxlr1324145294.ps tmp/6sxlr1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/72zji1324145294.ps tmp/72zji1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/8nhn11324145294.ps tmp/8nhn11324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/9beiz1324145294.ps tmp/9beiz1324145294.png",intern=TRUE))
character(0)
> try(system("convert tmp/10x1d31324145294.ps tmp/10x1d31324145294.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
9.270 0.695 9.978