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(210907
+ ,30
+ ,79
+ ,112285
+ ,120982
+ ,28
+ ,58
+ ,84786
+ ,176508
+ ,38
+ ,60
+ ,83123
+ ,179321
+ ,30
+ ,108
+ ,101193
+ ,123185
+ ,22
+ ,49
+ ,38361
+ ,52746
+ ,26
+ ,0
+ ,68504
+ ,385534
+ ,25
+ ,121
+ ,119182
+ ,33170
+ ,18
+ ,1
+ ,22807
+ ,101645
+ ,11
+ ,20
+ ,17140
+ ,149061
+ ,26
+ ,43
+ ,116174
+ ,165446
+ ,25
+ ,69
+ ,57635
+ ,237213
+ ,38
+ ,78
+ ,66198
+ ,173326
+ ,44
+ ,86
+ ,71701
+ ,133131
+ ,30
+ ,44
+ ,57793
+ ,258873
+ ,40
+ ,104
+ ,80444
+ ,180083
+ ,34
+ ,63
+ ,53855
+ ,324799
+ ,47
+ ,158
+ ,97668
+ ,230964
+ ,30
+ ,102
+ ,133824
+ ,236785
+ ,31
+ ,77
+ ,101481
+ ,135473
+ ,23
+ ,82
+ ,99645
+ ,202925
+ ,36
+ ,115
+ ,114789
+ ,215147
+ ,36
+ ,101
+ ,99052
+ ,344297
+ ,30
+ ,80
+ ,67654
+ ,153935
+ ,25
+ ,50
+ ,65553
+ ,132943
+ ,39
+ ,83
+ ,97500
+ ,174724
+ ,34
+ ,123
+ ,69112
+ ,174415
+ ,31
+ ,73
+ ,82753
+ ,225548
+ ,31
+ ,81
+ ,85323
+ ,223632
+ ,33
+ ,105
+ ,72654
+ ,124817
+ ,25
+ ,47
+ ,30727
+ ,221698
+ ,33
+ ,105
+ ,77873
+ ,210767
+ ,35
+ ,94
+ ,117478
+ ,170266
+ ,42
+ ,44
+ ,74007
+ ,260561
+ ,43
+ ,114
+ ,90183
+ ,84853
+ ,30
+ ,38
+ ,61542
+ ,294424
+ ,33
+ ,107
+ ,101494
+ ,101011
+ ,13
+ ,30
+ ,27570
+ ,215641
+ ,32
+ ,71
+ ,55813
+ ,325107
+ ,36
+ ,84
+ ,79215
+ ,7176
+ ,0
+ ,0
+ ,1423
+ ,167542
+ ,28
+ ,59
+ ,55461
+ ,106408
+ ,14
+ ,33
+ ,31081
+ ,96560
+ ,17
+ ,42
+ ,22996
+ ,265769
+ ,32
+ ,96
+ ,83122
+ ,269651
+ ,30
+ ,106
+ ,70106
+ ,149112
+ ,35
+ ,56
+ ,60578
+ ,175824
+ ,20
+ ,57
+ ,39992
+ ,152871
+ ,28
+ ,59
+ ,79892
+ ,111665
+ ,28
+ ,39
+ ,49810
+ ,116408
+ ,39
+ ,34
+ ,71570
+ ,362301
+ ,34
+ ,76
+ ,100708
+ ,78800
+ ,26
+ ,20
+ ,33032
+ ,183167
+ ,39
+ ,91
+ ,82875
+ ,277965
+ ,39
+ ,115
+ ,139077
+ ,150629
+ ,33
+ ,85
+ ,71595
+ ,168809
+ ,28
+ ,76
+ ,72260
+ ,24188
+ ,4
+ ,8
+ ,5950
+ ,329267
+ ,39
+ ,79
+ ,115762
+ ,65029
+ ,18
+ ,21
+ ,32551
+ ,101097
+ ,14
+ ,30
+ ,31701
+ ,218946
+ ,29
+ ,76
+ ,80670
+ ,244052
+ ,44
+ ,101
+ ,143558
+ ,341570
+ ,21
+ ,94
+ ,117105
+ ,103597
+ ,16
+ ,27
+ ,23789
+ ,233328
+ ,28
+ ,92
+ ,120733
+ ,256462
+ ,35
+ ,123
+ ,105195
+ ,206161
+ ,28
+ ,75
+ ,73107
+ ,311473
+ ,38
+ ,128
+ ,132068
+ ,235800
+ ,23
+ ,105
+ ,149193
+ ,177939
+ ,36
+ ,55
+ ,46821
+ ,207176
+ ,32
+ ,56
+ ,87011
+ ,196553
+ ,29
+ ,41
+ ,95260
+ ,174184
+ ,25
+ ,72
+ ,55183
+ ,143246
+ ,27
+ ,67
+ ,106671
+ ,187559
+ ,36
+ ,75
+ ,73511
+ ,187681
+ ,28
+ ,114
+ ,92945
+ ,119016
+ ,23
+ ,118
+ ,78664
+ ,182192
+ ,40
+ ,77
+ ,70054
+ ,73566
+ ,23
+ ,22
+ ,22618
+ ,194979
+ ,40
+ ,66
+ ,74011
+ ,167488
+ ,28
+ ,69
+ ,83737
+ ,143756
+ ,34
+ ,105
+ ,69094
+ ,275541
+ ,33
+ ,116
+ ,93133
+ ,243199
+ ,28
+ ,88
+ ,95536
+ ,182999
+ ,34
+ ,73
+ ,225920
+ ,135649
+ ,30
+ ,99
+ ,62133
+ ,152299
+ ,33
+ ,62
+ ,61370
+ ,120221
+ ,22
+ ,53
+ ,43836
+ ,346485
+ ,38
+ ,118
+ ,106117
+ ,145790
+ ,26
+ ,30
+ ,38692
+ ,193339
+ ,35
+ ,100
+ ,84651
+ ,80953
+ ,8
+ ,49
+ ,56622
+ ,122774
+ ,24
+ ,24
+ ,15986
+ ,130585
+ ,29
+ ,67
+ ,95364
+ ,112611
+ ,20
+ ,46
+ ,26706
+ ,286468
+ ,29
+ ,57
+ ,89691
+ ,241066
+ ,45
+ ,75
+ ,67267
+ ,148446
+ ,37
+ ,135
+ ,126846
+ ,204713
+ ,33
+ ,68
+ ,41140
+ ,182079
+ ,33
+ ,124
+ ,102860
+ ,140344
+ ,25
+ ,33
+ ,51715
+ ,220516
+ ,32
+ ,98
+ ,55801
+ ,243060
+ ,29
+ ,58
+ ,111813
+ ,162765
+ ,28
+ ,68
+ ,120293
+ ,182613
+ ,28
+ ,81
+ ,138599
+ ,232138
+ ,31
+ ,131
+ ,161647
+ ,265318
+ ,52
+ ,110
+ ,115929
+ ,85574
+ ,21
+ ,37
+ ,24266
+ ,310839
+ ,24
+ ,130
+ ,162901
+ ,225060
+ ,41
+ ,93
+ ,109825
+ ,232317
+ ,33
+ ,118
+ ,129838
+ ,144966
+ ,32
+ ,39
+ ,37510
+ ,43287
+ ,19
+ ,13
+ ,43750
+ ,155754
+ ,20
+ ,74
+ ,40652
+ ,164709
+ ,31
+ ,81
+ ,87771
+ ,201940
+ ,31
+ ,109
+ ,85872
+ ,235454
+ ,32
+ ,151
+ ,89275
+ ,220801
+ ,18
+ ,51
+ ,44418
+ ,99466
+ ,23
+ ,28
+ ,192565
+ ,92661
+ ,17
+ ,40
+ ,35232
+ ,133328
+ ,20
+ ,56
+ ,40909
+ ,61361
+ ,12
+ ,27
+ ,13294
+ ,125930
+ ,17
+ ,37
+ ,32387
+ ,100750
+ ,30
+ ,83
+ ,140867
+ ,224549
+ ,31
+ ,54
+ ,120662
+ ,82316
+ ,10
+ ,27
+ ,21233
+ ,102010
+ ,13
+ ,28
+ ,44332
+ ,101523
+ ,22
+ ,59
+ ,61056
+ ,243511
+ ,42
+ ,133
+ ,101338
+ ,22938
+ ,1
+ ,12
+ ,1168
+ ,41566
+ ,9
+ ,0
+ ,13497
+ ,152474
+ ,32
+ ,106
+ ,65567
+ ,61857
+ ,11
+ ,23
+ ,25162
+ ,99923
+ ,25
+ ,44
+ ,32334
+ ,132487
+ ,36
+ ,71
+ ,40735
+ ,317394
+ ,31
+ ,116
+ ,91413
+ ,21054
+ ,0
+ ,4
+ ,855
+ ,209641
+ ,24
+ ,62
+ ,97068
+ ,22648
+ ,13
+ ,12
+ ,44339
+ ,31414
+ ,8
+ ,18
+ ,14116
+ ,46698
+ ,13
+ ,14
+ ,10288
+ ,131698
+ ,19
+ ,60
+ ,65622
+ ,91735
+ ,18
+ ,7
+ ,16563
+ ,244749
+ ,33
+ ,98
+ ,76643
+ ,184510
+ ,40
+ ,64
+ ,110681
+ ,79863
+ ,22
+ ,29
+ ,29011
+ ,128423
+ ,38
+ ,32
+ ,92696
+ ,97839
+ ,24
+ ,25
+ ,94785
+ ,38214
+ ,8
+ ,16
+ ,8773
+ ,151101
+ ,35
+ ,48
+ ,83209
+ ,272458
+ ,43
+ ,100
+ ,93815
+ ,172494
+ ,43
+ ,46
+ ,86687
+ ,108043
+ ,14
+ ,45
+ ,34553
+ ,328107
+ ,41
+ ,129
+ ,105547
+ ,250579
+ ,38
+ ,130
+ ,103487
+ ,351067
+ ,45
+ ,136
+ ,213688
+ ,158015
+ ,31
+ ,59
+ ,71220
+ ,98866
+ ,13
+ ,25
+ ,23517
+ ,85439
+ ,28
+ ,32
+ ,56926
+ ,229242
+ ,31
+ ,63
+ ,91721
+ ,351619
+ ,40
+ ,95
+ ,115168
+ ,84207
+ ,30
+ ,14
+ ,111194
+ ,120445
+ ,16
+ ,36
+ ,51009
+ ,324598
+ ,37
+ ,113
+ ,135777
+ ,131069
+ ,30
+ ,47
+ ,51513
+ ,204271
+ ,35
+ ,92
+ ,74163
+ ,165543
+ ,32
+ ,70
+ ,51633
+ ,141722
+ ,27
+ ,19
+ ,75345
+ ,116048
+ ,20
+ ,50
+ ,33416
+ ,250047
+ ,18
+ ,41
+ ,83305
+ ,299775
+ ,31
+ ,91
+ ,98952
+ ,195838
+ ,31
+ ,111
+ ,102372
+ ,173260
+ ,21
+ ,41
+ ,37238
+ ,254488
+ ,39
+ ,120
+ ,103772
+ ,104389
+ ,41
+ ,135
+ ,123969
+ ,136084
+ ,13
+ ,27
+ ,27142
+ ,199476
+ ,32
+ ,87
+ ,135400
+ ,92499
+ ,18
+ ,25
+ ,21399
+ ,224330
+ ,39
+ ,131
+ ,130115
+ ,135781
+ ,14
+ ,45
+ ,24874
+ ,74408
+ ,7
+ ,29
+ ,34988
+ ,81240
+ ,17
+ ,58
+ ,45549
+ ,14688
+ ,0
+ ,4
+ ,6023
+ ,181633
+ ,30
+ ,47
+ ,64466
+ ,271856
+ ,37
+ ,109
+ ,54990
+ ,7199
+ ,0
+ ,7
+ ,1644
+ ,46660
+ ,5
+ ,12
+ ,6179
+ ,17547
+ ,1
+ ,0
+ ,3926
+ ,133368
+ ,16
+ ,37
+ ,32755
+ ,95227
+ ,32
+ ,37
+ ,34777
+ ,152601
+ ,24
+ ,46
+ ,73224
+ ,98146
+ ,17
+ ,15
+ ,27114
+ ,79619
+ ,11
+ ,42
+ ,20760
+ ,59194
+ ,24
+ ,7
+ ,37636
+ ,139942
+ ,22
+ ,54
+ ,65461
+ ,118612
+ ,12
+ ,54
+ ,30080
+ ,72880
+ ,19
+ ,14
+ ,24094
+ ,65475
+ ,13
+ ,16
+ ,69008
+ ,99643
+ ,17
+ ,33
+ ,54968
+ ,71965
+ ,15
+ ,32
+ ,46090
+ ,77272
+ ,16
+ ,21
+ ,27507
+ ,49289
+ ,24
+ ,15
+ ,10672
+ ,135131
+ ,15
+ ,38
+ ,34029
+ ,108446
+ ,17
+ ,22
+ ,46300
+ ,89746
+ ,18
+ ,28
+ ,24760
+ ,44296
+ ,20
+ ,10
+ ,18779
+ ,77648
+ ,16
+ ,31
+ ,21280
+ ,181528
+ ,16
+ ,32
+ ,40662
+ ,134019
+ ,18
+ ,32
+ ,28987
+ ,124064
+ ,22
+ ,43
+ ,22827
+ ,92630
+ ,8
+ ,27
+ ,18513
+ ,121848
+ ,17
+ ,37
+ ,30594
+ ,52915
+ ,18
+ ,20
+ ,24006
+ ,81872
+ ,16
+ ,32
+ ,27913
+ ,58981
+ ,23
+ ,0
+ ,42744
+ ,53515
+ ,22
+ ,5
+ ,12934
+ ,60812
+ ,13
+ ,26
+ ,22574
+ ,56375
+ ,13
+ ,10
+ ,41385
+ ,65490
+ ,16
+ ,27
+ ,18653
+ ,80949
+ ,16
+ ,11
+ ,18472
+ ,76302
+ ,20
+ ,29
+ ,30976
+ ,104011
+ ,22
+ ,25
+ ,63339
+ ,98104
+ ,17
+ ,55
+ ,25568
+ ,67989
+ ,18
+ ,23
+ ,33747
+ ,30989
+ ,17
+ ,5
+ ,4154
+ ,135458
+ ,12
+ ,43
+ ,19474
+ ,73504
+ ,7
+ ,23
+ ,35130
+ ,63123
+ ,17
+ ,34
+ ,39067
+ ,61254
+ ,14
+ ,36
+ ,13310
+ ,74914
+ ,23
+ ,35
+ ,65892
+ ,31774
+ ,17
+ ,0
+ ,4143
+ ,81437
+ ,14
+ ,37
+ ,28579
+ ,87186
+ ,15
+ ,28
+ ,51776
+ ,50090
+ ,17
+ ,16
+ ,21152
+ ,65745
+ ,21
+ ,26
+ ,38084
+ ,56653
+ ,18
+ ,38
+ ,27717
+ ,158399
+ ,18
+ ,23
+ ,32928
+ ,46455
+ ,17
+ ,22
+ ,11342
+ ,73624
+ ,17
+ ,30
+ ,19499
+ ,38395
+ ,16
+ ,16
+ ,16380
+ ,91899
+ ,15
+ ,18
+ ,36874
+ ,139526
+ ,21
+ ,28
+ ,48259
+ ,52164
+ ,16
+ ,32
+ ,16734
+ ,51567
+ ,14
+ ,21
+ ,28207
+ ,70551
+ ,15
+ ,23
+ ,30143
+ ,84856
+ ,17
+ ,29
+ ,41369
+ ,102538
+ ,15
+ ,50
+ ,45833
+ ,86678
+ ,15
+ ,12
+ ,29156
+ ,85709
+ ,10
+ ,21
+ ,35944
+ ,34662
+ ,6
+ ,18
+ ,36278
+ ,150580
+ ,22
+ ,27
+ ,45588
+ ,99611
+ ,21
+ ,41
+ ,45097
+ ,19349
+ ,1
+ ,13
+ ,3895
+ ,99373
+ ,18
+ ,12
+ ,28394
+ ,86230
+ ,17
+ ,21
+ ,18632
+ ,30837
+ ,4
+ ,8
+ ,2325
+ ,31706
+ ,10
+ ,26
+ ,25139
+ ,89806
+ ,16
+ ,27
+ ,27975
+ ,62088
+ ,16
+ ,13
+ ,14483
+ ,40151
+ ,9
+ ,16
+ ,13127
+ ,27634
+ ,16
+ ,2
+ ,5839
+ ,76990
+ ,17
+ ,42
+ ,24069
+ ,37460
+ ,7
+ ,5
+ ,3738
+ ,54157
+ ,15
+ ,37
+ ,18625
+ ,49862
+ ,14
+ ,17
+ ,36341
+ ,84337
+ ,14
+ ,38
+ ,24548
+ ,64175
+ ,18
+ ,37
+ ,21792
+ ,59382
+ ,12
+ ,29
+ ,26263
+ ,119308
+ ,16
+ ,32
+ ,23686
+ ,76702
+ ,21
+ ,35
+ ,49303
+ ,103425
+ ,19
+ ,17
+ ,25659
+ ,70344
+ ,16
+ ,20
+ ,28904
+ ,43410
+ ,1
+ ,7
+ ,2781
+ ,104838
+ ,16
+ ,46
+ ,29236
+ ,62215
+ ,10
+ ,24
+ ,19546
+ ,69304
+ ,19
+ ,40
+ ,22818
+ ,53117
+ ,12
+ ,3
+ ,32689
+ ,19764
+ ,2
+ ,10
+ ,5752
+ ,86680
+ ,14
+ ,37
+ ,22197
+ ,84105
+ ,17
+ ,17
+ ,20055
+ ,77945
+ ,19
+ ,28
+ ,25272
+ ,89113
+ ,14
+ ,19
+ ,82206
+ ,91005
+ ,11
+ ,29
+ ,32073
+ ,40248
+ ,4
+ ,8
+ ,5444
+ ,64187
+ ,16
+ ,10
+ ,20154
+ ,50857
+ ,20
+ ,15
+ ,36944
+ ,56613
+ ,12
+ ,15
+ ,8019
+ ,62792
+ ,15
+ ,28
+ ,30884
+ ,72535
+ ,16
+ ,17
+ ,19540)
+ ,dim=c(4
+ ,289)
+ ,dimnames=list(c('time'
+ ,'feedback'
+ ,'blogs'
+ ,'dgddg')
+ ,1:289))
> y <- array(NA,dim=c(4,289),dimnames=list(c('time','feedback','blogs','dgddg'),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 = '1'
> #'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
time feedback blogs dgddg
1 210907 30 79 112285
2 120982 28 58 84786
3 176508 38 60 83123
4 179321 30 108 101193
5 123185 22 49 38361
6 52746 26 0 68504
7 385534 25 121 119182
8 33170 18 1 22807
9 101645 11 20 17140
10 149061 26 43 116174
11 165446 25 69 57635
12 237213 38 78 66198
13 173326 44 86 71701
14 133131 30 44 57793
15 258873 40 104 80444
16 180083 34 63 53855
17 324799 47 158 97668
18 230964 30 102 133824
19 236785 31 77 101481
20 135473 23 82 99645
21 202925 36 115 114789
22 215147 36 101 99052
23 344297 30 80 67654
24 153935 25 50 65553
25 132943 39 83 97500
26 174724 34 123 69112
27 174415 31 73 82753
28 225548 31 81 85323
29 223632 33 105 72654
30 124817 25 47 30727
31 221698 33 105 77873
32 210767 35 94 117478
33 170266 42 44 74007
34 260561 43 114 90183
35 84853 30 38 61542
36 294424 33 107 101494
37 101011 13 30 27570
38 215641 32 71 55813
39 325107 36 84 79215
40 7176 0 0 1423
41 167542 28 59 55461
42 106408 14 33 31081
43 96560 17 42 22996
44 265769 32 96 83122
45 269651 30 106 70106
46 149112 35 56 60578
47 175824 20 57 39992
48 152871 28 59 79892
49 111665 28 39 49810
50 116408 39 34 71570
51 362301 34 76 100708
52 78800 26 20 33032
53 183167 39 91 82875
54 277965 39 115 139077
55 150629 33 85 71595
56 168809 28 76 72260
57 24188 4 8 5950
58 329267 39 79 115762
59 65029 18 21 32551
60 101097 14 30 31701
61 218946 29 76 80670
62 244052 44 101 143558
63 341570 21 94 117105
64 103597 16 27 23789
65 233328 28 92 120733
66 256462 35 123 105195
67 206161 28 75 73107
68 311473 38 128 132068
69 235800 23 105 149193
70 177939 36 55 46821
71 207176 32 56 87011
72 196553 29 41 95260
73 174184 25 72 55183
74 143246 27 67 106671
75 187559 36 75 73511
76 187681 28 114 92945
77 119016 23 118 78664
78 182192 40 77 70054
79 73566 23 22 22618
80 194979 40 66 74011
81 167488 28 69 83737
82 143756 34 105 69094
83 275541 33 116 93133
84 243199 28 88 95536
85 182999 34 73 225920
86 135649 30 99 62133
87 152299 33 62 61370
88 120221 22 53 43836
89 346485 38 118 106117
90 145790 26 30 38692
91 193339 35 100 84651
92 80953 8 49 56622
93 122774 24 24 15986
94 130585 29 67 95364
95 112611 20 46 26706
96 286468 29 57 89691
97 241066 45 75 67267
98 148446 37 135 126846
99 204713 33 68 41140
100 182079 33 124 102860
101 140344 25 33 51715
102 220516 32 98 55801
103 243060 29 58 111813
104 162765 28 68 120293
105 182613 28 81 138599
106 232138 31 131 161647
107 265318 52 110 115929
108 85574 21 37 24266
109 310839 24 130 162901
110 225060 41 93 109825
111 232317 33 118 129838
112 144966 32 39 37510
113 43287 19 13 43750
114 155754 20 74 40652
115 164709 31 81 87771
116 201940 31 109 85872
117 235454 32 151 89275
118 220801 18 51 44418
119 99466 23 28 192565
120 92661 17 40 35232
121 133328 20 56 40909
122 61361 12 27 13294
123 125930 17 37 32387
124 100750 30 83 140867
125 224549 31 54 120662
126 82316 10 27 21233
127 102010 13 28 44332
128 101523 22 59 61056
129 243511 42 133 101338
130 22938 1 12 1168
131 41566 9 0 13497
132 152474 32 106 65567
133 61857 11 23 25162
134 99923 25 44 32334
135 132487 36 71 40735
136 317394 31 116 91413
137 21054 0 4 855
138 209641 24 62 97068
139 22648 13 12 44339
140 31414 8 18 14116
141 46698 13 14 10288
142 131698 19 60 65622
143 91735 18 7 16563
144 244749 33 98 76643
145 184510 40 64 110681
146 79863 22 29 29011
147 128423 38 32 92696
148 97839 24 25 94785
149 38214 8 16 8773
150 151101 35 48 83209
151 272458 43 100 93815
152 172494 43 46 86687
153 108043 14 45 34553
154 328107 41 129 105547
155 250579 38 130 103487
156 351067 45 136 213688
157 158015 31 59 71220
158 98866 13 25 23517
159 85439 28 32 56926
160 229242 31 63 91721
161 351619 40 95 115168
162 84207 30 14 111194
163 120445 16 36 51009
164 324598 37 113 135777
165 131069 30 47 51513
166 204271 35 92 74163
167 165543 32 70 51633
168 141722 27 19 75345
169 116048 20 50 33416
170 250047 18 41 83305
171 299775 31 91 98952
172 195838 31 111 102372
173 173260 21 41 37238
174 254488 39 120 103772
175 104389 41 135 123969
176 136084 13 27 27142
177 199476 32 87 135400
178 92499 18 25 21399
179 224330 39 131 130115
180 135781 14 45 24874
181 74408 7 29 34988
182 81240 17 58 45549
183 14688 0 4 6023
184 181633 30 47 64466
185 271856 37 109 54990
186 7199 0 7 1644
187 46660 5 12 6179
188 17547 1 0 3926
189 133368 16 37 32755
190 95227 32 37 34777
191 152601 24 46 73224
192 98146 17 15 27114
193 79619 11 42 20760
194 59194 24 7 37636
195 139942 22 54 65461
196 118612 12 54 30080
197 72880 19 14 24094
198 65475 13 16 69008
199 99643 17 33 54968
200 71965 15 32 46090
201 77272 16 21 27507
202 49289 24 15 10672
203 135131 15 38 34029
204 108446 17 22 46300
205 89746 18 28 24760
206 44296 20 10 18779
207 77648 16 31 21280
208 181528 16 32 40662
209 134019 18 32 28987
210 124064 22 43 22827
211 92630 8 27 18513
212 121848 17 37 30594
213 52915 18 20 24006
214 81872 16 32 27913
215 58981 23 0 42744
216 53515 22 5 12934
217 60812 13 26 22574
218 56375 13 10 41385
219 65490 16 27 18653
220 80949 16 11 18472
221 76302 20 29 30976
222 104011 22 25 63339
223 98104 17 55 25568
224 67989 18 23 33747
225 30989 17 5 4154
226 135458 12 43 19474
227 73504 7 23 35130
228 63123 17 34 39067
229 61254 14 36 13310
230 74914 23 35 65892
231 31774 17 0 4143
232 81437 14 37 28579
233 87186 15 28 51776
234 50090 17 16 21152
235 65745 21 26 38084
236 56653 18 38 27717
237 158399 18 23 32928
238 46455 17 22 11342
239 73624 17 30 19499
240 38395 16 16 16380
241 91899 15 18 36874
242 139526 21 28 48259
243 52164 16 32 16734
244 51567 14 21 28207
245 70551 15 23 30143
246 84856 17 29 41369
247 102538 15 50 45833
248 86678 15 12 29156
249 85709 10 21 35944
250 34662 6 18 36278
251 150580 22 27 45588
252 99611 21 41 45097
253 19349 1 13 3895
254 99373 18 12 28394
255 86230 17 21 18632
256 30837 4 8 2325
257 31706 10 26 25139
258 89806 16 27 27975
259 62088 16 13 14483
260 40151 9 16 13127
261 27634 16 2 5839
262 76990 17 42 24069
263 37460 7 5 3738
264 54157 15 37 18625
265 49862 14 17 36341
266 84337 14 38 24548
267 64175 18 37 21792
268 59382 12 29 26263
269 119308 16 32 23686
270 76702 21 35 49303
271 103425 19 17 25659
272 70344 16 20 28904
273 43410 1 7 2781
274 104838 16 46 29236
275 62215 10 24 19546
276 69304 19 40 22818
277 53117 12 3 32689
278 19764 2 10 5752
279 86680 14 37 22197
280 84105 17 17 20055
281 77945 19 28 25272
282 89113 14 19 82206
283 91005 11 29 32073
284 40248 4 8 5444
285 64187 16 10 20154
286 50857 20 15 36944
287 56613 12 15 8019
288 62792 15 28 30884
289 72535 16 17 19540
> k <- length(x[1,])
> df <- as.data.frame(x)
> (mylm <- lm(df))
Call:
lm(formula = df)
Coefficients:
(Intercept) feedback blogs dgddg
1.470e+04 1.774e+03 1.148e+03 3.816e-01
> (mysum <- summary(mylm))
Call:
lm(formula = df)
Residuals:
Min 1Q Median 3Q Max
-185310 -22442 -5975 18293 161620
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 1.470e+04 6.091e+03 2.413 0.016437 *
feedback 1.774e+03 3.899e+02 4.549 7.97e-06 ***
blogs 1.148e+03 1.170e+02 9.811 < 2e-16 ***
dgddg 3.816e-01 1.019e-01 3.747 0.000217 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 42190 on 285 degrees of freedom
Multiple R-squared: 0.7402, Adjusted R-squared: 0.7375
F-statistic: 270.7 on 3 and 285 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.9633273 7.334541e-02 3.667271e-02
[2,] 0.9257279 1.485443e-01 7.427213e-02
[3,] 0.8768804 2.462392e-01 1.231196e-01
[4,] 0.8200441 3.599118e-01 1.799559e-01
[5,] 0.7383030 5.233941e-01 2.616970e-01
[6,] 0.8112905 3.774190e-01 1.887095e-01
[7,] 0.7594717 4.810566e-01 2.405283e-01
[8,] 0.6868835 6.262331e-01 3.131165e-01
[9,] 0.6200817 7.598367e-01 3.799183e-01
[10,] 0.5574093 8.851814e-01 4.425907e-01
[11,] 0.4975217 9.950434e-01 5.024783e-01
[12,] 0.4600968 9.201937e-01 5.399032e-01
[13,] 0.4376141 8.752281e-01 5.623859e-01
[14,] 0.6120534 7.758932e-01 3.879466e-01
[15,] 0.6763612 6.472777e-01 3.236388e-01
[16,] 0.6224309 7.551381e-01 3.775691e-01
[17,] 0.9700739 5.985224e-02 2.992612e-02
[18,] 0.9581349 8.373013e-02 4.186507e-02
[19,] 0.9710654 5.786914e-02 2.893457e-02
[20,] 0.9904115 1.917696e-02 9.588478e-03
[21,] 0.9861807 2.763861e-02 1.381930e-02
[22,] 0.9827731 3.445381e-02 1.722691e-02
[23,] 0.9761684 4.766321e-02 2.383160e-02
[24,] 0.9675975 6.480504e-02 3.240252e-02
[25,] 0.9570187 8.596255e-02 4.298127e-02
[26,] 0.9441225 1.117549e-01 5.587747e-02
[27,] 0.9387014 1.225971e-01 6.129856e-02
[28,] 0.9223384 1.553233e-01 7.766164e-02
[29,] 0.9193413 1.613174e-01 8.065870e-02
[30,] 0.9272756 1.454487e-01 7.272436e-02
[31,] 0.9087922 1.824156e-01 9.120778e-02
[32,] 0.9046413 1.907174e-01 9.535870e-02
[33,] 0.9805305 3.893892e-02 1.946946e-02
[34,] 0.9775209 4.495829e-02 2.247914e-02
[35,] 0.9707514 5.849715e-02 2.924858e-02
[36,] 0.9622262 7.554764e-02 3.777382e-02
[37,] 0.9535637 9.287251e-02 4.643625e-02
[38,] 0.9528895 9.422097e-02 4.711048e-02
[39,] 0.9493447 1.013107e-01 5.065533e-02
[40,] 0.9369971 1.260059e-01 6.300293e-02
[41,] 0.9299416 1.401169e-01 7.005844e-02
[42,] 0.9142877 1.714246e-01 8.571230e-02
[43,] 0.8971648 2.056705e-01 1.028352e-01
[44,] 0.8797056 2.405888e-01 1.202944e-01
[45,] 0.9940341 1.193187e-02 5.965935e-03
[46,] 0.9921493 1.570142e-02 7.850710e-03
[47,] 0.9916068 1.678635e-02 8.393177e-03
[48,] 0.9889503 2.209941e-02 1.104971e-02
[49,] 0.9902236 1.955285e-02 9.776423e-03
[50,] 0.9877638 2.447239e-02 1.223619e-02
[51,] 0.9849806 3.003879e-02 1.501939e-02
[52,] 0.9964427 7.114664e-03 3.557332e-03
[53,] 0.9954115 9.177044e-03 4.588522e-03
[54,] 0.9939692 1.206169e-02 6.030843e-03
[55,] 0.9928776 1.424489e-02 7.122443e-03
[56,] 0.9914190 1.716195e-02 8.580977e-03
[57,] 0.9984559 3.088292e-03 1.544146e-03
[58,] 0.9979859 4.028297e-03 2.014148e-03
[59,] 0.9974259 5.148169e-03 2.574085e-03
[60,] 0.9968412 6.317583e-03 3.158792e-03
[61,] 0.9960796 7.840718e-03 3.920359e-03
[62,] 0.9952936 9.412861e-03 4.706431e-03
[63,] 0.9949850 1.003003e-02 5.015017e-03
[64,] 0.9939767 1.204652e-02 6.023260e-03
[65,] 0.9935483 1.290335e-02 6.451673e-03
[66,] 0.9936621 1.267576e-02 6.337879e-03
[67,] 0.9918799 1.624014e-02 8.120070e-03
[68,] 0.9922839 1.543223e-02 7.716117e-03
[69,] 0.9900850 1.982993e-02 9.914966e-03
[70,] 0.9920044 1.599126e-02 7.995630e-03
[71,] 0.9983734 3.253165e-03 1.626582e-03
[72,] 0.9979053 4.189477e-03 2.094739e-03
[73,] 0.9972958 5.408389e-03 2.704194e-03
[74,] 0.9964460 7.107905e-03 3.553952e-03
[75,] 0.9954686 9.062800e-03 4.531400e-03
[76,] 0.9975266 4.946870e-03 2.473435e-03
[77,] 0.9971925 5.615033e-03 2.807516e-03
[78,] 0.9969937 6.012553e-03 3.006277e-03
[79,] 0.9982258 3.548397e-03 1.774198e-03
[80,] 0.9989451 2.109741e-03 1.054870e-03
[81,] 0.9986283 2.743328e-03 1.371664e-03
[82,] 0.9982163 3.567464e-03 1.783732e-03
[83,] 0.9992817 1.436545e-03 7.182723e-04
[84,] 0.9992102 1.579503e-03 7.897516e-04
[85,] 0.9990926 1.814731e-03 9.073656e-04
[86,] 0.9989594 2.081250e-03 1.040625e-03
[87,] 0.9988095 2.380928e-03 1.190464e-03
[88,] 0.9989193 2.161383e-03 1.080692e-03
[89,] 0.9985568 2.886433e-03 1.443216e-03
[90,] 0.9998211 3.577045e-04 1.788522e-04
[91,] 0.9998047 3.906925e-04 1.953463e-04
[92,] 0.9999918 1.648890e-05 8.244452e-06
[93,] 0.9999911 1.785267e-05 8.926334e-06
[94,] 0.9999956 8.884210e-06 4.442105e-06
[95,] 0.9999941 1.184290e-05 5.921452e-06
[96,] 0.9999917 1.652321e-05 8.261606e-06
[97,] 0.9999954 9.242968e-06 4.621484e-06
[98,] 0.9999941 1.175476e-05 5.877380e-06
[99,] 0.9999926 1.470804e-05 7.354018e-06
[100,] 0.9999934 1.328948e-05 6.644742e-06
[101,] 0.9999904 1.923762e-05 9.618808e-06
[102,] 0.9999871 2.579673e-05 1.289837e-05
[103,] 0.9999872 2.555532e-05 1.277766e-05
[104,] 0.9999817 3.652449e-05 1.826224e-05
[105,] 0.9999765 4.703084e-05 2.351542e-05
[106,] 0.9999677 6.456489e-05 3.228244e-05
[107,] 0.9999660 6.803971e-05 3.401985e-05
[108,] 0.9999517 9.654824e-05 4.827412e-05
[109,] 0.9999429 1.141246e-04 5.706231e-05
[110,] 0.9999281 1.438702e-04 7.193509e-05
[111,] 0.9999277 1.446899e-04 7.234493e-05
[112,] 0.9999843 3.132663e-05 1.566331e-05
[113,] 0.9999907 1.850725e-05 9.253624e-06
[114,] 0.9999871 2.586945e-05 1.293473e-05
[115,] 0.9999813 3.732450e-05 1.866225e-05
[116,] 0.9999742 5.162130e-05 2.581065e-05
[117,] 0.9999675 6.494151e-05 3.247075e-05
[118,] 0.9999981 3.764178e-06 1.882089e-06
[119,] 0.9999982 3.671061e-06 1.835531e-06
[120,] 0.9999973 5.403397e-06 2.701699e-06
[121,] 0.9999961 7.820873e-06 3.910437e-06
[122,] 0.9999964 7.256887e-06 3.628443e-06
[123,] 0.9999959 8.128161e-06 4.064080e-06
[124,] 0.9999942 1.160523e-05 5.802613e-06
[125,] 0.9999915 1.691928e-05 8.459638e-06
[126,] 0.9999953 9.371661e-06 4.685831e-06
[127,] 0.9999933 1.345890e-05 6.729449e-06
[128,] 0.9999911 1.781234e-05 8.906169e-06
[129,] 0.9999910 1.797880e-05 8.989398e-06
[130,] 0.9999967 6.635661e-06 3.317831e-06
[131,] 0.9999951 9.828142e-06 4.914071e-06
[132,] 0.9999951 9.811005e-06 4.905502e-06
[133,] 0.9999958 8.416433e-06 4.208217e-06
[134,] 0.9999946 1.079116e-05 5.395581e-06
[135,] 0.9999923 1.546216e-05 7.731081e-06
[136,] 0.9999891 2.186850e-05 1.093425e-05
[137,] 0.9999874 2.514843e-05 1.257421e-05
[138,] 0.9999852 2.962176e-05 1.481088e-05
[139,] 0.9999799 4.024987e-05 2.012494e-05
[140,] 0.9999728 5.443902e-05 2.721951e-05
[141,] 0.9999663 6.730285e-05 3.365143e-05
[142,] 0.9999603 7.940902e-05 3.970451e-05
[143,] 0.9999453 1.093132e-04 5.465661e-05
[144,] 0.9999256 1.488948e-04 7.444742e-05
[145,] 0.9999165 1.670369e-04 8.351843e-05
[146,] 0.9998829 2.341095e-04 1.170547e-04
[147,] 0.9998372 3.256524e-04 1.628262e-04
[148,] 0.9998808 2.383041e-04 1.191520e-04
[149,] 0.9998417 3.165333e-04 1.582666e-04
[150,] 0.9997895 4.210215e-04 2.105107e-04
[151,] 0.9997111 5.777521e-04 2.888760e-04
[152,] 0.9996440 7.119719e-04 3.559859e-04
[153,] 0.9996270 7.459820e-04 3.729910e-04
[154,] 0.9996886 6.227039e-04 3.113519e-04
[155,] 0.9999771 4.576605e-05 2.288302e-05
[156,] 0.9999818 3.640790e-05 1.820395e-05
[157,] 0.9999749 5.024912e-05 2.512456e-05
[158,] 0.9999891 2.189695e-05 1.094848e-05
[159,] 0.9999840 3.193549e-05 1.596774e-05
[160,] 0.9999776 4.479505e-05 2.239752e-05
[161,] 0.9999684 6.321353e-05 3.160677e-05
[162,] 0.9999612 7.769866e-05 3.884933e-05
[163,] 0.9999445 1.109972e-04 5.549862e-05
[164,] 0.9999987 2.506967e-06 1.253484e-06
[165,] 1.0000000 9.460124e-08 4.730062e-08
[166,] 0.9999999 1.331074e-07 6.655368e-08
[167,] 1.0000000 4.346300e-08 2.173150e-08
[168,] 1.0000000 4.908674e-08 2.454337e-08
[169,] 1.0000000 8.367739e-13 4.183870e-13
[170,] 1.0000000 2.237062e-13 1.118531e-13
[171,] 1.0000000 3.804879e-13 1.902439e-13
[172,] 1.0000000 6.934349e-13 3.467174e-13
[173,] 1.0000000 1.239825e-13 6.199126e-14
[174,] 1.0000000 1.141113e-13 5.705567e-14
[175,] 1.0000000 2.330708e-13 1.165354e-13
[176,] 1.0000000 8.937278e-14 4.468639e-14
[177,] 1.0000000 1.812949e-13 9.064746e-14
[178,] 1.0000000 1.750911e-13 8.754557e-14
[179,] 1.0000000 7.559650e-14 3.779825e-14
[180,] 1.0000000 1.258091e-13 6.290454e-14
[181,] 1.0000000 2.529540e-13 1.264770e-13
[182,] 1.0000000 5.195727e-13 2.597864e-13
[183,] 1.0000000 3.959381e-13 1.979690e-13
[184,] 1.0000000 6.174988e-13 3.087494e-13
[185,] 1.0000000 1.062105e-12 5.310523e-13
[186,] 1.0000000 1.201499e-12 6.007495e-13
[187,] 1.0000000 2.355531e-12 1.177766e-12
[188,] 1.0000000 4.160587e-12 2.080294e-12
[189,] 1.0000000 8.378632e-12 4.189316e-12
[190,] 1.0000000 1.508320e-11 7.541601e-12
[191,] 1.0000000 2.950234e-11 1.475117e-11
[192,] 1.0000000 4.255835e-11 2.127918e-11
[193,] 1.0000000 8.318434e-11 4.159217e-11
[194,] 1.0000000 1.160628e-10 5.803142e-11
[195,] 1.0000000 2.248085e-10 1.124043e-10
[196,] 1.0000000 3.532603e-10 1.766301e-10
[197,] 1.0000000 2.474879e-10 1.237440e-10
[198,] 1.0000000 3.657206e-10 1.828603e-10
[199,] 1.0000000 6.794805e-10 3.397403e-10
[200,] 1.0000000 1.051239e-09 5.256193e-10
[201,] 1.0000000 1.997308e-09 9.986542e-10
[202,] 1.0000000 1.834332e-11 9.171661e-12
[203,] 1.0000000 6.875563e-12 3.437781e-12
[204,] 1.0000000 8.130955e-12 4.065478e-12
[205,] 1.0000000 8.190545e-12 4.095272e-12
[206,] 1.0000000 6.823110e-12 3.411555e-12
[207,] 1.0000000 1.097848e-11 5.489242e-12
[208,] 1.0000000 2.335513e-11 1.167756e-11
[209,] 1.0000000 4.166897e-11 2.083448e-11
[210,] 1.0000000 8.416908e-11 4.208454e-11
[211,] 1.0000000 1.675959e-10 8.379795e-11
[212,] 1.0000000 3.183636e-10 1.591818e-10
[213,] 1.0000000 6.344778e-10 3.172389e-10
[214,] 1.0000000 9.473124e-10 4.736562e-10
[215,] 1.0000000 1.801723e-09 9.008616e-10
[216,] 1.0000000 3.666680e-09 1.833340e-09
[217,] 1.0000000 7.174948e-09 3.587474e-09
[218,] 1.0000000 1.286365e-08 6.431824e-09
[219,] 1.0000000 1.999768e-08 9.998842e-09
[220,] 1.0000000 2.686172e-09 1.343086e-09
[221,] 1.0000000 5.010112e-09 2.505056e-09
[222,] 1.0000000 6.400808e-09 3.200404e-09
[223,] 1.0000000 1.242744e-08 6.213718e-09
[224,] 1.0000000 5.774308e-09 2.887154e-09
[225,] 1.0000000 9.135014e-09 4.567507e-09
[226,] 1.0000000 1.864724e-08 9.323622e-09
[227,] 1.0000000 3.791390e-08 1.895695e-08
[228,] 1.0000000 5.610118e-08 2.805059e-08
[229,] 1.0000000 5.719475e-08 2.859738e-08
[230,] 1.0000000 5.508622e-08 2.754311e-08
[231,] 1.0000000 8.407517e-10 4.203759e-10
[232,] 1.0000000 1.225375e-09 6.126876e-10
[233,] 1.0000000 2.817013e-09 1.408507e-09
[234,] 1.0000000 2.782753e-09 1.391377e-09
[235,] 1.0000000 5.178761e-09 2.589380e-09
[236,] 1.0000000 2.294107e-09 1.147053e-09
[237,] 1.0000000 2.974851e-09 1.487426e-09
[238,] 1.0000000 4.743892e-09 2.371946e-09
[239,] 1.0000000 1.115982e-08 5.579910e-09
[240,] 1.0000000 2.623947e-08 1.311974e-08
[241,] 1.0000000 5.701166e-08 2.850583e-08
[242,] 1.0000000 9.470399e-08 4.735200e-08
[243,] 0.9999999 1.275306e-07 6.376529e-08
[244,] 0.9999999 2.144321e-07 1.072160e-07
[245,] 1.0000000 1.441051e-08 7.205254e-09
[246,] 1.0000000 3.671126e-08 1.835563e-08
[247,] 1.0000000 7.811165e-08 3.905582e-08
[248,] 1.0000000 5.583792e-08 2.791896e-08
[249,] 1.0000000 9.224648e-08 4.612324e-08
[250,] 0.9999999 2.340378e-07 1.170189e-07
[251,] 0.9999999 1.249486e-07 6.247432e-08
[252,] 0.9999999 2.321923e-07 1.160961e-07
[253,] 0.9999997 5.987794e-07 2.993897e-07
[254,] 0.9999994 1.245337e-06 6.226685e-07
[255,] 0.9999992 1.678587e-06 8.392934e-07
[256,] 0.9999980 3.949302e-06 1.974651e-06
[257,] 0.9999951 9.808251e-06 4.904125e-06
[258,] 0.9999948 1.038637e-05 5.193186e-06
[259,] 0.9999918 1.645555e-05 8.227776e-06
[260,] 0.9999791 4.173944e-05 2.086972e-05
[261,] 0.9999733 5.331013e-05 2.665507e-05
[262,] 0.9999506 9.883743e-05 4.941872e-05
[263,] 0.9999869 2.616793e-05 1.308396e-05
[264,] 0.9999812 3.769588e-05 1.884794e-05
[265,] 0.9999944 1.118467e-05 5.592336e-06
[266,] 0.9999810 3.790612e-05 1.895306e-05
[267,] 0.9999483 1.034654e-04 5.173269e-05
[268,] 0.9998816 2.368044e-04 1.184022e-04
[269,] 0.9996250 7.499315e-04 3.749657e-04
[270,] 0.9995306 9.387744e-04 4.693872e-04
[271,] 0.9985531 2.893857e-03 1.446928e-03
[272,] 0.9983878 3.224401e-03 1.612200e-03
[273,] 0.9946113 1.077735e-02 5.388674e-03
[274,] 0.9932989 1.340220e-02 6.701101e-03
[275,] 0.9769704 4.605921e-02 2.302960e-02
[276,] 0.9370958 1.258085e-01 6.290424e-02
> postscript(file="/var/wessaorg/rcomp/tmp/17m4g1324127192.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/2kdvx1324127192.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/3650n1324127192.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/469sk1324127192.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/5i1h81324127192.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
9459.55117 -42318.42267 -6191.83910 -51181.44191 -1424.77691
6 7 8 9 10
-34217.21900 142113.49876 -23311.33777 37933.80275 -5452.25711
11 12 13 14 15
5202.08209 40310.65109 -45502.31467 -7345.85894 23141.98373
16 17 18 19 20
12204.28474 8093.56625 -5104.05047 39982.50566 -52177.46897
21 22 23 24 25
-51444.03012 -17146.35953 158733.95290 12478.74058 -83418.08288
26 27 28 29 30
-67848.80908 -10648.99555 30320.34322 2142.80883 94.76229
31 32 33 34 35
-1782.88705 -18749.30410 2315.48642 4313.23727 -50167.40189
36 37 38 39 40
59633.04958 18293.18762 41379.84449 119897.55241 -8067.82469
41 42 43 44 45
14284.84145 17132.88928 -5281.92931 52389.55026 53307.82657
46 47 48 49 50
-15071.18073 44956.57391 -9709.61554 -16478.38248 -33812.38589
51 52 53 54 55
161619.86137 -17583.46955 -36795.71870 9005.58453 -47498.83271
56 57 58 59 60
-10372.69665 -9061.65831 110528.14257 -18128.10104 15028.86516
61 62 63 64 65
34781.01161 -19416.58746 137029.71115 20444.19031 17282.06624
66 67 68 69 70
-1654.78295 27803.92871 32040.02889 2840.02156 18379.84862
71 72 73 74 75
38226.85155 46995.24938 11432.24145 -36963.17969 -5142.91568
76 77 78 79 80
-43013.30532 -101950.60841 -18581.69736 -15817.45803 5321.68653
81 82 83 84 85
-8038.56928 -78148.44311 33610.05983 41360.29097 -62022.46549
86 87 88 89 90
-69616.45843 -15525.92548 -11069.61226 88434.16372 35767.92906
91 92 93 94 95
-30536.89723 -25791.94407 31851.91723 -48856.82370 -559.69883
96 97 98 99 100
120669.73879 34782.44303 -135255.32333 37721.16333 -72746.88468
101 102 103 104 105
23682.28990 15267.17887 67671.59079 -25564.35686 -27624.55868
106 107 108 109 110
-49609.77672 -12128.11343 -18108.62133 42177.36240 -11030.87558
111 112 113 114 115
-25917.17252 14421.25911 -36734.85479 5121.06623 -31452.87244
116 117 118 119 120
-25637.27442 -43405.93363 98679.33854 -61660.48481 -11554.75924
121 122 123 124 125
3258.48559 -10691.33122 26243.54396 -116196.47094 46827.37628
126 127 128 129 130
10781.62257 15191.12727 -43226.34386 -37029.30306 -7756.67525
131 132 133 134 135
5749.94481 -65684.64897 -8359.16811 -21968.92476 -43115.36233
136 137 138 139 140
79667.12025 1435.49405 44157.36125 -45807.76922 -23525.94368
141 142 143 144 145
-11058.81079 -10620.19957 30749.35717 29772.53703 -16845.74526
146 147 148 149 150
-18221.37351 -25790.02002 -24302.53665 -12391.20456 -12535.82700
151 152 153 154 155
30894.23177 -4365.06184 3668.55944 52325.71859 -20242.49620
156 157 158 159 160
18886.15233 -6577.67244 23434.21451 -37384.99266 52234.20934
161 162 163 164 165
112967.21883 -42213.13655 16573.63549 62741.33280 -10454.84274
166 167 168 169 170
-6419.53683 -5975.10608 28564.89972 -4274.83988 124563.73369
171 172 173 174 175
87867.57984 -40331.79298 60035.50787 -6737.48329 -185309.72510
176 177 178 179 180
56973.10550 -23523.24308 9006.32656 -59575.07568 35100.29855
181 182 183 184 185
650.16261 -47574.47123 -6902.73900 35165.98104 45421.03807
186 187 188 189 190
-16164.19011 6957.67253 -425.86320 35314.93984 -31979.04053
191 192 193 194 195
14582.57981 25724.78799 -10726.61674 -20476.60749 -749.09312
196 197 198 199 200
9161.48249 -788.51354 -16986.49645 -4069.46484 -23663.87858
201 202 203 204 205
-412.52230 -29274.38424 37217.72342 20667.92260 1527.10537
206 207 208 209 210
-24526.57156 -9138.75869 86196.74319 39595.53682 12269.53655
211 212 213 214 215
25681.30704 22845.79586 -25833.26267 -8593.93189 -12830.08305
216 217 218 219 220
-10885.34145 -15407.77494 -8657.72993 -15702.78861 18191.06023
221 222 223 224 225
-18984.59768 -2582.31993 -19641.65844 -17920.24526 -21191.51468
226 227 228 229 230
42681.45731 6579.13757 -35669.12186 -24682.85265 -45906.04895
231 232 233 234 235
-14663.01217 -11474.73109 -6021.34917 -21203.83039 -30584.43213
236 237 238 239 240
-44172.96623 72802.30483 -27982.26411 -13109.05887 -29303.88693
241 242 243 244 245
15857.22161 37017.82138 -34035.75672 -22836.99206 -8661.37160
246 247 248 249 250
-9075.31576 -13654.29794 20468.76156 15447.71686 -25187.82865
251 252 253 254 255
48465.16659 -16612.67549 -13534.22494 28133.05847 10158.55754
256 257 258 259 260
-1029.27115 -40171.14001 5055.71207 -1443.36333 -13889.62917
261 262 263 264 265
-19972.13487 -25261.41191 3176.57626 -36729.87890 -23054.67831
266 267 268 269 270
-8184.26550 -33241.98285 -19915.33529 30455.19294 -34239.62063
271 272 273 274 275
25715.66205 -6725.79016 17839.06982 -2202.87352 -5231.99488
276 277 278 279 280
-33721.94544 1211.73279 -12158.15154 -3796.20644 12081.95025
281 282 283 284 285
-12243.11968 -3602.58689 11264.26255 7191.44346 1935.02952
286 287 288 289
-30637.07695 348.06676 -22442.45967 2482.32017
> postscript(file="/var/wessaorg/rcomp/tmp/6su7a1324127192.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 9459.55117 NA
1 -42318.42267 9459.55117
2 -6191.83910 -42318.42267
3 -51181.44191 -6191.83910
4 -1424.77691 -51181.44191
5 -34217.21900 -1424.77691
6 142113.49876 -34217.21900
7 -23311.33777 142113.49876
8 37933.80275 -23311.33777
9 -5452.25711 37933.80275
10 5202.08209 -5452.25711
11 40310.65109 5202.08209
12 -45502.31467 40310.65109
13 -7345.85894 -45502.31467
14 23141.98373 -7345.85894
15 12204.28474 23141.98373
16 8093.56625 12204.28474
17 -5104.05047 8093.56625
18 39982.50566 -5104.05047
19 -52177.46897 39982.50566
20 -51444.03012 -52177.46897
21 -17146.35953 -51444.03012
22 158733.95290 -17146.35953
23 12478.74058 158733.95290
24 -83418.08288 12478.74058
25 -67848.80908 -83418.08288
26 -10648.99555 -67848.80908
27 30320.34322 -10648.99555
28 2142.80883 30320.34322
29 94.76229 2142.80883
30 -1782.88705 94.76229
31 -18749.30410 -1782.88705
32 2315.48642 -18749.30410
33 4313.23727 2315.48642
34 -50167.40189 4313.23727
35 59633.04958 -50167.40189
36 18293.18762 59633.04958
37 41379.84449 18293.18762
38 119897.55241 41379.84449
39 -8067.82469 119897.55241
40 14284.84145 -8067.82469
41 17132.88928 14284.84145
42 -5281.92931 17132.88928
43 52389.55026 -5281.92931
44 53307.82657 52389.55026
45 -15071.18073 53307.82657
46 44956.57391 -15071.18073
47 -9709.61554 44956.57391
48 -16478.38248 -9709.61554
49 -33812.38589 -16478.38248
50 161619.86137 -33812.38589
51 -17583.46955 161619.86137
52 -36795.71870 -17583.46955
53 9005.58453 -36795.71870
54 -47498.83271 9005.58453
55 -10372.69665 -47498.83271
56 -9061.65831 -10372.69665
57 110528.14257 -9061.65831
58 -18128.10104 110528.14257
59 15028.86516 -18128.10104
60 34781.01161 15028.86516
61 -19416.58746 34781.01161
62 137029.71115 -19416.58746
63 20444.19031 137029.71115
64 17282.06624 20444.19031
65 -1654.78295 17282.06624
66 27803.92871 -1654.78295
67 32040.02889 27803.92871
68 2840.02156 32040.02889
69 18379.84862 2840.02156
70 38226.85155 18379.84862
71 46995.24938 38226.85155
72 11432.24145 46995.24938
73 -36963.17969 11432.24145
74 -5142.91568 -36963.17969
75 -43013.30532 -5142.91568
76 -101950.60841 -43013.30532
77 -18581.69736 -101950.60841
78 -15817.45803 -18581.69736
79 5321.68653 -15817.45803
80 -8038.56928 5321.68653
81 -78148.44311 -8038.56928
82 33610.05983 -78148.44311
83 41360.29097 33610.05983
84 -62022.46549 41360.29097
85 -69616.45843 -62022.46549
86 -15525.92548 -69616.45843
87 -11069.61226 -15525.92548
88 88434.16372 -11069.61226
89 35767.92906 88434.16372
90 -30536.89723 35767.92906
91 -25791.94407 -30536.89723
92 31851.91723 -25791.94407
93 -48856.82370 31851.91723
94 -559.69883 -48856.82370
95 120669.73879 -559.69883
96 34782.44303 120669.73879
97 -135255.32333 34782.44303
98 37721.16333 -135255.32333
99 -72746.88468 37721.16333
100 23682.28990 -72746.88468
101 15267.17887 23682.28990
102 67671.59079 15267.17887
103 -25564.35686 67671.59079
104 -27624.55868 -25564.35686
105 -49609.77672 -27624.55868
106 -12128.11343 -49609.77672
107 -18108.62133 -12128.11343
108 42177.36240 -18108.62133
109 -11030.87558 42177.36240
110 -25917.17252 -11030.87558
111 14421.25911 -25917.17252
112 -36734.85479 14421.25911
113 5121.06623 -36734.85479
114 -31452.87244 5121.06623
115 -25637.27442 -31452.87244
116 -43405.93363 -25637.27442
117 98679.33854 -43405.93363
118 -61660.48481 98679.33854
119 -11554.75924 -61660.48481
120 3258.48559 -11554.75924
121 -10691.33122 3258.48559
122 26243.54396 -10691.33122
123 -116196.47094 26243.54396
124 46827.37628 -116196.47094
125 10781.62257 46827.37628
126 15191.12727 10781.62257
127 -43226.34386 15191.12727
128 -37029.30306 -43226.34386
129 -7756.67525 -37029.30306
130 5749.94481 -7756.67525
131 -65684.64897 5749.94481
132 -8359.16811 -65684.64897
133 -21968.92476 -8359.16811
134 -43115.36233 -21968.92476
135 79667.12025 -43115.36233
136 1435.49405 79667.12025
137 44157.36125 1435.49405
138 -45807.76922 44157.36125
139 -23525.94368 -45807.76922
140 -11058.81079 -23525.94368
141 -10620.19957 -11058.81079
142 30749.35717 -10620.19957
143 29772.53703 30749.35717
144 -16845.74526 29772.53703
145 -18221.37351 -16845.74526
146 -25790.02002 -18221.37351
147 -24302.53665 -25790.02002
148 -12391.20456 -24302.53665
149 -12535.82700 -12391.20456
150 30894.23177 -12535.82700
151 -4365.06184 30894.23177
152 3668.55944 -4365.06184
153 52325.71859 3668.55944
154 -20242.49620 52325.71859
155 18886.15233 -20242.49620
156 -6577.67244 18886.15233
157 23434.21451 -6577.67244
158 -37384.99266 23434.21451
159 52234.20934 -37384.99266
160 112967.21883 52234.20934
161 -42213.13655 112967.21883
162 16573.63549 -42213.13655
163 62741.33280 16573.63549
164 -10454.84274 62741.33280
165 -6419.53683 -10454.84274
166 -5975.10608 -6419.53683
167 28564.89972 -5975.10608
168 -4274.83988 28564.89972
169 124563.73369 -4274.83988
170 87867.57984 124563.73369
171 -40331.79298 87867.57984
172 60035.50787 -40331.79298
173 -6737.48329 60035.50787
174 -185309.72510 -6737.48329
175 56973.10550 -185309.72510
176 -23523.24308 56973.10550
177 9006.32656 -23523.24308
178 -59575.07568 9006.32656
179 35100.29855 -59575.07568
180 650.16261 35100.29855
181 -47574.47123 650.16261
182 -6902.73900 -47574.47123
183 35165.98104 -6902.73900
184 45421.03807 35165.98104
185 -16164.19011 45421.03807
186 6957.67253 -16164.19011
187 -425.86320 6957.67253
188 35314.93984 -425.86320
189 -31979.04053 35314.93984
190 14582.57981 -31979.04053
191 25724.78799 14582.57981
192 -10726.61674 25724.78799
193 -20476.60749 -10726.61674
194 -749.09312 -20476.60749
195 9161.48249 -749.09312
196 -788.51354 9161.48249
197 -16986.49645 -788.51354
198 -4069.46484 -16986.49645
199 -23663.87858 -4069.46484
200 -412.52230 -23663.87858
201 -29274.38424 -412.52230
202 37217.72342 -29274.38424
203 20667.92260 37217.72342
204 1527.10537 20667.92260
205 -24526.57156 1527.10537
206 -9138.75869 -24526.57156
207 86196.74319 -9138.75869
208 39595.53682 86196.74319
209 12269.53655 39595.53682
210 25681.30704 12269.53655
211 22845.79586 25681.30704
212 -25833.26267 22845.79586
213 -8593.93189 -25833.26267
214 -12830.08305 -8593.93189
215 -10885.34145 -12830.08305
216 -15407.77494 -10885.34145
217 -8657.72993 -15407.77494
218 -15702.78861 -8657.72993
219 18191.06023 -15702.78861
220 -18984.59768 18191.06023
221 -2582.31993 -18984.59768
222 -19641.65844 -2582.31993
223 -17920.24526 -19641.65844
224 -21191.51468 -17920.24526
225 42681.45731 -21191.51468
226 6579.13757 42681.45731
227 -35669.12186 6579.13757
228 -24682.85265 -35669.12186
229 -45906.04895 -24682.85265
230 -14663.01217 -45906.04895
231 -11474.73109 -14663.01217
232 -6021.34917 -11474.73109
233 -21203.83039 -6021.34917
234 -30584.43213 -21203.83039
235 -44172.96623 -30584.43213
236 72802.30483 -44172.96623
237 -27982.26411 72802.30483
238 -13109.05887 -27982.26411
239 -29303.88693 -13109.05887
240 15857.22161 -29303.88693
241 37017.82138 15857.22161
242 -34035.75672 37017.82138
243 -22836.99206 -34035.75672
244 -8661.37160 -22836.99206
245 -9075.31576 -8661.37160
246 -13654.29794 -9075.31576
247 20468.76156 -13654.29794
248 15447.71686 20468.76156
249 -25187.82865 15447.71686
250 48465.16659 -25187.82865
251 -16612.67549 48465.16659
252 -13534.22494 -16612.67549
253 28133.05847 -13534.22494
254 10158.55754 28133.05847
255 -1029.27115 10158.55754
256 -40171.14001 -1029.27115
257 5055.71207 -40171.14001
258 -1443.36333 5055.71207
259 -13889.62917 -1443.36333
260 -19972.13487 -13889.62917
261 -25261.41191 -19972.13487
262 3176.57626 -25261.41191
263 -36729.87890 3176.57626
264 -23054.67831 -36729.87890
265 -8184.26550 -23054.67831
266 -33241.98285 -8184.26550
267 -19915.33529 -33241.98285
268 30455.19294 -19915.33529
269 -34239.62063 30455.19294
270 25715.66205 -34239.62063
271 -6725.79016 25715.66205
272 17839.06982 -6725.79016
273 -2202.87352 17839.06982
274 -5231.99488 -2202.87352
275 -33721.94544 -5231.99488
276 1211.73279 -33721.94544
277 -12158.15154 1211.73279
278 -3796.20644 -12158.15154
279 12081.95025 -3796.20644
280 -12243.11968 12081.95025
281 -3602.58689 -12243.11968
282 11264.26255 -3602.58689
283 7191.44346 11264.26255
284 1935.02952 7191.44346
285 -30637.07695 1935.02952
286 348.06676 -30637.07695
287 -22442.45967 348.06676
288 2482.32017 -22442.45967
289 NA 2482.32017
> dum1 <- dum[2:length(myerror),]
> dum1
lag(myerror, k = 1) myerror
[1,] -42318.42267 9459.55117
[2,] -6191.83910 -42318.42267
[3,] -51181.44191 -6191.83910
[4,] -1424.77691 -51181.44191
[5,] -34217.21900 -1424.77691
[6,] 142113.49876 -34217.21900
[7,] -23311.33777 142113.49876
[8,] 37933.80275 -23311.33777
[9,] -5452.25711 37933.80275
[10,] 5202.08209 -5452.25711
[11,] 40310.65109 5202.08209
[12,] -45502.31467 40310.65109
[13,] -7345.85894 -45502.31467
[14,] 23141.98373 -7345.85894
[15,] 12204.28474 23141.98373
[16,] 8093.56625 12204.28474
[17,] -5104.05047 8093.56625
[18,] 39982.50566 -5104.05047
[19,] -52177.46897 39982.50566
[20,] -51444.03012 -52177.46897
[21,] -17146.35953 -51444.03012
[22,] 158733.95290 -17146.35953
[23,] 12478.74058 158733.95290
[24,] -83418.08288 12478.74058
[25,] -67848.80908 -83418.08288
[26,] -10648.99555 -67848.80908
[27,] 30320.34322 -10648.99555
[28,] 2142.80883 30320.34322
[29,] 94.76229 2142.80883
[30,] -1782.88705 94.76229
[31,] -18749.30410 -1782.88705
[32,] 2315.48642 -18749.30410
[33,] 4313.23727 2315.48642
[34,] -50167.40189 4313.23727
[35,] 59633.04958 -50167.40189
[36,] 18293.18762 59633.04958
[37,] 41379.84449 18293.18762
[38,] 119897.55241 41379.84449
[39,] -8067.82469 119897.55241
[40,] 14284.84145 -8067.82469
[41,] 17132.88928 14284.84145
[42,] -5281.92931 17132.88928
[43,] 52389.55026 -5281.92931
[44,] 53307.82657 52389.55026
[45,] -15071.18073 53307.82657
[46,] 44956.57391 -15071.18073
[47,] -9709.61554 44956.57391
[48,] -16478.38248 -9709.61554
[49,] -33812.38589 -16478.38248
[50,] 161619.86137 -33812.38589
[51,] -17583.46955 161619.86137
[52,] -36795.71870 -17583.46955
[53,] 9005.58453 -36795.71870
[54,] -47498.83271 9005.58453
[55,] -10372.69665 -47498.83271
[56,] -9061.65831 -10372.69665
[57,] 110528.14257 -9061.65831
[58,] -18128.10104 110528.14257
[59,] 15028.86516 -18128.10104
[60,] 34781.01161 15028.86516
[61,] -19416.58746 34781.01161
[62,] 137029.71115 -19416.58746
[63,] 20444.19031 137029.71115
[64,] 17282.06624 20444.19031
[65,] -1654.78295 17282.06624
[66,] 27803.92871 -1654.78295
[67,] 32040.02889 27803.92871
[68,] 2840.02156 32040.02889
[69,] 18379.84862 2840.02156
[70,] 38226.85155 18379.84862
[71,] 46995.24938 38226.85155
[72,] 11432.24145 46995.24938
[73,] -36963.17969 11432.24145
[74,] -5142.91568 -36963.17969
[75,] -43013.30532 -5142.91568
[76,] -101950.60841 -43013.30532
[77,] -18581.69736 -101950.60841
[78,] -15817.45803 -18581.69736
[79,] 5321.68653 -15817.45803
[80,] -8038.56928 5321.68653
[81,] -78148.44311 -8038.56928
[82,] 33610.05983 -78148.44311
[83,] 41360.29097 33610.05983
[84,] -62022.46549 41360.29097
[85,] -69616.45843 -62022.46549
[86,] -15525.92548 -69616.45843
[87,] -11069.61226 -15525.92548
[88,] 88434.16372 -11069.61226
[89,] 35767.92906 88434.16372
[90,] -30536.89723 35767.92906
[91,] -25791.94407 -30536.89723
[92,] 31851.91723 -25791.94407
[93,] -48856.82370 31851.91723
[94,] -559.69883 -48856.82370
[95,] 120669.73879 -559.69883
[96,] 34782.44303 120669.73879
[97,] -135255.32333 34782.44303
[98,] 37721.16333 -135255.32333
[99,] -72746.88468 37721.16333
[100,] 23682.28990 -72746.88468
[101,] 15267.17887 23682.28990
[102,] 67671.59079 15267.17887
[103,] -25564.35686 67671.59079
[104,] -27624.55868 -25564.35686
[105,] -49609.77672 -27624.55868
[106,] -12128.11343 -49609.77672
[107,] -18108.62133 -12128.11343
[108,] 42177.36240 -18108.62133
[109,] -11030.87558 42177.36240
[110,] -25917.17252 -11030.87558
[111,] 14421.25911 -25917.17252
[112,] -36734.85479 14421.25911
[113,] 5121.06623 -36734.85479
[114,] -31452.87244 5121.06623
[115,] -25637.27442 -31452.87244
[116,] -43405.93363 -25637.27442
[117,] 98679.33854 -43405.93363
[118,] -61660.48481 98679.33854
[119,] -11554.75924 -61660.48481
[120,] 3258.48559 -11554.75924
[121,] -10691.33122 3258.48559
[122,] 26243.54396 -10691.33122
[123,] -116196.47094 26243.54396
[124,] 46827.37628 -116196.47094
[125,] 10781.62257 46827.37628
[126,] 15191.12727 10781.62257
[127,] -43226.34386 15191.12727
[128,] -37029.30306 -43226.34386
[129,] -7756.67525 -37029.30306
[130,] 5749.94481 -7756.67525
[131,] -65684.64897 5749.94481
[132,] -8359.16811 -65684.64897
[133,] -21968.92476 -8359.16811
[134,] -43115.36233 -21968.92476
[135,] 79667.12025 -43115.36233
[136,] 1435.49405 79667.12025
[137,] 44157.36125 1435.49405
[138,] -45807.76922 44157.36125
[139,] -23525.94368 -45807.76922
[140,] -11058.81079 -23525.94368
[141,] -10620.19957 -11058.81079
[142,] 30749.35717 -10620.19957
[143,] 29772.53703 30749.35717
[144,] -16845.74526 29772.53703
[145,] -18221.37351 -16845.74526
[146,] -25790.02002 -18221.37351
[147,] -24302.53665 -25790.02002
[148,] -12391.20456 -24302.53665
[149,] -12535.82700 -12391.20456
[150,] 30894.23177 -12535.82700
[151,] -4365.06184 30894.23177
[152,] 3668.55944 -4365.06184
[153,] 52325.71859 3668.55944
[154,] -20242.49620 52325.71859
[155,] 18886.15233 -20242.49620
[156,] -6577.67244 18886.15233
[157,] 23434.21451 -6577.67244
[158,] -37384.99266 23434.21451
[159,] 52234.20934 -37384.99266
[160,] 112967.21883 52234.20934
[161,] -42213.13655 112967.21883
[162,] 16573.63549 -42213.13655
[163,] 62741.33280 16573.63549
[164,] -10454.84274 62741.33280
[165,] -6419.53683 -10454.84274
[166,] -5975.10608 -6419.53683
[167,] 28564.89972 -5975.10608
[168,] -4274.83988 28564.89972
[169,] 124563.73369 -4274.83988
[170,] 87867.57984 124563.73369
[171,] -40331.79298 87867.57984
[172,] 60035.50787 -40331.79298
[173,] -6737.48329 60035.50787
[174,] -185309.72510 -6737.48329
[175,] 56973.10550 -185309.72510
[176,] -23523.24308 56973.10550
[177,] 9006.32656 -23523.24308
[178,] -59575.07568 9006.32656
[179,] 35100.29855 -59575.07568
[180,] 650.16261 35100.29855
[181,] -47574.47123 650.16261
[182,] -6902.73900 -47574.47123
[183,] 35165.98104 -6902.73900
[184,] 45421.03807 35165.98104
[185,] -16164.19011 45421.03807
[186,] 6957.67253 -16164.19011
[187,] -425.86320 6957.67253
[188,] 35314.93984 -425.86320
[189,] -31979.04053 35314.93984
[190,] 14582.57981 -31979.04053
[191,] 25724.78799 14582.57981
[192,] -10726.61674 25724.78799
[193,] -20476.60749 -10726.61674
[194,] -749.09312 -20476.60749
[195,] 9161.48249 -749.09312
[196,] -788.51354 9161.48249
[197,] -16986.49645 -788.51354
[198,] -4069.46484 -16986.49645
[199,] -23663.87858 -4069.46484
[200,] -412.52230 -23663.87858
[201,] -29274.38424 -412.52230
[202,] 37217.72342 -29274.38424
[203,] 20667.92260 37217.72342
[204,] 1527.10537 20667.92260
[205,] -24526.57156 1527.10537
[206,] -9138.75869 -24526.57156
[207,] 86196.74319 -9138.75869
[208,] 39595.53682 86196.74319
[209,] 12269.53655 39595.53682
[210,] 25681.30704 12269.53655
[211,] 22845.79586 25681.30704
[212,] -25833.26267 22845.79586
[213,] -8593.93189 -25833.26267
[214,] -12830.08305 -8593.93189
[215,] -10885.34145 -12830.08305
[216,] -15407.77494 -10885.34145
[217,] -8657.72993 -15407.77494
[218,] -15702.78861 -8657.72993
[219,] 18191.06023 -15702.78861
[220,] -18984.59768 18191.06023
[221,] -2582.31993 -18984.59768
[222,] -19641.65844 -2582.31993
[223,] -17920.24526 -19641.65844
[224,] -21191.51468 -17920.24526
[225,] 42681.45731 -21191.51468
[226,] 6579.13757 42681.45731
[227,] -35669.12186 6579.13757
[228,] -24682.85265 -35669.12186
[229,] -45906.04895 -24682.85265
[230,] -14663.01217 -45906.04895
[231,] -11474.73109 -14663.01217
[232,] -6021.34917 -11474.73109
[233,] -21203.83039 -6021.34917
[234,] -30584.43213 -21203.83039
[235,] -44172.96623 -30584.43213
[236,] 72802.30483 -44172.96623
[237,] -27982.26411 72802.30483
[238,] -13109.05887 -27982.26411
[239,] -29303.88693 -13109.05887
[240,] 15857.22161 -29303.88693
[241,] 37017.82138 15857.22161
[242,] -34035.75672 37017.82138
[243,] -22836.99206 -34035.75672
[244,] -8661.37160 -22836.99206
[245,] -9075.31576 -8661.37160
[246,] -13654.29794 -9075.31576
[247,] 20468.76156 -13654.29794
[248,] 15447.71686 20468.76156
[249,] -25187.82865 15447.71686
[250,] 48465.16659 -25187.82865
[251,] -16612.67549 48465.16659
[252,] -13534.22494 -16612.67549
[253,] 28133.05847 -13534.22494
[254,] 10158.55754 28133.05847
[255,] -1029.27115 10158.55754
[256,] -40171.14001 -1029.27115
[257,] 5055.71207 -40171.14001
[258,] -1443.36333 5055.71207
[259,] -13889.62917 -1443.36333
[260,] -19972.13487 -13889.62917
[261,] -25261.41191 -19972.13487
[262,] 3176.57626 -25261.41191
[263,] -36729.87890 3176.57626
[264,] -23054.67831 -36729.87890
[265,] -8184.26550 -23054.67831
[266,] -33241.98285 -8184.26550
[267,] -19915.33529 -33241.98285
[268,] 30455.19294 -19915.33529
[269,] -34239.62063 30455.19294
[270,] 25715.66205 -34239.62063
[271,] -6725.79016 25715.66205
[272,] 17839.06982 -6725.79016
[273,] -2202.87352 17839.06982
[274,] -5231.99488 -2202.87352
[275,] -33721.94544 -5231.99488
[276,] 1211.73279 -33721.94544
[277,] -12158.15154 1211.73279
[278,] -3796.20644 -12158.15154
[279,] 12081.95025 -3796.20644
[280,] -12243.11968 12081.95025
[281,] -3602.58689 -12243.11968
[282,] 11264.26255 -3602.58689
[283,] 7191.44346 11264.26255
[284,] 1935.02952 7191.44346
[285,] -30637.07695 1935.02952
[286,] 348.06676 -30637.07695
[287,] -22442.45967 348.06676
[288,] 2482.32017 -22442.45967
> z <- as.data.frame(dum1)
> z
lag(myerror, k = 1) myerror
1 -42318.42267 9459.55117
2 -6191.83910 -42318.42267
3 -51181.44191 -6191.83910
4 -1424.77691 -51181.44191
5 -34217.21900 -1424.77691
6 142113.49876 -34217.21900
7 -23311.33777 142113.49876
8 37933.80275 -23311.33777
9 -5452.25711 37933.80275
10 5202.08209 -5452.25711
11 40310.65109 5202.08209
12 -45502.31467 40310.65109
13 -7345.85894 -45502.31467
14 23141.98373 -7345.85894
15 12204.28474 23141.98373
16 8093.56625 12204.28474
17 -5104.05047 8093.56625
18 39982.50566 -5104.05047
19 -52177.46897 39982.50566
20 -51444.03012 -52177.46897
21 -17146.35953 -51444.03012
22 158733.95290 -17146.35953
23 12478.74058 158733.95290
24 -83418.08288 12478.74058
25 -67848.80908 -83418.08288
26 -10648.99555 -67848.80908
27 30320.34322 -10648.99555
28 2142.80883 30320.34322
29 94.76229 2142.80883
30 -1782.88705 94.76229
31 -18749.30410 -1782.88705
32 2315.48642 -18749.30410
33 4313.23727 2315.48642
34 -50167.40189 4313.23727
35 59633.04958 -50167.40189
36 18293.18762 59633.04958
37 41379.84449 18293.18762
38 119897.55241 41379.84449
39 -8067.82469 119897.55241
40 14284.84145 -8067.82469
41 17132.88928 14284.84145
42 -5281.92931 17132.88928
43 52389.55026 -5281.92931
44 53307.82657 52389.55026
45 -15071.18073 53307.82657
46 44956.57391 -15071.18073
47 -9709.61554 44956.57391
48 -16478.38248 -9709.61554
49 -33812.38589 -16478.38248
50 161619.86137 -33812.38589
51 -17583.46955 161619.86137
52 -36795.71870 -17583.46955
53 9005.58453 -36795.71870
54 -47498.83271 9005.58453
55 -10372.69665 -47498.83271
56 -9061.65831 -10372.69665
57 110528.14257 -9061.65831
58 -18128.10104 110528.14257
59 15028.86516 -18128.10104
60 34781.01161 15028.86516
61 -19416.58746 34781.01161
62 137029.71115 -19416.58746
63 20444.19031 137029.71115
64 17282.06624 20444.19031
65 -1654.78295 17282.06624
66 27803.92871 -1654.78295
67 32040.02889 27803.92871
68 2840.02156 32040.02889
69 18379.84862 2840.02156
70 38226.85155 18379.84862
71 46995.24938 38226.85155
72 11432.24145 46995.24938
73 -36963.17969 11432.24145
74 -5142.91568 -36963.17969
75 -43013.30532 -5142.91568
76 -101950.60841 -43013.30532
77 -18581.69736 -101950.60841
78 -15817.45803 -18581.69736
79 5321.68653 -15817.45803
80 -8038.56928 5321.68653
81 -78148.44311 -8038.56928
82 33610.05983 -78148.44311
83 41360.29097 33610.05983
84 -62022.46549 41360.29097
85 -69616.45843 -62022.46549
86 -15525.92548 -69616.45843
87 -11069.61226 -15525.92548
88 88434.16372 -11069.61226
89 35767.92906 88434.16372
90 -30536.89723 35767.92906
91 -25791.94407 -30536.89723
92 31851.91723 -25791.94407
93 -48856.82370 31851.91723
94 -559.69883 -48856.82370
95 120669.73879 -559.69883
96 34782.44303 120669.73879
97 -135255.32333 34782.44303
98 37721.16333 -135255.32333
99 -72746.88468 37721.16333
100 23682.28990 -72746.88468
101 15267.17887 23682.28990
102 67671.59079 15267.17887
103 -25564.35686 67671.59079
104 -27624.55868 -25564.35686
105 -49609.77672 -27624.55868
106 -12128.11343 -49609.77672
107 -18108.62133 -12128.11343
108 42177.36240 -18108.62133
109 -11030.87558 42177.36240
110 -25917.17252 -11030.87558
111 14421.25911 -25917.17252
112 -36734.85479 14421.25911
113 5121.06623 -36734.85479
114 -31452.87244 5121.06623
115 -25637.27442 -31452.87244
116 -43405.93363 -25637.27442
117 98679.33854 -43405.93363
118 -61660.48481 98679.33854
119 -11554.75924 -61660.48481
120 3258.48559 -11554.75924
121 -10691.33122 3258.48559
122 26243.54396 -10691.33122
123 -116196.47094 26243.54396
124 46827.37628 -116196.47094
125 10781.62257 46827.37628
126 15191.12727 10781.62257
127 -43226.34386 15191.12727
128 -37029.30306 -43226.34386
129 -7756.67525 -37029.30306
130 5749.94481 -7756.67525
131 -65684.64897 5749.94481
132 -8359.16811 -65684.64897
133 -21968.92476 -8359.16811
134 -43115.36233 -21968.92476
135 79667.12025 -43115.36233
136 1435.49405 79667.12025
137 44157.36125 1435.49405
138 -45807.76922 44157.36125
139 -23525.94368 -45807.76922
140 -11058.81079 -23525.94368
141 -10620.19957 -11058.81079
142 30749.35717 -10620.19957
143 29772.53703 30749.35717
144 -16845.74526 29772.53703
145 -18221.37351 -16845.74526
146 -25790.02002 -18221.37351
147 -24302.53665 -25790.02002
148 -12391.20456 -24302.53665
149 -12535.82700 -12391.20456
150 30894.23177 -12535.82700
151 -4365.06184 30894.23177
152 3668.55944 -4365.06184
153 52325.71859 3668.55944
154 -20242.49620 52325.71859
155 18886.15233 -20242.49620
156 -6577.67244 18886.15233
157 23434.21451 -6577.67244
158 -37384.99266 23434.21451
159 52234.20934 -37384.99266
160 112967.21883 52234.20934
161 -42213.13655 112967.21883
162 16573.63549 -42213.13655
163 62741.33280 16573.63549
164 -10454.84274 62741.33280
165 -6419.53683 -10454.84274
166 -5975.10608 -6419.53683
167 28564.89972 -5975.10608
168 -4274.83988 28564.89972
169 124563.73369 -4274.83988
170 87867.57984 124563.73369
171 -40331.79298 87867.57984
172 60035.50787 -40331.79298
173 -6737.48329 60035.50787
174 -185309.72510 -6737.48329
175 56973.10550 -185309.72510
176 -23523.24308 56973.10550
177 9006.32656 -23523.24308
178 -59575.07568 9006.32656
179 35100.29855 -59575.07568
180 650.16261 35100.29855
181 -47574.47123 650.16261
182 -6902.73900 -47574.47123
183 35165.98104 -6902.73900
184 45421.03807 35165.98104
185 -16164.19011 45421.03807
186 6957.67253 -16164.19011
187 -425.86320 6957.67253
188 35314.93984 -425.86320
189 -31979.04053 35314.93984
190 14582.57981 -31979.04053
191 25724.78799 14582.57981
192 -10726.61674 25724.78799
193 -20476.60749 -10726.61674
194 -749.09312 -20476.60749
195 9161.48249 -749.09312
196 -788.51354 9161.48249
197 -16986.49645 -788.51354
198 -4069.46484 -16986.49645
199 -23663.87858 -4069.46484
200 -412.52230 -23663.87858
201 -29274.38424 -412.52230
202 37217.72342 -29274.38424
203 20667.92260 37217.72342
204 1527.10537 20667.92260
205 -24526.57156 1527.10537
206 -9138.75869 -24526.57156
207 86196.74319 -9138.75869
208 39595.53682 86196.74319
209 12269.53655 39595.53682
210 25681.30704 12269.53655
211 22845.79586 25681.30704
212 -25833.26267 22845.79586
213 -8593.93189 -25833.26267
214 -12830.08305 -8593.93189
215 -10885.34145 -12830.08305
216 -15407.77494 -10885.34145
217 -8657.72993 -15407.77494
218 -15702.78861 -8657.72993
219 18191.06023 -15702.78861
220 -18984.59768 18191.06023
221 -2582.31993 -18984.59768
222 -19641.65844 -2582.31993
223 -17920.24526 -19641.65844
224 -21191.51468 -17920.24526
225 42681.45731 -21191.51468
226 6579.13757 42681.45731
227 -35669.12186 6579.13757
228 -24682.85265 -35669.12186
229 -45906.04895 -24682.85265
230 -14663.01217 -45906.04895
231 -11474.73109 -14663.01217
232 -6021.34917 -11474.73109
233 -21203.83039 -6021.34917
234 -30584.43213 -21203.83039
235 -44172.96623 -30584.43213
236 72802.30483 -44172.96623
237 -27982.26411 72802.30483
238 -13109.05887 -27982.26411
239 -29303.88693 -13109.05887
240 15857.22161 -29303.88693
241 37017.82138 15857.22161
242 -34035.75672 37017.82138
243 -22836.99206 -34035.75672
244 -8661.37160 -22836.99206
245 -9075.31576 -8661.37160
246 -13654.29794 -9075.31576
247 20468.76156 -13654.29794
248 15447.71686 20468.76156
249 -25187.82865 15447.71686
250 48465.16659 -25187.82865
251 -16612.67549 48465.16659
252 -13534.22494 -16612.67549
253 28133.05847 -13534.22494
254 10158.55754 28133.05847
255 -1029.27115 10158.55754
256 -40171.14001 -1029.27115
257 5055.71207 -40171.14001
258 -1443.36333 5055.71207
259 -13889.62917 -1443.36333
260 -19972.13487 -13889.62917
261 -25261.41191 -19972.13487
262 3176.57626 -25261.41191
263 -36729.87890 3176.57626
264 -23054.67831 -36729.87890
265 -8184.26550 -23054.67831
266 -33241.98285 -8184.26550
267 -19915.33529 -33241.98285
268 30455.19294 -19915.33529
269 -34239.62063 30455.19294
270 25715.66205 -34239.62063
271 -6725.79016 25715.66205
272 17839.06982 -6725.79016
273 -2202.87352 17839.06982
274 -5231.99488 -2202.87352
275 -33721.94544 -5231.99488
276 1211.73279 -33721.94544
277 -12158.15154 1211.73279
278 -3796.20644 -12158.15154
279 12081.95025 -3796.20644
280 -12243.11968 12081.95025
281 -3602.58689 -12243.11968
282 11264.26255 -3602.58689
283 7191.44346 11264.26255
284 1935.02952 7191.44346
285 -30637.07695 1935.02952
286 348.06676 -30637.07695
287 -22442.45967 348.06676
288 2482.32017 -22442.45967
> 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/7opul1324127192.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/81c781324127192.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/9ycpj1324127192.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/10pkhv1324127192.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/11phjv1324127192.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/12rbfa1324127192.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/13dwi51324127192.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/14rpfg1324127192.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/15ea2g1324127192.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/16m9pl1324127192.tab")
+ }
>
> try(system("convert tmp/17m4g1324127192.ps tmp/17m4g1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kdvx1324127192.ps tmp/2kdvx1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/3650n1324127192.ps tmp/3650n1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/469sk1324127192.ps tmp/469sk1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/5i1h81324127192.ps tmp/5i1h81324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/6su7a1324127192.ps tmp/6su7a1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/7opul1324127192.ps tmp/7opul1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/81c781324127192.ps tmp/81c781324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/9ycpj1324127192.ps tmp/9ycpj1324127192.png",intern=TRUE))
character(0)
> try(system("convert tmp/10pkhv1324127192.ps tmp/10pkhv1324127192.png",intern=TRUE))
character(0)
>
>
> proc.time()
user system elapsed
7.526 0.637 8.255