R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> x <- array(list(58.58527778
+ ,145
+ ,30
+ ,94
+ ,112285
+ ,14.36009776
+ ,1
+ ,33.60611111
+ ,101
+ ,28
+ ,103
+ ,84786
+ ,11.30900293
+ ,1
+ ,49.03
+ ,98
+ ,38
+ ,93
+ ,83123
+ ,12.45810893
+ ,1
+ ,49.81138889
+ ,132
+ ,30
+ ,103
+ ,101193
+ ,13.50394283
+ ,1
+ ,34.21805556
+ ,60
+ ,22
+ ,51
+ ,38361
+ ,7.219358341
+ ,1
+ ,14.65166667
+ ,38
+ ,26
+ ,70
+ ,68504
+ ,7.381367182
+ ,1
+ ,107.0927778
+ ,144
+ ,25
+ ,91
+ ,119182
+ ,15.26096233
+ ,1
+ ,9.213888889
+ ,5
+ ,18
+ ,22
+ ,22807
+ ,3.165415029
+ ,1
+ ,41.40583333
+ ,84
+ ,26
+ ,93
+ ,116174
+ ,11.78143911
+ ,1
+ ,45.95722222
+ ,79
+ ,25
+ ,60
+ ,57635
+ ,9.277808364
+ ,1
+ ,65.8925
+ ,127
+ ,38
+ ,123
+ ,66198
+ ,14.22780485
+ ,1
+ ,48.14611111
+ ,78
+ ,44
+ ,148
+ ,71701
+ ,13.3753608
+ ,1
+ ,36.98083333
+ ,60
+ ,30
+ ,90
+ ,57793
+ ,9.487750499
+ ,1
+ ,71.90916667
+ ,131
+ ,40
+ ,124
+ ,80444
+ ,15.21511671
+ ,1
+ ,50.02305556
+ ,84
+ ,34
+ ,70
+ ,53855
+ ,10.42712012
+ ,1
+ ,90.22194444
+ ,133
+ ,47
+ ,168
+ ,97668
+ ,17.84978648
+ ,1
+ ,64.15666667
+ ,150
+ ,30
+ ,115
+ ,133824
+ ,15.89534882
+ ,1
+ ,65.77361111
+ ,91
+ ,31
+ ,71
+ ,101481
+ ,12.55115438
+ ,1
+ ,37.63138889
+ ,132
+ ,23
+ ,66
+ ,99645
+ ,11.43396085
+ ,1
+ ,56.36805556
+ ,136
+ ,36
+ ,134
+ ,114789
+ ,15.53507626
+ ,1
+ ,59.76305556
+ ,124
+ ,36
+ ,117
+ ,99052
+ ,14.52239325
+ ,1
+ ,95.63805556
+ ,118
+ ,30
+ ,108
+ ,67654
+ ,13.89331579
+ ,1
+ ,42.75972222
+ ,70
+ ,25
+ ,84
+ ,65553
+ ,9.712278591
+ ,1
+ ,36.92861111
+ ,107
+ ,39
+ ,156
+ ,97500
+ ,14.12942504
+ ,1
+ ,48.53444444
+ ,119
+ ,34
+ ,120
+ ,69112
+ ,12.89537363
+ ,1
+ ,48.44861111
+ ,89
+ ,31
+ ,114
+ ,82753
+ ,12.1602659
+ ,1
+ ,62.65222222
+ ,112
+ ,31
+ ,94
+ ,85323
+ ,13.01181689
+ ,1
+ ,62.12
+ ,108
+ ,33
+ ,120
+ ,72654
+ ,13.29873685
+ ,1
+ ,34.67138889
+ ,52
+ ,25
+ ,81
+ ,30727
+ ,7.775415068
+ ,1
+ ,61.58277778
+ ,112
+ ,33
+ ,110
+ ,77873
+ ,13.28429728
+ ,1
+ ,58.54638889
+ ,116
+ ,35
+ ,133
+ ,117478
+ ,15.11451039
+ ,1
+ ,47.29611111
+ ,123
+ ,42
+ ,122
+ ,74007
+ ,13.74797244
+ ,1
+ ,72.37805556
+ ,125
+ ,43
+ ,158
+ ,90183
+ ,16.42942616
+ ,1
+ ,23.57027778
+ ,27
+ ,30
+ ,109
+ ,61542
+ ,8.584377305
+ ,1
+ ,81.78444444
+ ,162
+ ,33
+ ,124
+ ,101494
+ ,16.47410685
+ ,1
+ ,59.90027778
+ ,64
+ ,32
+ ,92
+ ,55813
+ ,10.83539947
+ ,1
+ ,90.3075
+ ,92
+ ,36
+ ,126
+ ,79215
+ ,14.48035784
+ ,1
+ ,46.53944444
+ ,83
+ ,28
+ ,70
+ ,55461
+ ,9.817153316
+ ,1
+ ,29.55777778
+ ,41
+ ,14
+ ,37
+ ,31081
+ ,5.345232896
+ ,1
+ ,73.82472222
+ ,120
+ ,32
+ ,120
+ ,83122
+ ,14.39310875
+ ,1
+ ,74.90305556
+ ,105
+ ,30
+ ,93
+ ,70106
+ ,12.87630593
+ ,1
+ ,41.42
+ ,79
+ ,35
+ ,95
+ ,60578
+ ,10.76877965
+ ,1
+ ,42.46416667
+ ,70
+ ,28
+ ,90
+ ,79892
+ ,10.49694082
+ ,1
+ ,31.01805556
+ ,55
+ ,28
+ ,80
+ ,49810
+ ,8.443671407
+ ,1
+ ,32.33555556
+ ,39
+ ,39
+ ,31
+ ,71570
+ ,8.407782836
+ ,1
+ ,100.6391667
+ ,67
+ ,34
+ ,110
+ ,100708
+ ,13.94065255
+ ,1
+ ,21.88888889
+ ,21
+ ,26
+ ,66
+ ,33032
+ ,6.16326461
+ ,1
+ ,50.87972222
+ ,127
+ ,39
+ ,138
+ ,82875
+ ,14.43688293
+ ,1
+ ,77.2125
+ ,152
+ ,39
+ ,133
+ ,139077
+ ,17.77583855
+ ,1
+ ,41.84138889
+ ,113
+ ,33
+ ,113
+ ,71595
+ ,12.2486023
+ ,1
+ ,46.89138889
+ ,99
+ ,28
+ ,100
+ ,72260
+ ,11.45685566
+ ,1
+ ,6.718888889
+ ,7
+ ,4
+ ,7
+ ,5950
+ ,1.153093302
+ ,1
+ ,91.46305556
+ ,141
+ ,39
+ ,140
+ ,115762
+ ,17.35686844
+ ,1
+ ,18.06361111
+ ,21
+ ,18
+ ,61
+ ,32551
+ ,5.224652076
+ ,1
+ ,28.0825
+ ,35
+ ,14
+ ,41
+ ,31701
+ ,5.238696901
+ ,1
+ ,60.81833333
+ ,109
+ ,29
+ ,96
+ ,80670
+ ,12.60523152
+ ,1
+ ,67.79222222
+ ,133
+ ,44
+ ,164
+ ,143558
+ ,17.98404039
+ ,1
+ ,64.81333333
+ ,230
+ ,28
+ ,102
+ ,120733
+ ,15.30279334
+ ,1
+ ,71.23944444
+ ,166
+ ,35
+ ,124
+ ,105195
+ ,16.2306506
+ ,1
+ ,57.26694444
+ ,68
+ ,28
+ ,99
+ ,73107
+ ,11.17728596
+ ,1
+ ,86.52027778
+ ,147
+ ,38
+ ,129
+ ,132068
+ ,17.64250157
+ ,1
+ ,65.5
+ ,179
+ ,23
+ ,62
+ ,149193
+ ,14.43288916
+ ,1
+ ,49.4275
+ ,61
+ ,36
+ ,73
+ ,46821
+ ,9.841441509
+ ,1
+ ,57.54888889
+ ,101
+ ,32
+ ,114
+ ,87011
+ ,13.10185804
+ ,1
+ ,54.59805556
+ ,108
+ ,29
+ ,99
+ ,95260
+ ,12.78229154
+ ,1
+ ,48.38444444
+ ,90
+ ,25
+ ,70
+ ,55183
+ ,9.841353028
+ ,1
+ ,39.79055556
+ ,114
+ ,27
+ ,104
+ ,106671
+ ,12.52267132
+ ,1
+ ,52.09972222
+ ,103
+ ,36
+ ,116
+ ,73511
+ ,12.85177133
+ ,1
+ ,52.13361111
+ ,142
+ ,28
+ ,91
+ ,92945
+ ,13.17889823
+ ,1
+ ,33.06
+ ,79
+ ,23
+ ,74
+ ,78664
+ ,9.459926349
+ ,1
+ ,50.60888889
+ ,88
+ ,40
+ ,138
+ ,70054
+ ,13.14460233
+ ,1
+ ,20.435
+ ,25
+ ,23
+ ,67
+ ,22618
+ ,5.679993974
+ ,1
+ ,54.16083333
+ ,83
+ ,40
+ ,151
+ ,74011
+ ,13.621099
+ ,1
+ ,46.52444444
+ ,113
+ ,28
+ ,72
+ ,83737
+ ,11.4493559
+ ,1
+ ,39.93222222
+ ,118
+ ,34
+ ,120
+ ,69094
+ ,12.45612266
+ ,1
+ ,76.53916667
+ ,110
+ ,33
+ ,115
+ ,93133
+ ,14.52255588
+ ,1
+ ,67.55527778
+ ,129
+ ,28
+ ,105
+ ,95536
+ ,14.00916003
+ ,1
+ ,50.83305556
+ ,51
+ ,34
+ ,104
+ ,225920
+ ,12.87917848
+ ,1
+ ,37.68027778
+ ,93
+ ,30
+ ,108
+ ,62133
+ ,10.91453539
+ ,1
+ ,42.30527778
+ ,76
+ ,33
+ ,98
+ ,61370
+ ,10.67710167
+ ,1
+ ,33.39472222
+ ,49
+ ,22
+ ,69
+ ,43836
+ ,7.500506329
+ ,1
+ ,96.24583333
+ ,118
+ ,38
+ ,111
+ ,106117
+ ,15.71511893
+ ,1
+ ,40.49722222
+ ,38
+ ,26
+ ,99
+ ,38692
+ ,8.451414364
+ ,1
+ ,53.70527778
+ ,141
+ ,35
+ ,71
+ ,84651
+ ,13.04281114
+ ,1
+ ,22.48694444
+ ,58
+ ,8
+ ,27
+ ,56622
+ ,5.475735648
+ ,1
+ ,34.10388889
+ ,27
+ ,24
+ ,69
+ ,15986
+ ,6.318561991
+ ,1
+ ,36.27361111
+ ,91
+ ,29
+ ,107
+ ,95364
+ ,11.67102852
+ ,1
+ ,79.57444444
+ ,63
+ ,29
+ ,107
+ ,89691
+ ,12.88262502
+ ,1
+ ,66.96277778
+ ,56
+ ,45
+ ,93
+ ,67267
+ ,12.34103121
+ ,1
+ ,41.235
+ ,144
+ ,37
+ ,129
+ ,126846
+ ,15.31887313
+ ,1
+ ,56.86472222
+ ,73
+ ,33
+ ,69
+ ,41140
+ ,10.00450285
+ ,1
+ ,50.5775
+ ,168
+ ,33
+ ,118
+ ,102860
+ ,14.86835659
+ ,1
+ ,38.98444444
+ ,64
+ ,25
+ ,73
+ ,51715
+ ,8.706868119
+ ,1
+ ,61.25444444
+ ,97
+ ,32
+ ,119
+ ,55801
+ ,12.38401963
+ ,1
+ ,67.51666667
+ ,117
+ ,29
+ ,104
+ ,111813
+ ,14.23699191
+ ,1
+ ,45.2125
+ ,100
+ ,28
+ ,107
+ ,120293
+ ,12.98177075
+ ,1
+ ,50.72583333
+ ,149
+ ,28
+ ,99
+ ,138599
+ ,14.75112225
+ ,1
+ ,64.48277778
+ ,187
+ ,31
+ ,90
+ ,161647
+ ,15.68212424
+ ,1
+ ,73.69944444
+ ,127
+ ,52
+ ,197
+ ,115929
+ ,18.1400865
+ ,1
+ ,86.34416667
+ ,245
+ ,24
+ ,85
+ ,162901
+ ,15.91867058
+ ,1
+ ,62.51666667
+ ,87
+ ,41
+ ,139
+ ,109825
+ ,14.96179352
+ ,1
+ ,64.5325
+ ,177
+ ,33
+ ,106
+ ,129838
+ ,16.04078615
+ ,1
+ ,40.26833333
+ ,49
+ ,32
+ ,50
+ ,37510
+ ,7.958769343
+ ,1
+ ,12.02416667
+ ,49
+ ,19
+ ,64
+ ,43750
+ ,6.117022453
+ ,1
+ ,43.265
+ ,73
+ ,20
+ ,31
+ ,40652
+ ,7.408636579
+ ,1
+ ,45.7525
+ ,177
+ ,31
+ ,63
+ ,87771
+ ,12.70501195
+ ,1
+ ,56.09444444
+ ,94
+ ,31
+ ,92
+ ,85872
+ ,12.21180573
+ ,1
+ ,65.40388889
+ ,117
+ ,32
+ ,106
+ ,89275
+ ,13.7541755
+ ,1
+ ,27.62944444
+ ,55
+ ,23
+ ,69
+ ,192565
+ ,10.16398398
+ ,1
+ ,27.98611111
+ ,58
+ ,30
+ ,93
+ ,140867
+ ,11.38027871
+ ,1
+ ,62.37472222
+ ,95
+ ,31
+ ,114
+ ,120662
+ ,14.09406222
+ ,1
+ ,67.64194444
+ ,129
+ ,42
+ ,110
+ ,101338
+ ,15.38882422
+ ,1
+ ,6.371666667
+ ,11
+ ,1
+ ,0
+ ,1168
+ ,0.694512385
+ ,1
+ ,42.35388889
+ ,101
+ ,32
+ ,83
+ ,65567
+ ,10.98992814
+ ,1
+ ,17.1825
+ ,28
+ ,11
+ ,30
+ ,25162
+ ,3.847727223
+ ,1
+ ,36.80194444
+ ,89
+ ,36
+ ,98
+ ,40735
+ ,10.36537965
+ ,1
+ ,88.165
+ ,193
+ ,31
+ ,82
+ ,91413
+ ,15.07282766
+ ,1
+ ,5.848333333
+ ,4
+ ,0
+ ,0
+ ,855
+ ,0.406821823
+ ,1
+ ,58.23361111
+ ,84
+ ,24
+ ,60
+ ,97068
+ ,11.07449882
+ ,1
+ ,8.726111111
+ ,39
+ ,8
+ ,9
+ ,14116
+ ,2.651799226
+ ,1
+ ,67.98583333
+ ,101
+ ,33
+ ,115
+ ,76643
+ ,13.40079885
+ ,1
+ ,51.25277778
+ ,82
+ ,40
+ ,140
+ ,110681
+ ,14.266541
+ ,1
+ ,35.67305556
+ ,36
+ ,38
+ ,120
+ ,92696
+ ,11.1935173
+ ,1
+ ,27.1775
+ ,75
+ ,24
+ ,66
+ ,94785
+ ,9.43355457
+ ,1
+ ,10.615
+ ,16
+ ,8
+ ,21
+ ,8773
+ ,2.298582242
+ ,1
+ ,41.9725
+ ,55
+ ,35
+ ,124
+ ,83209
+ ,11.55994194
+ ,1
+ ,75.68277778
+ ,131
+ ,43
+ ,152
+ ,93815
+ ,16.70032444
+ ,1
+ ,47.915
+ ,131
+ ,43
+ ,139
+ ,86687
+ ,14.84074833
+ ,1
+ ,91.14083333
+ ,144
+ ,41
+ ,144
+ ,105547
+ ,17.38403831
+ ,1
+ ,69.60527778
+ ,139
+ ,38
+ ,120
+ ,103487
+ ,15.73039586
+ ,1
+ ,97.51861111
+ ,211
+ ,45
+ ,160
+ ,213688
+ ,19.36265885
+ ,1
+ ,43.89305556
+ ,78
+ ,31
+ ,114
+ ,71220
+ ,11.32605381
+ ,1
+ ,23.73305556
+ ,39
+ ,28
+ ,78
+ ,56926
+ ,7.851840967
+ ,1
+ ,63.67833333
+ ,90
+ ,31
+ ,119
+ ,91721
+ ,13.30225262
+ ,1
+ ,97.67194444
+ ,166
+ ,40
+ ,141
+ ,115168
+ ,17.89762797
+ ,1
+ ,23.39083333
+ ,12
+ ,30
+ ,101
+ ,111194
+ ,9.46320715
+ ,1
+ ,90.16611111
+ ,133
+ ,37
+ ,133
+ ,135777
+ ,17.41915883
+ ,1
+ ,36.40805556
+ ,69
+ ,30
+ ,83
+ ,51513
+ ,9.332238773
+ ,1
+ ,56.74194444
+ ,119
+ ,35
+ ,116
+ ,74163
+ ,13.4187927
+ ,1
+ ,45.98416667
+ ,119
+ ,32
+ ,90
+ ,51633
+ ,11.37725913
+ ,1
+ ,39.36722222
+ ,65
+ ,27
+ ,36
+ ,75345
+ ,8.702323456
+ ,1
+ ,83.27083333
+ ,101
+ ,31
+ ,97
+ ,98952
+ ,14.19611665
+ ,1
+ ,54.39944444
+ ,196
+ ,31
+ ,98
+ ,102372
+ ,14.3980936
+ ,1
+ ,48.12777778
+ ,15
+ ,21
+ ,78
+ ,37238
+ ,7.301383906
+ ,1
+ ,70.69111111
+ ,136
+ ,39
+ ,117
+ ,103772
+ ,15.72027262
+ ,1
+ ,28.99694444
+ ,89
+ ,41
+ ,148
+ ,123969
+ ,14.03405468
+ ,1
+ ,55.41
+ ,123
+ ,32
+ ,105
+ ,135400
+ ,14.75498272
+ ,1
+ ,62.31388889
+ ,163
+ ,39
+ ,132
+ ,130115
+ ,17.03707132
+ ,1
+ ,4.08
+ ,5
+ ,0
+ ,0
+ ,6023
+ ,0.498724543
+ ,1
+ ,50.45361111
+ ,96
+ ,30
+ ,73
+ ,64466
+ ,10.82410397
+ ,1
+ ,75.51555556
+ ,151
+ ,37
+ ,86
+ ,54990
+ ,13.99003257
+ ,1
+ ,1.999722222
+ ,6
+ ,0
+ ,0
+ ,1644
+ ,0.295218143
+ ,1
+ ,12.96111111
+ ,13
+ ,5
+ ,13
+ ,6179
+ ,1.833602184
+ ,1
+ ,4.874166667
+ ,3
+ ,1
+ ,4
+ ,3926
+ ,0.59944272
+ ,1
+ ,26.45194444
+ ,23
+ ,32
+ ,48
+ ,34777
+ ,6.512322257
+ ,1
+ ,42.38916667
+ ,57
+ ,24
+ ,46
+ ,73224
+ ,8.594439161
+ ,1
+ ,28.23472222
+ ,28
+ ,11
+ ,38
+ ,17140
+ ,10.17463007
+ ,0
+ ,28.05861111
+ ,32
+ ,13
+ ,39
+ ,27570
+ ,11.70685947
+ ,0
+ ,1.993333333
+ ,0
+ ,0
+ ,0
+ ,1423
+ ,0.336629956
+ ,0
+ ,26.82222222
+ ,47
+ ,17
+ ,38
+ ,22996
+ ,12.81545815
+ ,0
+ ,48.84
+ ,65
+ ,20
+ ,77
+ ,39992
+ ,18.48134929
+ ,0
+ ,94.88055556
+ ,123
+ ,21
+ ,78
+ ,117105
+ ,19.15
+ ,0
+ ,28.77694444
+ ,26
+ ,16
+ ,49
+ ,23789
+ ,11.98976509
+ ,0
+ ,31.28083333
+ ,48
+ ,20
+ ,73
+ ,26706
+ ,15.96311694
+ ,0
+ ,23.77055556
+ ,37
+ ,21
+ ,36
+ ,24266
+ ,12.36867636
+ ,0
+ ,61.33361111
+ ,60
+ ,18
+ ,63
+ ,44418
+ ,17.85359417
+ ,0
+ ,25.73916667
+ ,39
+ ,17
+ ,41
+ ,35232
+ ,13.33972078
+ ,0
+ ,37.03555556
+ ,64
+ ,20
+ ,56
+ ,40909
+ ,17.43586247
+ ,0
+ ,17.04472222
+ ,26
+ ,12
+ ,25
+ ,13294
+ ,7.971763963
+ ,0
+ ,34.98055556
+ ,64
+ ,17
+ ,65
+ ,32387
+ ,16.49423975
+ ,0
+ ,22.86555556
+ ,25
+ ,10
+ ,38
+ ,21233
+ ,9.589676646
+ ,0
+ ,28.33611111
+ ,26
+ ,13
+ ,44
+ ,44332
+ ,13.0201353
+ ,0
+ ,28.20083333
+ ,76
+ ,22
+ ,87
+ ,61056
+ ,18.30998406
+ ,0
+ ,11.54611111
+ ,2
+ ,9
+ ,27
+ ,13497
+ ,5.311440585
+ ,0
+ ,27.75638889
+ ,36
+ ,25
+ ,80
+ ,32334
+ ,16.14544025
+ ,0
+ ,6.291111111
+ ,23
+ ,13
+ ,28
+ ,44339
+ ,9.61483208
+ ,0
+ ,12.97166667
+ ,14
+ ,13
+ ,33
+ ,10288
+ ,6.97001928
+ ,0
+ ,36.58277778
+ ,78
+ ,19
+ ,59
+ ,65622
+ ,17.8328823
+ ,0
+ ,25.48194444
+ ,14
+ ,18
+ ,49
+ ,16563
+ ,10.43631816
+ ,0
+ ,22.18416667
+ ,24
+ ,22
+ ,49
+ ,29011
+ ,12.49521893
+ ,0
+ ,30.01194444
+ ,39
+ ,14
+ ,38
+ ,34553
+ ,13.13014142
+ ,0
+ ,27.46277778
+ ,50
+ ,13
+ ,39
+ ,23517
+ ,12.5995365
+ ,0
+ ,33.45694444
+ ,57
+ ,16
+ ,56
+ ,51009
+ ,16.89272827
+ ,0
+ ,32.23555556
+ ,61
+ ,20
+ ,50
+ ,33416
+ ,15.95822221
+ ,0
+ ,69.4575
+ ,49
+ ,18
+ ,61
+ ,83305
+ ,17.44624506
+ ,0
+ ,37.80111111
+ ,40
+ ,13
+ ,41
+ ,27142
+ ,13.36096345
+ ,0
+ ,25.69416667
+ ,21
+ ,18
+ ,55
+ ,21399
+ ,11.70208318
+ ,0
+ ,37.71694444
+ ,29
+ ,14
+ ,44
+ ,24874
+ ,12.6700532
+ ,0
+ ,20.66888889
+ ,35
+ ,7
+ ,21
+ ,34988
+ ,9.942922584
+ ,0
+ ,22.56666667
+ ,13
+ ,17
+ ,50
+ ,45549
+ ,12.47787034
+ ,0
+ ,37.04666667
+ ,56
+ ,16
+ ,57
+ ,32755
+ ,16.1794073
+ ,0
+ ,27.26277778
+ ,14
+ ,17
+ ,48
+ ,27114
+ ,11.34269609
+ ,0
+ ,22.11638889
+ ,43
+ ,11
+ ,32
+ ,20760
+ ,10.61458743
+ ,0
+ ,16.44277778
+ ,20
+ ,24
+ ,68
+ ,37636
+ ,13.63362964
+ ,0
+ ,38.87277778
+ ,72
+ ,22
+ ,87
+ ,65461
+ ,19.3
+ ,0
+ ,32.94777778
+ ,87
+ ,12
+ ,43
+ ,30080
+ ,14.17862044
+ ,0
+ ,20.24444444
+ ,21
+ ,19
+ ,67
+ ,24094
+ ,12.13110302
+ ,0
+ ,18.1875
+ ,56
+ ,13
+ ,46
+ ,69008
+ ,14.29122225
+ ,0
+ ,27.67861111
+ ,59
+ ,17
+ ,46
+ ,54968
+ ,15.90424527
+ ,0
+ ,19.99027778
+ ,82
+ ,15
+ ,56
+ ,46090
+ ,15.30537883
+ ,0
+ ,21.46444444
+ ,43
+ ,16
+ ,48
+ ,27507
+ ,12.71711954
+ ,0
+ ,13.69138889
+ ,25
+ ,24
+ ,44
+ ,10672
+ ,10.1041787
+ ,0
+ ,37.53638889
+ ,38
+ ,15
+ ,60
+ ,34029
+ ,15.10633448
+ ,0
+ ,30.12388889
+ ,25
+ ,17
+ ,65
+ ,46300
+ ,14.97472539
+ ,0
+ ,24.92944444
+ ,38
+ ,18
+ ,55
+ ,24760
+ ,13.14940394
+ ,0
+ ,12.30444444
+ ,12
+ ,20
+ ,38
+ ,18779
+ ,8.803370342
+ ,0
+ ,21.56888889
+ ,29
+ ,16
+ ,52
+ ,21280
+ ,11.37671367
+ ,0
+ ,50.42444444
+ ,47
+ ,16
+ ,60
+ ,40662
+ ,16.48829016
+ ,0
+ ,37.2275
+ ,45
+ ,18
+ ,54
+ ,28987
+ ,15.28690517
+ ,0
+ ,34.46222222
+ ,40
+ ,22
+ ,86
+ ,22827
+ ,15.87449367
+ ,0
+ ,25.73055556
+ ,30
+ ,8
+ ,24
+ ,18513
+ ,8.991891655
+ ,0
+ ,33.84666667
+ ,41
+ ,17
+ ,52
+ ,30594
+ ,14.52069002
+ ,0
+ ,14.69861111
+ ,25
+ ,18
+ ,49
+ ,24006
+ ,10.73329293
+ ,0
+ ,22.74222222
+ ,23
+ ,16
+ ,61
+ ,27913
+ ,12.11255989
+ ,0
+ ,16.38361111
+ ,14
+ ,23
+ ,61
+ ,42744
+ ,13.12039076
+ ,0
+ ,14.86527778
+ ,16
+ ,22
+ ,81
+ ,12934
+ ,11.17618446
+ ,0
+ ,16.89222222
+ ,26
+ ,13
+ ,43
+ ,22574
+ ,9.852452086
+ ,0
+ ,15.65972222
+ ,21
+ ,13
+ ,40
+ ,41385
+ ,10.83826353
+ ,0
+ ,18.19166667
+ ,27
+ ,16
+ ,40
+ ,18653
+ ,10.01602482
+ ,0
+ ,22.48583333
+ ,9
+ ,16
+ ,56
+ ,18472
+ ,9.98429329
+ ,0
+ ,21.195
+ ,33
+ ,20
+ ,68
+ ,30976
+ ,13.90654664
+ ,0
+ ,28.89194444
+ ,42
+ ,22
+ ,79
+ ,63339
+ ,17.43829448
+ ,0
+ ,27.25111111
+ ,68
+ ,17
+ ,47
+ ,25568
+ ,14.13651406
+ ,0
+ ,18.88583333
+ ,32
+ ,18
+ ,57
+ ,33747
+ ,12.95465875
+ ,0
+ ,8.608055556
+ ,6
+ ,17
+ ,41
+ ,4154
+ ,6.405875596
+ ,0
+ ,37.62722222
+ ,67
+ ,12
+ ,29
+ ,19474
+ ,13.00827743
+ ,0
+ ,20.41777778
+ ,33
+ ,7
+ ,3
+ ,35130
+ ,8.843896805
+ ,0
+ ,17.53416667
+ ,77
+ ,17
+ ,60
+ ,39067
+ ,14.95275123
+ ,0
+ ,17.015
+ ,46
+ ,14
+ ,30
+ ,13310
+ ,9.985399433
+ ,0
+ ,20.80944444
+ ,30
+ ,23
+ ,79
+ ,65892
+ ,15.85289071
+ ,0
+ ,8.826111111
+ ,0
+ ,17
+ ,47
+ ,4143
+ ,6.304871788
+ ,0
+ ,22.62138889
+ ,36
+ ,14
+ ,40
+ ,28579
+ ,11.70744131
+ ,0
+ ,24.21833333
+ ,46
+ ,15
+ ,48
+ ,51776
+ ,14.6847188
+ ,0
+ ,13.91388889
+ ,18
+ ,17
+ ,36
+ ,21152
+ ,9.063740747
+ ,0
+ ,18.2625
+ ,48
+ ,21
+ ,42
+ ,38084
+ ,14.09669644
+ ,0
+ ,15.73694444
+ ,29
+ ,18
+ ,49
+ ,27717
+ ,11.45807414
+ ,0
+ ,43.99972222
+ ,28
+ ,18
+ ,57
+ ,32928
+ ,14.5766967
+ ,0
+ ,12.90416667
+ ,34
+ ,17
+ ,12
+ ,11342
+ ,8.013459929
+ ,0
+ ,20.45111111
+ ,33
+ ,17
+ ,40
+ ,19499
+ ,10.91719297
+ ,0
+ ,10.66527778
+ ,34
+ ,16
+ ,43
+ ,16380
+ ,9.680449615
+ ,0
+ ,25.5275
+ ,33
+ ,15
+ ,33
+ ,36874
+ ,12.30631181
+ ,0
+ ,38.75722222
+ ,80
+ ,21
+ ,77
+ ,48259
+ ,19.15
+ ,0
+ ,14.49
+ ,32
+ ,16
+ ,43
+ ,16734
+ ,9.974038576
+ ,0
+ ,14.32416667
+ ,30
+ ,14
+ ,45
+ ,28207
+ ,10.61396958
+ ,0
+ ,19.5975
+ ,41
+ ,15
+ ,47
+ ,30143
+ ,12.39969169
+ ,0
+ ,23.57111111
+ ,41
+ ,17
+ ,43
+ ,41369
+ ,13.89235106
+ ,0
+ ,28.48277778
+ ,51
+ ,15
+ ,45
+ ,45833
+ ,15.3368091
+ ,0
+ ,24.07722222
+ ,18
+ ,15
+ ,50
+ ,29156
+ ,11.27570477
+ ,0
+ ,23.80805556
+ ,34
+ ,10
+ ,35
+ ,35944
+ ,11.46890639
+ ,0
+ ,9.628333333
+ ,31
+ ,6
+ ,7
+ ,36278
+ ,7.705473469
+ ,0
+ ,41.82777778
+ ,39
+ ,22
+ ,71
+ ,45588
+ ,17.80919895
+ ,0
+ ,27.66972222
+ ,54
+ ,21
+ ,67
+ ,45097
+ ,17.4519668
+ ,0
+ ,5.374722222
+ ,14
+ ,1
+ ,0
+ ,3895
+ ,2.080910919
+ ,0
+ ,27.60361111
+ ,24
+ ,18
+ ,62
+ ,28394
+ ,13.09820752
+ ,0
+ ,23.95277778
+ ,24
+ ,17
+ ,54
+ ,18632
+ ,11.29135539
+ ,0
+ ,8.565833333
+ ,8
+ ,4
+ ,4
+ ,2325
+ ,2.507172277
+ ,0
+ ,8.807222222
+ ,26
+ ,10
+ ,25
+ ,25139
+ ,7.823664497
+ ,0
+ ,24.94611111
+ ,19
+ ,16
+ ,40
+ ,27975
+ ,10.96662538
+ ,0
+ ,17.24666667
+ ,11
+ ,16
+ ,38
+ ,14483
+ ,8.284173862
+ ,0
+ ,11.15305556
+ ,14
+ ,9
+ ,19
+ ,13127
+ ,5.69261529
+ ,0
+ ,7.676111111
+ ,1
+ ,16
+ ,17
+ ,5839
+ ,4.687276457
+ ,0
+ ,21.38611111
+ ,39
+ ,17
+ ,67
+ ,24069
+ ,13.25987215
+ ,0
+ ,10.40555556
+ ,5
+ ,7
+ ,14
+ ,3738
+ ,3.580093526
+ ,0
+ ,15.04361111
+ ,37
+ ,15
+ ,30
+ ,18625
+ ,9.733116837
+ ,0
+ ,13.85055556
+ ,32
+ ,14
+ ,54
+ ,36341
+ ,11.88651299
+ ,0
+ ,23.42694444
+ ,38
+ ,14
+ ,35
+ ,24548
+ ,11.3271032
+ ,0
+ ,17.82638889
+ ,47
+ ,18
+ ,59
+ ,21792
+ ,12.99614065
+ ,0
+ ,16.495
+ ,47
+ ,12
+ ,24
+ ,26263
+ ,10.51715507
+ ,0
+ ,33.14111111
+ ,37
+ ,16
+ ,58
+ ,23686
+ ,13.71616844
+ ,0
+ ,21.30611111
+ ,51
+ ,21
+ ,42
+ ,49303
+ ,15.32447894
+ ,0
+ ,28.72916667
+ ,45
+ ,19
+ ,46
+ ,25659
+ ,13.82274766
+ ,0
+ ,19.54
+ ,21
+ ,16
+ ,61
+ ,28904
+ ,11.71158746
+ ,0
+ ,12.05833333
+ ,1
+ ,1
+ ,3
+ ,2781
+ ,1.908370816
+ ,0
+ ,29.12166667
+ ,42
+ ,16
+ ,52
+ ,29236
+ ,13.82088431
+ ,0
+ ,17.28194444
+ ,26
+ ,10
+ ,25
+ ,19546
+ ,8.241327352
+ ,0
+ ,19.25111111
+ ,21
+ ,19
+ ,40
+ ,22818
+ ,10.5053078
+ ,0
+ ,14.75472222
+ ,4
+ ,12
+ ,32
+ ,32689
+ ,8.180917703
+ ,0
+ ,5.49
+ ,10
+ ,2
+ ,4
+ ,5752
+ ,2.322655495
+ ,0
+ ,24.07777778
+ ,43
+ ,14
+ ,49
+ ,22197
+ ,12.28598307
+ ,0
+ ,23.3625
+ ,34
+ ,17
+ ,63
+ ,20055
+ ,12.54906448
+ ,0
+ ,21.65138889
+ ,31
+ ,19
+ ,67
+ ,25272
+ ,13.11109062
+ ,0
+ ,24.75361111
+ ,19
+ ,14
+ ,32
+ ,82206
+ ,11.7934321
+ ,0
+ ,25.27916667
+ ,34
+ ,11
+ ,23
+ ,32073
+ ,10.81286184
+ ,0
+ ,11.18
+ ,6
+ ,4
+ ,7
+ ,5444
+ ,3.068772109
+ ,0
+ ,17.82972222
+ ,11
+ ,16
+ ,54
+ ,20154
+ ,9.674855276
+ ,0
+ ,14.12694444
+ ,24
+ ,20
+ ,37
+ ,36944
+ ,11.39973018
+ ,0
+ ,15.72583333
+ ,16
+ ,12
+ ,35
+ ,8019
+ ,7.166265475
+ ,0
+ ,17.44222222
+ ,72
+ ,15
+ ,51
+ ,30884
+ ,13.46103292
+ ,0
+ ,20.14861111
+ ,21
+ ,16
+ ,39
+ ,19540
+ ,9.81357389
+ ,0)
+ ,dim=c(7
+ ,289)
+ ,dimnames=list(c('Tijd'
+ ,'totblogs'
+ ,'Reviews'
+ ,'LFM'
+ ,'totsize'
+ ,'Score'
+ ,'Student')
+ ,1:289))
> y <- array(NA,dim=c(7,289),dimnames=list(c('Tijd','totblogs','Reviews','LFM','totsize','Score','Student'),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 = '6'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(lattice)
> library(lmtest)
Loading required package: zoo
> n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test
> par1 <- as.numeric(par1)
> x <- t(y)
> k <- length(x[1,])
> n <- length(x[,1])
> x1 <- cbind(x[,par1], x[,1:k!=par1])
> mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1])
> colnames(x1) <- mycolnames #colnames(x)[par1]
> x <- x1
> if (par3 == 'First Differences'){
+ x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep='')))
+ for (i in 1:n-1) {
+ for (j in 1:k) {
+ x2[i,j] <- x[i+1,j] - x[i,j]
+ }
+ }
+ x <- x2
+ }
> if (par2 == 'Include Monthly Dummies'){
+ x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep ='')))
+ for (i in 1:11){
+ x2[seq(i,n,12),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> if (par2 == 'Include Quarterly Dummies'){
+ x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep ='')))
+ for (i in 1:3){
+ x2[seq(i,n,4),i] <- 1
+ }
+ x <- cbind(x, x2)
+ }
> k <- length(x[1,])
> if (par3 == 'Linear Trend'){
+ x <- cbind(x, c(1:n))
+ colnames(x)[k+1] <- 't'
+ }
> x
Score Tijd totblogs Reviews LFM totsize Student
1 14.3600978 58.585278 145 30 94 112285 1
2 11.3090029 33.606111 101 28 103 84786 1
3 12.4581089 49.030000 98 38 93 83123 1
4 13.5039428 49.811389 132 30 103 101193 1
5 7.2193583 34.218056 60 22 51 38361 1
6 7.3813672 14.651667 38 26 70 68504 1
7 15.2609623 107.092778 144 25 91 119182 1
8 3.1654150 9.213889 5 18 22 22807 1
9 11.7814391 41.405833 84 26 93 116174 1
10 9.2778084 45.957222 79 25 60 57635 1
11 14.2278049 65.892500 127 38 123 66198 1
12 13.3753608 48.146111 78 44 148 71701 1
13 9.4877505 36.980833 60 30 90 57793 1
14 15.2151167 71.909167 131 40 124 80444 1
15 10.4271201 50.023056 84 34 70 53855 1
16 17.8497865 90.221944 133 47 168 97668 1
17 15.8953488 64.156667 150 30 115 133824 1
18 12.5511544 65.773611 91 31 71 101481 1
19 11.4339609 37.631389 132 23 66 99645 1
20 15.5350763 56.368056 136 36 134 114789 1
21 14.5223933 59.763056 124 36 117 99052 1
22 13.8933158 95.638056 118 30 108 67654 1
23 9.7122786 42.759722 70 25 84 65553 1
24 14.1294250 36.928611 107 39 156 97500 1
25 12.8953736 48.534444 119 34 120 69112 1
26 12.1602659 48.448611 89 31 114 82753 1
27 13.0118169 62.652222 112 31 94 85323 1
28 13.2987368 62.120000 108 33 120 72654 1
29 7.7754151 34.671389 52 25 81 30727 1
30 13.2842973 61.582778 112 33 110 77873 1
31 15.1145104 58.546389 116 35 133 117478 1
32 13.7479724 47.296111 123 42 122 74007 1
33 16.4294262 72.378056 125 43 158 90183 1
34 8.5843773 23.570278 27 30 109 61542 1
35 16.4741068 81.784444 162 33 124 101494 1
36 10.8353995 59.900278 64 32 92 55813 1
37 14.4803578 90.307500 92 36 126 79215 1
38 9.8171533 46.539444 83 28 70 55461 1
39 5.3452329 29.557778 41 14 37 31081 1
40 14.3931087 73.824722 120 32 120 83122 1
41 12.8763059 74.903056 105 30 93 70106 1
42 10.7687797 41.420000 79 35 95 60578 1
43 10.4969408 42.464167 70 28 90 79892 1
44 8.4436714 31.018056 55 28 80 49810 1
45 8.4077828 32.335556 39 39 31 71570 1
46 13.9406525 100.639167 67 34 110 100708 1
47 6.1632646 21.888889 21 26 66 33032 1
48 14.4368829 50.879722 127 39 138 82875 1
49 17.7758385 77.212500 152 39 133 139077 1
50 12.2486023 41.841389 113 33 113 71595 1
51 11.4568557 46.891389 99 28 100 72260 1
52 1.1530933 6.718889 7 4 7 5950 1
53 17.3568684 91.463056 141 39 140 115762 1
54 5.2246521 18.063611 21 18 61 32551 1
55 5.2386969 28.082500 35 14 41 31701 1
56 12.6052315 60.818333 109 29 96 80670 1
57 17.9840404 67.792222 133 44 164 143558 1
58 15.3027933 64.813333 230 28 102 120733 1
59 16.2306506 71.239444 166 35 124 105195 1
60 11.1772860 57.266944 68 28 99 73107 1
61 17.6425016 86.520278 147 38 129 132068 1
62 14.4328892 65.500000 179 23 62 149193 1
63 9.8414415 49.427500 61 36 73 46821 1
64 13.1018580 57.548889 101 32 114 87011 1
65 12.7822915 54.598056 108 29 99 95260 1
66 9.8413530 48.384444 90 25 70 55183 1
67 12.5226713 39.790556 114 27 104 106671 1
68 12.8517713 52.099722 103 36 116 73511 1
69 13.1788982 52.133611 142 28 91 92945 1
70 9.4599263 33.060000 79 23 74 78664 1
71 13.1446023 50.608889 88 40 138 70054 1
72 5.6799940 20.435000 25 23 67 22618 1
73 13.6210990 54.160833 83 40 151 74011 1
74 11.4493559 46.524444 113 28 72 83737 1
75 12.4561227 39.932222 118 34 120 69094 1
76 14.5225559 76.539167 110 33 115 93133 1
77 14.0091600 67.555278 129 28 105 95536 1
78 12.8791785 50.833056 51 34 104 225920 1
79 10.9145354 37.680278 93 30 108 62133 1
80 10.6771017 42.305278 76 33 98 61370 1
81 7.5005063 33.394722 49 22 69 43836 1
82 15.7151189 96.245833 118 38 111 106117 1
83 8.4514144 40.497222 38 26 99 38692 1
84 13.0428111 53.705278 141 35 71 84651 1
85 5.4757356 22.486944 58 8 27 56622 1
86 6.3185620 34.103889 27 24 69 15986 1
87 11.6710285 36.273611 91 29 107 95364 1
88 12.8826250 79.574444 63 29 107 89691 1
89 12.3410312 66.962778 56 45 93 67267 1
90 15.3188731 41.235000 144 37 129 126846 1
91 10.0045028 56.864722 73 33 69 41140 1
92 14.8683566 50.577500 168 33 118 102860 1
93 8.7068681 38.984444 64 25 73 51715 1
94 12.3840196 61.254444 97 32 119 55801 1
95 14.2369919 67.516667 117 29 104 111813 1
96 12.9817708 45.212500 100 28 107 120293 1
97 14.7511222 50.725833 149 28 99 138599 1
98 15.6821242 64.482778 187 31 90 161647 1
99 18.1400865 73.699444 127 52 197 115929 1
100 15.9186706 86.344167 245 24 85 162901 1
101 14.9617935 62.516667 87 41 139 109825 1
102 16.0407861 64.532500 177 33 106 129838 1
103 7.9587693 40.268333 49 32 50 37510 1
104 6.1170225 12.024167 49 19 64 43750 1
105 7.4086366 43.265000 73 20 31 40652 1
106 12.7050119 45.752500 177 31 63 87771 1
107 12.2118057 56.094444 94 31 92 85872 1
108 13.7541755 65.403889 117 32 106 89275 1
109 10.1639840 27.629444 55 23 69 192565 1
110 11.3802787 27.986111 58 30 93 140867 1
111 14.0940622 62.374722 95 31 114 120662 1
112 15.3888242 67.641944 129 42 110 101338 1
113 0.6945124 6.371667 11 1 0 1168 1
114 10.9899281 42.353889 101 32 83 65567 1
115 3.8477272 17.182500 28 11 30 25162 1
116 10.3653796 36.801944 89 36 98 40735 1
117 15.0728277 88.165000 193 31 82 91413 1
118 0.4068218 5.848333 4 0 0 855 1
119 11.0744988 58.233611 84 24 60 97068 1
120 2.6517992 8.726111 39 8 9 14116 1
121 13.4007988 67.985833 101 33 115 76643 1
122 14.2665410 51.252778 82 40 140 110681 1
123 11.1935173 35.673056 36 38 120 92696 1
124 9.4335546 27.177500 75 24 66 94785 1
125 2.2985822 10.615000 16 8 21 8773 1
126 11.5599419 41.972500 55 35 124 83209 1
127 16.7003244 75.682778 131 43 152 93815 1
128 14.8407483 47.915000 131 43 139 86687 1
129 17.3840383 91.140833 144 41 144 105547 1
130 15.7303959 69.605278 139 38 120 103487 1
131 19.3626588 97.518611 211 45 160 213688 1
132 11.3260538 43.893056 78 31 114 71220 1
133 7.8518410 23.733056 39 28 78 56926 1
134 13.3022526 63.678333 90 31 119 91721 1
135 17.8976280 97.671944 166 40 141 115168 1
136 9.4632072 23.390833 12 30 101 111194 1
137 17.4191588 90.166111 133 37 133 135777 1
138 9.3322388 36.408056 69 30 83 51513 1
139 13.4187927 56.741944 119 35 116 74163 1
140 11.3772591 45.984167 119 32 90 51633 1
141 8.7023235 39.367222 65 27 36 75345 1
142 14.1961167 83.270833 101 31 97 98952 1
143 14.3980936 54.399444 196 31 98 102372 1
144 7.3013839 48.127778 15 21 78 37238 1
145 15.7202726 70.691111 136 39 117 103772 1
146 14.0340547 28.996944 89 41 148 123969 1
147 14.7549827 55.410000 123 32 105 135400 1
148 17.0370713 62.313889 163 39 132 130115 1
149 0.4987245 4.080000 5 0 0 6023 1
150 10.8241040 50.453611 96 30 73 64466 1
151 13.9900326 75.515556 151 37 86 54990 1
152 0.2952181 1.999722 6 0 0 1644 1
153 1.8336022 12.961111 13 5 13 6179 1
154 0.5994427 4.874167 3 1 4 3926 1
155 6.5123223 26.451944 23 32 48 34777 1
156 8.5944392 42.389167 57 24 46 73224 1
157 10.1746301 28.234722 28 11 38 17140 0
158 11.7068595 28.058611 32 13 39 27570 0
159 0.3366300 1.993333 0 0 0 1423 0
160 12.8154581 26.822222 47 17 38 22996 0
161 18.4813493 48.840000 65 20 77 39992 0
162 19.1500000 94.880556 123 21 78 117105 0
163 11.9897651 28.776944 26 16 49 23789 0
164 15.9631169 31.280833 48 20 73 26706 0
165 12.3686764 23.770556 37 21 36 24266 0
166 17.8535942 61.333611 60 18 63 44418 0
167 13.3397208 25.739167 39 17 41 35232 0
168 17.4358625 37.035556 64 20 56 40909 0
169 7.9717640 17.044722 26 12 25 13294 0
170 16.4942397 34.980556 64 17 65 32387 0
171 9.5896766 22.865556 25 10 38 21233 0
172 13.0201353 28.336111 26 13 44 44332 0
173 18.3099841 28.200833 76 22 87 61056 0
174 5.3114406 11.546111 2 9 27 13497 0
175 16.1454403 27.756389 36 25 80 32334 0
176 9.6148321 6.291111 23 13 28 44339 0
177 6.9700193 12.971667 14 13 33 10288 0
178 17.8328823 36.582778 78 19 59 65622 0
179 10.4363182 25.481944 14 18 49 16563 0
180 12.4952189 22.184167 24 22 49 29011 0
181 13.1301414 30.011944 39 14 38 34553 0
182 12.5995365 27.462778 50 13 39 23517 0
183 16.8927283 33.456944 57 16 56 51009 0
184 15.9582222 32.235556 61 20 50 33416 0
185 17.4462451 69.457500 49 18 61 83305 0
186 13.3609634 37.801111 40 13 41 27142 0
187 11.7020832 25.694167 21 18 55 21399 0
188 12.6700532 37.716944 29 14 44 24874 0
189 9.9429226 20.668889 35 7 21 34988 0
190 12.4778703 22.566667 13 17 50 45549 0
191 16.1794073 37.046667 56 16 57 32755 0
192 11.3426961 27.262778 14 17 48 27114 0
193 10.6145874 22.116389 43 11 32 20760 0
194 13.6336296 16.442778 20 24 68 37636 0
195 19.3000000 38.872778 72 22 87 65461 0
196 14.1786204 32.947778 87 12 43 30080 0
197 12.1311030 20.244444 21 19 67 24094 0
198 14.2912223 18.187500 56 13 46 69008 0
199 15.9042453 27.678611 59 17 46 54968 0
200 15.3053788 19.990278 82 15 56 46090 0
201 12.7171195 21.464444 43 16 48 27507 0
202 10.1041787 13.691389 25 24 44 10672 0
203 15.1063345 37.536389 38 15 60 34029 0
204 14.9747254 30.123889 25 17 65 46300 0
205 13.1494039 24.929444 38 18 55 24760 0
206 8.8033703 12.304444 12 20 38 18779 0
207 11.3767137 21.568889 29 16 52 21280 0
208 16.4882902 50.424444 47 16 60 40662 0
209 15.2869052 37.227500 45 18 54 28987 0
210 15.8744937 34.462222 40 22 86 22827 0
211 8.9918917 25.730556 30 8 24 18513 0
212 14.5206900 33.846667 41 17 52 30594 0
213 10.7332929 14.698611 25 18 49 24006 0
214 12.1125599 22.742222 23 16 61 27913 0
215 13.1203908 16.383611 14 23 61 42744 0
216 11.1761845 14.865278 16 22 81 12934 0
217 9.8524521 16.892222 26 13 43 22574 0
218 10.8382635 15.659722 21 13 40 41385 0
219 10.0160248 18.191667 27 16 40 18653 0
220 9.9842933 22.485833 9 16 56 18472 0
221 13.9065466 21.195000 33 20 68 30976 0
222 17.4382945 28.891944 42 22 79 63339 0
223 14.1365141 27.251111 68 17 47 25568 0
224 12.9546588 18.885833 32 18 57 33747 0
225 6.4058756 8.608056 6 17 41 4154 0
226 13.0082774 37.627222 67 12 29 19474 0
227 8.8438968 20.417778 33 7 3 35130 0
228 14.9527512 17.534167 77 17 60 39067 0
229 9.9853994 17.015000 46 14 30 13310 0
230 15.8528907 20.809444 30 23 79 65892 0
231 6.3048718 8.826111 0 17 47 4143 0
232 11.7074413 22.621389 36 14 40 28579 0
233 14.6847188 24.218333 46 15 48 51776 0
234 9.0637407 13.913889 18 17 36 21152 0
235 14.0966964 18.262500 48 21 42 38084 0
236 11.4580741 15.736944 29 18 49 27717 0
237 14.5766967 43.999722 28 18 57 32928 0
238 8.0134599 12.904167 34 17 12 11342 0
239 10.9171930 20.451111 33 17 40 19499 0
240 9.6804496 10.665278 34 16 43 16380 0
241 12.3063118 25.527500 33 15 33 36874 0
242 19.1500000 38.757222 80 21 77 48259 0
243 9.9740386 14.490000 32 16 43 16734 0
244 10.6139696 14.324167 30 14 45 28207 0
245 12.3996917 19.597500 41 15 47 30143 0
246 13.8923511 23.571111 41 17 43 41369 0
247 15.3368091 28.482778 51 15 45 45833 0
248 11.2757048 24.077222 18 15 50 29156 0
249 11.4689064 23.808056 34 10 35 35944 0
250 7.7054735 9.628333 31 6 7 36278 0
251 17.8091989 41.827778 39 22 71 45588 0
252 17.4519668 27.669722 54 21 67 45097 0
253 2.0809109 5.374722 14 1 0 3895 0
254 13.0982075 27.603611 24 18 62 28394 0
255 11.2913554 23.952778 24 17 54 18632 0
256 2.5071723 8.565833 8 4 4 2325 0
257 7.8236645 8.807222 26 10 25 25139 0
258 10.9666254 24.946111 19 16 40 27975 0
259 8.2841739 17.246667 11 16 38 14483 0
260 5.6926153 11.153056 14 9 19 13127 0
261 4.6872765 7.676111 1 16 17 5839 0
262 13.2598721 21.386111 39 17 67 24069 0
263 3.5800935 10.405556 5 7 14 3738 0
264 9.7331168 15.043611 37 15 30 18625 0
265 11.8865130 13.850556 32 14 54 36341 0
266 11.3271032 23.426944 38 14 35 24548 0
267 12.9961406 17.826389 47 18 59 21792 0
268 10.5171551 16.495000 47 12 24 26263 0
269 13.7161684 33.141111 37 16 58 23686 0
270 15.3244789 21.306111 51 21 42 49303 0
271 13.8227477 28.729167 45 19 46 25659 0
272 11.7115875 19.540000 21 16 61 28904 0
273 1.9083708 12.058333 1 1 3 2781 0
274 13.8208843 29.121667 42 16 52 29236 0
275 8.2413274 17.281944 26 10 25 19546 0
276 10.5053078 19.251111 21 19 40 22818 0
277 8.1809177 14.754722 4 12 32 32689 0
278 2.3226555 5.490000 10 2 4 5752 0
279 12.2859831 24.077778 43 14 49 22197 0
280 12.5490645 23.362500 34 17 63 20055 0
281 13.1110906 21.651389 31 19 67 25272 0
282 11.7934321 24.753611 19 14 32 82206 0
283 10.8128618 25.279167 34 11 23 32073 0
284 3.0687721 11.180000 6 4 7 5444 0
285 9.6748553 17.829722 11 16 54 20154 0
286 11.3997302 14.126944 24 20 37 36944 0
287 7.1662655 15.725833 16 12 35 8019 0
288 13.4610329 17.442222 72 15 51 30884 0
289 9.8135739 20.148611 21 16 39 19540 0
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) Tijd totblogs Reviews LFM totsize
6.031e+00 4.962e-02 2.594e-02 9.035e-02 3.768e-02 2.180e-05
Student
-7.392e+00
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-5.8242 -0.4574 0.0870 0.6519 3.2702
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 6.031e+00 2.118e-01 28.476 < 2e-16 ***
Tijd 4.962e-02 7.083e-03 7.006 1.80e-11 ***
totblogs 2.594e-02 3.425e-03 7.572 5.27e-13 ***
Reviews 9.035e-02 2.143e-02 4.215 3.36e-05 ***
LFM 3.768e-02 5.796e-03 6.502 3.60e-10 ***
totsize 2.180e-05 3.864e-06 5.641 4.10e-08 ***
Student -7.392e+00 2.396e-01 -30.854 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 1.415 on 282 degrees of freedom
Multiple R-squared: 0.8779, Adjusted R-squared: 0.8753
F-statistic: 337.9 on 6 and 282 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,] 2.834494e-03 5.668988e-03 9.971655e-01
[2,] 5.020765e-04 1.004153e-03 9.994979e-01
[3,] 5.426159e-05 1.085232e-04 9.999457e-01
[4,] 5.792339e-06 1.158468e-05 9.999942e-01
[5,] 7.638308e-07 1.527662e-06 9.999992e-01
[6,] 7.061596e-08 1.412319e-07 9.999999e-01
[7,] 6.267841e-09 1.253568e-08 1.000000e+00
[8,] 1.160981e-09 2.321962e-09 1.000000e+00
[9,] 9.739555e-10 1.947911e-09 1.000000e+00
[10,] 1.964185e-10 3.928371e-10 1.000000e+00
[11,] 2.016899e-11 4.033798e-11 1.000000e+00
[12,] 2.053987e-12 4.107973e-12 1.000000e+00
[13,] 3.217125e-13 6.434251e-13 1.000000e+00
[14,] 6.750273e-14 1.350055e-13 1.000000e+00
[15,] 1.102250e-14 2.204500e-14 1.000000e+00
[16,] 1.093925e-15 2.187849e-15 1.000000e+00
[17,] 1.669543e-16 3.339086e-16 1.000000e+00
[18,] 3.105274e-17 6.210547e-17 1.000000e+00
[19,] 7.409633e-18 1.481927e-17 1.000000e+00
[20,] 8.003154e-19 1.600631e-18 1.000000e+00
[21,] 1.153452e-19 2.306904e-19 1.000000e+00
[22,] 1.600159e-20 3.200318e-20 1.000000e+00
[23,] 5.168278e-21 1.033656e-20 1.000000e+00
[24,] 7.945064e-22 1.589013e-21 1.000000e+00
[25,] 8.032361e-23 1.606472e-22 1.000000e+00
[26,] 1.567671e-23 3.135343e-23 1.000000e+00
[27,] 4.260427e-24 8.520853e-24 1.000000e+00
[28,] 5.568148e-25 1.113630e-24 1.000000e+00
[29,] 5.504157e-26 1.100831e-25 1.000000e+00
[30,] 6.831013e-27 1.366203e-26 1.000000e+00
[31,] 2.329914e-27 4.659827e-27 1.000000e+00
[32,] 9.212886e-28 1.842577e-27 1.000000e+00
[33,] 1.071930e-28 2.143861e-28 1.000000e+00
[34,] 1.123406e-29 2.246812e-29 1.000000e+00
[35,] 1.198441e-30 2.396882e-30 1.000000e+00
[36,] 1.493929e-31 2.987858e-31 1.000000e+00
[37,] 2.901038e-32 5.802076e-32 1.000000e+00
[38,] 3.001925e-33 6.003851e-33 1.000000e+00
[39,] 3.979933e-34 7.959866e-34 1.000000e+00
[40,] 4.738698e-35 9.477397e-35 1.000000e+00
[41,] 6.655261e-36 1.331052e-35 1.000000e+00
[42,] 6.160461e-37 1.232092e-36 1.000000e+00
[43,] 6.197978e-38 1.239596e-37 1.000000e+00
[44,] 6.716263e-39 1.343253e-38 1.000000e+00
[45,] 6.125295e-40 1.225059e-39 1.000000e+00
[46,] 6.323203e-41 1.264641e-40 1.000000e+00
[47,] 8.481967e-42 1.696393e-41 1.000000e+00
[48,] 9.428982e-43 1.885796e-42 1.000000e+00
[49,] 8.650000e-30 1.730000e-29 1.000000e+00
[50,] 1.659529e-30 3.319058e-30 1.000000e+00
[51,] 2.910306e-31 5.820613e-31 1.000000e+00
[52,] 6.098341e-32 1.219668e-31 1.000000e+00
[53,] 1.929544e-32 3.859088e-32 1.000000e+00
[54,] 3.124362e-33 6.248723e-33 1.000000e+00
[55,] 5.395378e-34 1.079076e-33 1.000000e+00
[56,] 1.068516e-34 2.137031e-34 1.000000e+00
[57,] 2.037957e-35 4.075913e-35 1.000000e+00
[58,] 3.395894e-36 6.791788e-36 1.000000e+00
[59,] 5.232231e-37 1.046446e-36 1.000000e+00
[60,] 1.319209e-37 2.638418e-37 1.000000e+00
[61,] 2.071171e-38 4.142343e-38 1.000000e+00
[62,] 3.728500e-39 7.457000e-39 1.000000e+00
[63,] 7.480238e-40 1.496048e-39 1.000000e+00
[64,] 1.366588e-40 2.733176e-40 1.000000e+00
[65,] 2.586554e-41 5.173107e-41 1.000000e+00
[66,] 3.881532e-42 7.763065e-42 1.000000e+00
[67,] 9.653106e-43 1.930621e-42 1.000000e+00
[68,] 2.834983e-43 5.669966e-43 1.000000e+00
[69,] 5.538955e-36 1.107791e-35 1.000000e+00
[70,] 1.161097e-36 2.322194e-36 1.000000e+00
[71,] 2.229457e-37 4.458913e-37 1.000000e+00
[72,] 4.208044e-38 8.416087e-38 1.000000e+00
[73,] 9.533766e-39 1.906753e-38 1.000000e+00
[74,] 1.886349e-39 3.772697e-39 1.000000e+00
[75,] 3.560880e-40 7.121760e-40 1.000000e+00
[76,] 9.456538e-41 1.891308e-40 1.000000e+00
[77,] 2.325265e-41 4.650531e-41 1.000000e+00
[78,] 4.829880e-42 9.659760e-42 1.000000e+00
[79,] 2.723110e-42 5.446219e-42 1.000000e+00
[80,] 5.106534e-43 1.021307e-42 1.000000e+00
[81,] 1.152322e-43 2.304645e-43 1.000000e+00
[82,] 2.106674e-44 4.213348e-44 1.000000e+00
[83,] 4.454133e-45 8.908265e-45 1.000000e+00
[84,] 7.800833e-46 1.560167e-45 1.000000e+00
[85,] 1.505087e-46 3.010173e-46 1.000000e+00
[86,] 6.995579e-47 1.399116e-46 1.000000e+00
[87,] 2.843392e-47 5.686784e-47 1.000000e+00
[88,] 1.112337e-47 2.224675e-47 1.000000e+00
[89,] 5.610602e-47 1.122120e-46 1.000000e+00
[90,] 3.203611e-46 6.407221e-46 1.000000e+00
[91,] 8.725334e-39 1.745067e-38 1.000000e+00
[92,] 2.437023e-39 4.874047e-39 1.000000e+00
[93,] 6.056694e-40 1.211339e-39 1.000000e+00
[94,] 1.552563e-40 3.105125e-40 1.000000e+00
[95,] 3.982901e-41 7.965803e-41 1.000000e+00
[96,] 1.067697e-41 2.135394e-41 1.000000e+00
[97,] 2.519565e-42 5.039129e-42 1.000000e+00
[98,] 6.541308e-43 1.308262e-42 1.000000e+00
[99,] 1.763790e-43 3.527581e-43 1.000000e+00
[100,] 7.240886e-44 1.448177e-43 1.000000e+00
[101,] 3.387539e-44 6.775078e-44 1.000000e+00
[102,] 1.615944e-44 3.231888e-44 1.000000e+00
[103,] 4.370156e-45 8.740311e-45 1.000000e+00
[104,] 1.529064e-45 3.058128e-45 1.000000e+00
[105,] 3.264815e-46 6.529631e-46 1.000000e+00
[106,] 9.058789e-47 1.811758e-46 1.000000e+00
[107,] 2.927145e-47 5.854291e-47 1.000000e+00
[108,] 2.232116e-47 4.464233e-47 1.000000e+00
[109,] 9.181339e-48 1.836268e-47 1.000000e+00
[110,] 5.996891e-48 1.199378e-47 1.000000e+00
[111,] 2.477393e-48 4.954787e-48 1.000000e+00
[112,] 5.675356e-49 1.135071e-48 1.000000e+00
[113,] 1.354849e-49 2.709698e-49 1.000000e+00
[114,] 2.847223e-50 5.694446e-50 1.000000e+00
[115,] 1.006024e-50 2.012049e-50 1.000000e+00
[116,] 4.492292e-51 8.984584e-51 1.000000e+00
[117,] 9.670108e-52 1.934022e-51 1.000000e+00
[118,] 3.149513e-52 6.299027e-52 1.000000e+00
[119,] 8.731456e-53 1.746291e-52 1.000000e+00
[120,] 3.815769e-53 7.631537e-53 1.000000e+00
[121,] 1.206599e-53 2.413198e-53 1.000000e+00
[122,] 1.716217e-40 3.432434e-40 1.000000e+00
[123,] 5.201091e-41 1.040218e-40 1.000000e+00
[124,] 1.560831e-41 3.121662e-41 1.000000e+00
[125,] 6.067805e-42 1.213561e-41 1.000000e+00
[126,] 2.165435e-41 4.330870e-41 1.000000e+00
[127,] 6.461385e-42 1.292277e-41 1.000000e+00
[128,] 1.406565e-41 2.813130e-41 1.000000e+00
[129,] 4.022198e-42 8.044395e-42 1.000000e+00
[130,] 2.107335e-42 4.214670e-42 1.000000e+00
[131,] 6.942096e-43 1.388419e-42 1.000000e+00
[132,] 3.711178e-43 7.422356e-43 1.000000e+00
[133,] 3.196550e-43 6.393100e-43 1.000000e+00
[134,] 9.715433e-43 1.943087e-42 1.000000e+00
[135,] 2.921122e-43 5.842244e-43 1.000000e+00
[136,] 7.315927e-43 1.463185e-42 1.000000e+00
[137,] 3.485666e-42 6.971331e-42 1.000000e+00
[138,] 5.028470e-41 1.005694e-40 1.000000e+00
[139,] 6.486953e-36 1.297391e-35 1.000000e+00
[140,] 1.128727e-35 2.257454e-35 1.000000e+00
[141,] 6.446006e-35 1.289201e-34 1.000000e+00
[142,] 5.646320e-30 1.129264e-29 1.000000e+00
[143,] 1.234544e-29 2.469089e-29 1.000000e+00
[144,] 1.664289e-29 3.328579e-29 1.000000e+00
[145,] 8.904996e-27 1.780999e-26 1.000000e+00
[146,] 4.526858e-27 9.053716e-27 1.000000e+00
[147,] 2.008433e-27 4.016866e-27 1.000000e+00
[148,] 7.848857e-28 1.569771e-27 1.000000e+00
[149,] 1.085941e-27 2.171883e-27 1.000000e+00
[150,] 9.965455e-14 1.993091e-13 1.000000e+00
[151,] 9.794144e-12 1.958829e-11 1.000000e+00
[152,] 8.789869e-08 1.757974e-07 9.999999e-01
[153,] 7.627887e-01 4.744226e-01 2.372113e-01
[154,] 7.495674e-01 5.008653e-01 2.504326e-01
[155,] 8.671869e-01 2.656262e-01 1.328131e-01
[156,] 8.512746e-01 2.974509e-01 1.487254e-01
[157,] 9.416611e-01 1.166778e-01 5.833889e-02
[158,] 9.428183e-01 1.143634e-01 5.718171e-02
[159,] 9.783203e-01 4.335940e-02 2.167970e-02
[160,] 9.883682e-01 2.326358e-02 1.163179e-02
[161,] 9.928963e-01 1.420737e-02 7.103686e-03
[162,] 9.932287e-01 1.354261e-02 6.771303e-03
[163,] 9.939652e-01 1.206954e-02 6.034772e-03
[164,] 9.986034e-01 2.793276e-03 1.396638e-03
[165,] 9.998521e-01 2.958583e-04 1.479291e-04
[166,] 9.998644e-01 2.711657e-04 1.355829e-04
[167,] 9.998517e-01 2.965821e-04 1.482911e-04
[168,] 9.999544e-01 9.112626e-05 4.556313e-05
[169,] 9.999906e-01 1.887614e-05 9.438071e-06
[170,] 9.999892e-01 2.167859e-05 1.083930e-05
[171,] 9.999838e-01 3.245186e-05 1.622593e-05
[172,] 9.999834e-01 3.325325e-05 1.662662e-05
[173,] 9.999780e-01 4.393777e-05 2.196889e-05
[174,] 9.999924e-01 1.521221e-05 7.606106e-06
[175,] 9.999934e-01 1.323264e-05 6.616321e-06
[176,] 1.000000e+00 4.171323e-12 2.085662e-12
[177,] 1.000000e+00 7.060001e-12 3.530000e-12
[178,] 1.000000e+00 1.252441e-11 6.262204e-12
[179,] 1.000000e+00 2.447027e-11 1.223513e-11
[180,] 1.000000e+00 2.775217e-11 1.387609e-11
[181,] 1.000000e+00 4.369427e-11 2.184713e-11
[182,] 1.000000e+00 2.370129e-11 1.185065e-11
[183,] 1.000000e+00 4.159284e-11 2.079642e-11
[184,] 1.000000e+00 6.373541e-11 3.186770e-11
[185,] 1.000000e+00 1.155283e-10 5.776415e-11
[186,] 1.000000e+00 2.073352e-13 1.036676e-13
[187,] 1.000000e+00 1.348508e-14 6.742542e-15
[188,] 1.000000e+00 2.802891e-14 1.401445e-14
[189,] 1.000000e+00 2.539369e-14 1.269685e-14
[190,] 1.000000e+00 3.077104e-14 1.538552e-14
[191,] 1.000000e+00 9.685894e-15 4.842947e-15
[192,] 1.000000e+00 1.941479e-14 9.707397e-15
[193,] 1.000000e+00 1.910080e-14 9.550400e-15
[194,] 1.000000e+00 2.199147e-14 1.099573e-14
[195,] 1.000000e+00 1.864083e-14 9.320413e-15
[196,] 1.000000e+00 4.078413e-14 2.039207e-14
[197,] 1.000000e+00 3.551461e-14 1.775730e-14
[198,] 1.000000e+00 7.264262e-14 3.632131e-14
[199,] 1.000000e+00 4.168693e-14 2.084347e-14
[200,] 1.000000e+00 6.568378e-14 3.284189e-14
[201,] 1.000000e+00 1.133193e-13 5.665966e-14
[202,] 1.000000e+00 1.645905e-13 8.229525e-14
[203,] 1.000000e+00 2.831069e-13 1.415535e-13
[204,] 1.000000e+00 5.362664e-13 2.681332e-13
[205,] 1.000000e+00 1.078745e-12 5.393727e-13
[206,] 1.000000e+00 1.937576e-12 9.687882e-13
[207,] 1.000000e+00 2.355285e-12 1.177643e-12
[208,] 1.000000e+00 3.975618e-12 1.987809e-12
[209,] 1.000000e+00 5.057029e-12 2.528515e-12
[210,] 1.000000e+00 9.040301e-12 4.520150e-12
[211,] 1.000000e+00 1.474983e-11 7.374913e-12
[212,] 1.000000e+00 2.789232e-11 1.394616e-11
[213,] 1.000000e+00 2.219433e-11 1.109717e-11
[214,] 1.000000e+00 2.338357e-11 1.169179e-11
[215,] 1.000000e+00 4.311138e-11 2.155569e-11
[216,] 1.000000e+00 1.690301e-11 8.451504e-12
[217,] 1.000000e+00 1.106479e-11 5.532394e-12
[218,] 1.000000e+00 1.649015e-11 8.245074e-12
[219,] 1.000000e+00 4.288184e-12 2.144092e-12
[220,] 1.000000e+00 8.275272e-12 4.137636e-12
[221,] 1.000000e+00 4.679428e-12 2.339714e-12
[222,] 1.000000e+00 1.961122e-12 9.805611e-13
[223,] 1.000000e+00 4.540173e-12 2.270086e-12
[224,] 1.000000e+00 8.305946e-12 4.152973e-12
[225,] 1.000000e+00 1.637679e-11 8.188396e-12
[226,] 1.000000e+00 3.181269e-11 1.590635e-11
[227,] 1.000000e+00 7.929842e-11 3.964921e-11
[228,] 1.000000e+00 1.441307e-10 7.206533e-11
[229,] 1.000000e+00 2.182562e-10 1.091281e-10
[230,] 1.000000e+00 5.182398e-10 2.591199e-10
[231,] 1.000000e+00 1.065126e-09 5.325631e-10
[232,] 1.000000e+00 1.871971e-09 9.359856e-10
[233,] 1.000000e+00 3.552087e-12 1.776043e-12
[234,] 1.000000e+00 8.678309e-12 4.339155e-12
[235,] 1.000000e+00 2.238333e-11 1.119167e-11
[236,] 1.000000e+00 6.203180e-11 3.101590e-11
[237,] 1.000000e+00 1.029325e-10 5.146624e-11
[238,] 1.000000e+00 9.927649e-11 4.963824e-11
[239,] 1.000000e+00 2.480016e-10 1.240008e-10
[240,] 1.000000e+00 3.155904e-10 1.577952e-10
[241,] 1.000000e+00 1.033421e-10 5.167105e-11
[242,] 1.000000e+00 2.346622e-10 1.173311e-10
[243,] 1.000000e+00 2.369405e-10 1.184703e-10
[244,] 1.000000e+00 2.518494e-10 1.259247e-10
[245,] 1.000000e+00 8.104910e-10 4.052455e-10
[246,] 1.000000e+00 2.311288e-09 1.155644e-09
[247,] 1.000000e+00 2.942609e-09 1.471305e-09
[248,] 1.000000e+00 5.226916e-09 2.613458e-09
[249,] 1.000000e+00 1.565138e-08 7.825690e-09
[250,] 1.000000e+00 3.638920e-08 1.819460e-08
[251,] 1.000000e+00 8.834429e-08 4.417215e-08
[252,] 1.000000e+00 5.958033e-08 2.979016e-08
[253,] 9.999999e-01 1.992885e-07 9.964425e-08
[254,] 9.999998e-01 3.007415e-07 1.503707e-07
[255,] 9.999996e-01 8.759806e-07 4.379903e-07
[256,] 9.999998e-01 4.072972e-07 2.036486e-07
[257,] 9.999993e-01 1.459014e-06 7.295071e-07
[258,] 9.999974e-01 5.104011e-06 2.552006e-06
[259,] 9.999938e-01 1.239139e-05 6.195697e-06
[260,] 9.999793e-01 4.133988e-05 2.066994e-05
[261,] 9.999745e-01 5.091147e-05 2.545573e-05
[262,] 9.999221e-01 1.557077e-04 7.785383e-05
[263,] 9.998395e-01 3.210950e-04 1.605475e-04
[264,] 9.996272e-01 7.455128e-04 3.727564e-04
[265,] 9.989171e-01 2.165793e-03 1.082897e-03
[266,] 9.967926e-01 6.414824e-03 3.207412e-03
[267,] 9.932158e-01 1.356836e-02 6.784179e-03
[268,] 9.931096e-01 1.378073e-02 6.890363e-03
[269,] 9.879964e-01 2.400723e-02 1.200362e-02
[270,] 9.570889e-01 8.582222e-02 4.291111e-02
> postscript(file="/var/www/rcomp/tmp/1ctay1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index')
> points(x[,1]-mysum$resid)
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/2zjfc1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/3g24b1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/47m681323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals')
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/5xqre1323794347.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
0.353447776 0.124101316 0.095249427 0.172508056 0.581000130 0.550180223
7 8 9 10 11 12
-0.712364949 0.987571355 0.523754888 0.533811918 -0.485513448 -0.790472529
13 14 15 16 17 18
0.096375531 -0.429417208 0.244142492 -1.421196835 0.221735266 0.600255555
19 20 21 22 23 24
0.767306941 -0.232030576 -0.118316463 -0.806351258 0.283438142 -0.644035368
25 26 27 28 29 30
-0.338186557 -0.091125263 0.156656404 -0.310518351 0.086982731 -0.138997750
31 32 33 34 35 36
-0.172526852 -0.432723930 -0.847246355 -0.083265197 -0.290938990 -0.009937046
37 38 39 40 41 42
-0.752768921 0.340050109 0.839940259 -0.246043291 0.054517472 -0.036535583
43 44 45 46 47 48
0.273080682 0.209389089 0.901381854 -0.841620803 0.337821978 -0.550555661
49 50 51 52 53 54
-0.203398334 -0.197215551 0.050665788 1.244699734 -0.799757231 0.510692758
55 56 57 58 59 60
0.797995505 0.125640322 -0.752461235 -1.522498142 -0.376244684 0.079517918
61 62 63 64 65 66
-0.275068584 0.234961524 0.144166545 -0.095520714 0.206229094 0.368233803
67 68 69 70 71 72
0.269318098 -0.269562050 0.285512551 0.550688302 -0.628713118 0.283306502
73 74 75 76 77 78
-0.774896756 0.503157277 -0.324245780 -0.112074798 0.103685425 -1.519850763
79 80 81 82 83 84
-0.140352161 -0.043945583 0.390707890 -0.688930172 -0.105215498 0.399366152
85 86 87 88 89 90
1.242534598 0.170571206 0.141485639 0.054303828 -0.109120715 -0.069551408
91 92 93 94 95 96
0.172524453 -0.307441438 0.337123054 -0.401587378 0.237216275 0.322173494
97 98 99 100 101 102
0.449426534 -0.722147556 -2.097930133 -2.281288338 -0.371420974 -0.196810165
103 104 105 106 107 108
0.458260498 0.529001176 0.868544092 0.117281433 0.212327662 0.004115381
109 110 111 112 113 114
-0.147653850 0.563248674 0.169670702 -0.100680341 1.338645234 0.182114568
115 116 117 118 119 120
0.957541739 -0.241010432 -0.829841975 1.355655266 0.822429303 1.199054064
121 122 123 124 125 126
-0.216514063 -0.344039932 -0.124484509 0.779673182 1.012956291 -0.236355757
127 128 129 130 131 132
-0.749041618 -0.585484487 -0.943165420 -0.178179883 -4.340352904 -0.162574828
133 134 135 136 137 138
0.314306861 -0.114698750 -1.330614045 0.412681893 -0.457428101 0.136508727
139 140 141 142 143 144
-0.271763566 -0.037612159 0.986085844 0.193000369 -0.748648570 0.237407042
145 146 147 148 149 150
-0.147886858 -0.335149481 0.377478559 -0.255142071 1.396718590 0.325557728
151 152 153 154 155 156
-0.094435062 1.365955093 1.138403646 1.314540944 0.506724660 0.876204020
157 158 159 160 161 162
-0.782913969 0.208579713 -5.824205929 0.765466644 2.760821401 -2.168455292
163 164 165 166 167 168
0.046028164 1.995182434 0.415783805 2.254548464 1.171211530 3.098356918
169 170 171 172 173 174
-1.895285765 2.376375092 -1.022496186 1.109929608 2.311651802 -3.469004122
175 176 177 178 179 180
1.825466102 -0.520887752 -2.709940362 2.593334964 -1.055846337 0.274598560
181 182 183 184 185 186
1.148491785 0.752294925 3.055593113 2.326134965 0.957100322 1.105705348
187 188 189 190 191 192
-0.313676969 0.550330324 -0.207825290 0.577142269 2.550293702 -0.339797264
193 194 195 196 197 198
-0.281232354 0.717014454 2.779832415 0.896076309 -0.215497981 1.493245884
199 200 201 202 203 204
2.502139719 1.685623035 0.651945895 -1.313485601 1.869326349 1.806149511
205 206 207 208 209 210
0.657391754 -1.797522786 -0.345509288 2.143379090 1.948615463 1.370192978
211 212 213 214 215 216
-1.124613552 1.584595355 -0.671330696 0.004021190 0.605110161 -1.329119578
217 218 219 220 221 222
-0.977925331 -0.098259811 -0.977297995 -1.354198645 0.923487214 2.539229161
223 224 225 226 227 228
1.125356491 0.646910196 -3.379182231 0.771004777 -0.567345065 1.406227299
229 230 231 232 233 234
-0.768373983 1.520151084 -3.561230345 0.225197131 1.966415975 -1.477965045
235 236 237 238 239 240
1.604524376 -0.182715239 1.144341888 -1.774950939 -0.452660946 -1.184456984
241 242 243 244 245 246
0.750283451 3.270192015 -1.036500911 -0.481220976 0.549604221 1.570415221
247 248 249 250 251 252
2.519795574 -0.291644646 0.368905148 -1.203899498 3.034423324 3.242458975
253 254 255 256 257 258
-4.755072447 0.493638876 -0.527463208 -4.719079375 -1.712122635 -0.357566977
259 260 261 262 263 264
-2.080996866 -3.070058304 -3.963877220 0.570982487 -4.338291005 -0.895595933
265 266 267 268 269 270
0.246513618 0.029284312 0.537130953 -0.112316544 0.933671690 2.358914802
271 272 273 274 275 276
1.189812322 -0.207776252 -5.010842972 1.213265329 -1.593080158 -0.746789407
277 278 279 280 281 282
-1.688409350 -4.696848935 0.349887485 0.130011051 0.409621313 -0.221176021
283 284 285 286 287 288
0.086070227 -4.416356886 -1.445765682 0.038876400 -2.637794407 0.746955607
289
-1.102885265
> postscript(file="/var/www/rcomp/tmp/6zbt61323794347.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 0.353447776 NA
1 0.124101316 0.353447776
2 0.095249427 0.124101316
3 0.172508056 0.095249427
4 0.581000130 0.172508056
5 0.550180223 0.581000130
6 -0.712364949 0.550180223
7 0.987571355 -0.712364949
8 0.523754888 0.987571355
9 0.533811918 0.523754888
10 -0.485513448 0.533811918
11 -0.790472529 -0.485513448
12 0.096375531 -0.790472529
13 -0.429417208 0.096375531
14 0.244142492 -0.429417208
15 -1.421196835 0.244142492
16 0.221735266 -1.421196835
17 0.600255555 0.221735266
18 0.767306941 0.600255555
19 -0.232030576 0.767306941
20 -0.118316463 -0.232030576
21 -0.806351258 -0.118316463
22 0.283438142 -0.806351258
23 -0.644035368 0.283438142
24 -0.338186557 -0.644035368
25 -0.091125263 -0.338186557
26 0.156656404 -0.091125263
27 -0.310518351 0.156656404
28 0.086982731 -0.310518351
29 -0.138997750 0.086982731
30 -0.172526852 -0.138997750
31 -0.432723930 -0.172526852
32 -0.847246355 -0.432723930
33 -0.083265197 -0.847246355
34 -0.290938990 -0.083265197
35 -0.009937046 -0.290938990
36 -0.752768921 -0.009937046
37 0.340050109 -0.752768921
38 0.839940259 0.340050109
39 -0.246043291 0.839940259
40 0.054517472 -0.246043291
41 -0.036535583 0.054517472
42 0.273080682 -0.036535583
43 0.209389089 0.273080682
44 0.901381854 0.209389089
45 -0.841620803 0.901381854
46 0.337821978 -0.841620803
47 -0.550555661 0.337821978
48 -0.203398334 -0.550555661
49 -0.197215551 -0.203398334
50 0.050665788 -0.197215551
51 1.244699734 0.050665788
52 -0.799757231 1.244699734
53 0.510692758 -0.799757231
54 0.797995505 0.510692758
55 0.125640322 0.797995505
56 -0.752461235 0.125640322
57 -1.522498142 -0.752461235
58 -0.376244684 -1.522498142
59 0.079517918 -0.376244684
60 -0.275068584 0.079517918
61 0.234961524 -0.275068584
62 0.144166545 0.234961524
63 -0.095520714 0.144166545
64 0.206229094 -0.095520714
65 0.368233803 0.206229094
66 0.269318098 0.368233803
67 -0.269562050 0.269318098
68 0.285512551 -0.269562050
69 0.550688302 0.285512551
70 -0.628713118 0.550688302
71 0.283306502 -0.628713118
72 -0.774896756 0.283306502
73 0.503157277 -0.774896756
74 -0.324245780 0.503157277
75 -0.112074798 -0.324245780
76 0.103685425 -0.112074798
77 -1.519850763 0.103685425
78 -0.140352161 -1.519850763
79 -0.043945583 -0.140352161
80 0.390707890 -0.043945583
81 -0.688930172 0.390707890
82 -0.105215498 -0.688930172
83 0.399366152 -0.105215498
84 1.242534598 0.399366152
85 0.170571206 1.242534598
86 0.141485639 0.170571206
87 0.054303828 0.141485639
88 -0.109120715 0.054303828
89 -0.069551408 -0.109120715
90 0.172524453 -0.069551408
91 -0.307441438 0.172524453
92 0.337123054 -0.307441438
93 -0.401587378 0.337123054
94 0.237216275 -0.401587378
95 0.322173494 0.237216275
96 0.449426534 0.322173494
97 -0.722147556 0.449426534
98 -2.097930133 -0.722147556
99 -2.281288338 -2.097930133
100 -0.371420974 -2.281288338
101 -0.196810165 -0.371420974
102 0.458260498 -0.196810165
103 0.529001176 0.458260498
104 0.868544092 0.529001176
105 0.117281433 0.868544092
106 0.212327662 0.117281433
107 0.004115381 0.212327662
108 -0.147653850 0.004115381
109 0.563248674 -0.147653850
110 0.169670702 0.563248674
111 -0.100680341 0.169670702
112 1.338645234 -0.100680341
113 0.182114568 1.338645234
114 0.957541739 0.182114568
115 -0.241010432 0.957541739
116 -0.829841975 -0.241010432
117 1.355655266 -0.829841975
118 0.822429303 1.355655266
119 1.199054064 0.822429303
120 -0.216514063 1.199054064
121 -0.344039932 -0.216514063
122 -0.124484509 -0.344039932
123 0.779673182 -0.124484509
124 1.012956291 0.779673182
125 -0.236355757 1.012956291
126 -0.749041618 -0.236355757
127 -0.585484487 -0.749041618
128 -0.943165420 -0.585484487
129 -0.178179883 -0.943165420
130 -4.340352904 -0.178179883
131 -0.162574828 -4.340352904
132 0.314306861 -0.162574828
133 -0.114698750 0.314306861
134 -1.330614045 -0.114698750
135 0.412681893 -1.330614045
136 -0.457428101 0.412681893
137 0.136508727 -0.457428101
138 -0.271763566 0.136508727
139 -0.037612159 -0.271763566
140 0.986085844 -0.037612159
141 0.193000369 0.986085844
142 -0.748648570 0.193000369
143 0.237407042 -0.748648570
144 -0.147886858 0.237407042
145 -0.335149481 -0.147886858
146 0.377478559 -0.335149481
147 -0.255142071 0.377478559
148 1.396718590 -0.255142071
149 0.325557728 1.396718590
150 -0.094435062 0.325557728
151 1.365955093 -0.094435062
152 1.138403646 1.365955093
153 1.314540944 1.138403646
154 0.506724660 1.314540944
155 0.876204020 0.506724660
156 -0.782913969 0.876204020
157 0.208579713 -0.782913969
158 -5.824205929 0.208579713
159 0.765466644 -5.824205929
160 2.760821401 0.765466644
161 -2.168455292 2.760821401
162 0.046028164 -2.168455292
163 1.995182434 0.046028164
164 0.415783805 1.995182434
165 2.254548464 0.415783805
166 1.171211530 2.254548464
167 3.098356918 1.171211530
168 -1.895285765 3.098356918
169 2.376375092 -1.895285765
170 -1.022496186 2.376375092
171 1.109929608 -1.022496186
172 2.311651802 1.109929608
173 -3.469004122 2.311651802
174 1.825466102 -3.469004122
175 -0.520887752 1.825466102
176 -2.709940362 -0.520887752
177 2.593334964 -2.709940362
178 -1.055846337 2.593334964
179 0.274598560 -1.055846337
180 1.148491785 0.274598560
181 0.752294925 1.148491785
182 3.055593113 0.752294925
183 2.326134965 3.055593113
184 0.957100322 2.326134965
185 1.105705348 0.957100322
186 -0.313676969 1.105705348
187 0.550330324 -0.313676969
188 -0.207825290 0.550330324
189 0.577142269 -0.207825290
190 2.550293702 0.577142269
191 -0.339797264 2.550293702
192 -0.281232354 -0.339797264
193 0.717014454 -0.281232354
194 2.779832415 0.717014454
195 0.896076309 2.779832415
196 -0.215497981 0.896076309
197 1.493245884 -0.215497981
198 2.502139719 1.493245884
199 1.685623035 2.502139719
200 0.651945895 1.685623035
201 -1.313485601 0.651945895
202 1.869326349 -1.313485601
203 1.806149511 1.869326349
204 0.657391754 1.806149511
205 -1.797522786 0.657391754
206 -0.345509288 -1.797522786
207 2.143379090 -0.345509288
208 1.948615463 2.143379090
209 1.370192978 1.948615463
210 -1.124613552 1.370192978
211 1.584595355 -1.124613552
212 -0.671330696 1.584595355
213 0.004021190 -0.671330696
214 0.605110161 0.004021190
215 -1.329119578 0.605110161
216 -0.977925331 -1.329119578
217 -0.098259811 -0.977925331
218 -0.977297995 -0.098259811
219 -1.354198645 -0.977297995
220 0.923487214 -1.354198645
221 2.539229161 0.923487214
222 1.125356491 2.539229161
223 0.646910196 1.125356491
224 -3.379182231 0.646910196
225 0.771004777 -3.379182231
226 -0.567345065 0.771004777
227 1.406227299 -0.567345065
228 -0.768373983 1.406227299
229 1.520151084 -0.768373983
230 -3.561230345 1.520151084
231 0.225197131 -3.561230345
232 1.966415975 0.225197131
233 -1.477965045 1.966415975
234 1.604524376 -1.477965045
235 -0.182715239 1.604524376
236 1.144341888 -0.182715239
237 -1.774950939 1.144341888
238 -0.452660946 -1.774950939
239 -1.184456984 -0.452660946
240 0.750283451 -1.184456984
241 3.270192015 0.750283451
242 -1.036500911 3.270192015
243 -0.481220976 -1.036500911
244 0.549604221 -0.481220976
245 1.570415221 0.549604221
246 2.519795574 1.570415221
247 -0.291644646 2.519795574
248 0.368905148 -0.291644646
249 -1.203899498 0.368905148
250 3.034423324 -1.203899498
251 3.242458975 3.034423324
252 -4.755072447 3.242458975
253 0.493638876 -4.755072447
254 -0.527463208 0.493638876
255 -4.719079375 -0.527463208
256 -1.712122635 -4.719079375
257 -0.357566977 -1.712122635
258 -2.080996866 -0.357566977
259 -3.070058304 -2.080996866
260 -3.963877220 -3.070058304
261 0.570982487 -3.963877220
262 -4.338291005 0.570982487
263 -0.895595933 -4.338291005
264 0.246513618 -0.895595933
265 0.029284312 0.246513618
266 0.537130953 0.029284312
267 -0.112316544 0.537130953
268 0.933671690 -0.112316544
269 2.358914802 0.933671690
270 1.189812322 2.358914802
271 -0.207776252 1.189812322
272 -5.010842972 -0.207776252
273 1.213265329 -5.010842972
274 -1.593080158 1.213265329
275 -0.746789407 -1.593080158
276 -1.688409350 -0.746789407
277 -4.696848935 -1.688409350
278 0.349887485 -4.696848935
279 0.130011051 0.349887485
280 0.409621313 0.130011051
281 -0.221176021 0.409621313
282 0.086070227 -0.221176021
283 -4.416356886 0.086070227
284 -1.445765682 -4.416356886
285 0.038876400 -1.445765682
286 -2.637794407 0.038876400
287 0.746955607 -2.637794407
288 -1.102885265 0.746955607
289 NA -1.102885265
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] 0.124101316 0.353447776
[2,] 0.095249427 0.124101316
[3,] 0.172508056 0.095249427
[4,] 0.581000130 0.172508056
[5,] 0.550180223 0.581000130
[6,] -0.712364949 0.550180223
[7,] 0.987571355 -0.712364949
[8,] 0.523754888 0.987571355
[9,] 0.533811918 0.523754888
[10,] -0.485513448 0.533811918
[11,] -0.790472529 -0.485513448
[12,] 0.096375531 -0.790472529
[13,] -0.429417208 0.096375531
[14,] 0.244142492 -0.429417208
[15,] -1.421196835 0.244142492
[16,] 0.221735266 -1.421196835
[17,] 0.600255555 0.221735266
[18,] 0.767306941 0.600255555
[19,] -0.232030576 0.767306941
[20,] -0.118316463 -0.232030576
[21,] -0.806351258 -0.118316463
[22,] 0.283438142 -0.806351258
[23,] -0.644035368 0.283438142
[24,] -0.338186557 -0.644035368
[25,] -0.091125263 -0.338186557
[26,] 0.156656404 -0.091125263
[27,] -0.310518351 0.156656404
[28,] 0.086982731 -0.310518351
[29,] -0.138997750 0.086982731
[30,] -0.172526852 -0.138997750
[31,] -0.432723930 -0.172526852
[32,] -0.847246355 -0.432723930
[33,] -0.083265197 -0.847246355
[34,] -0.290938990 -0.083265197
[35,] -0.009937046 -0.290938990
[36,] -0.752768921 -0.009937046
[37,] 0.340050109 -0.752768921
[38,] 0.839940259 0.340050109
[39,] -0.246043291 0.839940259
[40,] 0.054517472 -0.246043291
[41,] -0.036535583 0.054517472
[42,] 0.273080682 -0.036535583
[43,] 0.209389089 0.273080682
[44,] 0.901381854 0.209389089
[45,] -0.841620803 0.901381854
[46,] 0.337821978 -0.841620803
[47,] -0.550555661 0.337821978
[48,] -0.203398334 -0.550555661
[49,] -0.197215551 -0.203398334
[50,] 0.050665788 -0.197215551
[51,] 1.244699734 0.050665788
[52,] -0.799757231 1.244699734
[53,] 0.510692758 -0.799757231
[54,] 0.797995505 0.510692758
[55,] 0.125640322 0.797995505
[56,] -0.752461235 0.125640322
[57,] -1.522498142 -0.752461235
[58,] -0.376244684 -1.522498142
[59,] 0.079517918 -0.376244684
[60,] -0.275068584 0.079517918
[61,] 0.234961524 -0.275068584
[62,] 0.144166545 0.234961524
[63,] -0.095520714 0.144166545
[64,] 0.206229094 -0.095520714
[65,] 0.368233803 0.206229094
[66,] 0.269318098 0.368233803
[67,] -0.269562050 0.269318098
[68,] 0.285512551 -0.269562050
[69,] 0.550688302 0.285512551
[70,] -0.628713118 0.550688302
[71,] 0.283306502 -0.628713118
[72,] -0.774896756 0.283306502
[73,] 0.503157277 -0.774896756
[74,] -0.324245780 0.503157277
[75,] -0.112074798 -0.324245780
[76,] 0.103685425 -0.112074798
[77,] -1.519850763 0.103685425
[78,] -0.140352161 -1.519850763
[79,] -0.043945583 -0.140352161
[80,] 0.390707890 -0.043945583
[81,] -0.688930172 0.390707890
[82,] -0.105215498 -0.688930172
[83,] 0.399366152 -0.105215498
[84,] 1.242534598 0.399366152
[85,] 0.170571206 1.242534598
[86,] 0.141485639 0.170571206
[87,] 0.054303828 0.141485639
[88,] -0.109120715 0.054303828
[89,] -0.069551408 -0.109120715
[90,] 0.172524453 -0.069551408
[91,] -0.307441438 0.172524453
[92,] 0.337123054 -0.307441438
[93,] -0.401587378 0.337123054
[94,] 0.237216275 -0.401587378
[95,] 0.322173494 0.237216275
[96,] 0.449426534 0.322173494
[97,] -0.722147556 0.449426534
[98,] -2.097930133 -0.722147556
[99,] -2.281288338 -2.097930133
[100,] -0.371420974 -2.281288338
[101,] -0.196810165 -0.371420974
[102,] 0.458260498 -0.196810165
[103,] 0.529001176 0.458260498
[104,] 0.868544092 0.529001176
[105,] 0.117281433 0.868544092
[106,] 0.212327662 0.117281433
[107,] 0.004115381 0.212327662
[108,] -0.147653850 0.004115381
[109,] 0.563248674 -0.147653850
[110,] 0.169670702 0.563248674
[111,] -0.100680341 0.169670702
[112,] 1.338645234 -0.100680341
[113,] 0.182114568 1.338645234
[114,] 0.957541739 0.182114568
[115,] -0.241010432 0.957541739
[116,] -0.829841975 -0.241010432
[117,] 1.355655266 -0.829841975
[118,] 0.822429303 1.355655266
[119,] 1.199054064 0.822429303
[120,] -0.216514063 1.199054064
[121,] -0.344039932 -0.216514063
[122,] -0.124484509 -0.344039932
[123,] 0.779673182 -0.124484509
[124,] 1.012956291 0.779673182
[125,] -0.236355757 1.012956291
[126,] -0.749041618 -0.236355757
[127,] -0.585484487 -0.749041618
[128,] -0.943165420 -0.585484487
[129,] -0.178179883 -0.943165420
[130,] -4.340352904 -0.178179883
[131,] -0.162574828 -4.340352904
[132,] 0.314306861 -0.162574828
[133,] -0.114698750 0.314306861
[134,] -1.330614045 -0.114698750
[135,] 0.412681893 -1.330614045
[136,] -0.457428101 0.412681893
[137,] 0.136508727 -0.457428101
[138,] -0.271763566 0.136508727
[139,] -0.037612159 -0.271763566
[140,] 0.986085844 -0.037612159
[141,] 0.193000369 0.986085844
[142,] -0.748648570 0.193000369
[143,] 0.237407042 -0.748648570
[144,] -0.147886858 0.237407042
[145,] -0.335149481 -0.147886858
[146,] 0.377478559 -0.335149481
[147,] -0.255142071 0.377478559
[148,] 1.396718590 -0.255142071
[149,] 0.325557728 1.396718590
[150,] -0.094435062 0.325557728
[151,] 1.365955093 -0.094435062
[152,] 1.138403646 1.365955093
[153,] 1.314540944 1.138403646
[154,] 0.506724660 1.314540944
[155,] 0.876204020 0.506724660
[156,] -0.782913969 0.876204020
[157,] 0.208579713 -0.782913969
[158,] -5.824205929 0.208579713
[159,] 0.765466644 -5.824205929
[160,] 2.760821401 0.765466644
[161,] -2.168455292 2.760821401
[162,] 0.046028164 -2.168455292
[163,] 1.995182434 0.046028164
[164,] 0.415783805 1.995182434
[165,] 2.254548464 0.415783805
[166,] 1.171211530 2.254548464
[167,] 3.098356918 1.171211530
[168,] -1.895285765 3.098356918
[169,] 2.376375092 -1.895285765
[170,] -1.022496186 2.376375092
[171,] 1.109929608 -1.022496186
[172,] 2.311651802 1.109929608
[173,] -3.469004122 2.311651802
[174,] 1.825466102 -3.469004122
[175,] -0.520887752 1.825466102
[176,] -2.709940362 -0.520887752
[177,] 2.593334964 -2.709940362
[178,] -1.055846337 2.593334964
[179,] 0.274598560 -1.055846337
[180,] 1.148491785 0.274598560
[181,] 0.752294925 1.148491785
[182,] 3.055593113 0.752294925
[183,] 2.326134965 3.055593113
[184,] 0.957100322 2.326134965
[185,] 1.105705348 0.957100322
[186,] -0.313676969 1.105705348
[187,] 0.550330324 -0.313676969
[188,] -0.207825290 0.550330324
[189,] 0.577142269 -0.207825290
[190,] 2.550293702 0.577142269
[191,] -0.339797264 2.550293702
[192,] -0.281232354 -0.339797264
[193,] 0.717014454 -0.281232354
[194,] 2.779832415 0.717014454
[195,] 0.896076309 2.779832415
[196,] -0.215497981 0.896076309
[197,] 1.493245884 -0.215497981
[198,] 2.502139719 1.493245884
[199,] 1.685623035 2.502139719
[200,] 0.651945895 1.685623035
[201,] -1.313485601 0.651945895
[202,] 1.869326349 -1.313485601
[203,] 1.806149511 1.869326349
[204,] 0.657391754 1.806149511
[205,] -1.797522786 0.657391754
[206,] -0.345509288 -1.797522786
[207,] 2.143379090 -0.345509288
[208,] 1.948615463 2.143379090
[209,] 1.370192978 1.948615463
[210,] -1.124613552 1.370192978
[211,] 1.584595355 -1.124613552
[212,] -0.671330696 1.584595355
[213,] 0.004021190 -0.671330696
[214,] 0.605110161 0.004021190
[215,] -1.329119578 0.605110161
[216,] -0.977925331 -1.329119578
[217,] -0.098259811 -0.977925331
[218,] -0.977297995 -0.098259811
[219,] -1.354198645 -0.977297995
[220,] 0.923487214 -1.354198645
[221,] 2.539229161 0.923487214
[222,] 1.125356491 2.539229161
[223,] 0.646910196 1.125356491
[224,] -3.379182231 0.646910196
[225,] 0.771004777 -3.379182231
[226,] -0.567345065 0.771004777
[227,] 1.406227299 -0.567345065
[228,] -0.768373983 1.406227299
[229,] 1.520151084 -0.768373983
[230,] -3.561230345 1.520151084
[231,] 0.225197131 -3.561230345
[232,] 1.966415975 0.225197131
[233,] -1.477965045 1.966415975
[234,] 1.604524376 -1.477965045
[235,] -0.182715239 1.604524376
[236,] 1.144341888 -0.182715239
[237,] -1.774950939 1.144341888
[238,] -0.452660946 -1.774950939
[239,] -1.184456984 -0.452660946
[240,] 0.750283451 -1.184456984
[241,] 3.270192015 0.750283451
[242,] -1.036500911 3.270192015
[243,] -0.481220976 -1.036500911
[244,] 0.549604221 -0.481220976
[245,] 1.570415221 0.549604221
[246,] 2.519795574 1.570415221
[247,] -0.291644646 2.519795574
[248,] 0.368905148 -0.291644646
[249,] -1.203899498 0.368905148
[250,] 3.034423324 -1.203899498
[251,] 3.242458975 3.034423324
[252,] -4.755072447 3.242458975
[253,] 0.493638876 -4.755072447
[254,] -0.527463208 0.493638876
[255,] -4.719079375 -0.527463208
[256,] -1.712122635 -4.719079375
[257,] -0.357566977 -1.712122635
[258,] -2.080996866 -0.357566977
[259,] -3.070058304 -2.080996866
[260,] -3.963877220 -3.070058304
[261,] 0.570982487 -3.963877220
[262,] -4.338291005 0.570982487
[263,] -0.895595933 -4.338291005
[264,] 0.246513618 -0.895595933
[265,] 0.029284312 0.246513618
[266,] 0.537130953 0.029284312
[267,] -0.112316544 0.537130953
[268,] 0.933671690 -0.112316544
[269,] 2.358914802 0.933671690
[270,] 1.189812322 2.358914802
[271,] -0.207776252 1.189812322
[272,] -5.010842972 -0.207776252
[273,] 1.213265329 -5.010842972
[274,] -1.593080158 1.213265329
[275,] -0.746789407 -1.593080158
[276,] -1.688409350 -0.746789407
[277,] -4.696848935 -1.688409350
[278,] 0.349887485 -4.696848935
[279,] 0.130011051 0.349887485
[280,] 0.409621313 0.130011051
[281,] -0.221176021 0.409621313
[282,] 0.086070227 -0.221176021
[283,] -4.416356886 0.086070227
[284,] -1.445765682 -4.416356886
[285,] 0.038876400 -1.445765682
[286,] -2.637794407 0.038876400
[287,] 0.746955607 -2.637794407
[288,] -1.102885265 0.746955607
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 0.124101316 0.353447776
2 0.095249427 0.124101316
3 0.172508056 0.095249427
4 0.581000130 0.172508056
5 0.550180223 0.581000130
6 -0.712364949 0.550180223
7 0.987571355 -0.712364949
8 0.523754888 0.987571355
9 0.533811918 0.523754888
10 -0.485513448 0.533811918
11 -0.790472529 -0.485513448
12 0.096375531 -0.790472529
13 -0.429417208 0.096375531
14 0.244142492 -0.429417208
15 -1.421196835 0.244142492
16 0.221735266 -1.421196835
17 0.600255555 0.221735266
18 0.767306941 0.600255555
19 -0.232030576 0.767306941
20 -0.118316463 -0.232030576
21 -0.806351258 -0.118316463
22 0.283438142 -0.806351258
23 -0.644035368 0.283438142
24 -0.338186557 -0.644035368
25 -0.091125263 -0.338186557
26 0.156656404 -0.091125263
27 -0.310518351 0.156656404
28 0.086982731 -0.310518351
29 -0.138997750 0.086982731
30 -0.172526852 -0.138997750
31 -0.432723930 -0.172526852
32 -0.847246355 -0.432723930
33 -0.083265197 -0.847246355
34 -0.290938990 -0.083265197
35 -0.009937046 -0.290938990
36 -0.752768921 -0.009937046
37 0.340050109 -0.752768921
38 0.839940259 0.340050109
39 -0.246043291 0.839940259
40 0.054517472 -0.246043291
41 -0.036535583 0.054517472
42 0.273080682 -0.036535583
43 0.209389089 0.273080682
44 0.901381854 0.209389089
45 -0.841620803 0.901381854
46 0.337821978 -0.841620803
47 -0.550555661 0.337821978
48 -0.203398334 -0.550555661
49 -0.197215551 -0.203398334
50 0.050665788 -0.197215551
51 1.244699734 0.050665788
52 -0.799757231 1.244699734
53 0.510692758 -0.799757231
54 0.797995505 0.510692758
55 0.125640322 0.797995505
56 -0.752461235 0.125640322
57 -1.522498142 -0.752461235
58 -0.376244684 -1.522498142
59 0.079517918 -0.376244684
60 -0.275068584 0.079517918
61 0.234961524 -0.275068584
62 0.144166545 0.234961524
63 -0.095520714 0.144166545
64 0.206229094 -0.095520714
65 0.368233803 0.206229094
66 0.269318098 0.368233803
67 -0.269562050 0.269318098
68 0.285512551 -0.269562050
69 0.550688302 0.285512551
70 -0.628713118 0.550688302
71 0.283306502 -0.628713118
72 -0.774896756 0.283306502
73 0.503157277 -0.774896756
74 -0.324245780 0.503157277
75 -0.112074798 -0.324245780
76 0.103685425 -0.112074798
77 -1.519850763 0.103685425
78 -0.140352161 -1.519850763
79 -0.043945583 -0.140352161
80 0.390707890 -0.043945583
81 -0.688930172 0.390707890
82 -0.105215498 -0.688930172
83 0.399366152 -0.105215498
84 1.242534598 0.399366152
85 0.170571206 1.242534598
86 0.141485639 0.170571206
87 0.054303828 0.141485639
88 -0.109120715 0.054303828
89 -0.069551408 -0.109120715
90 0.172524453 -0.069551408
91 -0.307441438 0.172524453
92 0.337123054 -0.307441438
93 -0.401587378 0.337123054
94 0.237216275 -0.401587378
95 0.322173494 0.237216275
96 0.449426534 0.322173494
97 -0.722147556 0.449426534
98 -2.097930133 -0.722147556
99 -2.281288338 -2.097930133
100 -0.371420974 -2.281288338
101 -0.196810165 -0.371420974
102 0.458260498 -0.196810165
103 0.529001176 0.458260498
104 0.868544092 0.529001176
105 0.117281433 0.868544092
106 0.212327662 0.117281433
107 0.004115381 0.212327662
108 -0.147653850 0.004115381
109 0.563248674 -0.147653850
110 0.169670702 0.563248674
111 -0.100680341 0.169670702
112 1.338645234 -0.100680341
113 0.182114568 1.338645234
114 0.957541739 0.182114568
115 -0.241010432 0.957541739
116 -0.829841975 -0.241010432
117 1.355655266 -0.829841975
118 0.822429303 1.355655266
119 1.199054064 0.822429303
120 -0.216514063 1.199054064
121 -0.344039932 -0.216514063
122 -0.124484509 -0.344039932
123 0.779673182 -0.124484509
124 1.012956291 0.779673182
125 -0.236355757 1.012956291
126 -0.749041618 -0.236355757
127 -0.585484487 -0.749041618
128 -0.943165420 -0.585484487
129 -0.178179883 -0.943165420
130 -4.340352904 -0.178179883
131 -0.162574828 -4.340352904
132 0.314306861 -0.162574828
133 -0.114698750 0.314306861
134 -1.330614045 -0.114698750
135 0.412681893 -1.330614045
136 -0.457428101 0.412681893
137 0.136508727 -0.457428101
138 -0.271763566 0.136508727
139 -0.037612159 -0.271763566
140 0.986085844 -0.037612159
141 0.193000369 0.986085844
142 -0.748648570 0.193000369
143 0.237407042 -0.748648570
144 -0.147886858 0.237407042
145 -0.335149481 -0.147886858
146 0.377478559 -0.335149481
147 -0.255142071 0.377478559
148 1.396718590 -0.255142071
149 0.325557728 1.396718590
150 -0.094435062 0.325557728
151 1.365955093 -0.094435062
152 1.138403646 1.365955093
153 1.314540944 1.138403646
154 0.506724660 1.314540944
155 0.876204020 0.506724660
156 -0.782913969 0.876204020
157 0.208579713 -0.782913969
158 -5.824205929 0.208579713
159 0.765466644 -5.824205929
160 2.760821401 0.765466644
161 -2.168455292 2.760821401
162 0.046028164 -2.168455292
163 1.995182434 0.046028164
164 0.415783805 1.995182434
165 2.254548464 0.415783805
166 1.171211530 2.254548464
167 3.098356918 1.171211530
168 -1.895285765 3.098356918
169 2.376375092 -1.895285765
170 -1.022496186 2.376375092
171 1.109929608 -1.022496186
172 2.311651802 1.109929608
173 -3.469004122 2.311651802
174 1.825466102 -3.469004122
175 -0.520887752 1.825466102
176 -2.709940362 -0.520887752
177 2.593334964 -2.709940362
178 -1.055846337 2.593334964
179 0.274598560 -1.055846337
180 1.148491785 0.274598560
181 0.752294925 1.148491785
182 3.055593113 0.752294925
183 2.326134965 3.055593113
184 0.957100322 2.326134965
185 1.105705348 0.957100322
186 -0.313676969 1.105705348
187 0.550330324 -0.313676969
188 -0.207825290 0.550330324
189 0.577142269 -0.207825290
190 2.550293702 0.577142269
191 -0.339797264 2.550293702
192 -0.281232354 -0.339797264
193 0.717014454 -0.281232354
194 2.779832415 0.717014454
195 0.896076309 2.779832415
196 -0.215497981 0.896076309
197 1.493245884 -0.215497981
198 2.502139719 1.493245884
199 1.685623035 2.502139719
200 0.651945895 1.685623035
201 -1.313485601 0.651945895
202 1.869326349 -1.313485601
203 1.806149511 1.869326349
204 0.657391754 1.806149511
205 -1.797522786 0.657391754
206 -0.345509288 -1.797522786
207 2.143379090 -0.345509288
208 1.948615463 2.143379090
209 1.370192978 1.948615463
210 -1.124613552 1.370192978
211 1.584595355 -1.124613552
212 -0.671330696 1.584595355
213 0.004021190 -0.671330696
214 0.605110161 0.004021190
215 -1.329119578 0.605110161
216 -0.977925331 -1.329119578
217 -0.098259811 -0.977925331
218 -0.977297995 -0.098259811
219 -1.354198645 -0.977297995
220 0.923487214 -1.354198645
221 2.539229161 0.923487214
222 1.125356491 2.539229161
223 0.646910196 1.125356491
224 -3.379182231 0.646910196
225 0.771004777 -3.379182231
226 -0.567345065 0.771004777
227 1.406227299 -0.567345065
228 -0.768373983 1.406227299
229 1.520151084 -0.768373983
230 -3.561230345 1.520151084
231 0.225197131 -3.561230345
232 1.966415975 0.225197131
233 -1.477965045 1.966415975
234 1.604524376 -1.477965045
235 -0.182715239 1.604524376
236 1.144341888 -0.182715239
237 -1.774950939 1.144341888
238 -0.452660946 -1.774950939
239 -1.184456984 -0.452660946
240 0.750283451 -1.184456984
241 3.270192015 0.750283451
242 -1.036500911 3.270192015
243 -0.481220976 -1.036500911
244 0.549604221 -0.481220976
245 1.570415221 0.549604221
246 2.519795574 1.570415221
247 -0.291644646 2.519795574
248 0.368905148 -0.291644646
249 -1.203899498 0.368905148
250 3.034423324 -1.203899498
251 3.242458975 3.034423324
252 -4.755072447 3.242458975
253 0.493638876 -4.755072447
254 -0.527463208 0.493638876
255 -4.719079375 -0.527463208
256 -1.712122635 -4.719079375
257 -0.357566977 -1.712122635
258 -2.080996866 -0.357566977
259 -3.070058304 -2.080996866
260 -3.963877220 -3.070058304
261 0.570982487 -3.963877220
262 -4.338291005 0.570982487
263 -0.895595933 -4.338291005
264 0.246513618 -0.895595933
265 0.029284312 0.246513618
266 0.537130953 0.029284312
267 -0.112316544 0.537130953
268 0.933671690 -0.112316544
269 2.358914802 0.933671690
270 1.189812322 2.358914802
271 -0.207776252 1.189812322
272 -5.010842972 -0.207776252
273 1.213265329 -5.010842972
274 -1.593080158 1.213265329
275 -0.746789407 -1.593080158
276 -1.688409350 -0.746789407
277 -4.696848935 -1.688409350
278 0.349887485 -4.696848935
279 0.130011051 0.349887485
280 0.409621313 0.130011051
281 -0.221176021 0.409621313
282 0.086070227 -0.221176021
283 -4.416356886 0.086070227
284 -1.445765682 -4.416356886
285 0.038876400 -1.445765682
286 -2.637794407 0.038876400
287 0.746955607 -2.637794407
288 -1.102885265 0.746955607
> plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals')
> lines(lowess(z))
> abline(lm(z))
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/7zshq1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/81pne1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function')
> grid()
> dev.off()
null device
1
> postscript(file="/var/www/rcomp/tmp/9jyci1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
> opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
> plot(mylm, las = 1, sub='Residual Diagnostics')
> par(opar)
> dev.off()
null device
1
> if (n > n25) {
+ postscript(file="/var/www/rcomp/tmp/10hsuy1323794347.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556)
+ plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint')
+ grid()
+ dev.off()
+ }
null device
1
>
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/rcomp/createtable")
>
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE)
> a<-table.row.end(a)
> myeq <- colnames(x)[1]
> myeq <- paste(myeq, '[t] = ', sep='')
> for (i in 1:k){
+ if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '')
+ myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ')
+ if (rownames(mysum$coefficients)[i] != '(Intercept)') {
+ myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='')
+ if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='')
+ }
+ }
> myeq <- paste(myeq, ' + e[t]')
> a<-table.row.start(a)
> a<-table.element(a, myeq)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/11uz911323794347.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Variable',header=TRUE)
> a<-table.element(a,'Parameter',header=TRUE)
> a<-table.element(a,'S.D.',header=TRUE)
> a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE)
> a<-table.element(a,'2-tail p-value',header=TRUE)
> a<-table.element(a,'1-tail p-value',header=TRUE)
> a<-table.row.end(a)
> for (i in 1:k){
+ a<-table.row.start(a)
+ a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE)
+ a<-table.element(a,mysum$coefficients[i,1])
+ a<-table.element(a, round(mysum$coefficients[i,2],6))
+ a<-table.element(a, round(mysum$coefficients[i,3],4))
+ a<-table.element(a, round(mysum$coefficients[i,4],6))
+ a<-table.element(a, round(mysum$coefficients[i,4]/2,6))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/12mf9j1323794347.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple R',1,TRUE)
> a<-table.element(a, sqrt(mysum$r.squared))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'R-squared',1,TRUE)
> a<-table.element(a, mysum$r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Adjusted R-squared',1,TRUE)
> a<-table.element(a, mysum$adj.r.squared)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (value)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[1])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[2])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE)
> a<-table.element(a, mysum$fstatistic[3])
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'p-value',1,TRUE)
> a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3]))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Residual Standard Deviation',1,TRUE)
> a<-table.element(a, mysum$sigma)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Sum Squared Residuals',1,TRUE)
> a<-table.element(a, sum(myerror*myerror))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/1360y51323794347.tab")
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a, 'Time or Index', 1, TRUE)
> a<-table.element(a, 'Actuals', 1, TRUE)
> a<-table.element(a, 'Interpolation
Forecast', 1, TRUE)
> a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE)
> a<-table.row.end(a)
> for (i in 1:n) {
+ a<-table.row.start(a)
+ a<-table.element(a,i, 1, TRUE)
+ a<-table.element(a,x[i])
+ a<-table.element(a,x[i]-mysum$resid[i])
+ a<-table.element(a,mysum$resid[i])
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/www/rcomp/tmp/14crcc1323794347.tab")
> if (n > n25) {
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'p-values',header=TRUE)
+ a<-table.element(a,'Alternative Hypothesis',3,header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'breakpoint index',header=TRUE)
+ a<-table.element(a,'greater',header=TRUE)
+ a<-table.element(a,'2-sided',header=TRUE)
+ a<-table.element(a,'less',header=TRUE)
+ a<-table.row.end(a)
+ for (mypoint in kp3:nmkm3) {
+ a<-table.row.start(a)
+ a<-table.element(a,mypoint,header=TRUE)
+ a<-table.element(a,gqarr[mypoint-kp3+1,1])
+ a<-table.element(a,gqarr[mypoint-kp3+1,2])
+ a<-table.element(a,gqarr[mypoint-kp3+1,3])
+ a<-table.row.end(a)
+ }
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/15ku3m1323794347.tab")
+ a<-table.start()
+ a<-table.row.start(a)
+ a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'Description',header=TRUE)
+ a<-table.element(a,'# significant tests',header=TRUE)
+ a<-table.element(a,'% significant tests',header=TRUE)
+ a<-table.element(a,'OK/NOK',header=TRUE)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'1% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant1)
+ a<-table.element(a,numsignificant1/numgqtests)
+ if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'5% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant5)
+ a<-table.element(a,numsignificant5/numgqtests)
+ if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,'10% type I error level',header=TRUE)
+ a<-table.element(a,numsignificant10)
+ a<-table.element(a,numsignificant10/numgqtests)
+ if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK'
+ a<-table.element(a,dum)
+ a<-table.row.end(a)
+ a<-table.end(a)
+ table.save(a,file="/var/www/rcomp/tmp/16jq021323794347.tab")
+ }
>
> try(system("convert tmp/1ctay1323794347.ps tmp/1ctay1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/2zjfc1323794347.ps tmp/2zjfc1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/3g24b1323794347.ps tmp/3g24b1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/47m681323794347.ps tmp/47m681323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/5xqre1323794347.ps tmp/5xqre1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/6zbt61323794347.ps tmp/6zbt61323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/7zshq1323794347.ps tmp/7zshq1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/81pne1323794347.ps tmp/81pne1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/9jyci1323794347.ps tmp/9jyci1323794347.png",intern=TRUE))
character(0)
> try(system("convert tmp/10hsuy1323794347.ps tmp/10hsuy1323794347.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
9.060 0.330 9.364