R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(0 + ,807 + ,6282154 + ,0 + ,444 + ,4321023 + ,0 + ,412 + ,4111912 + ,1 + ,428 + ,223193 + ,0 + ,315 + ,1491348 + ,0 + ,168 + ,1629616 + ,1 + ,263 + ,1398893 + ,0 + ,267 + ,1926517 + ,0 + ,228 + ,983660 + ,0 + ,129 + ,1443586 + ,1 + ,104 + ,1073089 + ,1 + ,122 + ,984885 + ,0 + ,393 + ,1405225 + ,1 + ,190 + ,227132 + ,0 + ,280 + ,929118 + ,1 + ,63 + ,1071292 + ,1 + ,102 + ,638830 + ,0 + ,265 + ,856956 + ,0 + ,234 + ,992426 + ,1 + ,277 + ,444477 + ,0 + ,73 + ,857217 + ,0 + ,67 + ,711969 + ,1 + ,103 + ,702380 + ,1 + ,290 + ,358589 + ,1 + ,83 + ,297978 + ,1 + ,56 + ,585715 + ,0 + ,236 + ,657954 + ,1 + ,73 + ,209458 + ,1 + ,34 + ,786690 + ,1 + ,139 + ,439798 + ,0 + ,26 + ,688779 + ,0 + ,70 + ,574339 + ,0 + ,40 + ,741409 + ,0 + ,42 + ,597793 + ,1 + ,12 + ,644190 + ,1 + ,211 + ,377934 + ,1 + ,74 + ,640273 + ,0 + ,80 + ,697458 + ,0 + ,83 + ,550608 + ,1 + ,131 + ,207393 + ,1 + ,203 + ,301607 + ,1 + ,56 + ,345783 + ,1 + ,89 + ,501749 + ,1 + ,88 + ,379983 + ,1 + ,39 + ,387475 + ,0 + ,25 + ,377305 + ,0 + ,49 + ,370837 + ,0 + ,149 + ,430866 + ,0 + ,58 + ,469107 + ,1 + ,41 + ,194493 + ,0 + ,90 + ,530670 + ,0 + ,136 + ,518365 + ,0 + ,97 + ,491303 + ,0 + ,63 + ,527021 + ,0 + ,114 + ,233773 + ,1 + ,77 + ,405972 + ,1 + ,6 + ,652925 + ,1 + ,47 + ,446211 + ,1 + ,51 + ,341340 + ,0 + ,85 + ,387699 + ,0 + ,43 + ,493408 + ,1 + ,32 + ,146494 + ,0 + ,25 + ,414462 + ,0 + ,77 + ,364304 + ,1 + ,54 + ,355178 + ,1 + ,251 + ,357760 + ,1 + ,15 + ,261216 + ,0 + ,44 + ,397144 + ,1 + ,73 + ,374943 + ,0 + ,85 + ,424898 + ,0 + ,49 + ,202055 + ,1 + ,38 + ,378525 + ,1 + ,35 + ,310768 + ,1 + ,9 + ,325738 + ,1 + ,34 + ,394510 + ,0 + ,20 + ,247060 + ,1 + ,29 + ,368078 + ,1 + ,11 + ,236761 + ,1 + ,52 + ,312378 + ,0 + ,13 + ,339836 + ,1 + ,29 + ,347385 + ,0 + ,66 + ,426280 + ,1 + ,33 + ,352850 + ,1 + ,15 + ,301881 + ,1 + ,15 + ,377516 + ,0 + ,68 + ,357312 + ,1 + ,100 + ,458343 + ,1 + ,13 + ,354228 + ,1 + ,45 + ,308636 + ,1 + ,14 + ,386212 + ,1 + ,36 + ,393343 + ,0 + ,40 + ,378509 + ,1 + ,68 + ,452469 + ,0 + ,29 + ,364839 + ,1 + ,43 + ,358649 + ,0 + ,30 + ,376641 + ,1 + ,9 + ,429112 + ,0 + ,22 + ,330546 + ,1 + ,19 + ,403560 + ,0 + ,9 + ,317892 + ,0 + ,31 + ,307528 + ,1 + ,19 + ,235133 + ,1 + ,55 + ,299243 + ,1 + ,8 + ,314073 + ,1 + ,28 + ,368186 + ,1 + ,29 + ,269661 + ,0 + ,48 + ,125390 + ,0 + ,16 + ,510834 + ,1 + ,47 + ,321896 + ,0 + ,20 + ,249898 + ,1 + ,22 + ,408881 + ,0 + ,33 + ,158492 + ,1 + ,44 + ,292154 + ,0 + ,13 + ,289513 + ,0 + ,6 + ,378049 + ,1 + ,35 + ,343466 + ,1 + ,8 + ,332743 + ,1 + ,17 + ,442882 + ,0 + ,11 + ,214215 + ,1 + ,21 + ,315688 + ,1 + ,92 + ,375195 + ,1 + ,12 + ,334280 + ,1 + ,112 + ,355864 + ,0 + ,25 + ,480382 + ,0 + ,17 + ,353058 + ,0 + ,23 + ,217193 + ,1 + ,0 + ,315380 + ,1 + ,10 + ,314533 + ,1 + ,23 + ,318056 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,314353 + ,1 + ,25 + ,369448 + ,1 + ,0 + ,315380 + ,1 + ,20 + ,312846 + ,1 + ,4 + ,312075 + ,1 + ,4 + ,315009 + ,1 + ,10 + ,318903 + ,1 + ,1 + ,314887 + ,1 + ,4 + ,314913 + ,1 + ,0 + ,315380 + ,1 + ,8 + ,325506 + ,1 + ,0 + ,315380 + ,1 + ,11 + ,298568 + ,1 + ,4 + ,315834 + ,1 + ,15 + ,329784 + ,1 + ,9 + ,312878 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,314987 + ,1 + ,2 + ,325249 + ,1 + ,0 + ,315877 + ,1 + ,7 + ,291650 + ,1 + ,46 + ,305959 + ,1 + ,5 + ,315380 + ,1 + ,7 + ,297765 + ,1 + ,2 + ,315245 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,2 + ,315236 + ,1 + ,5 + ,336425 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,306268 + ,1 + ,24 + ,302187 + ,1 + ,1 + ,314882 + ,1 + ,0 + ,315380 + ,1 + ,18 + ,382712 + ,1 + ,55 + ,341570 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,3 + ,312412 + ,1 + ,0 + ,315380 + ,1 + ,9 + ,309596 + ,1 + ,0 + ,315380 + ,1 + ,8 + ,315547 + ,1 + ,113 + ,313267 + ,0 + ,0 + ,316176 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,19 + ,359335 + ,0 + ,11 + ,330068 + ,1 + ,25 + ,314289 + ,0 + ,16 + ,297413 + ,0 + ,5 + ,314806 + ,0 + ,11 + ,333210 + ,0 + ,23 + ,352108 + ,1 + ,6 + ,313332 + ,1 + ,5 + ,291787 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,318745 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,315366 + ,1 + ,0 + ,315380 + ,1 + ,3 + ,315688 + ,0 + ,0 + ,315380 + ,0 + ,89 + ,409642 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,19 + ,269587 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,12 + ,300962 + ,0 + ,12 + ,325479 + ,0 + ,5 + ,316155 + ,0 + ,2 + ,318574 + ,0 + ,0 + ,315380 + ,0 + ,26 + ,343613 + ,0 + ,3 + ,306948 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,11 + ,330059 + ,1 + ,10 + ,288985 + ,1 + ,5 + ,304485 + ,1 + ,2 + ,315380 + ,1 + ,6 + ,315688 + ,1 + ,7 + ,317736 + ,1 + ,2 + ,315380 + ,1 + ,28 + ,322331 + ,1 + ,3 + ,296656 + ,1 + ,0 + ,315380 + ,1 + ,1 + ,315354 + ,1 + ,20 + ,312161 + ,1 + ,1 + ,315576 + ,1 + ,22 + ,314922 + ,1 + ,9 + ,314551 + ,1 + ,0 + ,315380 + ,1 + ,2 + ,312339 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,298700 + ,1 + ,9 + ,321376 + ,1 + ,0 + ,315380 + ,1 + ,13 + ,303230 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315487 + ,1 + ,0 + ,315380 + ,1 + ,6 + ,315793 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,3 + ,312887 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,315637 + ,1 + ,2 + ,324385 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,15 + ,308989 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,9 + ,296702 + ,0 + ,0 + ,315380 + ,1 + ,1 + ,307322 + ,1 + ,38 + ,304376 + ,1 + ,57 + ,253588 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,309560 + ,1 + ,26 + ,298466 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,13 + ,343929 + ,1 + ,10 + ,331955 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,9 + ,381180 + ,0 + ,0 + ,315380 + ,1 + ,26 + ,331420 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,19 + ,310201 + ,1 + ,0 + ,315380 + ,1 + ,12 + ,320016 + ,1 + ,23 + ,320398 + ,1 + ,0 + ,315380 + ,0 + ,29 + ,291841 + ,1 + ,8 + ,310670 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,26 + ,313491 + ,1 + ,0 + ,315380 + ,1 + ,9 + ,331323 + ,1 + ,0 + ,315380 + ,0 + ,5 + ,319210 + ,1 + ,3 + ,318098 + ,0 + ,0 + ,315380 + ,1 + ,13 + ,292754 + ,0 + ,0 + ,315380 + ,1 + ,12 + ,325176 + ,1 + ,19 + ,365959 + ,1 + ,0 + ,315380 + ,1 + ,10 + ,302409 + ,0 + ,9 + ,340968 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,9 + ,315380 + ,0 + ,4 + ,313164 + ,1 + ,1 + ,301164 + ,1 + ,1 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,14 + ,344425 + ,1 + ,12 + ,315394 + ,0 + ,0 + ,315380 + ,0 + ,19 + ,316647 + ,1 + ,17 + ,309836 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,32 + ,346611 + ,0 + ,0 + ,315380 + ,0 + ,14 + ,322031 + ,1 + ,8 + ,315656 + ,1 + ,4 + ,339445 + ,1 + ,0 + ,314964 + ,1 + ,20 + ,297141 + ,1 + ,5 + ,315372 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,1 + ,312502 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,4 + ,313729 + ,1 + ,1 + ,315388 + ,1 + ,4 + ,315371 + ,1 + ,20 + ,296139 + ,1 + ,0 + ,315380 + ,1 + ,1 + ,313880 + ,1 + ,10 + ,317698 + ,1 + ,12 + ,295580 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,13 + ,308256 + ,1 + ,0 + ,315380 + ,1 + ,3 + ,303677 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,10 + ,319369 + ,1 + ,3 + ,318690 + ,1 + ,7 + ,314049 + ,1 + ,10 + ,325699 + ,1 + ,1 + ,314210 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,15 + ,322378 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,4 + ,315398 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,0 + ,28 + ,308336 + ,1 + ,9 + ,316386 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,315553 + ,1 + ,0 + ,315380 + ,1 + ,7 + ,323361 + ,1 + ,7 + ,336639 + ,1 + ,3 + ,307424 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,11 + ,295370 + ,1 + ,7 + ,322340 + ,1 + ,10 + ,319864 + ,1 + ,0 + ,315380 + ,1 + ,0 + ,315380 + ,1 + ,18 + ,317291 + ,1 + ,14 + ,280398 + ,1 + ,0 + ,315380 + ,1 + ,12 + ,317330 + ,1 + ,29 + ,238125 + ,1 + ,3 + ,327071 + ,1 + ,6 + ,309038 + ,1 + ,3 + ,314210 + ,1 + ,8 + ,307930 + ,1 + ,10 + ,322327 + ,1 + ,6 + ,292136 + ,1 + ,8 + ,263276 + ,1 + ,6 + ,367655 + ,0 + ,9 + ,283910 + ,1 + ,8 + ,283587 + ,0 + ,26 + ,243650 + ,0 + ,239 + ,438493 + ,1 + ,7 + ,296261 + ,0 + ,41 + ,230621 + ,1 + ,3 + ,304252 + ,1 + ,8 + ,333505 + ,1 + ,6 + ,296919 + ,0 + ,21 + ,278990 + ,1 + ,7 + ,276898 + ,1 + ,11 + ,327007 + ,1 + ,11 + ,317046 + ,1 + ,12 + ,304555 + ,1 + ,9 + ,298096 + ,1 + ,3 + ,231861 + ,1 + ,57 + ,309422 + ,0 + ,21 + ,286963 + ,0 + ,15 + ,269753 + ,0 + ,32 + ,448243 + ,1 + ,11 + ,165404 + ,1 + ,2 + ,204325 + ,1 + ,23 + ,407159 + ,0 + ,20 + ,290476 + ,1 + ,24 + ,275311 + ,1 + ,1 + ,246541 + ,1 + ,1 + ,253468 + ,1 + ,74 + ,240897 + ,0 + ,68 + ,-83265 + ,1 + ,20 + ,-42143 + ,1 + ,20 + ,272713 + ,0 + ,82 + ,215362 + ,1 + ,21 + ,42754 + ,1 + ,244 + ,306275 + ,1 + ,32 + ,253537 + ,1 + ,86 + ,372631 + ,1 + ,69 + ,-7170) + ,dim=c(3 + ,431) + ,dimnames=list(c('group' + ,'orders' + ,'wealth ') + ,1:431)) > y <- array(NA,dim=c(3,431),dimnames=list(c('group','orders','wealth '),1:431)) > 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 = '3' > #'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 Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > 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 wealth\r group orders 1 6282154 0 807 2 4321023 0 444 3 4111912 0 412 4 223193 1 428 5 1491348 0 315 6 1629616 0 168 7 1398893 1 263 8 1926517 0 267 9 983660 0 228 10 1443586 0 129 11 1073089 1 104 12 984885 1 122 13 1405225 0 393 14 227132 1 190 15 929118 0 280 16 1071292 1 63 17 638830 1 102 18 856956 0 265 19 992426 0 234 20 444477 1 277 21 857217 0 73 22 711969 0 67 23 702380 1 103 24 358589 1 290 25 297978 1 83 26 585715 1 56 27 657954 0 236 28 209458 1 73 29 786690 1 34 30 439798 1 139 31 688779 0 26 32 574339 0 70 33 741409 0 40 34 597793 0 42 35 644190 1 12 36 377934 1 211 37 640273 1 74 38 697458 0 80 39 550608 0 83 40 207393 1 131 41 301607 1 203 42 345783 1 56 43 501749 1 89 44 379983 1 88 45 387475 1 39 46 377305 0 25 47 370837 0 49 48 430866 0 149 49 469107 0 58 50 194493 1 41 51 530670 0 90 52 518365 0 136 53 491303 0 97 54 527021 0 63 55 233773 0 114 56 405972 1 77 57 652925 1 6 58 446211 1 47 59 341340 1 51 60 387699 0 85 61 493408 0 43 62 146494 1 32 63 414462 0 25 64 364304 0 77 65 355178 1 54 66 357760 1 251 67 261216 1 15 68 397144 0 44 69 374943 1 73 70 424898 0 85 71 202055 0 49 72 378525 1 38 73 310768 1 35 74 325738 1 9 75 394510 1 34 76 247060 0 20 77 368078 1 29 78 236761 1 11 79 312378 1 52 80 339836 0 13 81 347385 1 29 82 426280 0 66 83 352850 1 33 84 301881 1 15 85 377516 1 15 86 357312 0 68 87 458343 1 100 88 354228 1 13 89 308636 1 45 90 386212 1 14 91 393343 1 36 92 378509 0 40 93 452469 1 68 94 364839 0 29 95 358649 1 43 96 376641 0 30 97 429112 1 9 98 330546 0 22 99 403560 1 19 100 317892 0 9 101 307528 0 31 102 235133 1 19 103 299243 1 55 104 314073 1 8 105 368186 1 28 106 269661 1 29 107 125390 0 48 108 510834 0 16 109 321896 1 47 110 249898 0 20 111 408881 1 22 112 158492 0 33 113 292154 1 44 114 289513 0 13 115 378049 0 6 116 343466 1 35 117 332743 1 8 118 442882 1 17 119 214215 0 11 120 315688 1 21 121 375195 1 92 122 334280 1 12 123 355864 1 112 124 480382 0 25 125 353058 0 17 126 217193 0 23 127 315380 1 0 128 314533 1 10 129 318056 1 23 130 315380 1 0 131 314353 1 7 132 369448 1 25 133 315380 1 0 134 312846 1 20 135 312075 1 4 136 315009 1 4 137 318903 1 10 138 314887 1 1 139 314913 1 4 140 315380 1 0 141 325506 1 8 142 315380 1 0 143 298568 1 11 144 315834 1 4 145 329784 1 15 146 312878 1 9 147 315380 1 0 148 314987 1 7 149 325249 1 2 150 315877 1 0 151 291650 1 7 152 305959 1 46 153 315380 1 5 154 297765 1 7 155 315245 1 2 156 315380 1 0 157 315380 1 0 158 315236 1 2 159 336425 1 5 160 315380 1 0 161 315380 1 0 162 315380 1 0 163 315380 1 0 164 306268 1 7 165 302187 1 24 166 314882 1 1 167 315380 1 0 168 382712 1 18 169 341570 1 55 170 315380 1 0 171 315380 1 0 172 312412 1 3 173 315380 1 0 174 309596 1 9 175 315380 1 0 176 315547 1 8 177 313267 1 113 178 316176 0 0 179 315380 1 0 180 315380 1 0 181 359335 1 19 182 330068 0 11 183 314289 1 25 184 297413 0 16 185 314806 0 5 186 333210 0 11 187 352108 0 23 188 313332 1 6 189 291787 1 5 190 315380 1 0 191 318745 1 7 192 315380 1 0 193 315366 1 7 194 315380 1 0 195 315688 1 3 196 315380 0 0 197 409642 0 89 198 315380 0 0 199 315380 0 0 200 269587 0 19 201 315380 0 0 202 315380 0 0 203 315380 0 0 204 300962 0 12 205 325479 0 12 206 316155 0 5 207 318574 0 2 208 315380 0 0 209 343613 0 26 210 306948 0 3 211 315380 1 0 212 315380 1 0 213 330059 1 11 214 288985 1 10 215 304485 1 5 216 315380 1 2 217 315688 1 6 218 317736 1 7 219 315380 1 2 220 322331 1 28 221 296656 1 3 222 315380 1 0 223 315354 1 1 224 312161 1 20 225 315576 1 1 226 314922 1 22 227 314551 1 9 228 315380 1 0 229 312339 1 2 230 315380 1 0 231 298700 1 7 232 321376 1 9 233 315380 1 0 234 303230 1 13 235 315380 1 0 236 315487 1 0 237 315380 1 0 238 315793 1 6 239 315380 1 0 240 315380 1 0 241 315380 1 0 242 312887 1 3 243 315380 1 0 244 315637 1 7 245 324385 1 2 246 315380 0 0 247 315380 0 0 248 308989 1 15 249 315380 0 0 250 315380 0 0 251 296702 1 9 252 315380 0 0 253 307322 1 1 254 304376 1 38 255 253588 1 57 256 315380 1 0 257 309560 1 7 258 298466 1 26 259 315380 0 0 260 315380 1 0 261 315380 0 0 262 315380 0 0 263 343929 1 13 264 331955 1 10 265 315380 0 0 266 315380 1 0 267 315380 0 0 268 381180 1 9 269 315380 0 0 270 331420 1 26 271 315380 0 0 272 315380 0 0 273 315380 0 0 274 310201 1 19 275 315380 1 0 276 320016 1 12 277 320398 1 23 278 315380 1 0 279 291841 0 29 280 310670 1 8 281 315380 0 0 282 315380 1 0 283 313491 1 26 284 315380 1 0 285 331323 1 9 286 315380 1 0 287 319210 0 5 288 318098 1 3 289 315380 0 0 290 292754 1 13 291 315380 0 0 292 325176 1 12 293 365959 1 19 294 315380 1 0 295 302409 1 10 296 340968 0 9 297 315380 1 0 298 315380 0 0 299 315380 1 0 300 315380 1 9 301 313164 0 4 302 301164 1 1 303 315380 1 1 304 315380 0 0 305 344425 1 14 306 315394 1 12 307 315380 0 0 308 316647 0 19 309 309836 1 17 310 315380 0 0 311 315380 0 0 312 346611 1 32 313 315380 0 0 314 322031 0 14 315 315656 1 8 316 339445 1 4 317 314964 1 0 318 297141 1 20 319 315372 1 5 320 315380 1 0 321 315380 1 0 322 315380 1 0 323 315380 1 0 324 315380 0 0 325 315380 1 0 326 312502 1 1 327 315380 0 0 328 315380 0 0 329 315380 1 0 330 315380 1 0 331 315380 1 0 332 315380 1 0 333 315380 1 0 334 313729 1 4 335 315388 1 1 336 315371 1 4 337 296139 1 20 338 315380 1 0 339 313880 1 1 340 317698 1 10 341 295580 1 12 342 315380 1 0 343 315380 1 0 344 315380 1 0 345 308256 1 13 346 315380 1 0 347 303677 1 3 348 315380 1 0 349 315380 1 0 350 319369 1 10 351 318690 1 3 352 314049 1 7 353 325699 1 10 354 314210 1 1 355 315380 1 0 356 315380 1 0 357 322378 1 15 358 315380 1 0 359 315380 1 0 360 315380 1 0 361 315398 1 4 362 315380 1 0 363 315380 1 0 364 308336 0 28 365 316386 1 9 366 315380 1 0 367 315380 1 0 368 315380 1 0 369 315380 1 0 370 315553 1 7 371 315380 1 0 372 323361 1 7 373 336639 1 7 374 307424 1 3 375 315380 1 0 376 315380 1 0 377 295370 1 11 378 322340 1 7 379 319864 1 10 380 315380 1 0 381 315380 1 0 382 317291 1 18 383 280398 1 14 384 315380 1 0 385 317330 1 12 386 238125 1 29 387 327071 1 3 388 309038 1 6 389 314210 1 3 390 307930 1 8 391 322327 1 10 392 292136 1 6 393 263276 1 8 394 367655 1 6 395 283910 0 9 396 283587 1 8 397 243650 0 26 398 438493 0 239 399 296261 1 7 400 230621 0 41 401 304252 1 3 402 333505 1 8 403 296919 1 6 404 278990 0 21 405 276898 1 7 406 327007 1 11 407 317046 1 11 408 304555 1 12 409 298096 1 9 410 231861 1 3 411 309422 1 57 412 286963 0 21 413 269753 0 15 414 448243 0 32 415 165404 1 11 416 204325 1 2 417 407159 1 23 418 290476 0 20 419 275311 1 24 420 246541 1 1 421 253468 1 1 422 240897 1 74 423 -83265 0 68 424 -42143 1 20 425 272713 1 20 426 215362 0 82 427 42754 1 21 428 306275 1 244 429 253537 1 32 430 372631 1 86 431 -7170 1 69 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) group orders 282194 -50104 4263 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -1833510 -30828 34083 83290 2559626 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 282193.7 28180.1 10.014 <2e-16 *** group -50104.1 31230.8 -1.604 0.109 orders 4263.1 189.1 22.544 <2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 285000 on 428 degrees of freedom Multiple R-squared: 0.5637, Adjusted R-squared: 0.5616 F-statistic: 276.5 on 2 and 428 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,] 1.0000000 6.726328e-23 3.363164e-23 [2,] 1.0000000 3.082179e-46 1.541090e-46 [3,] 1.0000000 1.377808e-54 6.889039e-55 [4,] 1.0000000 1.009301e-60 5.046506e-61 [5,] 1.0000000 3.857151e-68 1.928575e-68 [6,] 1.0000000 5.449678e-90 2.724839e-90 [7,] 1.0000000 1.344880e-97 6.724398e-98 [8,] 1.0000000 3.642576e-117 1.821288e-117 [9,] 1.0000000 4.090349e-118 2.045175e-118 [10,] 1.0000000 4.188647e-126 2.094324e-126 [11,] 1.0000000 2.376011e-141 1.188005e-141 [12,] 1.0000000 4.164841e-142 2.082421e-142 [13,] 1.0000000 2.348058e-148 1.174029e-148 [14,] 1.0000000 6.209040e-155 3.104520e-155 [15,] 1.0000000 2.594458e-157 1.297229e-157 [16,] 1.0000000 2.757674e-162 1.378837e-162 [17,] 1.0000000 2.027162e-164 1.013581e-164 [18,] 1.0000000 1.741313e-167 8.706567e-168 [19,] 1.0000000 1.289633e-170 6.448166e-171 [20,] 1.0000000 7.914985e-170 3.957492e-170 [21,] 1.0000000 4.878305e-172 2.439153e-172 [22,] 1.0000000 4.309772e-176 2.154886e-176 [23,] 1.0000000 1.187706e-175 5.938531e-176 [24,] 1.0000000 2.671289e-184 1.335644e-184 [25,] 1.0000000 2.097498e-183 1.048749e-183 [26,] 1.0000000 9.444548e-187 4.722274e-187 [27,] 1.0000000 1.803520e-187 9.017600e-188 [28,] 1.0000000 1.049531e-192 5.247657e-193 [29,] 1.0000000 2.317962e-194 1.158981e-194 [30,] 1.0000000 4.402765e-200 2.201383e-200 [31,] 1.0000000 3.347369e-200 1.673685e-200 [32,] 1.0000000 6.966684e-204 3.483342e-204 [33,] 1.0000000 8.854423e-209 4.427212e-209 [34,] 1.0000000 4.370387e-210 2.185194e-210 [35,] 1.0000000 1.656524e-210 8.282618e-211 [36,] 1.0000000 5.495587e-211 2.747794e-211 [37,] 1.0000000 4.317012e-210 2.158506e-210 [38,] 1.0000000 1.729027e-210 8.645136e-211 [39,] 1.0000000 1.693906e-209 8.469529e-210 [40,] 1.0000000 5.109832e-209 2.554916e-209 [41,] 1.0000000 3.853302e-208 1.926651e-208 [42,] 1.0000000 2.332572e-207 1.166286e-207 [43,] 1.0000000 4.025577e-208 2.012789e-208 [44,] 1.0000000 4.224760e-208 2.112380e-208 [45,] 1.0000000 5.132186e-208 2.566093e-208 [46,] 1.0000000 3.952471e-209 1.976236e-209 [47,] 1.0000000 1.776731e-210 8.883655e-211 [48,] 1.0000000 4.066384e-211 2.033192e-211 [49,] 1.0000000 1.630917e-212 8.154586e-213 [50,] 1.0000000 2.673706e-213 1.336853e-213 [51,] 1.0000000 1.115930e-212 5.579648e-213 [52,] 1.0000000 6.899460e-221 3.449730e-221 [53,] 1.0000000 3.472489e-221 1.736245e-221 [54,] 1.0000000 3.208432e-220 1.604216e-220 [55,] 1.0000000 1.140296e-219 5.701479e-220 [56,] 1.0000000 1.046703e-220 5.233516e-221 [57,] 1.0000000 1.032100e-221 5.160502e-222 [58,] 1.0000000 2.361317e-221 1.180658e-221 [59,] 1.0000000 1.229232e-220 6.146158e-221 [60,] 1.0000000 9.967937e-220 4.983968e-220 [61,] 1.0000000 8.192262e-221 4.096131e-221 [62,] 1.0000000 3.526720e-220 1.763360e-220 [63,] 1.0000000 1.321497e-219 6.607486e-220 [64,] 1.0000000 9.152588e-219 4.576294e-219 [65,] 1.0000000 9.860747e-219 4.930374e-219 [66,] 1.0000000 8.071744e-219 4.035872e-219 [67,] 1.0000000 2.956260e-218 1.478130e-218 [68,] 1.0000000 2.541529e-217 1.270765e-217 [69,] 1.0000000 1.262817e-216 6.314086e-217 [70,] 1.0000000 2.962382e-216 1.481191e-216 [71,] 1.0000000 1.562874e-215 7.814370e-216 [72,] 1.0000000 6.717016e-215 3.358508e-215 [73,] 1.0000000 2.002139e-214 1.001070e-214 [74,] 1.0000000 2.151002e-213 1.075501e-213 [75,] 1.0000000 2.108402e-212 1.054201e-212 [76,] 1.0000000 1.370947e-211 6.854733e-212 [77,] 1.0000000 1.543117e-211 7.715586e-212 [78,] 1.0000000 9.842808e-211 4.921404e-211 [79,] 1.0000000 7.137662e-210 3.568831e-210 [80,] 1.0000000 2.037030e-209 1.018515e-209 [81,] 1.0000000 1.142408e-208 5.712042e-209 [82,] 1.0000000 3.221133e-209 1.610567e-209 [83,] 1.0000000 1.486053e-208 7.430263e-209 [84,] 1.0000000 1.536360e-207 7.681802e-208 [85,] 1.0000000 3.434063e-207 1.717031e-207 [86,] 1.0000000 8.682390e-207 4.341195e-207 [87,] 1.0000000 3.838682e-206 1.919341e-206 [88,] 1.0000000 9.262027e-207 4.631014e-207 [89,] 1.0000000 5.494880e-206 2.747440e-206 [90,] 1.0000000 3.339792e-205 1.669896e-205 [91,] 1.0000000 1.459655e-204 7.298277e-205 [92,] 1.0000000 5.429811e-205 2.714906e-205 [93,] 1.0000000 5.179814e-204 2.589907e-204 [94,] 1.0000000 6.676329e-204 3.338164e-204 [95,] 1.0000000 6.885725e-203 3.442863e-203 [96,] 1.0000000 6.857358e-202 3.428679e-202 [97,] 1.0000000 2.398643e-201 1.199322e-201 [98,] 1.0000000 2.466999e-200 1.233499e-200 [99,] 1.0000000 1.808342e-199 9.041711e-200 [100,] 1.0000000 7.830099e-199 3.915050e-199 [101,] 1.0000000 5.810222e-198 2.905111e-198 [102,] 1.0000000 6.730912e-200 3.365456e-200 [103,] 1.0000000 1.887291e-202 9.436453e-203 [104,] 1.0000000 1.885190e-201 9.425951e-202 [105,] 1.0000000 9.856601e-201 4.928300e-201 [106,] 1.0000000 9.718754e-201 4.859377e-201 [107,] 1.0000000 7.827697e-202 3.913848e-202 [108,] 1.0000000 7.796084e-201 3.898042e-201 [109,] 1.0000000 7.491519e-200 3.745760e-200 [110,] 1.0000000 2.804295e-199 1.402147e-199 [111,] 1.0000000 2.081273e-198 1.040636e-198 [112,] 1.0000000 1.397426e-197 6.987130e-198 [113,] 1.0000000 2.093055e-198 1.046527e-198 [114,] 1.0000000 3.015376e-198 1.507688e-198 [115,] 1.0000000 2.725664e-197 1.362832e-197 [116,] 1.0000000 8.266059e-197 4.133030e-197 [117,] 1.0000000 5.780207e-196 2.890103e-196 [118,] 1.0000000 2.018862e-195 1.009431e-195 [119,] 1.0000000 2.959550e-197 1.479775e-197 [120,] 1.0000000 1.923554e-196 9.617768e-197 [121,] 1.0000000 3.018466e-196 1.509233e-196 [122,] 1.0000000 2.212589e-195 1.106294e-195 [123,] 1.0000000 1.835511e-194 9.177556e-195 [124,] 1.0000000 1.659612e-193 8.298060e-194 [125,] 1.0000000 1.235918e-192 6.179591e-193 [126,] 1.0000000 1.004367e-191 5.021837e-192 [127,] 1.0000000 3.623713e-191 1.811856e-191 [128,] 1.0000000 2.744110e-190 1.372055e-190 [129,] 1.0000000 2.504589e-189 1.252295e-189 [130,] 1.0000000 2.015884e-188 1.007942e-188 [131,] 1.0000000 1.607680e-187 8.038400e-188 [132,] 1.0000000 1.324437e-186 6.622186e-187 [133,] 1.0000000 1.030988e-185 5.154942e-186 [134,] 1.0000000 8.278354e-185 4.139177e-185 [135,] 1.0000000 6.402712e-184 3.201356e-184 [136,] 1.0000000 4.954088e-183 2.477044e-183 [137,] 1.0000000 3.848546e-182 1.924273e-182 [138,] 1.0000000 3.328470e-181 1.664235e-181 [139,] 1.0000000 2.674965e-180 1.337482e-180 [140,] 1.0000000 2.070347e-179 1.035173e-179 [141,] 1.0000000 1.745965e-178 8.729823e-179 [142,] 1.0000000 1.359607e-177 6.798034e-178 [143,] 1.0000000 1.119825e-176 5.599125e-177 [144,] 1.0000000 8.307593e-176 4.153796e-176 [145,] 1.0000000 6.464785e-175 3.232392e-175 [146,] 1.0000000 5.239891e-174 2.619945e-174 [147,] 1.0000000 4.797956e-173 2.398978e-173 [148,] 1.0000000 3.858196e-172 1.929098e-172 [149,] 1.0000000 3.179787e-171 1.589893e-171 [150,] 1.0000000 2.499233e-170 1.249616e-170 [151,] 1.0000000 1.933679e-169 9.668393e-170 [152,] 1.0000000 1.495589e-168 7.477945e-169 [153,] 1.0000000 1.173045e-167 5.865223e-168 [154,] 1.0000000 7.802621e-167 3.901311e-167 [155,] 1.0000000 6.022576e-166 3.011288e-166 [156,] 1.0000000 4.644313e-165 2.322156e-165 [157,] 1.0000000 3.577677e-164 1.788838e-164 [158,] 1.0000000 2.752752e-163 1.376376e-163 [159,] 1.0000000 2.250704e-162 1.125352e-162 [160,] 1.0000000 1.943666e-161 9.718332e-162 [161,] 1.0000000 1.496661e-160 7.483306e-161 [162,] 1.0000000 1.141159e-159 5.705795e-160 [163,] 1.0000000 2.817966e-159 1.408983e-159 [164,] 1.0000000 1.692629e-158 8.463147e-159 [165,] 1.0000000 1.284836e-157 6.424178e-158 [166,] 1.0000000 9.735294e-157 4.867647e-157 [167,] 1.0000000 7.566151e-156 3.783076e-156 [168,] 1.0000000 5.708948e-155 2.854474e-155 [169,] 1.0000000 4.568908e-154 2.284454e-154 [170,] 1.0000000 3.430162e-153 1.715081e-153 [171,] 1.0000000 2.665102e-152 1.332551e-152 [172,] 1.0000000 1.121869e-151 5.609345e-152 [173,] 1.0000000 9.169723e-151 4.584862e-151 [174,] 1.0000000 6.791610e-150 3.395805e-150 [175,] 1.0000000 5.018933e-149 2.509466e-149 [176,] 1.0000000 2.223902e-148 1.111951e-148 [177,] 1.0000000 1.664524e-147 8.322621e-148 [178,] 1.0000000 1.314104e-146 6.570518e-147 [179,] 1.0000000 1.027910e-145 5.139552e-146 [180,] 1.0000000 8.210924e-145 4.105462e-145 [181,] 1.0000000 5.866128e-144 2.933064e-144 [182,] 1.0000000 3.143643e-143 1.571821e-143 [183,] 1.0000000 2.338737e-142 1.169368e-142 [184,] 1.0000000 1.705457e-141 8.527283e-142 [185,] 1.0000000 1.219372e-140 6.096860e-141 [186,] 1.0000000 8.803511e-140 4.401756e-140 [187,] 1.0000000 6.259172e-139 3.129586e-139 [188,] 1.0000000 4.563136e-138 2.281568e-138 [189,] 1.0000000 3.225440e-137 1.612720e-137 [190,] 1.0000000 2.299282e-136 1.149641e-136 [191,] 1.0000000 1.740652e-135 8.703262e-136 [192,] 1.0000000 8.219029e-136 4.109515e-136 [193,] 1.0000000 6.213304e-135 3.106652e-135 [194,] 1.0000000 4.674899e-134 2.337450e-134 [195,] 1.0000000 2.794588e-133 1.397294e-133 [196,] 1.0000000 2.087543e-132 1.043772e-132 [197,] 1.0000000 1.552039e-131 7.760193e-132 [198,] 1.0000000 1.148432e-130 5.742158e-131 [199,] 1.0000000 8.534374e-130 4.267187e-130 [200,] 1.0000000 5.910929e-129 2.955464e-129 [201,] 1.0000000 4.304301e-128 2.152151e-128 [202,] 1.0000000 3.087786e-127 1.543893e-127 [203,] 1.0000000 2.228433e-126 1.114216e-126 [204,] 1.0000000 1.177517e-125 5.887585e-126 [205,] 1.0000000 8.568856e-125 4.284428e-125 [206,] 1.0000000 5.743945e-124 2.871972e-124 [207,] 1.0000000 3.839176e-123 1.919588e-123 [208,] 1.0000000 2.396006e-122 1.198003e-122 [209,] 1.0000000 1.617302e-121 8.086509e-122 [210,] 1.0000000 1.111701e-120 5.558507e-121 [211,] 1.0000000 7.386867e-120 3.693433e-120 [212,] 1.0000000 4.944522e-119 2.472261e-119 [213,] 1.0000000 3.276635e-118 1.638317e-118 [214,] 1.0000000 2.154521e-117 1.077260e-117 [215,] 1.0000000 1.395665e-116 6.978324e-117 [216,] 1.0000000 9.231025e-116 4.615512e-116 [217,] 1.0000000 5.964796e-115 2.982398e-115 [218,] 1.0000000 3.854450e-114 1.927225e-114 [219,] 1.0000000 2.593910e-113 1.296955e-113 [220,] 1.0000000 1.662066e-112 8.310331e-113 [221,] 1.0000000 1.097416e-111 5.487082e-112 [222,] 1.0000000 7.140217e-111 3.570109e-111 [223,] 1.0000000 4.511684e-110 2.255842e-110 [224,] 1.0000000 2.884669e-109 1.442334e-109 [225,] 1.0000000 1.809976e-108 9.049880e-109 [226,] 1.0000000 1.172815e-107 5.864074e-108 [227,] 1.0000000 7.247830e-107 3.623915e-107 [228,] 1.0000000 4.497404e-106 2.248702e-106 [229,] 1.0000000 2.923441e-105 1.461721e-105 [230,] 1.0000000 1.800147e-104 9.000733e-105 [231,] 1.0000000 1.104032e-103 5.520159e-104 [232,] 1.0000000 6.748767e-103 3.374384e-103 [233,] 1.0000000 4.160306e-102 2.080153e-102 [234,] 1.0000000 2.523734e-101 1.261867e-101 [235,] 1.0000000 1.525225e-100 7.626124e-101 [236,] 1.0000000 9.182982e-100 4.591491e-100 [237,] 1.0000000 5.590293e-99 2.795147e-99 [238,] 1.0000000 3.339761e-98 1.669881e-98 [239,] 1.0000000 2.013517e-97 1.006759e-97 [240,] 1.0000000 1.147796e-96 5.738981e-97 [241,] 1.0000000 7.010722e-96 3.505361e-96 [242,] 1.0000000 4.259977e-95 2.129988e-95 [243,] 1.0000000 2.594147e-94 1.297073e-94 [244,] 1.0000000 1.560807e-93 7.804037e-94 [245,] 1.0000000 9.341138e-93 4.670569e-93 [246,] 1.0000000 5.580283e-92 2.790142e-92 [247,] 1.0000000 3.306510e-91 1.653255e-91 [248,] 1.0000000 1.931361e-90 9.656807e-91 [249,] 1.0000000 1.133052e-89 5.665262e-90 [250,] 1.0000000 4.781360e-89 2.390680e-89 [251,] 1.0000000 2.706072e-88 1.353036e-88 [252,] 1.0000000 1.568709e-87 7.843546e-88 [253,] 1.0000000 9.206815e-87 4.603408e-87 [254,] 1.0000000 5.284544e-86 2.642272e-86 [255,] 1.0000000 2.939047e-85 1.469524e-85 [256,] 1.0000000 1.670017e-84 8.350085e-85 [257,] 1.0000000 9.435649e-84 4.717824e-84 [258,] 1.0000000 4.280041e-83 2.140021e-83 [259,] 1.0000000 2.162712e-82 1.081356e-82 [260,] 1.0000000 1.206026e-81 6.030128e-82 [261,] 1.0000000 6.537767e-81 3.268884e-81 [262,] 1.0000000 3.607863e-80 1.803931e-80 [263,] 1.0000000 8.882770e-80 4.441385e-80 [264,] 1.0000000 4.874070e-79 2.437035e-79 [265,] 1.0000000 2.395643e-78 1.197822e-78 [266,] 1.0000000 1.301515e-77 6.507575e-78 [267,] 1.0000000 7.026577e-77 3.513289e-77 [268,] 1.0000000 3.769073e-76 1.884537e-76 [269,] 1.0000000 2.043665e-75 1.021832e-75 [270,] 1.0000000 1.067336e-74 5.336678e-75 [271,] 1.0000000 5.521608e-74 2.760804e-74 [272,] 1.0000000 2.833298e-73 1.416649e-73 [273,] 1.0000000 1.459212e-72 7.296059e-73 [274,] 1.0000000 7.516732e-72 3.758366e-72 [275,] 1.0000000 3.931302e-71 1.965651e-71 [276,] 1.0000000 2.028213e-70 1.014107e-70 [277,] 1.0000000 1.025376e-69 5.126879e-70 [278,] 1.0000000 5.248852e-69 2.624426e-69 [279,] 1.0000000 2.628153e-68 1.314076e-68 [280,] 1.0000000 1.214814e-67 6.074072e-68 [281,] 1.0000000 6.026698e-67 3.013349e-67 [282,] 1.0000000 2.963398e-66 1.481699e-66 [283,] 1.0000000 1.447834e-65 7.239170e-66 [284,] 1.0000000 7.157983e-65 3.578991e-65 [285,] 1.0000000 3.606523e-64 1.803261e-64 [286,] 1.0000000 1.760832e-63 8.804158e-64 [287,] 1.0000000 8.210000e-63 4.105000e-63 [288,] 1.0000000 2.489541e-62 1.244770e-62 [289,] 1.0000000 1.190921e-61 5.954603e-62 [290,] 1.0000000 5.886735e-61 2.943368e-61 [291,] 1.0000000 2.331152e-60 1.165576e-60 [292,] 1.0000000 1.100140e-59 5.500702e-60 [293,] 1.0000000 5.159553e-59 2.579777e-59 [294,] 1.0000000 2.410856e-58 1.205428e-58 [295,] 1.0000000 1.130542e-57 5.652709e-58 [296,] 1.0000000 5.244498e-57 2.622249e-57 [297,] 1.0000000 2.482107e-56 1.241053e-56 [298,] 1.0000000 1.137868e-55 5.689340e-56 [299,] 1.0000000 5.135263e-55 2.567632e-55 [300,] 1.0000000 1.946900e-54 9.734502e-55 [301,] 1.0000000 8.862554e-54 4.431277e-54 [302,] 1.0000000 3.925449e-53 1.962724e-53 [303,] 1.0000000 1.675474e-52 8.377371e-53 [304,] 1.0000000 7.640250e-52 3.820125e-52 [305,] 1.0000000 3.287417e-51 1.643708e-51 [306,] 1.0000000 1.395196e-50 6.975980e-51 [307,] 1.0000000 4.961187e-50 2.480594e-50 [308,] 1.0000000 2.069994e-49 1.034997e-49 [309,] 1.0000000 8.001725e-49 4.000862e-49 [310,] 1.0000000 3.471617e-48 1.735808e-48 [311,] 1.0000000 1.302932e-47 6.514661e-48 [312,] 1.0000000 5.564422e-47 2.782211e-47 [313,] 1.0000000 2.461785e-46 1.230893e-46 [314,] 1.0000000 1.042985e-45 5.214923e-46 [315,] 1.0000000 4.372211e-45 2.186106e-45 [316,] 1.0000000 1.822279e-44 9.111393e-45 [317,] 1.0000000 7.550888e-44 3.775444e-44 [318,] 1.0000000 3.110496e-43 1.555248e-43 [319,] 1.0000000 1.191362e-42 5.956812e-43 [320,] 1.0000000 4.853245e-42 2.426623e-42 [321,] 1.0000000 1.985223e-41 9.926113e-42 [322,] 1.0000000 7.342479e-41 3.671240e-41 [323,] 1.0000000 2.634627e-40 1.317313e-40 [324,] 1.0000000 1.049485e-39 5.247427e-40 [325,] 1.0000000 4.154705e-39 2.077353e-39 [326,] 1.0000000 1.634505e-38 8.172524e-39 [327,] 1.0000000 6.389832e-38 3.194916e-38 [328,] 1.0000000 2.482131e-37 1.241066e-37 [329,] 1.0000000 9.667222e-37 4.833611e-37 [330,] 1.0000000 3.710400e-36 1.855200e-36 [331,] 1.0000000 1.418496e-35 7.092482e-36 [332,] 1.0000000 5.608775e-35 2.804387e-35 [333,] 1.0000000 2.108687e-34 1.054344e-34 [334,] 1.0000000 7.921023e-34 3.960511e-34 [335,] 1.0000000 2.934168e-33 1.467084e-33 [336,] 1.0000000 1.129957e-32 5.649785e-33 [337,] 1.0000000 4.134485e-32 2.067243e-32 [338,] 1.0000000 1.502245e-31 7.511226e-32 [339,] 1.0000000 5.419857e-31 2.709929e-31 [340,] 1.0000000 1.995358e-30 9.976789e-31 [341,] 1.0000000 7.096139e-30 3.548070e-30 [342,] 1.0000000 2.587570e-29 1.293785e-29 [343,] 1.0000000 9.069535e-29 4.534767e-29 [344,] 1.0000000 3.155273e-28 1.577637e-28 [345,] 1.0000000 1.080671e-27 5.403355e-28 [346,] 1.0000000 3.667457e-27 1.833728e-27 [347,] 1.0000000 1.258212e-26 6.291059e-27 [348,] 1.0000000 4.097393e-26 2.048696e-26 [349,] 1.0000000 1.377795e-25 6.888977e-26 [350,] 1.0000000 4.574410e-25 2.287205e-25 [351,] 1.0000000 1.506466e-24 7.532332e-25 [352,] 1.0000000 4.827474e-24 2.413737e-24 [353,] 1.0000000 1.563724e-23 7.818622e-24 [354,] 1.0000000 5.022637e-23 2.511318e-23 [355,] 1.0000000 1.599494e-22 7.997471e-23 [356,] 1.0000000 5.063385e-22 2.531692e-22 [357,] 1.0000000 1.584486e-21 7.922432e-22 [358,] 1.0000000 4.914154e-21 2.457077e-21 [359,] 1.0000000 1.436833e-20 7.184164e-21 [360,] 1.0000000 4.395458e-20 2.197729e-20 [361,] 1.0000000 1.329370e-19 6.646848e-20 [362,] 1.0000000 3.982698e-19 1.991349e-19 [363,] 1.0000000 1.181761e-18 5.908806e-19 [364,] 1.0000000 3.472403e-18 1.736201e-18 [365,] 1.0000000 1.014466e-17 5.072328e-18 [366,] 1.0000000 2.921684e-17 1.460842e-17 [367,] 1.0000000 8.109839e-17 4.054920e-17 [368,] 1.0000000 2.083822e-16 1.041911e-16 [369,] 1.0000000 5.986175e-16 2.993087e-16 [370,] 1.0000000 1.652587e-15 8.262936e-16 [371,] 1.0000000 4.511336e-15 2.255668e-15 [372,] 1.0000000 1.296956e-14 6.484782e-15 [373,] 1.0000000 3.386516e-14 1.693258e-14 [374,] 1.0000000 8.847481e-14 4.423741e-14 [375,] 1.0000000 2.306212e-13 1.153106e-13 [376,] 1.0000000 5.936228e-13 2.968114e-13 [377,] 1.0000000 1.517267e-12 7.586334e-13 [378,] 1.0000000 4.179049e-12 2.089525e-12 [379,] 1.0000000 1.037201e-11 5.186005e-12 [380,] 1.0000000 2.545341e-11 1.272671e-11 [381,] 1.0000000 6.573680e-11 3.286840e-11 [382,] 1.0000000 1.505304e-10 7.526522e-11 [383,] 1.0000000 3.660334e-10 1.830167e-10 [384,] 1.0000000 8.585822e-10 4.292911e-10 [385,] 1.0000000 2.037454e-09 1.018727e-09 [386,] 1.0000000 4.511531e-09 2.255765e-09 [387,] 1.0000000 1.085894e-08 5.429469e-09 [388,] 1.0000000 2.695118e-08 1.347559e-08 [389,] 1.0000000 4.361734e-08 2.180867e-08 [390,] 0.9999999 1.037101e-07 5.185506e-08 [391,] 0.9999999 2.411496e-07 1.205748e-07 [392,] 0.9999997 5.663044e-07 2.831522e-07 [393,] 0.9999998 4.405580e-07 2.202790e-07 [394,] 0.9999995 1.010151e-06 5.050757e-07 [395,] 0.9999988 2.341153e-06 1.170577e-06 [396,] 0.9999975 5.089456e-06 2.544728e-06 [397,] 0.9999951 9.744057e-06 4.872029e-06 [398,] 0.9999896 2.070270e-05 1.035135e-05 [399,] 0.9999776 4.471630e-05 2.235815e-05 [400,] 0.9999523 9.537199e-05 4.768599e-05 [401,] 0.9999134 1.732003e-04 8.660015e-05 [402,] 0.9998414 3.172405e-04 1.586202e-04 [403,] 0.9997054 5.891871e-04 2.945936e-04 [404,] 0.9994584 1.083114e-03 5.415568e-04 [405,] 0.9989082 2.183586e-03 1.091793e-03 [406,] 0.9981264 3.747132e-03 1.873566e-03 [407,] 0.9965397 6.920645e-03 3.460322e-03 [408,] 0.9936916 1.261690e-02 6.308449e-03 [409,] 0.9962002 7.599505e-03 3.799753e-03 [410,] 0.9926660 1.466805e-02 7.334027e-03 [411,] 0.9858898 2.822038e-02 1.411019e-02 [412,] 0.9858736 2.825282e-02 1.412641e-02 [413,] 0.9861077 2.778468e-02 1.389234e-02 [414,] 0.9761110 4.777808e-02 2.388904e-02 [415,] 0.9594301 8.113985e-02 4.056993e-02 [416,] 0.9381104 1.237791e-01 6.188955e-02 [417,] 0.8907581 2.184838e-01 1.092419e-01 [418,] 0.8826826 2.346348e-01 1.173174e-01 [419,] 0.8624195 2.751610e-01 1.375805e-01 [420,] 0.7740979 4.518042e-01 2.259021e-01 > postscript(file="/var/www/html/rcomp/tmp/1hxin1291232858.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/html/rcomp/tmp/2hxin1291232858.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/html/rcomp/tmp/3r6hq1291232858.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/html/rcomp/tmp/4r6hq1291232858.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/html/rcomp/tmp/5r6hq1291232858.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 = 431 Frequency = 1 1 2 3 4 5 2559625.7921 2146005.8721 2073314.5816 -1833510.1393 -133727.1739 6 7 8 9 10 631218.8668 45603.9879 506071.3904 -270524.0886 611450.3878 11 12 13 14 15 397635.4196 232695.3330 -552373.2158 -814949.5497 -546748.1166 16 17 18 19 20 570626.1725 -28097.3485 -554963.3777 -287336.7841 -968495.6350 21 22 23 24 25 263815.8794 144146.5750 31189.5356 -1109804.1420 -287950.1460 26 27 28 29 30 114890.9839 -630335.0160 -333838.9868 409654.5342 -384864.6377 31 32 33 34 35 295744.3278 -6272.7728 288690.7049 136548.4730 360943.0845 36 37 38 39 40 -753672.9841 92712.8973 74215.0680 -85424.2798 -583164.7103 41 42 43 44 45 -795895.0567 -125041.0161 -109757.8415 -227260.7256 -10876.0454 46 47 48 49 50 -11466.5563 -120249.3384 -486531.9307 -60347.3817 -212384.2772 51 52 53 54 55 -135204.0912 -343612.4237 -204412.9027 -23748.9613 -534415.8734 56 57 58 59 60 -154377.4505 395256.7801 13755.0272 -108168.4365 -256859.5116 61 62 63 64 65 27900.3571 -222015.2339 25690.4437 -246149.5843 -107119.7842 66 67 68 69 70 -944371.6210 -34820.2632 -72626.7588 -168353.9868 -219660.5116 71 72 73 74 75 -289031.3384 -15562.9295 -70530.5817 55280.4323 17474.5342 76 77 78 79 80 -120395.9767 12358.1138 -42222.7996 -141393.5524 2221.8348 81 82 83 84 85 -8334.8862 -137279.3091 -19922.3499 5844.7368 81479.7368 86 87 88 89 90 -214773.5409 -200058.1167 66717.9686 -115293.7409 94438.8527 91 92 93 94 95 7781.3024 -74209.2951 -69512.4071 -40985.0200 -56754.5091 96 97 98 99 100 -33446.1359 158654.4323 -45436.2085 90471.2731 -2669.7015 101 102 103 104 105 -106822.2518 -77955.7269 -167317.9002 47878.5482 16729.2298 106 107 108 109 110 -86058.8862 -361433.2225 160430.4870 -110559.9728 -117557.9767 111 112 113 114 115 83002.9253 -264384.4836 -127512.6250 -48101.1652 70276.6463 116 117 118 119 120 -37832.5817 66548.5482 138319.5049 -114872.9333 -5926.9588 121 122 123 124 125 -249101.1893 51033.0845 -353694.5078 91610.4437 -1608.6289 126 127 128 129 130 -163052.3244 83290.4756 39812.3164 -12085.1906 83290.4756 131 132 133 134 135 52421.6641 30780.5775 83290.4756 -4505.8429 62933.0119 136 137 138 139 140 65867.0119 44182.3164 78534.3597 65771.0119 83290.4756 141 142 143 144 145 59311.5482 83290.4756 19584.2004 66692.0119 33747.7368 146 147 148 149 150 42420.4323 83290.4756 53055.6641 84633.2438 83787.4756 151 152 153 154 155 29718.6641 -122233.8568 61974.8960 35833.6641 74629.2438 156 157 158 159 160 83290.4756 83290.4756 74620.2438 83019.8960 83290.4756 161 162 163 164 165 83290.4756 83290.4756 83290.4756 44336.6641 -32217.3065 166 167 168 169 170 78529.3597 83290.4756 73886.3890 -124990.9002 83290.4756 171 172 173 174 175 83290.4756 67533.1278 83290.4756 39138.4323 83290.4756 176 177 178 179 180 49352.5482 -400554.6237 33982.3418 83290.4756 83290.4756 181 182 183 184 185 46246.2731 980.0667 -24378.4225 -52990.5130 11296.7622 186 187 188 189 190 4122.0667 -28137.3244 55663.7801 38381.8960 83290.4756 191 192 193 194 195 56813.6641 83290.4756 53434.6641 83290.4756 70809.1278 196 197 198 199 200 33186.3418 -251968.9753 33186.3418 33186.3418 -93605.8607 201 202 203 204 205 33186.3418 33186.3418 33186.3418 -32389.0493 -7872.0493 206 207 208 209 210 12645.7622 27854.1100 33186.3418 -49421.6722 11964.9940 211 212 213 214 215 83290.4756 83290.4756 51075.2004 14264.3164 51079.8960 216 217 218 219 220 74764.2438 58019.7801 55804.6641 74764.2438 -29125.7702 221 222 223 224 225 51777.1278 83290.4756 79001.3597 -5190.8429 79223.3597 226 227 228 229 230 -10956.0747 44093.4323 83290.4756 71723.2438 83290.4756 231 232 233 234 235 36768.6641 50918.4323 83290.4756 15719.9686 83290.4756 236 237 238 239 240 83397.4756 83290.4756 58124.7801 83290.4756 83290.4756 241 242 243 244 245 83290.4756 68008.1278 83290.4756 53705.6641 83769.2438 246 247 248 249 250 33186.3418 33186.3418 12952.7368 33186.3418 33186.3418 251 252 253 254 255 26244.4323 33186.3418 70969.3597 -89711.9295 -221499.1320 256 257 258 259 260 83290.4756 47628.6641 -44464.5384 33186.3418 83290.4756 261 262 263 264 265 33186.3418 33186.3418 56418.9686 57234.3164 33186.3418 266 267 268 269 270 83290.4756 33186.3418 110722.4323 33186.3418 -11510.5384 271 272 273 274 275 33186.3418 33186.3418 33186.3418 -2887.7269 83290.4756 276 277 278 279 280 36769.0845 -9743.1906 83290.4756 -113983.0200 44475.5482 281 282 283 284 285 33186.3418 83290.4756 -29439.5384 83290.4756 60865.4323 286 287 288 289 290 83290.4756 15700.7622 73219.1278 33186.3418 5243.9686 291 292 293 294 295 33186.3418 41929.0845 52870.2731 83290.4756 27688.3164 296 297 298 299 300 20406.2985 83290.4756 33186.3418 83290.4756 44922.4323 301 302 303 304 305 13917.8781 64811.3597 79027.3597 33186.3418 52651.8527 306 307 308 309 310 32147.0845 33186.3418 -46545.8607 5273.5049 33186.3418 311 312 313 314 315 33186.3418 -21898.2339 33186.3418 -19846.2811 49461.5482 316 317 318 319 320 90303.0119 82874.4756 -20210.8429 61966.8960 83290.4756 321 322 323 324 325 83290.4756 83290.4756 83290.4756 33186.3418 83290.4756 326 327 328 329 330 76149.3597 33186.3418 33186.3418 83290.4756 83290.4756 331 332 333 334 335 83290.4756 83290.4756 83290.4756 64587.0119 79035.3597 336 337 338 339 340 66229.0119 -21212.8429 83290.4756 77527.3597 42977.3164 341 342 343 344 345 12333.0845 83290.4756 83290.4756 83290.4756 20745.9686 346 347 348 349 350 83290.4756 58798.1278 83290.4756 83290.4756 44648.3164 351 352 353 354 355 73811.1278 52117.6641 50978.3164 77857.3597 83290.4756 356 357 358 359 360 83290.4756 26341.7368 83290.4756 83290.4756 83290.4756 361 362 363 364 365 66256.0119 83290.4756 83290.4756 -93224.9040 45928.4323 366 367 368 369 370 83290.4756 83290.4756 83290.4756 83290.4756 53621.6641 371 372 373 374 375 83290.4756 61429.6641 74707.6641 62545.1278 83290.4756 376 377 378 379 380 83290.4756 16386.2004 60408.6641 45143.3164 83290.4756 381 382 383 384 385 83290.4756 8465.3890 -11375.1473 83290.4756 34083.0845 386 387 388 389 390 -117594.8862 82192.1278 51369.7801 69331.1278 41735.5482 391 392 393 394 395 47606.3164 34467.7801 -2918.4518 109986.7801 -36651.7015 396 397 398 399 400 17392.5482 -149384.6722 -862585.3637 34329.6641 -226360.4110 401 402 403 404 405 59373.1278 67310.5482 39250.7801 -92729.0926 14966.6641 406 407 408 409 410 48023.2004 38062.2004 21308.0845 27638.4323 -13017.8722 411 412 413 414 415 -165665.1320 -84756.0926 -76387.3970 29629.6323 -113579.7996 416 417 418 419 420 -36290.7562 77017.8094 -76979.9767 -59093.3065 10188.3597 421 422 423 424 425 17115.3597 -306663.1027 -655350.5409 -359494.8429 -44638.8429 426 427 428 429 430 -416407.1639 -278860.9588 -966014.8096 -114972.2339 -226086.4938 431 -533414.5231 > postscript(file="/var/www/html/rcomp/tmp/62fzt1291232858.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 = 431 Frequency = 1 lag(myerror, k = 1) myerror 0 2559625.7921 NA 1 2146005.8721 2559625.7921 2 2073314.5816 2146005.8721 3 -1833510.1393 2073314.5816 4 -133727.1739 -1833510.1393 5 631218.8668 -133727.1739 6 45603.9879 631218.8668 7 506071.3904 45603.9879 8 -270524.0886 506071.3904 9 611450.3878 -270524.0886 10 397635.4196 611450.3878 11 232695.3330 397635.4196 12 -552373.2158 232695.3330 13 -814949.5497 -552373.2158 14 -546748.1166 -814949.5497 15 570626.1725 -546748.1166 16 -28097.3485 570626.1725 17 -554963.3777 -28097.3485 18 -287336.7841 -554963.3777 19 -968495.6350 -287336.7841 20 263815.8794 -968495.6350 21 144146.5750 263815.8794 22 31189.5356 144146.5750 23 -1109804.1420 31189.5356 24 -287950.1460 -1109804.1420 25 114890.9839 -287950.1460 26 -630335.0160 114890.9839 27 -333838.9868 -630335.0160 28 409654.5342 -333838.9868 29 -384864.6377 409654.5342 30 295744.3278 -384864.6377 31 -6272.7728 295744.3278 32 288690.7049 -6272.7728 33 136548.4730 288690.7049 34 360943.0845 136548.4730 35 -753672.9841 360943.0845 36 92712.8973 -753672.9841 37 74215.0680 92712.8973 38 -85424.2798 74215.0680 39 -583164.7103 -85424.2798 40 -795895.0567 -583164.7103 41 -125041.0161 -795895.0567 42 -109757.8415 -125041.0161 43 -227260.7256 -109757.8415 44 -10876.0454 -227260.7256 45 -11466.5563 -10876.0454 46 -120249.3384 -11466.5563 47 -486531.9307 -120249.3384 48 -60347.3817 -486531.9307 49 -212384.2772 -60347.3817 50 -135204.0912 -212384.2772 51 -343612.4237 -135204.0912 52 -204412.9027 -343612.4237 53 -23748.9613 -204412.9027 54 -534415.8734 -23748.9613 55 -154377.4505 -534415.8734 56 395256.7801 -154377.4505 57 13755.0272 395256.7801 58 -108168.4365 13755.0272 59 -256859.5116 -108168.4365 60 27900.3571 -256859.5116 61 -222015.2339 27900.3571 62 25690.4437 -222015.2339 63 -246149.5843 25690.4437 64 -107119.7842 -246149.5843 65 -944371.6210 -107119.7842 66 -34820.2632 -944371.6210 67 -72626.7588 -34820.2632 68 -168353.9868 -72626.7588 69 -219660.5116 -168353.9868 70 -289031.3384 -219660.5116 71 -15562.9295 -289031.3384 72 -70530.5817 -15562.9295 73 55280.4323 -70530.5817 74 17474.5342 55280.4323 75 -120395.9767 17474.5342 76 12358.1138 -120395.9767 77 -42222.7996 12358.1138 78 -141393.5524 -42222.7996 79 2221.8348 -141393.5524 80 -8334.8862 2221.8348 81 -137279.3091 -8334.8862 82 -19922.3499 -137279.3091 83 5844.7368 -19922.3499 84 81479.7368 5844.7368 85 -214773.5409 81479.7368 86 -200058.1167 -214773.5409 87 66717.9686 -200058.1167 88 -115293.7409 66717.9686 89 94438.8527 -115293.7409 90 7781.3024 94438.8527 91 -74209.2951 7781.3024 92 -69512.4071 -74209.2951 93 -40985.0200 -69512.4071 94 -56754.5091 -40985.0200 95 -33446.1359 -56754.5091 96 158654.4323 -33446.1359 97 -45436.2085 158654.4323 98 90471.2731 -45436.2085 99 -2669.7015 90471.2731 100 -106822.2518 -2669.7015 101 -77955.7269 -106822.2518 102 -167317.9002 -77955.7269 103 47878.5482 -167317.9002 104 16729.2298 47878.5482 105 -86058.8862 16729.2298 106 -361433.2225 -86058.8862 107 160430.4870 -361433.2225 108 -110559.9728 160430.4870 109 -117557.9767 -110559.9728 110 83002.9253 -117557.9767 111 -264384.4836 83002.9253 112 -127512.6250 -264384.4836 113 -48101.1652 -127512.6250 114 70276.6463 -48101.1652 115 -37832.5817 70276.6463 116 66548.5482 -37832.5817 117 138319.5049 66548.5482 118 -114872.9333 138319.5049 119 -5926.9588 -114872.9333 120 -249101.1893 -5926.9588 121 51033.0845 -249101.1893 122 -353694.5078 51033.0845 123 91610.4437 -353694.5078 124 -1608.6289 91610.4437 125 -163052.3244 -1608.6289 126 83290.4756 -163052.3244 127 39812.3164 83290.4756 128 -12085.1906 39812.3164 129 83290.4756 -12085.1906 130 52421.6641 83290.4756 131 30780.5775 52421.6641 132 83290.4756 30780.5775 133 -4505.8429 83290.4756 134 62933.0119 -4505.8429 135 65867.0119 62933.0119 136 44182.3164 65867.0119 137 78534.3597 44182.3164 138 65771.0119 78534.3597 139 83290.4756 65771.0119 140 59311.5482 83290.4756 141 83290.4756 59311.5482 142 19584.2004 83290.4756 143 66692.0119 19584.2004 144 33747.7368 66692.0119 145 42420.4323 33747.7368 146 83290.4756 42420.4323 147 53055.6641 83290.4756 148 84633.2438 53055.6641 149 83787.4756 84633.2438 150 29718.6641 83787.4756 151 -122233.8568 29718.6641 152 61974.8960 -122233.8568 153 35833.6641 61974.8960 154 74629.2438 35833.6641 155 83290.4756 74629.2438 156 83290.4756 83290.4756 157 74620.2438 83290.4756 158 83019.8960 74620.2438 159 83290.4756 83019.8960 160 83290.4756 83290.4756 161 83290.4756 83290.4756 162 83290.4756 83290.4756 163 44336.6641 83290.4756 164 -32217.3065 44336.6641 165 78529.3597 -32217.3065 166 83290.4756 78529.3597 167 73886.3890 83290.4756 168 -124990.9002 73886.3890 169 83290.4756 -124990.9002 170 83290.4756 83290.4756 171 67533.1278 83290.4756 172 83290.4756 67533.1278 173 39138.4323 83290.4756 174 83290.4756 39138.4323 175 49352.5482 83290.4756 176 -400554.6237 49352.5482 177 33982.3418 -400554.6237 178 83290.4756 33982.3418 179 83290.4756 83290.4756 180 46246.2731 83290.4756 181 980.0667 46246.2731 182 -24378.4225 980.0667 183 -52990.5130 -24378.4225 184 11296.7622 -52990.5130 185 4122.0667 11296.7622 186 -28137.3244 4122.0667 187 55663.7801 -28137.3244 188 38381.8960 55663.7801 189 83290.4756 38381.8960 190 56813.6641 83290.4756 191 83290.4756 56813.6641 192 53434.6641 83290.4756 193 83290.4756 53434.6641 194 70809.1278 83290.4756 195 33186.3418 70809.1278 196 -251968.9753 33186.3418 197 33186.3418 -251968.9753 198 33186.3418 33186.3418 199 -93605.8607 33186.3418 200 33186.3418 -93605.8607 201 33186.3418 33186.3418 202 33186.3418 33186.3418 203 -32389.0493 33186.3418 204 -7872.0493 -32389.0493 205 12645.7622 -7872.0493 206 27854.1100 12645.7622 207 33186.3418 27854.1100 208 -49421.6722 33186.3418 209 11964.9940 -49421.6722 210 83290.4756 11964.9940 211 83290.4756 83290.4756 212 51075.2004 83290.4756 213 14264.3164 51075.2004 214 51079.8960 14264.3164 215 74764.2438 51079.8960 216 58019.7801 74764.2438 217 55804.6641 58019.7801 218 74764.2438 55804.6641 219 -29125.7702 74764.2438 220 51777.1278 -29125.7702 221 83290.4756 51777.1278 222 79001.3597 83290.4756 223 -5190.8429 79001.3597 224 79223.3597 -5190.8429 225 -10956.0747 79223.3597 226 44093.4323 -10956.0747 227 83290.4756 44093.4323 228 71723.2438 83290.4756 229 83290.4756 71723.2438 230 36768.6641 83290.4756 231 50918.4323 36768.6641 232 83290.4756 50918.4323 233 15719.9686 83290.4756 234 83290.4756 15719.9686 235 83397.4756 83290.4756 236 83290.4756 83397.4756 237 58124.7801 83290.4756 238 83290.4756 58124.7801 239 83290.4756 83290.4756 240 83290.4756 83290.4756 241 68008.1278 83290.4756 242 83290.4756 68008.1278 243 53705.6641 83290.4756 244 83769.2438 53705.6641 245 33186.3418 83769.2438 246 33186.3418 33186.3418 247 12952.7368 33186.3418 248 33186.3418 12952.7368 249 33186.3418 33186.3418 250 26244.4323 33186.3418 251 33186.3418 26244.4323 252 70969.3597 33186.3418 253 -89711.9295 70969.3597 254 -221499.1320 -89711.9295 255 83290.4756 -221499.1320 256 47628.6641 83290.4756 257 -44464.5384 47628.6641 258 33186.3418 -44464.5384 259 83290.4756 33186.3418 260 33186.3418 83290.4756 261 33186.3418 33186.3418 262 56418.9686 33186.3418 263 57234.3164 56418.9686 264 33186.3418 57234.3164 265 83290.4756 33186.3418 266 33186.3418 83290.4756 267 110722.4323 33186.3418 268 33186.3418 110722.4323 269 -11510.5384 33186.3418 270 33186.3418 -11510.5384 271 33186.3418 33186.3418 272 33186.3418 33186.3418 273 -2887.7269 33186.3418 274 83290.4756 -2887.7269 275 36769.0845 83290.4756 276 -9743.1906 36769.0845 277 83290.4756 -9743.1906 278 -113983.0200 83290.4756 279 44475.5482 -113983.0200 280 33186.3418 44475.5482 281 83290.4756 33186.3418 282 -29439.5384 83290.4756 283 83290.4756 -29439.5384 284 60865.4323 83290.4756 285 83290.4756 60865.4323 286 15700.7622 83290.4756 287 73219.1278 15700.7622 288 33186.3418 73219.1278 289 5243.9686 33186.3418 290 33186.3418 5243.9686 291 41929.0845 33186.3418 292 52870.2731 41929.0845 293 83290.4756 52870.2731 294 27688.3164 83290.4756 295 20406.2985 27688.3164 296 83290.4756 20406.2985 297 33186.3418 83290.4756 298 83290.4756 33186.3418 299 44922.4323 83290.4756 300 13917.8781 44922.4323 301 64811.3597 13917.8781 302 79027.3597 64811.3597 303 33186.3418 79027.3597 304 52651.8527 33186.3418 305 32147.0845 52651.8527 306 33186.3418 32147.0845 307 -46545.8607 33186.3418 308 5273.5049 -46545.8607 309 33186.3418 5273.5049 310 33186.3418 33186.3418 311 -21898.2339 33186.3418 312 33186.3418 -21898.2339 313 -19846.2811 33186.3418 314 49461.5482 -19846.2811 315 90303.0119 49461.5482 316 82874.4756 90303.0119 317 -20210.8429 82874.4756 318 61966.8960 -20210.8429 319 83290.4756 61966.8960 320 83290.4756 83290.4756 321 83290.4756 83290.4756 322 83290.4756 83290.4756 323 33186.3418 83290.4756 324 83290.4756 33186.3418 325 76149.3597 83290.4756 326 33186.3418 76149.3597 327 33186.3418 33186.3418 328 83290.4756 33186.3418 329 83290.4756 83290.4756 330 83290.4756 83290.4756 331 83290.4756 83290.4756 332 83290.4756 83290.4756 333 64587.0119 83290.4756 334 79035.3597 64587.0119 335 66229.0119 79035.3597 336 -21212.8429 66229.0119 337 83290.4756 -21212.8429 338 77527.3597 83290.4756 339 42977.3164 77527.3597 340 12333.0845 42977.3164 341 83290.4756 12333.0845 342 83290.4756 83290.4756 343 83290.4756 83290.4756 344 20745.9686 83290.4756 345 83290.4756 20745.9686 346 58798.1278 83290.4756 347 83290.4756 58798.1278 348 83290.4756 83290.4756 349 44648.3164 83290.4756 350 73811.1278 44648.3164 351 52117.6641 73811.1278 352 50978.3164 52117.6641 353 77857.3597 50978.3164 354 83290.4756 77857.3597 355 83290.4756 83290.4756 356 26341.7368 83290.4756 357 83290.4756 26341.7368 358 83290.4756 83290.4756 359 83290.4756 83290.4756 360 66256.0119 83290.4756 361 83290.4756 66256.0119 362 83290.4756 83290.4756 363 -93224.9040 83290.4756 364 45928.4323 -93224.9040 365 83290.4756 45928.4323 366 83290.4756 83290.4756 367 83290.4756 83290.4756 368 83290.4756 83290.4756 369 53621.6641 83290.4756 370 83290.4756 53621.6641 371 61429.6641 83290.4756 372 74707.6641 61429.6641 373 62545.1278 74707.6641 374 83290.4756 62545.1278 375 83290.4756 83290.4756 376 16386.2004 83290.4756 377 60408.6641 16386.2004 378 45143.3164 60408.6641 379 83290.4756 45143.3164 380 83290.4756 83290.4756 381 8465.3890 83290.4756 382 -11375.1473 8465.3890 383 83290.4756 -11375.1473 384 34083.0845 83290.4756 385 -117594.8862 34083.0845 386 82192.1278 -117594.8862 387 51369.7801 82192.1278 388 69331.1278 51369.7801 389 41735.5482 69331.1278 390 47606.3164 41735.5482 391 34467.7801 47606.3164 392 -2918.4518 34467.7801 393 109986.7801 -2918.4518 394 -36651.7015 109986.7801 395 17392.5482 -36651.7015 396 -149384.6722 17392.5482 397 -862585.3637 -149384.6722 398 34329.6641 -862585.3637 399 -226360.4110 34329.6641 400 59373.1278 -226360.4110 401 67310.5482 59373.1278 402 39250.7801 67310.5482 403 -92729.0926 39250.7801 404 14966.6641 -92729.0926 405 48023.2004 14966.6641 406 38062.2004 48023.2004 407 21308.0845 38062.2004 408 27638.4323 21308.0845 409 -13017.8722 27638.4323 410 -165665.1320 -13017.8722 411 -84756.0926 -165665.1320 412 -76387.3970 -84756.0926 413 29629.6323 -76387.3970 414 -113579.7996 29629.6323 415 -36290.7562 -113579.7996 416 77017.8094 -36290.7562 417 -76979.9767 77017.8094 418 -59093.3065 -76979.9767 419 10188.3597 -59093.3065 420 17115.3597 10188.3597 421 -306663.1027 17115.3597 422 -655350.5409 -306663.1027 423 -359494.8429 -655350.5409 424 -44638.8429 -359494.8429 425 -416407.1639 -44638.8429 426 -278860.9588 -416407.1639 427 -966014.8096 -278860.9588 428 -114972.2339 -966014.8096 429 -226086.4938 -114972.2339 430 -533414.5231 -226086.4938 431 NA -533414.5231 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2146005.8721 2559625.7921 [2,] 2073314.5816 2146005.8721 [3,] -1833510.1393 2073314.5816 [4,] -133727.1739 -1833510.1393 [5,] 631218.8668 -133727.1739 [6,] 45603.9879 631218.8668 [7,] 506071.3904 45603.9879 [8,] -270524.0886 506071.3904 [9,] 611450.3878 -270524.0886 [10,] 397635.4196 611450.3878 [11,] 232695.3330 397635.4196 [12,] -552373.2158 232695.3330 [13,] -814949.5497 -552373.2158 [14,] -546748.1166 -814949.5497 [15,] 570626.1725 -546748.1166 [16,] -28097.3485 570626.1725 [17,] -554963.3777 -28097.3485 [18,] -287336.7841 -554963.3777 [19,] -968495.6350 -287336.7841 [20,] 263815.8794 -968495.6350 [21,] 144146.5750 263815.8794 [22,] 31189.5356 144146.5750 [23,] -1109804.1420 31189.5356 [24,] -287950.1460 -1109804.1420 [25,] 114890.9839 -287950.1460 [26,] -630335.0160 114890.9839 [27,] -333838.9868 -630335.0160 [28,] 409654.5342 -333838.9868 [29,] -384864.6377 409654.5342 [30,] 295744.3278 -384864.6377 [31,] -6272.7728 295744.3278 [32,] 288690.7049 -6272.7728 [33,] 136548.4730 288690.7049 [34,] 360943.0845 136548.4730 [35,] -753672.9841 360943.0845 [36,] 92712.8973 -753672.9841 [37,] 74215.0680 92712.8973 [38,] -85424.2798 74215.0680 [39,] -583164.7103 -85424.2798 [40,] -795895.0567 -583164.7103 [41,] -125041.0161 -795895.0567 [42,] -109757.8415 -125041.0161 [43,] -227260.7256 -109757.8415 [44,] -10876.0454 -227260.7256 [45,] -11466.5563 -10876.0454 [46,] -120249.3384 -11466.5563 [47,] -486531.9307 -120249.3384 [48,] -60347.3817 -486531.9307 [49,] -212384.2772 -60347.3817 [50,] -135204.0912 -212384.2772 [51,] -343612.4237 -135204.0912 [52,] -204412.9027 -343612.4237 [53,] -23748.9613 -204412.9027 [54,] -534415.8734 -23748.9613 [55,] -154377.4505 -534415.8734 [56,] 395256.7801 -154377.4505 [57,] 13755.0272 395256.7801 [58,] -108168.4365 13755.0272 [59,] -256859.5116 -108168.4365 [60,] 27900.3571 -256859.5116 [61,] -222015.2339 27900.3571 [62,] 25690.4437 -222015.2339 [63,] -246149.5843 25690.4437 [64,] -107119.7842 -246149.5843 [65,] -944371.6210 -107119.7842 [66,] -34820.2632 -944371.6210 [67,] -72626.7588 -34820.2632 [68,] -168353.9868 -72626.7588 [69,] -219660.5116 -168353.9868 [70,] -289031.3384 -219660.5116 [71,] -15562.9295 -289031.3384 [72,] -70530.5817 -15562.9295 [73,] 55280.4323 -70530.5817 [74,] 17474.5342 55280.4323 [75,] -120395.9767 17474.5342 [76,] 12358.1138 -120395.9767 [77,] -42222.7996 12358.1138 [78,] -141393.5524 -42222.7996 [79,] 2221.8348 -141393.5524 [80,] -8334.8862 2221.8348 [81,] -137279.3091 -8334.8862 [82,] -19922.3499 -137279.3091 [83,] 5844.7368 -19922.3499 [84,] 81479.7368 5844.7368 [85,] -214773.5409 81479.7368 [86,] -200058.1167 -214773.5409 [87,] 66717.9686 -200058.1167 [88,] -115293.7409 66717.9686 [89,] 94438.8527 -115293.7409 [90,] 7781.3024 94438.8527 [91,] -74209.2951 7781.3024 [92,] -69512.4071 -74209.2951 [93,] -40985.0200 -69512.4071 [94,] -56754.5091 -40985.0200 [95,] -33446.1359 -56754.5091 [96,] 158654.4323 -33446.1359 [97,] -45436.2085 158654.4323 [98,] 90471.2731 -45436.2085 [99,] -2669.7015 90471.2731 [100,] -106822.2518 -2669.7015 [101,] -77955.7269 -106822.2518 [102,] -167317.9002 -77955.7269 [103,] 47878.5482 -167317.9002 [104,] 16729.2298 47878.5482 [105,] -86058.8862 16729.2298 [106,] -361433.2225 -86058.8862 [107,] 160430.4870 -361433.2225 [108,] -110559.9728 160430.4870 [109,] -117557.9767 -110559.9728 [110,] 83002.9253 -117557.9767 [111,] -264384.4836 83002.9253 [112,] -127512.6250 -264384.4836 [113,] -48101.1652 -127512.6250 [114,] 70276.6463 -48101.1652 [115,] -37832.5817 70276.6463 [116,] 66548.5482 -37832.5817 [117,] 138319.5049 66548.5482 [118,] -114872.9333 138319.5049 [119,] -5926.9588 -114872.9333 [120,] -249101.1893 -5926.9588 [121,] 51033.0845 -249101.1893 [122,] -353694.5078 51033.0845 [123,] 91610.4437 -353694.5078 [124,] -1608.6289 91610.4437 [125,] -163052.3244 -1608.6289 [126,] 83290.4756 -163052.3244 [127,] 39812.3164 83290.4756 [128,] -12085.1906 39812.3164 [129,] 83290.4756 -12085.1906 [130,] 52421.6641 83290.4756 [131,] 30780.5775 52421.6641 [132,] 83290.4756 30780.5775 [133,] -4505.8429 83290.4756 [134,] 62933.0119 -4505.8429 [135,] 65867.0119 62933.0119 [136,] 44182.3164 65867.0119 [137,] 78534.3597 44182.3164 [138,] 65771.0119 78534.3597 [139,] 83290.4756 65771.0119 [140,] 59311.5482 83290.4756 [141,] 83290.4756 59311.5482 [142,] 19584.2004 83290.4756 [143,] 66692.0119 19584.2004 [144,] 33747.7368 66692.0119 [145,] 42420.4323 33747.7368 [146,] 83290.4756 42420.4323 [147,] 53055.6641 83290.4756 [148,] 84633.2438 53055.6641 [149,] 83787.4756 84633.2438 [150,] 29718.6641 83787.4756 [151,] -122233.8568 29718.6641 [152,] 61974.8960 -122233.8568 [153,] 35833.6641 61974.8960 [154,] 74629.2438 35833.6641 [155,] 83290.4756 74629.2438 [156,] 83290.4756 83290.4756 [157,] 74620.2438 83290.4756 [158,] 83019.8960 74620.2438 [159,] 83290.4756 83019.8960 [160,] 83290.4756 83290.4756 [161,] 83290.4756 83290.4756 [162,] 83290.4756 83290.4756 [163,] 44336.6641 83290.4756 [164,] -32217.3065 44336.6641 [165,] 78529.3597 -32217.3065 [166,] 83290.4756 78529.3597 [167,] 73886.3890 83290.4756 [168,] -124990.9002 73886.3890 [169,] 83290.4756 -124990.9002 [170,] 83290.4756 83290.4756 [171,] 67533.1278 83290.4756 [172,] 83290.4756 67533.1278 [173,] 39138.4323 83290.4756 [174,] 83290.4756 39138.4323 [175,] 49352.5482 83290.4756 [176,] -400554.6237 49352.5482 [177,] 33982.3418 -400554.6237 [178,] 83290.4756 33982.3418 [179,] 83290.4756 83290.4756 [180,] 46246.2731 83290.4756 [181,] 980.0667 46246.2731 [182,] -24378.4225 980.0667 [183,] -52990.5130 -24378.4225 [184,] 11296.7622 -52990.5130 [185,] 4122.0667 11296.7622 [186,] -28137.3244 4122.0667 [187,] 55663.7801 -28137.3244 [188,] 38381.8960 55663.7801 [189,] 83290.4756 38381.8960 [190,] 56813.6641 83290.4756 [191,] 83290.4756 56813.6641 [192,] 53434.6641 83290.4756 [193,] 83290.4756 53434.6641 [194,] 70809.1278 83290.4756 [195,] 33186.3418 70809.1278 [196,] -251968.9753 33186.3418 [197,] 33186.3418 -251968.9753 [198,] 33186.3418 33186.3418 [199,] -93605.8607 33186.3418 [200,] 33186.3418 -93605.8607 [201,] 33186.3418 33186.3418 [202,] 33186.3418 33186.3418 [203,] -32389.0493 33186.3418 [204,] -7872.0493 -32389.0493 [205,] 12645.7622 -7872.0493 [206,] 27854.1100 12645.7622 [207,] 33186.3418 27854.1100 [208,] -49421.6722 33186.3418 [209,] 11964.9940 -49421.6722 [210,] 83290.4756 11964.9940 [211,] 83290.4756 83290.4756 [212,] 51075.2004 83290.4756 [213,] 14264.3164 51075.2004 [214,] 51079.8960 14264.3164 [215,] 74764.2438 51079.8960 [216,] 58019.7801 74764.2438 [217,] 55804.6641 58019.7801 [218,] 74764.2438 55804.6641 [219,] -29125.7702 74764.2438 [220,] 51777.1278 -29125.7702 [221,] 83290.4756 51777.1278 [222,] 79001.3597 83290.4756 [223,] -5190.8429 79001.3597 [224,] 79223.3597 -5190.8429 [225,] -10956.0747 79223.3597 [226,] 44093.4323 -10956.0747 [227,] 83290.4756 44093.4323 [228,] 71723.2438 83290.4756 [229,] 83290.4756 71723.2438 [230,] 36768.6641 83290.4756 [231,] 50918.4323 36768.6641 [232,] 83290.4756 50918.4323 [233,] 15719.9686 83290.4756 [234,] 83290.4756 15719.9686 [235,] 83397.4756 83290.4756 [236,] 83290.4756 83397.4756 [237,] 58124.7801 83290.4756 [238,] 83290.4756 58124.7801 [239,] 83290.4756 83290.4756 [240,] 83290.4756 83290.4756 [241,] 68008.1278 83290.4756 [242,] 83290.4756 68008.1278 [243,] 53705.6641 83290.4756 [244,] 83769.2438 53705.6641 [245,] 33186.3418 83769.2438 [246,] 33186.3418 33186.3418 [247,] 12952.7368 33186.3418 [248,] 33186.3418 12952.7368 [249,] 33186.3418 33186.3418 [250,] 26244.4323 33186.3418 [251,] 33186.3418 26244.4323 [252,] 70969.3597 33186.3418 [253,] -89711.9295 70969.3597 [254,] -221499.1320 -89711.9295 [255,] 83290.4756 -221499.1320 [256,] 47628.6641 83290.4756 [257,] -44464.5384 47628.6641 [258,] 33186.3418 -44464.5384 [259,] 83290.4756 33186.3418 [260,] 33186.3418 83290.4756 [261,] 33186.3418 33186.3418 [262,] 56418.9686 33186.3418 [263,] 57234.3164 56418.9686 [264,] 33186.3418 57234.3164 [265,] 83290.4756 33186.3418 [266,] 33186.3418 83290.4756 [267,] 110722.4323 33186.3418 [268,] 33186.3418 110722.4323 [269,] -11510.5384 33186.3418 [270,] 33186.3418 -11510.5384 [271,] 33186.3418 33186.3418 [272,] 33186.3418 33186.3418 [273,] -2887.7269 33186.3418 [274,] 83290.4756 -2887.7269 [275,] 36769.0845 83290.4756 [276,] -9743.1906 36769.0845 [277,] 83290.4756 -9743.1906 [278,] -113983.0200 83290.4756 [279,] 44475.5482 -113983.0200 [280,] 33186.3418 44475.5482 [281,] 83290.4756 33186.3418 [282,] -29439.5384 83290.4756 [283,] 83290.4756 -29439.5384 [284,] 60865.4323 83290.4756 [285,] 83290.4756 60865.4323 [286,] 15700.7622 83290.4756 [287,] 73219.1278 15700.7622 [288,] 33186.3418 73219.1278 [289,] 5243.9686 33186.3418 [290,] 33186.3418 5243.9686 [291,] 41929.0845 33186.3418 [292,] 52870.2731 41929.0845 [293,] 83290.4756 52870.2731 [294,] 27688.3164 83290.4756 [295,] 20406.2985 27688.3164 [296,] 83290.4756 20406.2985 [297,] 33186.3418 83290.4756 [298,] 83290.4756 33186.3418 [299,] 44922.4323 83290.4756 [300,] 13917.8781 44922.4323 [301,] 64811.3597 13917.8781 [302,] 79027.3597 64811.3597 [303,] 33186.3418 79027.3597 [304,] 52651.8527 33186.3418 [305,] 32147.0845 52651.8527 [306,] 33186.3418 32147.0845 [307,] -46545.8607 33186.3418 [308,] 5273.5049 -46545.8607 [309,] 33186.3418 5273.5049 [310,] 33186.3418 33186.3418 [311,] -21898.2339 33186.3418 [312,] 33186.3418 -21898.2339 [313,] -19846.2811 33186.3418 [314,] 49461.5482 -19846.2811 [315,] 90303.0119 49461.5482 [316,] 82874.4756 90303.0119 [317,] -20210.8429 82874.4756 [318,] 61966.8960 -20210.8429 [319,] 83290.4756 61966.8960 [320,] 83290.4756 83290.4756 [321,] 83290.4756 83290.4756 [322,] 83290.4756 83290.4756 [323,] 33186.3418 83290.4756 [324,] 83290.4756 33186.3418 [325,] 76149.3597 83290.4756 [326,] 33186.3418 76149.3597 [327,] 33186.3418 33186.3418 [328,] 83290.4756 33186.3418 [329,] 83290.4756 83290.4756 [330,] 83290.4756 83290.4756 [331,] 83290.4756 83290.4756 [332,] 83290.4756 83290.4756 [333,] 64587.0119 83290.4756 [334,] 79035.3597 64587.0119 [335,] 66229.0119 79035.3597 [336,] -21212.8429 66229.0119 [337,] 83290.4756 -21212.8429 [338,] 77527.3597 83290.4756 [339,] 42977.3164 77527.3597 [340,] 12333.0845 42977.3164 [341,] 83290.4756 12333.0845 [342,] 83290.4756 83290.4756 [343,] 83290.4756 83290.4756 [344,] 20745.9686 83290.4756 [345,] 83290.4756 20745.9686 [346,] 58798.1278 83290.4756 [347,] 83290.4756 58798.1278 [348,] 83290.4756 83290.4756 [349,] 44648.3164 83290.4756 [350,] 73811.1278 44648.3164 [351,] 52117.6641 73811.1278 [352,] 50978.3164 52117.6641 [353,] 77857.3597 50978.3164 [354,] 83290.4756 77857.3597 [355,] 83290.4756 83290.4756 [356,] 26341.7368 83290.4756 [357,] 83290.4756 26341.7368 [358,] 83290.4756 83290.4756 [359,] 83290.4756 83290.4756 [360,] 66256.0119 83290.4756 [361,] 83290.4756 66256.0119 [362,] 83290.4756 83290.4756 [363,] -93224.9040 83290.4756 [364,] 45928.4323 -93224.9040 [365,] 83290.4756 45928.4323 [366,] 83290.4756 83290.4756 [367,] 83290.4756 83290.4756 [368,] 83290.4756 83290.4756 [369,] 53621.6641 83290.4756 [370,] 83290.4756 53621.6641 [371,] 61429.6641 83290.4756 [372,] 74707.6641 61429.6641 [373,] 62545.1278 74707.6641 [374,] 83290.4756 62545.1278 [375,] 83290.4756 83290.4756 [376,] 16386.2004 83290.4756 [377,] 60408.6641 16386.2004 [378,] 45143.3164 60408.6641 [379,] 83290.4756 45143.3164 [380,] 83290.4756 83290.4756 [381,] 8465.3890 83290.4756 [382,] -11375.1473 8465.3890 [383,] 83290.4756 -11375.1473 [384,] 34083.0845 83290.4756 [385,] -117594.8862 34083.0845 [386,] 82192.1278 -117594.8862 [387,] 51369.7801 82192.1278 [388,] 69331.1278 51369.7801 [389,] 41735.5482 69331.1278 [390,] 47606.3164 41735.5482 [391,] 34467.7801 47606.3164 [392,] -2918.4518 34467.7801 [393,] 109986.7801 -2918.4518 [394,] -36651.7015 109986.7801 [395,] 17392.5482 -36651.7015 [396,] -149384.6722 17392.5482 [397,] -862585.3637 -149384.6722 [398,] 34329.6641 -862585.3637 [399,] -226360.4110 34329.6641 [400,] 59373.1278 -226360.4110 [401,] 67310.5482 59373.1278 [402,] 39250.7801 67310.5482 [403,] -92729.0926 39250.7801 [404,] 14966.6641 -92729.0926 [405,] 48023.2004 14966.6641 [406,] 38062.2004 48023.2004 [407,] 21308.0845 38062.2004 [408,] 27638.4323 21308.0845 [409,] -13017.8722 27638.4323 [410,] -165665.1320 -13017.8722 [411,] -84756.0926 -165665.1320 [412,] -76387.3970 -84756.0926 [413,] 29629.6323 -76387.3970 [414,] -113579.7996 29629.6323 [415,] -36290.7562 -113579.7996 [416,] 77017.8094 -36290.7562 [417,] -76979.9767 77017.8094 [418,] -59093.3065 -76979.9767 [419,] 10188.3597 -59093.3065 [420,] 17115.3597 10188.3597 [421,] -306663.1027 17115.3597 [422,] -655350.5409 -306663.1027 [423,] -359494.8429 -655350.5409 [424,] -44638.8429 -359494.8429 [425,] -416407.1639 -44638.8429 [426,] -278860.9588 -416407.1639 [427,] -966014.8096 -278860.9588 [428,] -114972.2339 -966014.8096 [429,] -226086.4938 -114972.2339 [430,] -533414.5231 -226086.4938 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2146005.8721 2559625.7921 2 2073314.5816 2146005.8721 3 -1833510.1393 2073314.5816 4 -133727.1739 -1833510.1393 5 631218.8668 -133727.1739 6 45603.9879 631218.8668 7 506071.3904 45603.9879 8 -270524.0886 506071.3904 9 611450.3878 -270524.0886 10 397635.4196 611450.3878 11 232695.3330 397635.4196 12 -552373.2158 232695.3330 13 -814949.5497 -552373.2158 14 -546748.1166 -814949.5497 15 570626.1725 -546748.1166 16 -28097.3485 570626.1725 17 -554963.3777 -28097.3485 18 -287336.7841 -554963.3777 19 -968495.6350 -287336.7841 20 263815.8794 -968495.6350 21 144146.5750 263815.8794 22 31189.5356 144146.5750 23 -1109804.1420 31189.5356 24 -287950.1460 -1109804.1420 25 114890.9839 -287950.1460 26 -630335.0160 114890.9839 27 -333838.9868 -630335.0160 28 409654.5342 -333838.9868 29 -384864.6377 409654.5342 30 295744.3278 -384864.6377 31 -6272.7728 295744.3278 32 288690.7049 -6272.7728 33 136548.4730 288690.7049 34 360943.0845 136548.4730 35 -753672.9841 360943.0845 36 92712.8973 -753672.9841 37 74215.0680 92712.8973 38 -85424.2798 74215.0680 39 -583164.7103 -85424.2798 40 -795895.0567 -583164.7103 41 -125041.0161 -795895.0567 42 -109757.8415 -125041.0161 43 -227260.7256 -109757.8415 44 -10876.0454 -227260.7256 45 -11466.5563 -10876.0454 46 -120249.3384 -11466.5563 47 -486531.9307 -120249.3384 48 -60347.3817 -486531.9307 49 -212384.2772 -60347.3817 50 -135204.0912 -212384.2772 51 -343612.4237 -135204.0912 52 -204412.9027 -343612.4237 53 -23748.9613 -204412.9027 54 -534415.8734 -23748.9613 55 -154377.4505 -534415.8734 56 395256.7801 -154377.4505 57 13755.0272 395256.7801 58 -108168.4365 13755.0272 59 -256859.5116 -108168.4365 60 27900.3571 -256859.5116 61 -222015.2339 27900.3571 62 25690.4437 -222015.2339 63 -246149.5843 25690.4437 64 -107119.7842 -246149.5843 65 -944371.6210 -107119.7842 66 -34820.2632 -944371.6210 67 -72626.7588 -34820.2632 68 -168353.9868 -72626.7588 69 -219660.5116 -168353.9868 70 -289031.3384 -219660.5116 71 -15562.9295 -289031.3384 72 -70530.5817 -15562.9295 73 55280.4323 -70530.5817 74 17474.5342 55280.4323 75 -120395.9767 17474.5342 76 12358.1138 -120395.9767 77 -42222.7996 12358.1138 78 -141393.5524 -42222.7996 79 2221.8348 -141393.5524 80 -8334.8862 2221.8348 81 -137279.3091 -8334.8862 82 -19922.3499 -137279.3091 83 5844.7368 -19922.3499 84 81479.7368 5844.7368 85 -214773.5409 81479.7368 86 -200058.1167 -214773.5409 87 66717.9686 -200058.1167 88 -115293.7409 66717.9686 89 94438.8527 -115293.7409 90 7781.3024 94438.8527 91 -74209.2951 7781.3024 92 -69512.4071 -74209.2951 93 -40985.0200 -69512.4071 94 -56754.5091 -40985.0200 95 -33446.1359 -56754.5091 96 158654.4323 -33446.1359 97 -45436.2085 158654.4323 98 90471.2731 -45436.2085 99 -2669.7015 90471.2731 100 -106822.2518 -2669.7015 101 -77955.7269 -106822.2518 102 -167317.9002 -77955.7269 103 47878.5482 -167317.9002 104 16729.2298 47878.5482 105 -86058.8862 16729.2298 106 -361433.2225 -86058.8862 107 160430.4870 -361433.2225 108 -110559.9728 160430.4870 109 -117557.9767 -110559.9728 110 83002.9253 -117557.9767 111 -264384.4836 83002.9253 112 -127512.6250 -264384.4836 113 -48101.1652 -127512.6250 114 70276.6463 -48101.1652 115 -37832.5817 70276.6463 116 66548.5482 -37832.5817 117 138319.5049 66548.5482 118 -114872.9333 138319.5049 119 -5926.9588 -114872.9333 120 -249101.1893 -5926.9588 121 51033.0845 -249101.1893 122 -353694.5078 51033.0845 123 91610.4437 -353694.5078 124 -1608.6289 91610.4437 125 -163052.3244 -1608.6289 126 83290.4756 -163052.3244 127 39812.3164 83290.4756 128 -12085.1906 39812.3164 129 83290.4756 -12085.1906 130 52421.6641 83290.4756 131 30780.5775 52421.6641 132 83290.4756 30780.5775 133 -4505.8429 83290.4756 134 62933.0119 -4505.8429 135 65867.0119 62933.0119 136 44182.3164 65867.0119 137 78534.3597 44182.3164 138 65771.0119 78534.3597 139 83290.4756 65771.0119 140 59311.5482 83290.4756 141 83290.4756 59311.5482 142 19584.2004 83290.4756 143 66692.0119 19584.2004 144 33747.7368 66692.0119 145 42420.4323 33747.7368 146 83290.4756 42420.4323 147 53055.6641 83290.4756 148 84633.2438 53055.6641 149 83787.4756 84633.2438 150 29718.6641 83787.4756 151 -122233.8568 29718.6641 152 61974.8960 -122233.8568 153 35833.6641 61974.8960 154 74629.2438 35833.6641 155 83290.4756 74629.2438 156 83290.4756 83290.4756 157 74620.2438 83290.4756 158 83019.8960 74620.2438 159 83290.4756 83019.8960 160 83290.4756 83290.4756 161 83290.4756 83290.4756 162 83290.4756 83290.4756 163 44336.6641 83290.4756 164 -32217.3065 44336.6641 165 78529.3597 -32217.3065 166 83290.4756 78529.3597 167 73886.3890 83290.4756 168 -124990.9002 73886.3890 169 83290.4756 -124990.9002 170 83290.4756 83290.4756 171 67533.1278 83290.4756 172 83290.4756 67533.1278 173 39138.4323 83290.4756 174 83290.4756 39138.4323 175 49352.5482 83290.4756 176 -400554.6237 49352.5482 177 33982.3418 -400554.6237 178 83290.4756 33982.3418 179 83290.4756 83290.4756 180 46246.2731 83290.4756 181 980.0667 46246.2731 182 -24378.4225 980.0667 183 -52990.5130 -24378.4225 184 11296.7622 -52990.5130 185 4122.0667 11296.7622 186 -28137.3244 4122.0667 187 55663.7801 -28137.3244 188 38381.8960 55663.7801 189 83290.4756 38381.8960 190 56813.6641 83290.4756 191 83290.4756 56813.6641 192 53434.6641 83290.4756 193 83290.4756 53434.6641 194 70809.1278 83290.4756 195 33186.3418 70809.1278 196 -251968.9753 33186.3418 197 33186.3418 -251968.9753 198 33186.3418 33186.3418 199 -93605.8607 33186.3418 200 33186.3418 -93605.8607 201 33186.3418 33186.3418 202 33186.3418 33186.3418 203 -32389.0493 33186.3418 204 -7872.0493 -32389.0493 205 12645.7622 -7872.0493 206 27854.1100 12645.7622 207 33186.3418 27854.1100 208 -49421.6722 33186.3418 209 11964.9940 -49421.6722 210 83290.4756 11964.9940 211 83290.4756 83290.4756 212 51075.2004 83290.4756 213 14264.3164 51075.2004 214 51079.8960 14264.3164 215 74764.2438 51079.8960 216 58019.7801 74764.2438 217 55804.6641 58019.7801 218 74764.2438 55804.6641 219 -29125.7702 74764.2438 220 51777.1278 -29125.7702 221 83290.4756 51777.1278 222 79001.3597 83290.4756 223 -5190.8429 79001.3597 224 79223.3597 -5190.8429 225 -10956.0747 79223.3597 226 44093.4323 -10956.0747 227 83290.4756 44093.4323 228 71723.2438 83290.4756 229 83290.4756 71723.2438 230 36768.6641 83290.4756 231 50918.4323 36768.6641 232 83290.4756 50918.4323 233 15719.9686 83290.4756 234 83290.4756 15719.9686 235 83397.4756 83290.4756 236 83290.4756 83397.4756 237 58124.7801 83290.4756 238 83290.4756 58124.7801 239 83290.4756 83290.4756 240 83290.4756 83290.4756 241 68008.1278 83290.4756 242 83290.4756 68008.1278 243 53705.6641 83290.4756 244 83769.2438 53705.6641 245 33186.3418 83769.2438 246 33186.3418 33186.3418 247 12952.7368 33186.3418 248 33186.3418 12952.7368 249 33186.3418 33186.3418 250 26244.4323 33186.3418 251 33186.3418 26244.4323 252 70969.3597 33186.3418 253 -89711.9295 70969.3597 254 -221499.1320 -89711.9295 255 83290.4756 -221499.1320 256 47628.6641 83290.4756 257 -44464.5384 47628.6641 258 33186.3418 -44464.5384 259 83290.4756 33186.3418 260 33186.3418 83290.4756 261 33186.3418 33186.3418 262 56418.9686 33186.3418 263 57234.3164 56418.9686 264 33186.3418 57234.3164 265 83290.4756 33186.3418 266 33186.3418 83290.4756 267 110722.4323 33186.3418 268 33186.3418 110722.4323 269 -11510.5384 33186.3418 270 33186.3418 -11510.5384 271 33186.3418 33186.3418 272 33186.3418 33186.3418 273 -2887.7269 33186.3418 274 83290.4756 -2887.7269 275 36769.0845 83290.4756 276 -9743.1906 36769.0845 277 83290.4756 -9743.1906 278 -113983.0200 83290.4756 279 44475.5482 -113983.0200 280 33186.3418 44475.5482 281 83290.4756 33186.3418 282 -29439.5384 83290.4756 283 83290.4756 -29439.5384 284 60865.4323 83290.4756 285 83290.4756 60865.4323 286 15700.7622 83290.4756 287 73219.1278 15700.7622 288 33186.3418 73219.1278 289 5243.9686 33186.3418 290 33186.3418 5243.9686 291 41929.0845 33186.3418 292 52870.2731 41929.0845 293 83290.4756 52870.2731 294 27688.3164 83290.4756 295 20406.2985 27688.3164 296 83290.4756 20406.2985 297 33186.3418 83290.4756 298 83290.4756 33186.3418 299 44922.4323 83290.4756 300 13917.8781 44922.4323 301 64811.3597 13917.8781 302 79027.3597 64811.3597 303 33186.3418 79027.3597 304 52651.8527 33186.3418 305 32147.0845 52651.8527 306 33186.3418 32147.0845 307 -46545.8607 33186.3418 308 5273.5049 -46545.8607 309 33186.3418 5273.5049 310 33186.3418 33186.3418 311 -21898.2339 33186.3418 312 33186.3418 -21898.2339 313 -19846.2811 33186.3418 314 49461.5482 -19846.2811 315 90303.0119 49461.5482 316 82874.4756 90303.0119 317 -20210.8429 82874.4756 318 61966.8960 -20210.8429 319 83290.4756 61966.8960 320 83290.4756 83290.4756 321 83290.4756 83290.4756 322 83290.4756 83290.4756 323 33186.3418 83290.4756 324 83290.4756 33186.3418 325 76149.3597 83290.4756 326 33186.3418 76149.3597 327 33186.3418 33186.3418 328 83290.4756 33186.3418 329 83290.4756 83290.4756 330 83290.4756 83290.4756 331 83290.4756 83290.4756 332 83290.4756 83290.4756 333 64587.0119 83290.4756 334 79035.3597 64587.0119 335 66229.0119 79035.3597 336 -21212.8429 66229.0119 337 83290.4756 -21212.8429 338 77527.3597 83290.4756 339 42977.3164 77527.3597 340 12333.0845 42977.3164 341 83290.4756 12333.0845 342 83290.4756 83290.4756 343 83290.4756 83290.4756 344 20745.9686 83290.4756 345 83290.4756 20745.9686 346 58798.1278 83290.4756 347 83290.4756 58798.1278 348 83290.4756 83290.4756 349 44648.3164 83290.4756 350 73811.1278 44648.3164 351 52117.6641 73811.1278 352 50978.3164 52117.6641 353 77857.3597 50978.3164 354 83290.4756 77857.3597 355 83290.4756 83290.4756 356 26341.7368 83290.4756 357 83290.4756 26341.7368 358 83290.4756 83290.4756 359 83290.4756 83290.4756 360 66256.0119 83290.4756 361 83290.4756 66256.0119 362 83290.4756 83290.4756 363 -93224.9040 83290.4756 364 45928.4323 -93224.9040 365 83290.4756 45928.4323 366 83290.4756 83290.4756 367 83290.4756 83290.4756 368 83290.4756 83290.4756 369 53621.6641 83290.4756 370 83290.4756 53621.6641 371 61429.6641 83290.4756 372 74707.6641 61429.6641 373 62545.1278 74707.6641 374 83290.4756 62545.1278 375 83290.4756 83290.4756 376 16386.2004 83290.4756 377 60408.6641 16386.2004 378 45143.3164 60408.6641 379 83290.4756 45143.3164 380 83290.4756 83290.4756 381 8465.3890 83290.4756 382 -11375.1473 8465.3890 383 83290.4756 -11375.1473 384 34083.0845 83290.4756 385 -117594.8862 34083.0845 386 82192.1278 -117594.8862 387 51369.7801 82192.1278 388 69331.1278 51369.7801 389 41735.5482 69331.1278 390 47606.3164 41735.5482 391 34467.7801 47606.3164 392 -2918.4518 34467.7801 393 109986.7801 -2918.4518 394 -36651.7015 109986.7801 395 17392.5482 -36651.7015 396 -149384.6722 17392.5482 397 -862585.3637 -149384.6722 398 34329.6641 -862585.3637 399 -226360.4110 34329.6641 400 59373.1278 -226360.4110 401 67310.5482 59373.1278 402 39250.7801 67310.5482 403 -92729.0926 39250.7801 404 14966.6641 -92729.0926 405 48023.2004 14966.6641 406 38062.2004 48023.2004 407 21308.0845 38062.2004 408 27638.4323 21308.0845 409 -13017.8722 27638.4323 410 -165665.1320 -13017.8722 411 -84756.0926 -165665.1320 412 -76387.3970 -84756.0926 413 29629.6323 -76387.3970 414 -113579.7996 29629.6323 415 -36290.7562 -113579.7996 416 77017.8094 -36290.7562 417 -76979.9767 77017.8094 418 -59093.3065 -76979.9767 419 10188.3597 -59093.3065 420 17115.3597 10188.3597 421 -306663.1027 17115.3597 422 -655350.5409 -306663.1027 423 -359494.8429 -655350.5409 424 -44638.8429 -359494.8429 425 -416407.1639 -44638.8429 426 -278860.9588 -416407.1639 427 -966014.8096 -278860.9588 428 -114972.2339 -966014.8096 429 -226086.4938 -114972.2339 430 -533414.5231 -226086.4938 > 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/html/rcomp/tmp/7d6gw1291232858.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/html/rcomp/tmp/8d6gw1291232858.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/html/rcomp/tmp/9d6gw1291232858.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/html/rcomp/tmp/106yfh1291232858.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/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/html/rcomp/tmp/11rgd51291232858.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/html/rcomp/tmp/12chub1291232858.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/html/rcomp/tmp/138qa11291232858.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/html/rcomp/tmp/14c9871291232858.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/html/rcomp/tmp/15fapv1291232858.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/html/rcomp/tmp/161a5j1291232858.tab") + } > > try(system("convert tmp/1hxin1291232858.ps tmp/1hxin1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/2hxin1291232858.ps tmp/2hxin1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/3r6hq1291232858.ps tmp/3r6hq1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/4r6hq1291232858.ps tmp/4r6hq1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/5r6hq1291232858.ps tmp/5r6hq1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/62fzt1291232858.ps tmp/62fzt1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/7d6gw1291232858.ps tmp/7d6gw1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/8d6gw1291232858.ps tmp/8d6gw1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/9d6gw1291232858.ps tmp/9d6gw1291232858.png",intern=TRUE)) character(0) > try(system("convert tmp/106yfh1291232858.ps tmp/106yfh1291232858.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.085 2.220 23.406