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(162556 + ,807 + ,6282154 + ,29790 + ,444 + ,4321023 + ,87550 + ,412 + ,4111912 + ,84738 + ,428 + ,223193 + ,54660 + ,315 + ,1491348 + ,42634 + ,168 + ,1629616 + ,40949 + ,263 + ,1398893 + ,45187 + ,267 + ,1926517 + ,37704 + ,228 + ,983660 + ,16275 + ,129 + ,1443586 + ,25830 + ,104 + ,1073089 + ,12679 + ,122 + ,984885 + ,18014 + ,393 + ,1405225 + ,43556 + ,190 + ,227132 + ,24811 + ,280 + ,929118 + ,6575 + ,63 + ,1071292 + ,7123 + ,102 + ,638830 + ,21950 + ,265 + ,856956 + ,37597 + ,234 + ,992426 + ,17821 + ,277 + ,444477 + ,12988 + ,73 + ,857217 + ,22330 + ,67 + ,711969 + ,13326 + ,103 + ,702380 + ,16189 + ,290 + ,358589 + ,7146 + ,83 + ,297978 + ,15824 + ,56 + ,585715 + ,27664 + ,236 + ,657954 + ,11920 + ,73 + ,209458 + ,8568 + ,34 + ,786690 + ,14416 + ,139 + ,439798 + ,3369 + ,26 + ,688779 + ,11819 + ,70 + ,574339 + ,6984 + ,40 + ,741409 + ,4519 + ,42 + ,597793 + ,2220 + ,12 + ,644190 + ,18562 + ,211 + ,377934 + ,10327 + ,74 + ,640273 + ,5336 + ,80 + ,697458 + ,2365 + ,83 + ,550608 + ,4069 + ,131 + ,207393 + ,8636 + ,203 + ,301607 + ,13718 + ,56 + ,345783 + ,4525 + ,89 + ,501749 + ,6869 + ,88 + ,379983 + ,4628 + ,39 + ,387475 + ,3689 + ,25 + ,377305 + ,4891 + ,49 + ,370837 + ,7489 + ,149 + ,430866 + ,4901 + ,58 + ,469107 + ,2284 + ,41 + ,194493 + ,3160 + ,90 + ,530670 + ,4150 + ,136 + ,518365 + ,7285 + ,97 + ,491303 + ,1134 + ,63 + ,527021 + ,4658 + ,114 + ,233773 + ,2384 + ,77 + ,405972 + ,3748 + ,6 + ,652925 + ,5371 + ,47 + ,446211 + ,1285 + ,51 + ,341340 + ,9327 + ,85 + ,387699 + ,5565 + ,43 + ,493408 + ,1528 + ,32 + ,146494 + ,3122 + ,25 + ,414462 + ,7561 + ,77 + ,364304 + ,2675 + ,54 + ,355178 + ,13253 + ,251 + ,357760 + ,880 + ,15 + ,261216 + ,2053 + ,44 + ,397144 + ,1424 + ,73 + ,374943 + ,4036 + ,85 + ,424898 + ,3045 + ,49 + ,202055 + ,5119 + ,38 + ,378525 + ,1431 + ,35 + ,310768 + ,554 + ,9 + ,325738 + ,1975 + ,34 + ,394510 + ,1765 + ,20 + ,247060 + ,1012 + ,29 + ,368078 + ,810 + ,11 + ,236761 + ,1280 + ,52 + ,312378 + ,666 + ,13 + ,339836 + ,1380 + ,29 + ,347385 + ,4677 + ,66 + ,426280 + ,876 + ,33 + ,352850 + ,814 + ,15 + ,301881 + ,514 + ,15 + ,377516 + ,5692 + ,68 + ,357312 + ,3642 + ,100 + ,458343 + ,540 + ,13 + ,354228 + ,2099 + ,45 + ,308636 + ,567 + ,14 + ,386212 + ,2001 + ,36 + ,393343 + ,2949 + ,40 + ,378509 + ,2253 + ,68 + ,452469 + ,6533 + ,29 + ,364839 + ,1889 + ,43 + ,358649 + ,3055 + ,30 + ,376641 + ,272 + ,9 + ,429112 + ,1414 + ,22 + ,330546 + ,2564 + ,19 + ,403560 + ,1383 + ,9 + ,317892 + ,1261 + ,31 + ,307528 + ,975 + ,19 + ,235133 + ,3366 + ,55 + ,299243 + ,576 + ,8 + ,314073 + ,1686 + ,28 + ,368186 + ,746 + ,29 + ,269661 + ,3192 + ,48 + ,125390 + ,2045 + ,16 + ,510834 + ,5702 + ,47 + ,321896 + ,1932 + ,20 + ,249898 + ,936 + ,22 + ,408881 + ,3437 + ,33 + ,158492 + ,5131 + ,44 + ,292154 + ,2397 + ,13 + ,289513 + ,1389 + ,6 + ,378049 + ,1503 + ,35 + ,343466 + ,402 + ,8 + ,332743 + ,2239 + ,17 + ,442882 + ,2234 + ,11 + ,214215 + ,837 + ,21 + ,315688 + ,10579 + ,92 + ,375195 + ,875 + ,12 + ,334280 + ,1585 + ,112 + ,355864 + ,1659 + ,25 + ,480382 + ,2647 + ,17 + ,353058 + ,3294 + ,23 + ,217193 + ,0 + ,0 + ,315380 + ,94 + ,10 + ,314533 + ,422 + ,23 + ,318056 + ,0 + ,0 + ,315380 + ,34 + ,7 + ,314353 + ,1558 + ,25 + ,369448 + ,0 + ,0 + ,315380 + ,43 + ,20 + ,312846 + ,645 + ,4 + ,312075 + ,316 + ,4 + ,315009 + ,115 + ,10 + ,318903 + ,5 + ,1 + ,314887 + ,897 + ,4 + ,314913 + ,0 + ,0 + ,315380 + ,389 + ,8 + ,325506 + ,0 + ,0 + ,315380 + ,1002 + ,11 + ,298568 + ,36 + ,4 + ,315834 + ,460 + ,15 + ,329784 + ,309 + ,9 + ,312878 + ,0 + ,0 + ,315380 + ,9 + ,7 + ,314987 + ,271 + ,2 + ,325249 + ,14 + ,0 + ,315877 + ,520 + ,7 + ,291650 + ,1766 + ,46 + ,305959 + ,0 + ,5 + ,315380 + ,458 + ,7 + ,297765 + ,20 + ,2 + ,315245 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,98 + ,2 + ,315236 + ,405 + ,5 + ,336425 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,483 + ,7 + ,306268 + ,454 + ,24 + ,302187 + ,47 + ,1 + ,314882 + ,0 + ,0 + ,315380 + ,757 + ,18 + ,382712 + ,4655 + ,55 + ,341570 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,36 + ,3 + ,312412 + ,0 + ,0 + ,315380 + ,203 + ,9 + ,309596 + ,0 + ,0 + ,315380 + ,126 + ,8 + ,315547 + ,400 + ,113 + ,313267 + ,71 + ,0 + ,316176 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,972 + ,19 + ,359335 + ,531 + ,11 + ,330068 + ,2461 + ,25 + ,314289 + ,378 + ,16 + ,297413 + ,23 + ,5 + ,314806 + ,638 + ,11 + ,333210 + ,2300 + ,23 + ,352108 + ,149 + ,6 + ,313332 + ,226 + ,5 + ,291787 + ,0 + ,0 + ,315380 + ,275 + ,7 + ,318745 + ,0 + ,0 + ,315380 + ,141 + ,7 + ,315366 + ,0 + ,0 + ,315380 + ,28 + ,3 + ,315688 + ,0 + ,0 + ,315380 + ,4980 + ,89 + ,409642 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,472 + ,19 + ,269587 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,203 + ,12 + ,300962 + ,496 + ,12 + ,325479 + ,10 + ,5 + ,316155 + ,63 + ,2 + ,318574 + ,0 + ,0 + ,315380 + ,1136 + ,26 + ,343613 + ,265 + ,3 + ,306948 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,267 + ,11 + ,330059 + ,474 + ,10 + ,288985 + ,534 + ,5 + ,304485 + ,0 + ,2 + ,315380 + ,15 + ,6 + ,315688 + ,397 + ,7 + ,317736 + ,0 + ,2 + ,315380 + ,1866 + ,28 + ,322331 + ,288 + ,3 + ,296656 + ,0 + ,0 + ,315380 + ,3 + ,1 + ,315354 + ,468 + ,20 + ,312161 + ,20 + ,1 + ,315576 + ,278 + ,22 + ,314922 + ,61 + ,9 + ,314551 + ,0 + ,0 + ,315380 + ,192 + ,2 + ,312339 + ,0 + ,0 + ,315380 + ,317 + ,7 + ,298700 + ,738 + ,9 + ,321376 + ,0 + ,0 + ,315380 + ,368 + ,13 + ,303230 + ,0 + ,0 + ,315380 + ,2 + ,0 + ,315487 + ,0 + ,0 + ,315380 + ,53 + ,6 + ,315793 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,94 + ,3 + ,312887 + ,0 + ,0 + ,315380 + ,24 + ,7 + ,315637 + ,2332 + ,2 + ,324385 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,131 + ,15 + ,308989 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,206 + ,9 + ,296702 + ,0 + ,0 + ,315380 + ,167 + ,1 + ,307322 + ,622 + ,38 + ,304376 + ,2328 + ,57 + ,253588 + ,0 + ,0 + ,315380 + ,365 + ,7 + ,309560 + ,364 + ,26 + ,298466 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,226 + ,13 + ,343929 + ,307 + ,10 + ,331955 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,188 + ,9 + ,381180 + ,0 + ,0 + ,315380 + ,138 + ,26 + ,331420 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,125 + ,19 + ,310201 + ,0 + ,0 + ,315380 + ,282 + ,12 + ,320016 + ,335 + ,23 + ,320398 + ,0 + ,0 + ,315380 + ,1324 + ,29 + ,291841 + ,176 + ,8 + ,310670 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,249 + ,26 + ,313491 + ,0 + ,0 + ,315380 + ,333 + ,9 + ,331323 + ,0 + ,0 + ,315380 + ,601 + ,5 + ,319210 + ,30 + ,3 + ,318098 + ,0 + ,0 + ,315380 + ,249 + ,13 + ,292754 + ,0 + ,0 + ,315380 + ,165 + ,12 + ,325176 + ,453 + ,19 + ,365959 + ,0 + ,0 + ,315380 + ,53 + ,10 + ,302409 + ,382 + ,9 + ,340968 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,9 + ,315380 + ,30 + ,4 + ,313164 + ,290 + ,1 + ,301164 + ,0 + ,1 + ,315380 + ,0 + ,0 + ,315380 + ,366 + ,14 + ,344425 + ,2 + ,12 + ,315394 + ,0 + ,0 + ,315380 + ,209 + ,19 + ,316647 + ,384 + ,17 + ,309836 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,365 + ,32 + ,346611 + ,0 + ,0 + ,315380 + ,49 + ,14 + ,322031 + ,3 + ,8 + ,315656 + ,133 + ,4 + ,339445 + ,32 + ,0 + ,314964 + ,368 + ,20 + ,297141 + ,1 + ,5 + ,315372 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,22 + ,1 + ,312502 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,96 + ,4 + ,313729 + ,1 + ,1 + ,315388 + ,314 + ,4 + ,315371 + ,844 + ,20 + ,296139 + ,0 + ,0 + ,315380 + ,26 + ,1 + ,313880 + ,125 + ,10 + ,317698 + ,304 + ,12 + ,295580 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,621 + ,13 + ,308256 + ,0 + ,0 + ,315380 + ,119 + ,3 + ,303677 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1595 + ,10 + ,319369 + ,312 + ,3 + ,318690 + ,60 + ,7 + ,314049 + ,587 + ,10 + ,325699 + ,135 + ,1 + ,314210 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,514 + ,15 + ,322378 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1 + ,4 + ,315398 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,1763 + ,28 + ,308336 + ,180 + ,9 + ,316386 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,218 + ,7 + ,315553 + ,0 + ,0 + ,315380 + ,448 + ,7 + ,323361 + ,227 + ,7 + ,336639 + ,174 + ,3 + ,307424 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,121 + ,11 + ,295370 + ,607 + ,7 + ,322340 + ,2212 + ,10 + ,319864 + ,0 + ,0 + ,315380 + ,0 + ,0 + ,315380 + ,530 + ,18 + ,317291 + ,571 + ,14 + ,280398 + ,0 + ,0 + ,315380 + ,78 + ,12 + ,317330 + ,2489 + ,29 + ,238125 + ,131 + ,3 + ,327071 + ,923 + ,6 + ,309038 + ,72 + ,3 + ,314210 + ,572 + ,8 + ,307930 + ,397 + ,10 + ,322327 + ,450 + ,6 + ,292136 + ,622 + ,8 + ,263276 + ,694 + ,6 + ,367655 + ,3425 + ,9 + ,283910 + ,562 + ,8 + ,283587 + ,4917 + ,26 + ,243650 + ,1442 + ,239 + ,438493 + ,529 + ,7 + ,296261 + ,2126 + ,41 + ,230621 + ,1061 + ,3 + ,304252 + ,776 + ,8 + ,333505 + ,611 + ,6 + ,296919 + ,1526 + ,21 + ,278990 + ,592 + ,7 + ,276898 + ,1182 + ,11 + ,327007 + ,621 + ,11 + ,317046 + ,989 + ,12 + ,304555 + ,438 + ,9 + ,298096 + ,726 + ,3 + ,231861 + ,1303 + ,57 + ,309422 + ,7419 + ,21 + ,286963 + ,1164 + ,15 + ,269753 + ,3310 + ,32 + ,448243 + ,1920 + ,11 + ,165404 + ,965 + ,2 + ,204325 + ,3256 + ,23 + ,407159 + ,1135 + ,20 + ,290476 + ,1270 + ,24 + ,275311 + ,661 + ,1 + ,246541 + ,1013 + ,1 + ,253468 + ,2844 + ,74 + ,240897 + ,11528 + ,68 + ,-83265 + ,6526 + ,20 + ,-42143 + ,2264 + ,20 + ,272713 + ,5109 + ,82 + ,215362 + ,3999 + ,21 + ,42754 + ,35624 + ,244 + ,306275 + ,9252 + ,32 + ,253537 + ,15236 + ,86 + ,372631 + ,18073 + ,69 + ,-7170) + ,dim=c(3 + ,431) + ,dimnames=list(c('kosten' + ,'orders' + ,'rijkdom') + ,1:431)) > y <- array(NA,dim=c(3,431),dimnames=list(c('kosten','orders','rijkdom'),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 rijkdom kosten orders 1 6282154 162556 807 2 4321023 29790 444 3 4111912 87550 412 4 223193 84738 428 5 1491348 54660 315 6 1629616 42634 168 7 1398893 40949 263 8 1926517 45187 267 9 983660 37704 228 10 1443586 16275 129 11 1073089 25830 104 12 984885 12679 122 13 1405225 18014 393 14 227132 43556 190 15 929118 24811 280 16 1071292 6575 63 17 638830 7123 102 18 856956 21950 265 19 992426 37597 234 20 444477 17821 277 21 857217 12988 73 22 711969 22330 67 23 702380 13326 103 24 358589 16189 290 25 297978 7146 83 26 585715 15824 56 27 657954 27664 236 28 209458 11920 73 29 786690 8568 34 30 439798 14416 139 31 688779 3369 26 32 574339 11819 70 33 741409 6984 40 34 597793 4519 42 35 644190 2220 12 36 377934 18562 211 37 640273 10327 74 38 697458 5336 80 39 550608 2365 83 40 207393 4069 131 41 301607 8636 203 42 345783 13718 56 43 501749 4525 89 44 379983 6869 88 45 387475 4628 39 46 377305 3689 25 47 370837 4891 49 48 430866 7489 149 49 469107 4901 58 50 194493 2284 41 51 530670 3160 90 52 518365 4150 136 53 491303 7285 97 54 527021 1134 63 55 233773 4658 114 56 405972 2384 77 57 652925 3748 6 58 446211 5371 47 59 341340 1285 51 60 387699 9327 85 61 493408 5565 43 62 146494 1528 32 63 414462 3122 25 64 364304 7561 77 65 355178 2675 54 66 357760 13253 251 67 261216 880 15 68 397144 2053 44 69 374943 1424 73 70 424898 4036 85 71 202055 3045 49 72 378525 5119 38 73 310768 1431 35 74 325738 554 9 75 394510 1975 34 76 247060 1765 20 77 368078 1012 29 78 236761 810 11 79 312378 1280 52 80 339836 666 13 81 347385 1380 29 82 426280 4677 66 83 352850 876 33 84 301881 814 15 85 377516 514 15 86 357312 5692 68 87 458343 3642 100 88 354228 540 13 89 308636 2099 45 90 386212 567 14 91 393343 2001 36 92 378509 2949 40 93 452469 2253 68 94 364839 6533 29 95 358649 1889 43 96 376641 3055 30 97 429112 272 9 98 330546 1414 22 99 403560 2564 19 100 317892 1383 9 101 307528 1261 31 102 235133 975 19 103 299243 3366 55 104 314073 576 8 105 368186 1686 28 106 269661 746 29 107 125390 3192 48 108 510834 2045 16 109 321896 5702 47 110 249898 1932 20 111 408881 936 22 112 158492 3437 33 113 292154 5131 44 114 289513 2397 13 115 378049 1389 6 116 343466 1503 35 117 332743 402 8 118 442882 2239 17 119 214215 2234 11 120 315688 837 21 121 375195 10579 92 122 334280 875 12 123 355864 1585 112 124 480382 1659 25 125 353058 2647 17 126 217193 3294 23 127 315380 0 0 128 314533 94 10 129 318056 422 23 130 315380 0 0 131 314353 34 7 132 369448 1558 25 133 315380 0 0 134 312846 43 20 135 312075 645 4 136 315009 316 4 137 318903 115 10 138 314887 5 1 139 314913 897 4 140 315380 0 0 141 325506 389 8 142 315380 0 0 143 298568 1002 11 144 315834 36 4 145 329784 460 15 146 312878 309 9 147 315380 0 0 148 314987 9 7 149 325249 271 2 150 315877 14 0 151 291650 520 7 152 305959 1766 46 153 315380 0 5 154 297765 458 7 155 315245 20 2 156 315380 0 0 157 315380 0 0 158 315236 98 2 159 336425 405 5 160 315380 0 0 161 315380 0 0 162 315380 0 0 163 315380 0 0 164 306268 483 7 165 302187 454 24 166 314882 47 1 167 315380 0 0 168 382712 757 18 169 341570 4655 55 170 315380 0 0 171 315380 0 0 172 312412 36 3 173 315380 0 0 174 309596 203 9 175 315380 0 0 176 315547 126 8 177 313267 400 113 178 316176 71 0 179 315380 0 0 180 315380 0 0 181 359335 972 19 182 330068 531 11 183 314289 2461 25 184 297413 378 16 185 314806 23 5 186 333210 638 11 187 352108 2300 23 188 313332 149 6 189 291787 226 5 190 315380 0 0 191 318745 275 7 192 315380 0 0 193 315366 141 7 194 315380 0 0 195 315688 28 3 196 315380 0 0 197 409642 4980 89 198 315380 0 0 199 315380 0 0 200 269587 472 19 201 315380 0 0 202 315380 0 0 203 315380 0 0 204 300962 203 12 205 325479 496 12 206 316155 10 5 207 318574 63 2 208 315380 0 0 209 343613 1136 26 210 306948 265 3 211 315380 0 0 212 315380 0 0 213 330059 267 11 214 288985 474 10 215 304485 534 5 216 315380 0 2 217 315688 15 6 218 317736 397 7 219 315380 0 2 220 322331 1866 28 221 296656 288 3 222 315380 0 0 223 315354 3 1 224 312161 468 20 225 315576 20 1 226 314922 278 22 227 314551 61 9 228 315380 0 0 229 312339 192 2 230 315380 0 0 231 298700 317 7 232 321376 738 9 233 315380 0 0 234 303230 368 13 235 315380 0 0 236 315487 2 0 237 315380 0 0 238 315793 53 6 239 315380 0 0 240 315380 0 0 241 315380 0 0 242 312887 94 3 243 315380 0 0 244 315637 24 7 245 324385 2332 2 246 315380 0 0 247 315380 0 0 248 308989 131 15 249 315380 0 0 250 315380 0 0 251 296702 206 9 252 315380 0 0 253 307322 167 1 254 304376 622 38 255 253588 2328 57 256 315380 0 0 257 309560 365 7 258 298466 364 26 259 315380 0 0 260 315380 0 0 261 315380 0 0 262 315380 0 0 263 343929 226 13 264 331955 307 10 265 315380 0 0 266 315380 0 0 267 315380 0 0 268 381180 188 9 269 315380 0 0 270 331420 138 26 271 315380 0 0 272 315380 0 0 273 315380 0 0 274 310201 125 19 275 315380 0 0 276 320016 282 12 277 320398 335 23 278 315380 0 0 279 291841 1324 29 280 310670 176 8 281 315380 0 0 282 315380 0 0 283 313491 249 26 284 315380 0 0 285 331323 333 9 286 315380 0 0 287 319210 601 5 288 318098 30 3 289 315380 0 0 290 292754 249 13 291 315380 0 0 292 325176 165 12 293 365959 453 19 294 315380 0 0 295 302409 53 10 296 340968 382 9 297 315380 0 0 298 315380 0 0 299 315380 0 0 300 315380 0 9 301 313164 30 4 302 301164 290 1 303 315380 0 1 304 315380 0 0 305 344425 366 14 306 315394 2 12 307 315380 0 0 308 316647 209 19 309 309836 384 17 310 315380 0 0 311 315380 0 0 312 346611 365 32 313 315380 0 0 314 322031 49 14 315 315656 3 8 316 339445 133 4 317 314964 32 0 318 297141 368 20 319 315372 1 5 320 315380 0 0 321 315380 0 0 322 315380 0 0 323 315380 0 0 324 315380 0 0 325 315380 0 0 326 312502 22 1 327 315380 0 0 328 315380 0 0 329 315380 0 0 330 315380 0 0 331 315380 0 0 332 315380 0 0 333 315380 0 0 334 313729 96 4 335 315388 1 1 336 315371 314 4 337 296139 844 20 338 315380 0 0 339 313880 26 1 340 317698 125 10 341 295580 304 12 342 315380 0 0 343 315380 0 0 344 315380 0 0 345 308256 621 13 346 315380 0 0 347 303677 119 3 348 315380 0 0 349 315380 0 0 350 319369 1595 10 351 318690 312 3 352 314049 60 7 353 325699 587 10 354 314210 135 1 355 315380 0 0 356 315380 0 0 357 322378 514 15 358 315380 0 0 359 315380 0 0 360 315380 0 0 361 315398 1 4 362 315380 0 0 363 315380 0 0 364 308336 1763 28 365 316386 180 9 366 315380 0 0 367 315380 0 0 368 315380 0 0 369 315380 0 0 370 315553 218 7 371 315380 0 0 372 323361 448 7 373 336639 227 7 374 307424 174 3 375 315380 0 0 376 315380 0 0 377 295370 121 11 378 322340 607 7 379 319864 2212 10 380 315380 0 0 381 315380 0 0 382 317291 530 18 383 280398 571 14 384 315380 0 0 385 317330 78 12 386 238125 2489 29 387 327071 131 3 388 309038 923 6 389 314210 72 3 390 307930 572 8 391 322327 397 10 392 292136 450 6 393 263276 622 8 394 367655 694 6 395 283910 3425 9 396 283587 562 8 397 243650 4917 26 398 438493 1442 239 399 296261 529 7 400 230621 2126 41 401 304252 1061 3 402 333505 776 8 403 296919 611 6 404 278990 1526 21 405 276898 592 7 406 327007 1182 11 407 317046 621 11 408 304555 989 12 409 298096 438 9 410 231861 726 3 411 309422 1303 57 412 286963 7419 21 413 269753 1164 15 414 448243 3310 32 415 165404 1920 11 416 204325 965 2 417 407159 3256 23 418 290476 1135 20 419 275311 1270 24 420 246541 661 1 421 253468 1013 1 422 240897 2844 74 423 -83265 11528 68 424 -42143 6526 20 425 272713 2264 20 426 215362 5109 82 427 42754 3999 21 428 306275 35624 244 429 253537 9252 32 430 372631 15236 86 431 -7170 18073 69 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) kosten orders 269503.75 21.83 1290.12 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2448703 -8351 34265 45876 2828255 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 2.695e+05 1.388e+04 19.423 < 2e-16 *** kosten 2.184e+01 2.188e+00 9.979 < 2e-16 *** orders 1.290e+03 3.472e+02 3.716 0.000229 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 257500 on 428 degrees of freedom Multiple R-squared: 0.6439, Adjusted R-squared: 0.6422 F-statistic: 387 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 3.031121e-83 1.515560e-83 [2,] 1.0000000 7.131567e-84 3.565783e-84 [3,] 1.0000000 2.374157e-91 1.187079e-91 [4,] 1.0000000 2.005009e-91 1.002504e-91 [5,] 1.0000000 8.553319e-103 4.276660e-103 [6,] 1.0000000 9.376522e-108 4.688261e-108 [7,] 1.0000000 2.438239e-110 1.219119e-110 [8,] 1.0000000 1.251548e-128 6.257742e-129 [9,] 1.0000000 4.637866e-136 2.318933e-136 [10,] 1.0000000 1.224253e-139 6.121264e-140 [11,] 1.0000000 2.282291e-151 1.141145e-151 [12,] 1.0000000 1.302023e-151 6.510116e-152 [13,] 1.0000000 4.153631e-155 2.076816e-155 [14,] 1.0000000 3.358898e-158 1.679449e-158 [15,] 1.0000000 8.154764e-162 4.077382e-162 [16,] 1.0000000 3.041757e-167 1.520879e-167 [17,] 1.0000000 6.439966e-169 3.219983e-169 [18,] 1.0000000 8.009691e-171 4.004845e-171 [19,] 1.0000000 1.523329e-174 7.616643e-175 [20,] 1.0000000 1.307145e-173 6.535727e-174 [21,] 1.0000000 2.520764e-174 1.260382e-174 [22,] 1.0000000 3.349988e-176 1.674994e-176 [23,] 1.0000000 8.574509e-176 4.287255e-176 [24,] 1.0000000 4.261454e-183 2.130727e-183 [25,] 1.0000000 1.961060e-182 9.805302e-183 [26,] 1.0000000 2.904256e-187 1.452128e-187 [27,] 1.0000000 1.191407e-188 5.957033e-189 [28,] 1.0000000 1.056506e-195 5.282530e-196 [29,] 1.0000000 3.178940e-198 1.589470e-198 [30,] 1.0000000 5.465698e-203 2.732849e-203 [31,] 1.0000000 1.492737e-203 7.463685e-204 [32,] 1.0000000 7.124162e-208 3.562081e-208 [33,] 1.0000000 1.312931e-213 6.564656e-214 [34,] 1.0000000 1.606901e-214 8.034503e-215 [35,] 1.0000000 9.667649e-216 4.833824e-216 [36,] 1.0000000 1.631686e-216 8.158431e-217 [37,] 1.0000000 5.859008e-216 2.929504e-216 [38,] 1.0000000 4.555031e-216 2.277516e-216 [39,] 1.0000000 4.292153e-215 2.146076e-215 [40,] 1.0000000 2.224175e-214 1.112087e-214 [41,] 1.0000000 1.162643e-213 5.813215e-214 [42,] 1.0000000 9.458867e-213 4.729434e-213 [43,] 1.0000000 5.942340e-212 2.971170e-212 [44,] 1.0000000 3.438009e-212 1.719005e-212 [45,] 1.0000000 3.270241e-212 1.635121e-212 [46,] 1.0000000 7.167938e-213 3.583969e-213 [47,] 1.0000000 7.790233e-213 3.895117e-213 [48,] 1.0000000 1.981194e-213 9.905972e-214 [49,] 1.0000000 1.464028e-214 7.320138e-215 [50,] 1.0000000 9.801019e-215 4.900510e-215 [51,] 1.0000000 7.250524e-214 3.625262e-214 [52,] 1.0000000 6.733777e-223 3.366889e-223 [53,] 1.0000000 2.130061e-223 1.065031e-223 [54,] 1.0000000 2.616428e-222 1.308214e-222 [55,] 1.0000000 6.340337e-222 3.170169e-222 [56,] 1.0000000 5.267335e-224 2.633668e-224 [57,] 1.0000000 1.847421e-225 9.237107e-226 [58,] 1.0000000 1.739899e-225 8.699497e-226 [59,] 1.0000000 9.231271e-225 4.615635e-225 [60,] 1.0000000 1.001221e-223 5.006107e-224 [61,] 1.0000000 5.188583e-224 2.594291e-224 [62,] 1.0000000 3.527322e-223 1.763661e-223 [63,] 1.0000000 1.567839e-222 7.839196e-223 [64,] 1.0000000 1.781478e-221 8.907390e-222 [65,] 1.0000000 5.434636e-221 2.717318e-221 [66,] 1.0000000 3.807592e-221 1.903796e-221 [67,] 1.0000000 1.051930e-220 5.259649e-221 [68,] 1.0000000 1.173747e-219 5.868737e-220 [69,] 1.0000000 1.157295e-218 5.786474e-219 [70,] 1.0000000 4.012713e-218 2.006356e-218 [71,] 1.0000000 2.139418e-217 1.069709e-217 [72,] 1.0000000 1.529183e-216 7.645913e-217 [73,] 1.0000000 5.646611e-216 2.823305e-216 [74,] 1.0000000 5.536921e-215 2.768461e-215 [75,] 1.0000000 4.917277e-214 2.458639e-214 [76,] 1.0000000 4.578538e-213 2.289269e-213 [77,] 1.0000000 5.020193e-213 2.510097e-213 [78,] 1.0000000 4.679761e-212 2.339881e-212 [79,] 1.0000000 4.780762e-211 2.390381e-211 [80,] 1.0000000 2.180355e-210 1.090177e-210 [81,] 1.0000000 1.426621e-209 7.133106e-210 [82,] 1.0000000 1.430879e-209 7.154395e-210 [83,] 1.0000000 1.007994e-208 5.039971e-209 [84,] 1.0000000 1.053723e-207 5.268613e-208 [85,] 1.0000000 3.484210e-207 1.742105e-207 [86,] 1.0000000 1.095539e-206 5.477695e-207 [87,] 1.0000000 4.543387e-206 2.271694e-206 [88,] 1.0000000 2.669290e-206 1.334645e-206 [89,] 1.0000000 4.322523e-206 2.161261e-206 [90,] 1.0000000 3.352014e-205 1.676007e-205 [91,] 1.0000000 1.054615e-204 5.273073e-205 [92,] 1.0000000 5.049718e-205 2.524859e-205 [93,] 1.0000000 4.789167e-204 2.394584e-204 [94,] 1.0000000 4.358908e-204 2.179454e-204 [95,] 1.0000000 4.201895e-203 2.100947e-203 [96,] 1.0000000 4.297084e-202 2.148542e-202 [97,] 1.0000000 1.243226e-201 6.216129e-202 [98,] 1.0000000 1.201700e-200 6.008500e-201 [99,] 1.0000000 1.167627e-199 5.838133e-200 [100,] 1.0000000 5.869277e-199 2.934638e-199 [101,] 1.0000000 3.405227e-198 1.702614e-198 [102,] 1.0000000 1.425923e-200 7.129617e-201 [103,] 1.0000000 4.431270e-204 2.215635e-204 [104,] 1.0000000 3.331433e-203 1.665716e-203 [105,] 1.0000000 1.738070e-202 8.690348e-203 [106,] 1.0000000 2.015645e-202 1.007823e-202 [107,] 1.0000000 1.503619e-203 7.518093e-204 [108,] 1.0000000 1.518862e-202 7.594310e-203 [109,] 1.0000000 1.591812e-201 7.959060e-202 [110,] 1.0000000 3.404485e-201 1.702242e-201 [111,] 1.0000000 3.055910e-200 1.527955e-200 [112,] 1.0000000 2.646932e-199 1.323466e-199 [113,] 1.0000000 1.391162e-200 6.955808e-201 [114,] 1.0000000 2.738460e-200 1.369230e-200 [115,] 1.0000000 2.827946e-199 1.413973e-199 [116,] 1.0000000 1.056223e-199 5.281113e-200 [117,] 1.0000000 8.913536e-199 4.456768e-199 [118,] 1.0000000 8.572013e-198 4.286006e-198 [119,] 1.0000000 3.105122e-200 1.552561e-200 [120,] 1.0000000 1.233619e-199 6.168094e-200 [121,] 1.0000000 2.829858e-199 1.414929e-199 [122,] 1.0000000 2.751382e-198 1.375691e-198 [123,] 1.0000000 2.766446e-197 1.383223e-197 [124,] 1.0000000 2.838145e-196 1.419072e-196 [125,] 1.0000000 2.729189e-195 1.364594e-195 [126,] 1.0000000 2.689797e-194 1.344898e-194 [127,] 1.0000000 1.021587e-193 5.107936e-194 [128,] 1.0000000 9.745014e-193 4.872507e-193 [129,] 1.0000000 9.557103e-192 4.778552e-192 [130,] 1.0000000 9.268271e-191 4.634135e-191 [131,] 1.0000000 8.875615e-190 4.437807e-190 [132,] 1.0000000 8.540957e-189 4.270479e-189 [133,] 1.0000000 8.042221e-188 4.021110e-188 [134,] 1.0000000 7.477459e-187 3.738729e-187 [135,] 1.0000000 6.944076e-186 3.472038e-186 [136,] 1.0000000 6.189332e-185 3.094666e-185 [137,] 1.0000000 5.704780e-184 2.852390e-184 [138,] 1.0000000 5.458611e-183 2.729306e-183 [139,] 1.0000000 5.055158e-182 2.527579e-182 [140,] 1.0000000 4.453346e-181 2.226673e-181 [141,] 1.0000000 4.176469e-180 2.088234e-180 [142,] 1.0000000 3.771655e-179 1.885828e-179 [143,] 1.0000000 3.471104e-178 1.735552e-178 [144,] 1.0000000 2.895904e-177 1.447952e-177 [145,] 1.0000000 2.578447e-176 1.289223e-176 [146,] 1.0000000 2.242139e-175 1.121070e-175 [147,] 1.0000000 1.967866e-174 9.839331e-175 [148,] 1.0000000 1.762457e-173 8.812283e-174 [149,] 1.0000000 1.563083e-172 7.815413e-173 [150,] 1.0000000 1.376205e-171 6.881023e-172 [151,] 1.0000000 1.197838e-170 5.989192e-171 [152,] 1.0000000 1.038451e-169 5.192256e-170 [153,] 1.0000000 9.025131e-169 4.512565e-169 [154,] 1.0000000 6.395145e-168 3.197573e-168 [155,] 1.0000000 5.479831e-167 2.739915e-167 [156,] 1.0000000 4.676536e-166 2.338268e-166 [157,] 1.0000000 3.974795e-165 1.987397e-165 [158,] 1.0000000 3.364585e-164 1.682292e-164 [159,] 1.0000000 2.928573e-163 1.464287e-163 [160,] 1.0000000 2.385561e-162 1.192781e-162 [161,] 1.0000000 2.003256e-161 1.001628e-161 [162,] 1.0000000 1.667338e-160 8.336690e-161 [163,] 1.0000000 4.228547e-160 2.114274e-160 [164,] 1.0000000 2.193356e-159 1.096678e-159 [165,] 1.0000000 1.811629e-158 9.058143e-159 [166,] 1.0000000 1.490160e-157 7.450801e-158 [167,] 1.0000000 1.241878e-156 6.209390e-157 [168,] 1.0000000 1.012994e-155 5.064972e-156 [169,] 1.0000000 8.456779e-155 4.228390e-155 [170,] 1.0000000 6.839812e-154 3.419906e-154 [171,] 1.0000000 5.620298e-153 2.810149e-153 [172,] 1.0000000 1.473969e-152 7.369844e-153 [173,] 1.0000000 1.174411e-151 5.872056e-152 [174,] 1.0000000 9.369974e-151 4.684987e-151 [175,] 1.0000000 7.444461e-150 3.722230e-150 [176,] 1.0000000 3.510510e-149 1.755255e-149 [177,] 1.0000000 2.580623e-148 1.290312e-148 [178,] 1.0000000 2.002368e-147 1.001184e-147 [179,] 1.0000000 1.513582e-146 7.567909e-147 [180,] 1.0000000 1.198769e-145 5.993847e-146 [181,] 1.0000000 8.291007e-145 4.145504e-145 [182,] 1.0000000 3.549216e-144 1.774608e-144 [183,] 1.0000000 2.794352e-143 1.397176e-143 [184,] 1.0000000 2.081553e-142 1.040776e-142 [185,] 1.0000000 1.590676e-141 7.953378e-142 [186,] 1.0000000 1.214343e-140 6.071716e-141 [187,] 1.0000000 9.200386e-140 4.600193e-140 [188,] 1.0000000 7.070252e-139 3.535126e-139 [189,] 1.0000000 5.310867e-138 2.655434e-138 [190,] 1.0000000 4.005299e-137 2.002650e-137 [191,] 1.0000000 2.982943e-136 1.491471e-136 [192,] 1.0000000 3.311066e-136 1.655533e-136 [193,] 1.0000000 2.465500e-135 1.232750e-135 [194,] 1.0000000 1.828035e-134 9.140175e-135 [195,] 1.0000000 8.884029e-134 4.442015e-134 [196,] 1.0000000 6.545488e-133 3.272744e-133 [197,] 1.0000000 4.801884e-132 2.400942e-132 [198,] 1.0000000 3.507646e-131 1.753823e-131 [199,] 1.0000000 2.547654e-130 1.273827e-130 [200,] 1.0000000 1.794754e-129 8.973768e-130 [201,] 1.0000000 1.309192e-128 6.545959e-129 [202,] 1.0000000 9.310178e-128 4.655089e-128 [203,] 1.0000000 6.656406e-127 3.328203e-127 [204,] 1.0000000 3.972948e-126 1.986474e-126 [205,] 1.0000000 2.885767e-125 1.442883e-125 [206,] 1.0000000 2.037376e-124 1.018688e-124 [207,] 1.0000000 1.432171e-123 7.160856e-124 [208,] 1.0000000 9.495477e-123 4.747739e-123 [209,] 1.0000000 6.313217e-122 3.156609e-122 [210,] 1.0000000 4.503926e-121 2.251963e-121 [211,] 1.0000000 3.131531e-120 1.565766e-120 [212,] 1.0000000 2.186161e-119 1.093081e-119 [213,] 1.0000000 1.497999e-118 7.489996e-119 [214,] 1.0000000 1.027789e-117 5.138946e-118 [215,] 1.0000000 6.805576e-117 3.402788e-117 [216,] 1.0000000 4.651250e-116 2.325625e-116 [217,] 1.0000000 3.131166e-115 1.565583e-115 [218,] 1.0000000 2.104950e-114 1.052475e-114 [219,] 1.0000000 1.443321e-113 7.216607e-114 [220,] 1.0000000 9.604622e-113 4.802311e-113 [221,] 1.0000000 6.517694e-112 3.258847e-112 [222,] 1.0000000 4.380233e-111 2.190117e-111 [223,] 1.0000000 2.871270e-110 1.435635e-110 [224,] 1.0000000 1.895349e-109 9.476744e-110 [225,] 1.0000000 1.231377e-108 6.156884e-109 [226,] 1.0000000 8.052550e-108 4.026275e-108 [227,] 1.0000000 5.036219e-107 2.518110e-107 [228,] 1.0000000 3.228984e-106 1.614492e-106 [229,] 1.0000000 2.092817e-105 1.046409e-105 [230,] 1.0000000 1.330031e-104 6.650154e-105 [231,] 1.0000000 8.410791e-104 4.205396e-104 [232,] 1.0000000 5.297394e-103 2.648697e-103 [233,] 1.0000000 3.363310e-102 1.681655e-102 [234,] 1.0000000 2.099244e-101 1.049622e-101 [235,] 1.0000000 1.304346e-100 6.521729e-101 [236,] 1.0000000 8.067742e-100 4.033871e-100 [237,] 1.0000000 5.033656e-99 2.516828e-99 [238,] 1.0000000 3.085202e-98 1.542601e-98 [239,] 1.0000000 1.910990e-97 9.554949e-98 [240,] 1.0000000 8.739258e-97 4.369629e-97 [241,] 1.0000000 5.287413e-96 2.643706e-96 [242,] 1.0000000 3.184315e-95 1.592158e-95 [243,] 1.0000000 1.949304e-94 9.746520e-95 [244,] 1.0000000 1.163345e-93 5.816726e-94 [245,] 1.0000000 6.910807e-93 3.455403e-93 [246,] 1.0000000 4.082010e-92 2.041005e-92 [247,] 1.0000000 2.403511e-91 1.201756e-91 [248,] 1.0000000 1.437097e-90 7.185487e-91 [249,] 1.0000000 7.976304e-90 3.988152e-90 [250,] 1.0000000 2.028640e-89 1.014320e-89 [251,] 1.0000000 1.179327e-88 5.896637e-89 [252,] 1.0000000 6.974840e-88 3.487420e-88 [253,] 1.0000000 3.792069e-87 1.896034e-87 [254,] 1.0000000 2.175337e-86 1.087668e-86 [255,] 1.0000000 1.242010e-85 6.210051e-86 [256,] 1.0000000 7.057753e-85 3.528877e-85 [257,] 1.0000000 3.991577e-84 1.995788e-84 [258,] 1.0000000 1.940624e-83 9.703122e-84 [259,] 1.0000000 1.018520e-82 5.092598e-83 [260,] 1.0000000 5.682157e-82 2.841079e-82 [261,] 1.0000000 3.154782e-81 1.577391e-81 [262,] 1.0000000 1.743136e-80 8.715682e-81 [263,] 1.0000000 4.189782e-80 2.094891e-80 [264,] 1.0000000 2.301701e-79 1.150850e-79 [265,] 1.0000000 1.273005e-78 6.365025e-79 [266,] 1.0000000 6.925541e-78 3.462770e-78 [267,] 1.0000000 3.749287e-77 1.874643e-77 [268,] 1.0000000 2.019798e-76 1.009899e-76 [269,] 1.0000000 1.105614e-75 5.528068e-76 [270,] 1.0000000 5.899069e-75 2.949534e-75 [271,] 1.0000000 3.159306e-74 1.579653e-74 [272,] 1.0000000 1.712204e-73 8.561021e-74 [273,] 1.0000000 9.000298e-73 4.500149e-73 [274,] 1.0000000 4.548468e-72 2.274234e-72 [275,] 1.0000000 2.425400e-71 1.212700e-71 [276,] 1.0000000 1.256916e-70 6.284579e-71 [277,] 1.0000000 6.480823e-70 3.240412e-70 [278,] 1.0000000 3.391832e-69 1.695916e-69 [279,] 1.0000000 1.731629e-68 8.658144e-69 [280,] 1.0000000 8.235681e-68 4.117841e-68 [281,] 1.0000000 4.162352e-67 2.081176e-67 [282,] 1.0000000 2.041252e-66 1.020626e-66 [283,] 1.0000000 1.018799e-65 5.093996e-66 [284,] 1.0000000 5.069690e-65 2.534845e-65 [285,] 1.0000000 2.476666e-64 1.238333e-64 [286,] 1.0000000 1.220625e-63 6.103125e-64 [287,] 1.0000000 5.943787e-63 2.971893e-63 [288,] 1.0000000 1.995089e-62 9.975446e-63 [289,] 1.0000000 9.696541e-62 4.848271e-62 [290,] 1.0000000 4.791553e-61 2.395776e-61 [291,] 1.0000000 1.981954e-60 9.909772e-61 [292,] 1.0000000 9.488062e-60 4.744031e-60 [293,] 1.0000000 4.517650e-59 2.258825e-59 [294,] 1.0000000 2.139378e-58 1.069689e-58 [295,] 1.0000000 1.028776e-57 5.143880e-58 [296,] 1.0000000 4.897574e-57 2.448787e-57 [297,] 1.0000000 2.346128e-56 1.173064e-56 [298,] 1.0000000 1.090289e-55 5.451445e-56 [299,] 1.0000000 5.024048e-55 2.512024e-55 [300,] 1.0000000 1.999610e-54 9.998050e-55 [301,] 1.0000000 9.349052e-54 4.674526e-54 [302,] 1.0000000 4.238831e-53 2.119415e-53 [303,] 1.0000000 1.966225e-52 9.831123e-53 [304,] 1.0000000 9.086901e-52 4.543450e-52 [305,] 1.0000000 4.051376e-51 2.025688e-51 [306,] 1.0000000 1.795901e-50 8.979505e-51 [307,] 1.0000000 7.537517e-50 3.768758e-50 [308,] 1.0000000 3.303520e-49 1.651760e-49 [309,] 1.0000000 1.457987e-48 7.289935e-49 [310,] 1.0000000 6.431942e-48 3.215971e-48 [311,] 1.0000000 2.407674e-47 1.203837e-47 [312,] 1.0000000 1.031654e-46 5.158269e-47 [313,] 1.0000000 4.446723e-46 2.223362e-46 [314,] 1.0000000 1.907777e-45 9.538883e-46 [315,] 1.0000000 8.030875e-45 4.015437e-45 [316,] 1.0000000 3.359996e-44 1.679998e-44 [317,] 1.0000000 1.397132e-43 6.985662e-44 [318,] 1.0000000 5.773528e-43 2.886764e-43 [319,] 1.0000000 2.371000e-42 1.185500e-42 [320,] 1.0000000 9.675861e-42 4.837930e-42 [321,] 1.0000000 3.969942e-41 1.984971e-41 [322,] 1.0000000 1.599673e-40 7.998363e-41 [323,] 1.0000000 6.404527e-40 3.202263e-40 [324,] 1.0000000 2.547601e-39 1.273800e-39 [325,] 1.0000000 1.006799e-38 5.033997e-39 [326,] 1.0000000 3.952759e-38 1.976379e-38 [327,] 1.0000000 1.541634e-37 7.708169e-38 [328,] 1.0000000 5.972596e-37 2.986298e-37 [329,] 1.0000000 2.329532e-36 1.164766e-36 [330,] 1.0000000 8.930177e-36 4.465088e-36 [331,] 1.0000000 3.393832e-35 1.696916e-35 [332,] 1.0000000 1.322450e-34 6.612248e-35 [333,] 1.0000000 4.955462e-34 2.477731e-34 [334,] 1.0000000 1.857433e-33 9.287165e-34 [335,] 1.0000000 6.967922e-33 3.483961e-33 [336,] 1.0000000 2.648756e-32 1.324378e-32 [337,] 1.0000000 9.656592e-32 4.828296e-32 [338,] 1.0000000 3.495095e-31 1.747548e-31 [339,] 1.0000000 1.255792e-30 6.278962e-31 [340,] 1.0000000 4.620686e-30 2.310343e-30 [341,] 1.0000000 1.636147e-29 8.180736e-30 [342,] 1.0000000 5.945144e-29 2.972572e-29 [343,] 1.0000000 2.074084e-28 1.037042e-28 [344,] 1.0000000 7.180569e-28 3.590284e-28 [345,] 1.0000000 2.339259e-27 1.169629e-27 [346,] 1.0000000 7.855837e-27 3.927919e-27 [347,] 1.0000000 2.712840e-26 1.356420e-26 [348,] 1.0000000 8.808039e-26 4.404019e-26 [349,] 1.0000000 2.940289e-25 1.470144e-25 [350,] 1.0000000 9.709432e-25 4.854716e-25 [351,] 1.0000000 3.179831e-24 1.589915e-24 [352,] 1.0000000 1.034631e-23 5.173156e-24 [353,] 1.0000000 3.332345e-23 1.666173e-23 [354,] 1.0000000 1.064121e-22 5.320605e-23 [355,] 1.0000000 3.368691e-22 1.684345e-22 [356,] 1.0000000 1.069586e-21 5.347929e-22 [357,] 1.0000000 3.326869e-21 1.663435e-21 [358,] 1.0000000 1.025496e-20 5.127480e-21 [359,] 1.0000000 3.245582e-20 1.622791e-20 [360,] 1.0000000 9.994554e-20 4.997277e-20 [361,] 1.0000000 2.998708e-19 1.499354e-19 [362,] 1.0000000 8.911874e-19 4.455937e-19 [363,] 1.0000000 2.623042e-18 1.311521e-18 [364,] 1.0000000 7.644992e-18 3.822496e-18 [365,] 1.0000000 2.236332e-17 1.118166e-17 [366,] 1.0000000 6.389299e-17 3.194649e-17 [367,] 1.0000000 1.761165e-16 8.805823e-17 [368,] 1.0000000 4.532034e-16 2.266017e-16 [369,] 1.0000000 1.294156e-15 6.470781e-16 [370,] 1.0000000 3.546293e-15 1.773146e-15 [371,] 1.0000000 9.610652e-15 4.805326e-15 [372,] 1.0000000 2.760638e-14 1.380319e-14 [373,] 1.0000000 7.128828e-14 3.564414e-14 [374,] 1.0000000 1.726107e-13 8.630536e-14 [375,] 1.0000000 4.475471e-13 2.237735e-13 [376,] 1.0000000 1.146299e-12 5.731497e-13 [377,] 1.0000000 3.001393e-12 1.500697e-12 [378,] 1.0000000 8.171533e-12 4.085766e-12 [379,] 1.0000000 2.022368e-11 1.011184e-11 [380,] 1.0000000 5.091401e-11 2.545701e-11 [381,] 1.0000000 1.254559e-10 6.272795e-11 [382,] 1.0000000 2.880587e-10 1.440293e-10 [383,] 1.0000000 6.897764e-10 3.448882e-10 [384,] 1.0000000 1.628101e-09 8.140504e-10 [385,] 1.0000000 3.871869e-09 1.935935e-09 [386,] 1.0000000 8.730066e-09 4.365033e-09 [387,] 1.0000000 2.097880e-08 1.048940e-08 [388,] 1.0000000 5.106215e-08 2.553108e-08 [389,] 1.0000000 8.174199e-08 4.087100e-08 [390,] 0.9999999 1.802571e-07 9.012855e-08 [391,] 0.9999998 4.186767e-07 2.093384e-07 [392,] 0.9999995 9.596505e-07 4.798252e-07 [393,] 0.9999992 1.654665e-06 8.273326e-07 [394,] 0.9999982 3.603466e-06 1.801733e-06 [395,] 0.9999965 7.082044e-06 3.541022e-06 [396,] 0.9999930 1.407388e-05 7.036939e-06 [397,] 0.9999872 2.563926e-05 1.281963e-05 [398,] 0.9999739 5.210188e-05 2.605094e-05 [399,] 0.9999449 1.102990e-04 5.514949e-05 [400,] 0.9998871 2.258447e-04 1.129223e-04 [401,] 0.9998032 3.935247e-04 1.967623e-04 [402,] 0.9996429 7.141098e-04 3.570549e-04 [403,] 0.9993487 1.302538e-03 6.512691e-04 [404,] 0.9988154 2.369239e-03 1.184619e-03 [405,] 0.9977313 4.537322e-03 2.268661e-03 [406,] 0.9957575 8.484940e-03 4.242470e-03 [407,] 0.9942121 1.157583e-02 5.787914e-03 [408,] 0.9897312 2.053756e-02 1.026878e-02 [409,] 0.9918348 1.633031e-02 8.165155e-03 [410,] 0.9851772 2.964564e-02 1.482282e-02 [411,] 0.9735342 5.293160e-02 2.646580e-02 [412,] 0.9773280 4.534402e-02 2.267201e-02 [413,] 0.9638779 7.224420e-02 3.612210e-02 [414,] 0.9416571 1.166859e-01 5.834294e-02 [415,] 0.9137935 1.724131e-01 8.620654e-02 [416,] 0.8869231 2.261538e-01 1.130769e-01 [417,] 0.8118077 3.763847e-01 1.881923e-01 [418,] 0.8517170 2.965659e-01 1.482830e-01 [419,] 0.8304385 3.391229e-01 1.695615e-01 [420,] 0.7214623 5.570754e-01 2.785377e-01 > postscript(file="/var/www/html/rcomp/tmp/17eoc1292875566.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/20o6x1292875566.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/30o6x1292875566.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/40o6x1292875566.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/5axnh1292875566.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 1422177.7855 2828254.5535 1399259.2669 -2448702.6109 -378022.7100 6 7 8 9 10 212475.4498 -104017.1613 325911.2340 -403242.7649 652298.8444 11 12 13 14 15 105424.8794 281145.9654 235376.9863 -1238522.6120 -243356.9332 16 17 18 19 20 576948.3241 82206.4518 -233698.3394 -399881.1519 -571503.5540 21 22 23 24 25 209946.7167 -131539.6641 9026.1537 -638528.9441 -234635.5535 26 27 28 29 30 -101546.3212 -520050.0461 -414492.9031 286243.2595 -323799.7397 31 32 33 34 35 312171.4007 -43536.2599 267807.6146 175433.7361 310732.0000 36 37 38 39 40 -569078.4006 49814.5402 208235.4649 122385.7917 -319960.9457 41 42 43 44 45 -418353.9839 -295494.5996 18623.3084 -153032.9388 -33393.8921 46 47 48 49 50 -4999.5642 -68675.5518 -194384.4411 17765.0640 -177775.7603 51 52 53 54 55 76058.4571 -17208.1270 -62407.1809 151479.5222 -284509.5798 56 57 58 59 60 -14924.3743 293844.3840 -1201.9415 -22017.1224 -195116.1026 61 62 63 64 65 46919.6024 -197656.7454 44537.6686 -169630.2312 -42399.5980 66 67 68 69 70 -524936.9858 -46853.9476 26048.6923 -19831.6724 -42390.0984 71 72 73 74 75 -197150.8329 -51774.5780 -35135.1327 32526.8323 38018.9463 76 77 78 79 80 -86784.1679 39064.2674 -64620.0625 -52160.0646 39018.8929 81 82 83 84 85 10336.1251 -30491.9024 21645.3152 -4747.8623 77437.5254 86 87 88 89 90 -124202.2780 -19693.9857 56162.0558 -64753.8158 86266.4056 91 92 93 94 95 33704.0154 -6989.6710 46043.9995 -84723.7008 -7575.3138 96 97 98 99 100 1729.0119 142058.1967 1785.5551 53560.0816 6579.9276 101 102 103 104 105 -29502.7850 -80171.6983 -114712.4397 21671.5858 25745.8450 106 107 108 109 110 -53544.7222 -275735.4076 176036.5982 -132744.2026 -87592.5503 111 112 113 114 115 90557.5062 -228631.1611 -146148.2854 -49099.8440 70476.2658 116 117 118 119 120 -4009.2257 44140.8107 102558.5655 -118258.5694 816.2494 121 122 123 124 125 -243987.8612 30189.5715 -92740.5446 142401.7259 3826.0383 126 127 128 129 130 -153906.6564 45876.2528 30075.6448 9665.3884 45876.2528 131 132 133 134 135 35076.0683 33673.0231 45876.2528 16601.0575 23327.4580 136 137 138 139 140 33445.0498 33987.1176 43983.9643 20663.1323 45876.2528 141 142 143 144 145 37187.6608 45876.2528 -7005.3106 40383.7450 30884.5952 146 147 148 149 150 25016.3156 45876.2528 36255.9339 47247.8386 46067.5680 151 152 153 154 155 1761.4402 -61450.0008 39425.6762 9230.1870 42724.3296 156 157 158 159 160 45876.2528 45876.2528 41012.2288 51627.6528 45876.2528 161 162 163 164 165 45876.2528 45876.2528 45876.2528 17187.3214 -8192.4349 166 167 168 169 170 43061.9100 45876.2528 73457.3654 -100530.2721 45876.2528 171 172 173 174 175 45876.2528 38251.8603 45876.2528 24048.7859 45876.2528 176 177 178 179 180 32971.1674 -110753.6285 45121.9943 45876.2528 45876.2528 181 182 183 184 185 44095.8056 34778.7980 -41202.6439 -986.0808 38349.4798 186 187 188 189 190 35584.4931 2711.9615 32834.2016 10898.0508 45876.2528 191 192 193 194 195 34205.9235 45876.2528 33752.7633 45876.2528 41702.5373 196 197 198 199 200 45876.2528 -83418.4463 45876.2528 45876.2528 -34734.8816 201 202 203 204 205 45876.2528 45876.2528 45876.2528 11544.4399 29663.8946 206 207 208 209 210 39982.3299 45114.4407 45876.2528 15762.1198 27787.7310 211 212 213 214 215 45876.2528 45876.2528 40534.1392 -3769.5130 16870.9861 216 217 218 219 220 43296.0221 38116.0415 30533.0992 43296.0221 -24039.3876 221 222 223 224 225 16993.5346 45876.2528 44494.6336 6636.3416 44345.4449 226 227 228 229 230 10965.6898 32104.3027 45876.2528 36062.7740 45876.2528 231 232 233 234 235 13243.8692 24147.2612 45876.2528 8919.6114 45876.2528 236 237 238 239 240 45939.5835 45876.2528 37391.3257 45876.2528 45876.2528 241 242 243 244 245 45876.2528 37460.4520 45876.2528 36578.4145 1382.6751 246 247 248 249 250 45876.2528 45876.2528 17273.1870 45876.2528 45876.2528 251 252 253 254 255 11089.2820 45876.2528 32881.7550 -27733.2665 -140283.3289 256 257 258 259 260 45876.2528 23055.8072 -12528.5493 45876.2528 45876.2528 261 262 263 264 265 45876.2528 45876.2528 52719.1282 42846.8695 45876.2528 266 267 268 269 270 45876.2528 45876.2528 95960.3053 45876.2528 25360.0761 271 272 273 274 275 45876.2528 45876.2528 45876.2528 13455.7335 45876.2528 276 277 278 279 280 28873.5045 13907.0008 45876.2528 -43985.1358 27002.4361 281 282 283 284 285 45876.2528 45876.2528 5007.4327 45876.2528 42937.2846 286 287 288 289 290 45876.2528 30133.0662 44068.8680 45876.2528 1041.9318 291 292 293 294 295 45876.2528 36588.1557 62051.9763 45876.2528 18846.8644 296 297 298 299 300 51512.3879 45876.2528 45876.2528 45876.2528 34265.2149 301 302 303 304 305 37844.7527 24038.0960 44586.1374 45876.2528 48868.1653 306 307 308 309 310 30365.1997 45876.2528 18067.6250 10015.7961 45876.2528 311 312 313 314 315 45876.2528 27853.9242 45876.2528 33395.7416 35765.8263 316 317 318 319 320 61876.7863 44761.5447 -6200.1958 39395.8415 45876.2528 321 322 323 324 325 45876.2528 45876.2528 45876.2528 45876.2528 45876.2528 326 327 328 329 330 41227.7757 45876.2528 45876.2528 45876.2528 45876.2528 331 332 333 334 335 45876.2528 45876.2528 45876.2528 36968.6674 44572.3028 336 337 338 339 340 33850.7190 -17595.4776 45876.2528 42518.4372 32563.7714 341 342 343 344 345 3957.1427 45876.2528 45876.2528 45876.2528 8421.4511 346 347 348 349 350 45876.2528 27704.5864 45876.2528 45876.2528 2137.8717 351 352 353 354 355 38503.5036 34204.3680 30477.1743 40468.4630 45876.2528 356 357 358 359 360 45876.2528 22299.5254 45876.2528 45876.2528 45876.2528 361 362 363 364 365 40711.9569 45876.2528 45876.2528 -35785.4212 31340.9823 366 367 368 369 370 45876.2528 45876.2528 45876.2528 45876.2528 32258.4971 371 372 373 374 375 45876.2528 35044.5333 53147.9855 30250.6819 45876.2528 376 377 378 379 380 45876.2528 9032.9946 30551.8278 -10839.0923 45876.2528 381 382 383 384 385 45876.2528 12992.8254 -19634.9329 45876.2528 30641.7681 386 387 388 389 390 -123138.4747 50836.5709 11640.2014 39263.8138 15615.9244 391 392 393 394 395 31253.7532 5065.9793 -30129.8069 75257.3307 -71988.3779 396 397 398 399 400 -8508.7294 -166757.5998 -170833.8386 6175.9286 -138197.8894 401 402 403 404 405 7711.3690 36736.6607 6333.6046 -50925.8076 -14562.6528 406 407 408 409 410 17503.4568 19791.6817 -2024.5758 7417.6489 -57365.0314 411 412 413 414 415 -62068.8376 -171624.2565 -44517.9813 65182.9517 -160213.4970 416 417 418 419 420 -88829.3916 36889.0594 -29612.3537 -52885.4894 -38685.5501 421 422 423 424 425 -39444.3383 -186172.9561 -692206.1533 -479941.8205 -72026.6461 426 427 428 429 430 -271484.3057 -341158.8368 -1055853.5887 -259264.3938 -340495.0209 431 -760308.8934 > postscript(file="/var/www/html/rcomp/tmp/6axnh1292875566.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 1422177.7855 NA 1 2828254.5535 1422177.7855 2 1399259.2669 2828254.5535 3 -2448702.6109 1399259.2669 4 -378022.7100 -2448702.6109 5 212475.4498 -378022.7100 6 -104017.1613 212475.4498 7 325911.2340 -104017.1613 8 -403242.7649 325911.2340 9 652298.8444 -403242.7649 10 105424.8794 652298.8444 11 281145.9654 105424.8794 12 235376.9863 281145.9654 13 -1238522.6120 235376.9863 14 -243356.9332 -1238522.6120 15 576948.3241 -243356.9332 16 82206.4518 576948.3241 17 -233698.3394 82206.4518 18 -399881.1519 -233698.3394 19 -571503.5540 -399881.1519 20 209946.7167 -571503.5540 21 -131539.6641 209946.7167 22 9026.1537 -131539.6641 23 -638528.9441 9026.1537 24 -234635.5535 -638528.9441 25 -101546.3212 -234635.5535 26 -520050.0461 -101546.3212 27 -414492.9031 -520050.0461 28 286243.2595 -414492.9031 29 -323799.7397 286243.2595 30 312171.4007 -323799.7397 31 -43536.2599 312171.4007 32 267807.6146 -43536.2599 33 175433.7361 267807.6146 34 310732.0000 175433.7361 35 -569078.4006 310732.0000 36 49814.5402 -569078.4006 37 208235.4649 49814.5402 38 122385.7917 208235.4649 39 -319960.9457 122385.7917 40 -418353.9839 -319960.9457 41 -295494.5996 -418353.9839 42 18623.3084 -295494.5996 43 -153032.9388 18623.3084 44 -33393.8921 -153032.9388 45 -4999.5642 -33393.8921 46 -68675.5518 -4999.5642 47 -194384.4411 -68675.5518 48 17765.0640 -194384.4411 49 -177775.7603 17765.0640 50 76058.4571 -177775.7603 51 -17208.1270 76058.4571 52 -62407.1809 -17208.1270 53 151479.5222 -62407.1809 54 -284509.5798 151479.5222 55 -14924.3743 -284509.5798 56 293844.3840 -14924.3743 57 -1201.9415 293844.3840 58 -22017.1224 -1201.9415 59 -195116.1026 -22017.1224 60 46919.6024 -195116.1026 61 -197656.7454 46919.6024 62 44537.6686 -197656.7454 63 -169630.2312 44537.6686 64 -42399.5980 -169630.2312 65 -524936.9858 -42399.5980 66 -46853.9476 -524936.9858 67 26048.6923 -46853.9476 68 -19831.6724 26048.6923 69 -42390.0984 -19831.6724 70 -197150.8329 -42390.0984 71 -51774.5780 -197150.8329 72 -35135.1327 -51774.5780 73 32526.8323 -35135.1327 74 38018.9463 32526.8323 75 -86784.1679 38018.9463 76 39064.2674 -86784.1679 77 -64620.0625 39064.2674 78 -52160.0646 -64620.0625 79 39018.8929 -52160.0646 80 10336.1251 39018.8929 81 -30491.9024 10336.1251 82 21645.3152 -30491.9024 83 -4747.8623 21645.3152 84 77437.5254 -4747.8623 85 -124202.2780 77437.5254 86 -19693.9857 -124202.2780 87 56162.0558 -19693.9857 88 -64753.8158 56162.0558 89 86266.4056 -64753.8158 90 33704.0154 86266.4056 91 -6989.6710 33704.0154 92 46043.9995 -6989.6710 93 -84723.7008 46043.9995 94 -7575.3138 -84723.7008 95 1729.0119 -7575.3138 96 142058.1967 1729.0119 97 1785.5551 142058.1967 98 53560.0816 1785.5551 99 6579.9276 53560.0816 100 -29502.7850 6579.9276 101 -80171.6983 -29502.7850 102 -114712.4397 -80171.6983 103 21671.5858 -114712.4397 104 25745.8450 21671.5858 105 -53544.7222 25745.8450 106 -275735.4076 -53544.7222 107 176036.5982 -275735.4076 108 -132744.2026 176036.5982 109 -87592.5503 -132744.2026 110 90557.5062 -87592.5503 111 -228631.1611 90557.5062 112 -146148.2854 -228631.1611 113 -49099.8440 -146148.2854 114 70476.2658 -49099.8440 115 -4009.2257 70476.2658 116 44140.8107 -4009.2257 117 102558.5655 44140.8107 118 -118258.5694 102558.5655 119 816.2494 -118258.5694 120 -243987.8612 816.2494 121 30189.5715 -243987.8612 122 -92740.5446 30189.5715 123 142401.7259 -92740.5446 124 3826.0383 142401.7259 125 -153906.6564 3826.0383 126 45876.2528 -153906.6564 127 30075.6448 45876.2528 128 9665.3884 30075.6448 129 45876.2528 9665.3884 130 35076.0683 45876.2528 131 33673.0231 35076.0683 132 45876.2528 33673.0231 133 16601.0575 45876.2528 134 23327.4580 16601.0575 135 33445.0498 23327.4580 136 33987.1176 33445.0498 137 43983.9643 33987.1176 138 20663.1323 43983.9643 139 45876.2528 20663.1323 140 37187.6608 45876.2528 141 45876.2528 37187.6608 142 -7005.3106 45876.2528 143 40383.7450 -7005.3106 144 30884.5952 40383.7450 145 25016.3156 30884.5952 146 45876.2528 25016.3156 147 36255.9339 45876.2528 148 47247.8386 36255.9339 149 46067.5680 47247.8386 150 1761.4402 46067.5680 151 -61450.0008 1761.4402 152 39425.6762 -61450.0008 153 9230.1870 39425.6762 154 42724.3296 9230.1870 155 45876.2528 42724.3296 156 45876.2528 45876.2528 157 41012.2288 45876.2528 158 51627.6528 41012.2288 159 45876.2528 51627.6528 160 45876.2528 45876.2528 161 45876.2528 45876.2528 162 45876.2528 45876.2528 163 17187.3214 45876.2528 164 -8192.4349 17187.3214 165 43061.9100 -8192.4349 166 45876.2528 43061.9100 167 73457.3654 45876.2528 168 -100530.2721 73457.3654 169 45876.2528 -100530.2721 170 45876.2528 45876.2528 171 38251.8603 45876.2528 172 45876.2528 38251.8603 173 24048.7859 45876.2528 174 45876.2528 24048.7859 175 32971.1674 45876.2528 176 -110753.6285 32971.1674 177 45121.9943 -110753.6285 178 45876.2528 45121.9943 179 45876.2528 45876.2528 180 44095.8056 45876.2528 181 34778.7980 44095.8056 182 -41202.6439 34778.7980 183 -986.0808 -41202.6439 184 38349.4798 -986.0808 185 35584.4931 38349.4798 186 2711.9615 35584.4931 187 32834.2016 2711.9615 188 10898.0508 32834.2016 189 45876.2528 10898.0508 190 34205.9235 45876.2528 191 45876.2528 34205.9235 192 33752.7633 45876.2528 193 45876.2528 33752.7633 194 41702.5373 45876.2528 195 45876.2528 41702.5373 196 -83418.4463 45876.2528 197 45876.2528 -83418.4463 198 45876.2528 45876.2528 199 -34734.8816 45876.2528 200 45876.2528 -34734.8816 201 45876.2528 45876.2528 202 45876.2528 45876.2528 203 11544.4399 45876.2528 204 29663.8946 11544.4399 205 39982.3299 29663.8946 206 45114.4407 39982.3299 207 45876.2528 45114.4407 208 15762.1198 45876.2528 209 27787.7310 15762.1198 210 45876.2528 27787.7310 211 45876.2528 45876.2528 212 40534.1392 45876.2528 213 -3769.5130 40534.1392 214 16870.9861 -3769.5130 215 43296.0221 16870.9861 216 38116.0415 43296.0221 217 30533.0992 38116.0415 218 43296.0221 30533.0992 219 -24039.3876 43296.0221 220 16993.5346 -24039.3876 221 45876.2528 16993.5346 222 44494.6336 45876.2528 223 6636.3416 44494.6336 224 44345.4449 6636.3416 225 10965.6898 44345.4449 226 32104.3027 10965.6898 227 45876.2528 32104.3027 228 36062.7740 45876.2528 229 45876.2528 36062.7740 230 13243.8692 45876.2528 231 24147.2612 13243.8692 232 45876.2528 24147.2612 233 8919.6114 45876.2528 234 45876.2528 8919.6114 235 45939.5835 45876.2528 236 45876.2528 45939.5835 237 37391.3257 45876.2528 238 45876.2528 37391.3257 239 45876.2528 45876.2528 240 45876.2528 45876.2528 241 37460.4520 45876.2528 242 45876.2528 37460.4520 243 36578.4145 45876.2528 244 1382.6751 36578.4145 245 45876.2528 1382.6751 246 45876.2528 45876.2528 247 17273.1870 45876.2528 248 45876.2528 17273.1870 249 45876.2528 45876.2528 250 11089.2820 45876.2528 251 45876.2528 11089.2820 252 32881.7550 45876.2528 253 -27733.2665 32881.7550 254 -140283.3289 -27733.2665 255 45876.2528 -140283.3289 256 23055.8072 45876.2528 257 -12528.5493 23055.8072 258 45876.2528 -12528.5493 259 45876.2528 45876.2528 260 45876.2528 45876.2528 261 45876.2528 45876.2528 262 52719.1282 45876.2528 263 42846.8695 52719.1282 264 45876.2528 42846.8695 265 45876.2528 45876.2528 266 45876.2528 45876.2528 267 95960.3053 45876.2528 268 45876.2528 95960.3053 269 25360.0761 45876.2528 270 45876.2528 25360.0761 271 45876.2528 45876.2528 272 45876.2528 45876.2528 273 13455.7335 45876.2528 274 45876.2528 13455.7335 275 28873.5045 45876.2528 276 13907.0008 28873.5045 277 45876.2528 13907.0008 278 -43985.1358 45876.2528 279 27002.4361 -43985.1358 280 45876.2528 27002.4361 281 45876.2528 45876.2528 282 5007.4327 45876.2528 283 45876.2528 5007.4327 284 42937.2846 45876.2528 285 45876.2528 42937.2846 286 30133.0662 45876.2528 287 44068.8680 30133.0662 288 45876.2528 44068.8680 289 1041.9318 45876.2528 290 45876.2528 1041.9318 291 36588.1557 45876.2528 292 62051.9763 36588.1557 293 45876.2528 62051.9763 294 18846.8644 45876.2528 295 51512.3879 18846.8644 296 45876.2528 51512.3879 297 45876.2528 45876.2528 298 45876.2528 45876.2528 299 34265.2149 45876.2528 300 37844.7527 34265.2149 301 24038.0960 37844.7527 302 44586.1374 24038.0960 303 45876.2528 44586.1374 304 48868.1653 45876.2528 305 30365.1997 48868.1653 306 45876.2528 30365.1997 307 18067.6250 45876.2528 308 10015.7961 18067.6250 309 45876.2528 10015.7961 310 45876.2528 45876.2528 311 27853.9242 45876.2528 312 45876.2528 27853.9242 313 33395.7416 45876.2528 314 35765.8263 33395.7416 315 61876.7863 35765.8263 316 44761.5447 61876.7863 317 -6200.1958 44761.5447 318 39395.8415 -6200.1958 319 45876.2528 39395.8415 320 45876.2528 45876.2528 321 45876.2528 45876.2528 322 45876.2528 45876.2528 323 45876.2528 45876.2528 324 45876.2528 45876.2528 325 41227.7757 45876.2528 326 45876.2528 41227.7757 327 45876.2528 45876.2528 328 45876.2528 45876.2528 329 45876.2528 45876.2528 330 45876.2528 45876.2528 331 45876.2528 45876.2528 332 45876.2528 45876.2528 333 36968.6674 45876.2528 334 44572.3028 36968.6674 335 33850.7190 44572.3028 336 -17595.4776 33850.7190 337 45876.2528 -17595.4776 338 42518.4372 45876.2528 339 32563.7714 42518.4372 340 3957.1427 32563.7714 341 45876.2528 3957.1427 342 45876.2528 45876.2528 343 45876.2528 45876.2528 344 8421.4511 45876.2528 345 45876.2528 8421.4511 346 27704.5864 45876.2528 347 45876.2528 27704.5864 348 45876.2528 45876.2528 349 2137.8717 45876.2528 350 38503.5036 2137.8717 351 34204.3680 38503.5036 352 30477.1743 34204.3680 353 40468.4630 30477.1743 354 45876.2528 40468.4630 355 45876.2528 45876.2528 356 22299.5254 45876.2528 357 45876.2528 22299.5254 358 45876.2528 45876.2528 359 45876.2528 45876.2528 360 40711.9569 45876.2528 361 45876.2528 40711.9569 362 45876.2528 45876.2528 363 -35785.4212 45876.2528 364 31340.9823 -35785.4212 365 45876.2528 31340.9823 366 45876.2528 45876.2528 367 45876.2528 45876.2528 368 45876.2528 45876.2528 369 32258.4971 45876.2528 370 45876.2528 32258.4971 371 35044.5333 45876.2528 372 53147.9855 35044.5333 373 30250.6819 53147.9855 374 45876.2528 30250.6819 375 45876.2528 45876.2528 376 9032.9946 45876.2528 377 30551.8278 9032.9946 378 -10839.0923 30551.8278 379 45876.2528 -10839.0923 380 45876.2528 45876.2528 381 12992.8254 45876.2528 382 -19634.9329 12992.8254 383 45876.2528 -19634.9329 384 30641.7681 45876.2528 385 -123138.4747 30641.7681 386 50836.5709 -123138.4747 387 11640.2014 50836.5709 388 39263.8138 11640.2014 389 15615.9244 39263.8138 390 31253.7532 15615.9244 391 5065.9793 31253.7532 392 -30129.8069 5065.9793 393 75257.3307 -30129.8069 394 -71988.3779 75257.3307 395 -8508.7294 -71988.3779 396 -166757.5998 -8508.7294 397 -170833.8386 -166757.5998 398 6175.9286 -170833.8386 399 -138197.8894 6175.9286 400 7711.3690 -138197.8894 401 36736.6607 7711.3690 402 6333.6046 36736.6607 403 -50925.8076 6333.6046 404 -14562.6528 -50925.8076 405 17503.4568 -14562.6528 406 19791.6817 17503.4568 407 -2024.5758 19791.6817 408 7417.6489 -2024.5758 409 -57365.0314 7417.6489 410 -62068.8376 -57365.0314 411 -171624.2565 -62068.8376 412 -44517.9813 -171624.2565 413 65182.9517 -44517.9813 414 -160213.4970 65182.9517 415 -88829.3916 -160213.4970 416 36889.0594 -88829.3916 417 -29612.3537 36889.0594 418 -52885.4894 -29612.3537 419 -38685.5501 -52885.4894 420 -39444.3383 -38685.5501 421 -186172.9561 -39444.3383 422 -692206.1533 -186172.9561 423 -479941.8205 -692206.1533 424 -72026.6461 -479941.8205 425 -271484.3057 -72026.6461 426 -341158.8368 -271484.3057 427 -1055853.5887 -341158.8368 428 -259264.3938 -1055853.5887 429 -340495.0209 -259264.3938 430 -760308.8934 -340495.0209 431 NA -760308.8934 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2828254.5535 1422177.7855 [2,] 1399259.2669 2828254.5535 [3,] -2448702.6109 1399259.2669 [4,] -378022.7100 -2448702.6109 [5,] 212475.4498 -378022.7100 [6,] -104017.1613 212475.4498 [7,] 325911.2340 -104017.1613 [8,] -403242.7649 325911.2340 [9,] 652298.8444 -403242.7649 [10,] 105424.8794 652298.8444 [11,] 281145.9654 105424.8794 [12,] 235376.9863 281145.9654 [13,] -1238522.6120 235376.9863 [14,] -243356.9332 -1238522.6120 [15,] 576948.3241 -243356.9332 [16,] 82206.4518 576948.3241 [17,] -233698.3394 82206.4518 [18,] -399881.1519 -233698.3394 [19,] -571503.5540 -399881.1519 [20,] 209946.7167 -571503.5540 [21,] -131539.6641 209946.7167 [22,] 9026.1537 -131539.6641 [23,] -638528.9441 9026.1537 [24,] -234635.5535 -638528.9441 [25,] -101546.3212 -234635.5535 [26,] -520050.0461 -101546.3212 [27,] -414492.9031 -520050.0461 [28,] 286243.2595 -414492.9031 [29,] -323799.7397 286243.2595 [30,] 312171.4007 -323799.7397 [31,] -43536.2599 312171.4007 [32,] 267807.6146 -43536.2599 [33,] 175433.7361 267807.6146 [34,] 310732.0000 175433.7361 [35,] -569078.4006 310732.0000 [36,] 49814.5402 -569078.4006 [37,] 208235.4649 49814.5402 [38,] 122385.7917 208235.4649 [39,] -319960.9457 122385.7917 [40,] -418353.9839 -319960.9457 [41,] -295494.5996 -418353.9839 [42,] 18623.3084 -295494.5996 [43,] -153032.9388 18623.3084 [44,] -33393.8921 -153032.9388 [45,] -4999.5642 -33393.8921 [46,] -68675.5518 -4999.5642 [47,] -194384.4411 -68675.5518 [48,] 17765.0640 -194384.4411 [49,] -177775.7603 17765.0640 [50,] 76058.4571 -177775.7603 [51,] -17208.1270 76058.4571 [52,] -62407.1809 -17208.1270 [53,] 151479.5222 -62407.1809 [54,] -284509.5798 151479.5222 [55,] -14924.3743 -284509.5798 [56,] 293844.3840 -14924.3743 [57,] -1201.9415 293844.3840 [58,] -22017.1224 -1201.9415 [59,] -195116.1026 -22017.1224 [60,] 46919.6024 -195116.1026 [61,] -197656.7454 46919.6024 [62,] 44537.6686 -197656.7454 [63,] -169630.2312 44537.6686 [64,] -42399.5980 -169630.2312 [65,] -524936.9858 -42399.5980 [66,] -46853.9476 -524936.9858 [67,] 26048.6923 -46853.9476 [68,] -19831.6724 26048.6923 [69,] -42390.0984 -19831.6724 [70,] -197150.8329 -42390.0984 [71,] -51774.5780 -197150.8329 [72,] -35135.1327 -51774.5780 [73,] 32526.8323 -35135.1327 [74,] 38018.9463 32526.8323 [75,] -86784.1679 38018.9463 [76,] 39064.2674 -86784.1679 [77,] -64620.0625 39064.2674 [78,] -52160.0646 -64620.0625 [79,] 39018.8929 -52160.0646 [80,] 10336.1251 39018.8929 [81,] -30491.9024 10336.1251 [82,] 21645.3152 -30491.9024 [83,] -4747.8623 21645.3152 [84,] 77437.5254 -4747.8623 [85,] -124202.2780 77437.5254 [86,] -19693.9857 -124202.2780 [87,] 56162.0558 -19693.9857 [88,] -64753.8158 56162.0558 [89,] 86266.4056 -64753.8158 [90,] 33704.0154 86266.4056 [91,] -6989.6710 33704.0154 [92,] 46043.9995 -6989.6710 [93,] -84723.7008 46043.9995 [94,] -7575.3138 -84723.7008 [95,] 1729.0119 -7575.3138 [96,] 142058.1967 1729.0119 [97,] 1785.5551 142058.1967 [98,] 53560.0816 1785.5551 [99,] 6579.9276 53560.0816 [100,] -29502.7850 6579.9276 [101,] -80171.6983 -29502.7850 [102,] -114712.4397 -80171.6983 [103,] 21671.5858 -114712.4397 [104,] 25745.8450 21671.5858 [105,] -53544.7222 25745.8450 [106,] -275735.4076 -53544.7222 [107,] 176036.5982 -275735.4076 [108,] -132744.2026 176036.5982 [109,] -87592.5503 -132744.2026 [110,] 90557.5062 -87592.5503 [111,] -228631.1611 90557.5062 [112,] -146148.2854 -228631.1611 [113,] -49099.8440 -146148.2854 [114,] 70476.2658 -49099.8440 [115,] -4009.2257 70476.2658 [116,] 44140.8107 -4009.2257 [117,] 102558.5655 44140.8107 [118,] -118258.5694 102558.5655 [119,] 816.2494 -118258.5694 [120,] -243987.8612 816.2494 [121,] 30189.5715 -243987.8612 [122,] -92740.5446 30189.5715 [123,] 142401.7259 -92740.5446 [124,] 3826.0383 142401.7259 [125,] -153906.6564 3826.0383 [126,] 45876.2528 -153906.6564 [127,] 30075.6448 45876.2528 [128,] 9665.3884 30075.6448 [129,] 45876.2528 9665.3884 [130,] 35076.0683 45876.2528 [131,] 33673.0231 35076.0683 [132,] 45876.2528 33673.0231 [133,] 16601.0575 45876.2528 [134,] 23327.4580 16601.0575 [135,] 33445.0498 23327.4580 [136,] 33987.1176 33445.0498 [137,] 43983.9643 33987.1176 [138,] 20663.1323 43983.9643 [139,] 45876.2528 20663.1323 [140,] 37187.6608 45876.2528 [141,] 45876.2528 37187.6608 [142,] -7005.3106 45876.2528 [143,] 40383.7450 -7005.3106 [144,] 30884.5952 40383.7450 [145,] 25016.3156 30884.5952 [146,] 45876.2528 25016.3156 [147,] 36255.9339 45876.2528 [148,] 47247.8386 36255.9339 [149,] 46067.5680 47247.8386 [150,] 1761.4402 46067.5680 [151,] -61450.0008 1761.4402 [152,] 39425.6762 -61450.0008 [153,] 9230.1870 39425.6762 [154,] 42724.3296 9230.1870 [155,] 45876.2528 42724.3296 [156,] 45876.2528 45876.2528 [157,] 41012.2288 45876.2528 [158,] 51627.6528 41012.2288 [159,] 45876.2528 51627.6528 [160,] 45876.2528 45876.2528 [161,] 45876.2528 45876.2528 [162,] 45876.2528 45876.2528 [163,] 17187.3214 45876.2528 [164,] -8192.4349 17187.3214 [165,] 43061.9100 -8192.4349 [166,] 45876.2528 43061.9100 [167,] 73457.3654 45876.2528 [168,] -100530.2721 73457.3654 [169,] 45876.2528 -100530.2721 [170,] 45876.2528 45876.2528 [171,] 38251.8603 45876.2528 [172,] 45876.2528 38251.8603 [173,] 24048.7859 45876.2528 [174,] 45876.2528 24048.7859 [175,] 32971.1674 45876.2528 [176,] -110753.6285 32971.1674 [177,] 45121.9943 -110753.6285 [178,] 45876.2528 45121.9943 [179,] 45876.2528 45876.2528 [180,] 44095.8056 45876.2528 [181,] 34778.7980 44095.8056 [182,] -41202.6439 34778.7980 [183,] -986.0808 -41202.6439 [184,] 38349.4798 -986.0808 [185,] 35584.4931 38349.4798 [186,] 2711.9615 35584.4931 [187,] 32834.2016 2711.9615 [188,] 10898.0508 32834.2016 [189,] 45876.2528 10898.0508 [190,] 34205.9235 45876.2528 [191,] 45876.2528 34205.9235 [192,] 33752.7633 45876.2528 [193,] 45876.2528 33752.7633 [194,] 41702.5373 45876.2528 [195,] 45876.2528 41702.5373 [196,] -83418.4463 45876.2528 [197,] 45876.2528 -83418.4463 [198,] 45876.2528 45876.2528 [199,] -34734.8816 45876.2528 [200,] 45876.2528 -34734.8816 [201,] 45876.2528 45876.2528 [202,] 45876.2528 45876.2528 [203,] 11544.4399 45876.2528 [204,] 29663.8946 11544.4399 [205,] 39982.3299 29663.8946 [206,] 45114.4407 39982.3299 [207,] 45876.2528 45114.4407 [208,] 15762.1198 45876.2528 [209,] 27787.7310 15762.1198 [210,] 45876.2528 27787.7310 [211,] 45876.2528 45876.2528 [212,] 40534.1392 45876.2528 [213,] -3769.5130 40534.1392 [214,] 16870.9861 -3769.5130 [215,] 43296.0221 16870.9861 [216,] 38116.0415 43296.0221 [217,] 30533.0992 38116.0415 [218,] 43296.0221 30533.0992 [219,] -24039.3876 43296.0221 [220,] 16993.5346 -24039.3876 [221,] 45876.2528 16993.5346 [222,] 44494.6336 45876.2528 [223,] 6636.3416 44494.6336 [224,] 44345.4449 6636.3416 [225,] 10965.6898 44345.4449 [226,] 32104.3027 10965.6898 [227,] 45876.2528 32104.3027 [228,] 36062.7740 45876.2528 [229,] 45876.2528 36062.7740 [230,] 13243.8692 45876.2528 [231,] 24147.2612 13243.8692 [232,] 45876.2528 24147.2612 [233,] 8919.6114 45876.2528 [234,] 45876.2528 8919.6114 [235,] 45939.5835 45876.2528 [236,] 45876.2528 45939.5835 [237,] 37391.3257 45876.2528 [238,] 45876.2528 37391.3257 [239,] 45876.2528 45876.2528 [240,] 45876.2528 45876.2528 [241,] 37460.4520 45876.2528 [242,] 45876.2528 37460.4520 [243,] 36578.4145 45876.2528 [244,] 1382.6751 36578.4145 [245,] 45876.2528 1382.6751 [246,] 45876.2528 45876.2528 [247,] 17273.1870 45876.2528 [248,] 45876.2528 17273.1870 [249,] 45876.2528 45876.2528 [250,] 11089.2820 45876.2528 [251,] 45876.2528 11089.2820 [252,] 32881.7550 45876.2528 [253,] -27733.2665 32881.7550 [254,] -140283.3289 -27733.2665 [255,] 45876.2528 -140283.3289 [256,] 23055.8072 45876.2528 [257,] -12528.5493 23055.8072 [258,] 45876.2528 -12528.5493 [259,] 45876.2528 45876.2528 [260,] 45876.2528 45876.2528 [261,] 45876.2528 45876.2528 [262,] 52719.1282 45876.2528 [263,] 42846.8695 52719.1282 [264,] 45876.2528 42846.8695 [265,] 45876.2528 45876.2528 [266,] 45876.2528 45876.2528 [267,] 95960.3053 45876.2528 [268,] 45876.2528 95960.3053 [269,] 25360.0761 45876.2528 [270,] 45876.2528 25360.0761 [271,] 45876.2528 45876.2528 [272,] 45876.2528 45876.2528 [273,] 13455.7335 45876.2528 [274,] 45876.2528 13455.7335 [275,] 28873.5045 45876.2528 [276,] 13907.0008 28873.5045 [277,] 45876.2528 13907.0008 [278,] -43985.1358 45876.2528 [279,] 27002.4361 -43985.1358 [280,] 45876.2528 27002.4361 [281,] 45876.2528 45876.2528 [282,] 5007.4327 45876.2528 [283,] 45876.2528 5007.4327 [284,] 42937.2846 45876.2528 [285,] 45876.2528 42937.2846 [286,] 30133.0662 45876.2528 [287,] 44068.8680 30133.0662 [288,] 45876.2528 44068.8680 [289,] 1041.9318 45876.2528 [290,] 45876.2528 1041.9318 [291,] 36588.1557 45876.2528 [292,] 62051.9763 36588.1557 [293,] 45876.2528 62051.9763 [294,] 18846.8644 45876.2528 [295,] 51512.3879 18846.8644 [296,] 45876.2528 51512.3879 [297,] 45876.2528 45876.2528 [298,] 45876.2528 45876.2528 [299,] 34265.2149 45876.2528 [300,] 37844.7527 34265.2149 [301,] 24038.0960 37844.7527 [302,] 44586.1374 24038.0960 [303,] 45876.2528 44586.1374 [304,] 48868.1653 45876.2528 [305,] 30365.1997 48868.1653 [306,] 45876.2528 30365.1997 [307,] 18067.6250 45876.2528 [308,] 10015.7961 18067.6250 [309,] 45876.2528 10015.7961 [310,] 45876.2528 45876.2528 [311,] 27853.9242 45876.2528 [312,] 45876.2528 27853.9242 [313,] 33395.7416 45876.2528 [314,] 35765.8263 33395.7416 [315,] 61876.7863 35765.8263 [316,] 44761.5447 61876.7863 [317,] -6200.1958 44761.5447 [318,] 39395.8415 -6200.1958 [319,] 45876.2528 39395.8415 [320,] 45876.2528 45876.2528 [321,] 45876.2528 45876.2528 [322,] 45876.2528 45876.2528 [323,] 45876.2528 45876.2528 [324,] 45876.2528 45876.2528 [325,] 41227.7757 45876.2528 [326,] 45876.2528 41227.7757 [327,] 45876.2528 45876.2528 [328,] 45876.2528 45876.2528 [329,] 45876.2528 45876.2528 [330,] 45876.2528 45876.2528 [331,] 45876.2528 45876.2528 [332,] 45876.2528 45876.2528 [333,] 36968.6674 45876.2528 [334,] 44572.3028 36968.6674 [335,] 33850.7190 44572.3028 [336,] -17595.4776 33850.7190 [337,] 45876.2528 -17595.4776 [338,] 42518.4372 45876.2528 [339,] 32563.7714 42518.4372 [340,] 3957.1427 32563.7714 [341,] 45876.2528 3957.1427 [342,] 45876.2528 45876.2528 [343,] 45876.2528 45876.2528 [344,] 8421.4511 45876.2528 [345,] 45876.2528 8421.4511 [346,] 27704.5864 45876.2528 [347,] 45876.2528 27704.5864 [348,] 45876.2528 45876.2528 [349,] 2137.8717 45876.2528 [350,] 38503.5036 2137.8717 [351,] 34204.3680 38503.5036 [352,] 30477.1743 34204.3680 [353,] 40468.4630 30477.1743 [354,] 45876.2528 40468.4630 [355,] 45876.2528 45876.2528 [356,] 22299.5254 45876.2528 [357,] 45876.2528 22299.5254 [358,] 45876.2528 45876.2528 [359,] 45876.2528 45876.2528 [360,] 40711.9569 45876.2528 [361,] 45876.2528 40711.9569 [362,] 45876.2528 45876.2528 [363,] -35785.4212 45876.2528 [364,] 31340.9823 -35785.4212 [365,] 45876.2528 31340.9823 [366,] 45876.2528 45876.2528 [367,] 45876.2528 45876.2528 [368,] 45876.2528 45876.2528 [369,] 32258.4971 45876.2528 [370,] 45876.2528 32258.4971 [371,] 35044.5333 45876.2528 [372,] 53147.9855 35044.5333 [373,] 30250.6819 53147.9855 [374,] 45876.2528 30250.6819 [375,] 45876.2528 45876.2528 [376,] 9032.9946 45876.2528 [377,] 30551.8278 9032.9946 [378,] -10839.0923 30551.8278 [379,] 45876.2528 -10839.0923 [380,] 45876.2528 45876.2528 [381,] 12992.8254 45876.2528 [382,] -19634.9329 12992.8254 [383,] 45876.2528 -19634.9329 [384,] 30641.7681 45876.2528 [385,] -123138.4747 30641.7681 [386,] 50836.5709 -123138.4747 [387,] 11640.2014 50836.5709 [388,] 39263.8138 11640.2014 [389,] 15615.9244 39263.8138 [390,] 31253.7532 15615.9244 [391,] 5065.9793 31253.7532 [392,] -30129.8069 5065.9793 [393,] 75257.3307 -30129.8069 [394,] -71988.3779 75257.3307 [395,] -8508.7294 -71988.3779 [396,] -166757.5998 -8508.7294 [397,] -170833.8386 -166757.5998 [398,] 6175.9286 -170833.8386 [399,] -138197.8894 6175.9286 [400,] 7711.3690 -138197.8894 [401,] 36736.6607 7711.3690 [402,] 6333.6046 36736.6607 [403,] -50925.8076 6333.6046 [404,] -14562.6528 -50925.8076 [405,] 17503.4568 -14562.6528 [406,] 19791.6817 17503.4568 [407,] -2024.5758 19791.6817 [408,] 7417.6489 -2024.5758 [409,] -57365.0314 7417.6489 [410,] -62068.8376 -57365.0314 [411,] -171624.2565 -62068.8376 [412,] -44517.9813 -171624.2565 [413,] 65182.9517 -44517.9813 [414,] -160213.4970 65182.9517 [415,] -88829.3916 -160213.4970 [416,] 36889.0594 -88829.3916 [417,] -29612.3537 36889.0594 [418,] -52885.4894 -29612.3537 [419,] -38685.5501 -52885.4894 [420,] -39444.3383 -38685.5501 [421,] -186172.9561 -39444.3383 [422,] -692206.1533 -186172.9561 [423,] -479941.8205 -692206.1533 [424,] -72026.6461 -479941.8205 [425,] -271484.3057 -72026.6461 [426,] -341158.8368 -271484.3057 [427,] -1055853.5887 -341158.8368 [428,] -259264.3938 -1055853.5887 [429,] -340495.0209 -259264.3938 [430,] -760308.8934 -340495.0209 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2828254.5535 1422177.7855 2 1399259.2669 2828254.5535 3 -2448702.6109 1399259.2669 4 -378022.7100 -2448702.6109 5 212475.4498 -378022.7100 6 -104017.1613 212475.4498 7 325911.2340 -104017.1613 8 -403242.7649 325911.2340 9 652298.8444 -403242.7649 10 105424.8794 652298.8444 11 281145.9654 105424.8794 12 235376.9863 281145.9654 13 -1238522.6120 235376.9863 14 -243356.9332 -1238522.6120 15 576948.3241 -243356.9332 16 82206.4518 576948.3241 17 -233698.3394 82206.4518 18 -399881.1519 -233698.3394 19 -571503.5540 -399881.1519 20 209946.7167 -571503.5540 21 -131539.6641 209946.7167 22 9026.1537 -131539.6641 23 -638528.9441 9026.1537 24 -234635.5535 -638528.9441 25 -101546.3212 -234635.5535 26 -520050.0461 -101546.3212 27 -414492.9031 -520050.0461 28 286243.2595 -414492.9031 29 -323799.7397 286243.2595 30 312171.4007 -323799.7397 31 -43536.2599 312171.4007 32 267807.6146 -43536.2599 33 175433.7361 267807.6146 34 310732.0000 175433.7361 35 -569078.4006 310732.0000 36 49814.5402 -569078.4006 37 208235.4649 49814.5402 38 122385.7917 208235.4649 39 -319960.9457 122385.7917 40 -418353.9839 -319960.9457 41 -295494.5996 -418353.9839 42 18623.3084 -295494.5996 43 -153032.9388 18623.3084 44 -33393.8921 -153032.9388 45 -4999.5642 -33393.8921 46 -68675.5518 -4999.5642 47 -194384.4411 -68675.5518 48 17765.0640 -194384.4411 49 -177775.7603 17765.0640 50 76058.4571 -177775.7603 51 -17208.1270 76058.4571 52 -62407.1809 -17208.1270 53 151479.5222 -62407.1809 54 -284509.5798 151479.5222 55 -14924.3743 -284509.5798 56 293844.3840 -14924.3743 57 -1201.9415 293844.3840 58 -22017.1224 -1201.9415 59 -195116.1026 -22017.1224 60 46919.6024 -195116.1026 61 -197656.7454 46919.6024 62 44537.6686 -197656.7454 63 -169630.2312 44537.6686 64 -42399.5980 -169630.2312 65 -524936.9858 -42399.5980 66 -46853.9476 -524936.9858 67 26048.6923 -46853.9476 68 -19831.6724 26048.6923 69 -42390.0984 -19831.6724 70 -197150.8329 -42390.0984 71 -51774.5780 -197150.8329 72 -35135.1327 -51774.5780 73 32526.8323 -35135.1327 74 38018.9463 32526.8323 75 -86784.1679 38018.9463 76 39064.2674 -86784.1679 77 -64620.0625 39064.2674 78 -52160.0646 -64620.0625 79 39018.8929 -52160.0646 80 10336.1251 39018.8929 81 -30491.9024 10336.1251 82 21645.3152 -30491.9024 83 -4747.8623 21645.3152 84 77437.5254 -4747.8623 85 -124202.2780 77437.5254 86 -19693.9857 -124202.2780 87 56162.0558 -19693.9857 88 -64753.8158 56162.0558 89 86266.4056 -64753.8158 90 33704.0154 86266.4056 91 -6989.6710 33704.0154 92 46043.9995 -6989.6710 93 -84723.7008 46043.9995 94 -7575.3138 -84723.7008 95 1729.0119 -7575.3138 96 142058.1967 1729.0119 97 1785.5551 142058.1967 98 53560.0816 1785.5551 99 6579.9276 53560.0816 100 -29502.7850 6579.9276 101 -80171.6983 -29502.7850 102 -114712.4397 -80171.6983 103 21671.5858 -114712.4397 104 25745.8450 21671.5858 105 -53544.7222 25745.8450 106 -275735.4076 -53544.7222 107 176036.5982 -275735.4076 108 -132744.2026 176036.5982 109 -87592.5503 -132744.2026 110 90557.5062 -87592.5503 111 -228631.1611 90557.5062 112 -146148.2854 -228631.1611 113 -49099.8440 -146148.2854 114 70476.2658 -49099.8440 115 -4009.2257 70476.2658 116 44140.8107 -4009.2257 117 102558.5655 44140.8107 118 -118258.5694 102558.5655 119 816.2494 -118258.5694 120 -243987.8612 816.2494 121 30189.5715 -243987.8612 122 -92740.5446 30189.5715 123 142401.7259 -92740.5446 124 3826.0383 142401.7259 125 -153906.6564 3826.0383 126 45876.2528 -153906.6564 127 30075.6448 45876.2528 128 9665.3884 30075.6448 129 45876.2528 9665.3884 130 35076.0683 45876.2528 131 33673.0231 35076.0683 132 45876.2528 33673.0231 133 16601.0575 45876.2528 134 23327.4580 16601.0575 135 33445.0498 23327.4580 136 33987.1176 33445.0498 137 43983.9643 33987.1176 138 20663.1323 43983.9643 139 45876.2528 20663.1323 140 37187.6608 45876.2528 141 45876.2528 37187.6608 142 -7005.3106 45876.2528 143 40383.7450 -7005.3106 144 30884.5952 40383.7450 145 25016.3156 30884.5952 146 45876.2528 25016.3156 147 36255.9339 45876.2528 148 47247.8386 36255.9339 149 46067.5680 47247.8386 150 1761.4402 46067.5680 151 -61450.0008 1761.4402 152 39425.6762 -61450.0008 153 9230.1870 39425.6762 154 42724.3296 9230.1870 155 45876.2528 42724.3296 156 45876.2528 45876.2528 157 41012.2288 45876.2528 158 51627.6528 41012.2288 159 45876.2528 51627.6528 160 45876.2528 45876.2528 161 45876.2528 45876.2528 162 45876.2528 45876.2528 163 17187.3214 45876.2528 164 -8192.4349 17187.3214 165 43061.9100 -8192.4349 166 45876.2528 43061.9100 167 73457.3654 45876.2528 168 -100530.2721 73457.3654 169 45876.2528 -100530.2721 170 45876.2528 45876.2528 171 38251.8603 45876.2528 172 45876.2528 38251.8603 173 24048.7859 45876.2528 174 45876.2528 24048.7859 175 32971.1674 45876.2528 176 -110753.6285 32971.1674 177 45121.9943 -110753.6285 178 45876.2528 45121.9943 179 45876.2528 45876.2528 180 44095.8056 45876.2528 181 34778.7980 44095.8056 182 -41202.6439 34778.7980 183 -986.0808 -41202.6439 184 38349.4798 -986.0808 185 35584.4931 38349.4798 186 2711.9615 35584.4931 187 32834.2016 2711.9615 188 10898.0508 32834.2016 189 45876.2528 10898.0508 190 34205.9235 45876.2528 191 45876.2528 34205.9235 192 33752.7633 45876.2528 193 45876.2528 33752.7633 194 41702.5373 45876.2528 195 45876.2528 41702.5373 196 -83418.4463 45876.2528 197 45876.2528 -83418.4463 198 45876.2528 45876.2528 199 -34734.8816 45876.2528 200 45876.2528 -34734.8816 201 45876.2528 45876.2528 202 45876.2528 45876.2528 203 11544.4399 45876.2528 204 29663.8946 11544.4399 205 39982.3299 29663.8946 206 45114.4407 39982.3299 207 45876.2528 45114.4407 208 15762.1198 45876.2528 209 27787.7310 15762.1198 210 45876.2528 27787.7310 211 45876.2528 45876.2528 212 40534.1392 45876.2528 213 -3769.5130 40534.1392 214 16870.9861 -3769.5130 215 43296.0221 16870.9861 216 38116.0415 43296.0221 217 30533.0992 38116.0415 218 43296.0221 30533.0992 219 -24039.3876 43296.0221 220 16993.5346 -24039.3876 221 45876.2528 16993.5346 222 44494.6336 45876.2528 223 6636.3416 44494.6336 224 44345.4449 6636.3416 225 10965.6898 44345.4449 226 32104.3027 10965.6898 227 45876.2528 32104.3027 228 36062.7740 45876.2528 229 45876.2528 36062.7740 230 13243.8692 45876.2528 231 24147.2612 13243.8692 232 45876.2528 24147.2612 233 8919.6114 45876.2528 234 45876.2528 8919.6114 235 45939.5835 45876.2528 236 45876.2528 45939.5835 237 37391.3257 45876.2528 238 45876.2528 37391.3257 239 45876.2528 45876.2528 240 45876.2528 45876.2528 241 37460.4520 45876.2528 242 45876.2528 37460.4520 243 36578.4145 45876.2528 244 1382.6751 36578.4145 245 45876.2528 1382.6751 246 45876.2528 45876.2528 247 17273.1870 45876.2528 248 45876.2528 17273.1870 249 45876.2528 45876.2528 250 11089.2820 45876.2528 251 45876.2528 11089.2820 252 32881.7550 45876.2528 253 -27733.2665 32881.7550 254 -140283.3289 -27733.2665 255 45876.2528 -140283.3289 256 23055.8072 45876.2528 257 -12528.5493 23055.8072 258 45876.2528 -12528.5493 259 45876.2528 45876.2528 260 45876.2528 45876.2528 261 45876.2528 45876.2528 262 52719.1282 45876.2528 263 42846.8695 52719.1282 264 45876.2528 42846.8695 265 45876.2528 45876.2528 266 45876.2528 45876.2528 267 95960.3053 45876.2528 268 45876.2528 95960.3053 269 25360.0761 45876.2528 270 45876.2528 25360.0761 271 45876.2528 45876.2528 272 45876.2528 45876.2528 273 13455.7335 45876.2528 274 45876.2528 13455.7335 275 28873.5045 45876.2528 276 13907.0008 28873.5045 277 45876.2528 13907.0008 278 -43985.1358 45876.2528 279 27002.4361 -43985.1358 280 45876.2528 27002.4361 281 45876.2528 45876.2528 282 5007.4327 45876.2528 283 45876.2528 5007.4327 284 42937.2846 45876.2528 285 45876.2528 42937.2846 286 30133.0662 45876.2528 287 44068.8680 30133.0662 288 45876.2528 44068.8680 289 1041.9318 45876.2528 290 45876.2528 1041.9318 291 36588.1557 45876.2528 292 62051.9763 36588.1557 293 45876.2528 62051.9763 294 18846.8644 45876.2528 295 51512.3879 18846.8644 296 45876.2528 51512.3879 297 45876.2528 45876.2528 298 45876.2528 45876.2528 299 34265.2149 45876.2528 300 37844.7527 34265.2149 301 24038.0960 37844.7527 302 44586.1374 24038.0960 303 45876.2528 44586.1374 304 48868.1653 45876.2528 305 30365.1997 48868.1653 306 45876.2528 30365.1997 307 18067.6250 45876.2528 308 10015.7961 18067.6250 309 45876.2528 10015.7961 310 45876.2528 45876.2528 311 27853.9242 45876.2528 312 45876.2528 27853.9242 313 33395.7416 45876.2528 314 35765.8263 33395.7416 315 61876.7863 35765.8263 316 44761.5447 61876.7863 317 -6200.1958 44761.5447 318 39395.8415 -6200.1958 319 45876.2528 39395.8415 320 45876.2528 45876.2528 321 45876.2528 45876.2528 322 45876.2528 45876.2528 323 45876.2528 45876.2528 324 45876.2528 45876.2528 325 41227.7757 45876.2528 326 45876.2528 41227.7757 327 45876.2528 45876.2528 328 45876.2528 45876.2528 329 45876.2528 45876.2528 330 45876.2528 45876.2528 331 45876.2528 45876.2528 332 45876.2528 45876.2528 333 36968.6674 45876.2528 334 44572.3028 36968.6674 335 33850.7190 44572.3028 336 -17595.4776 33850.7190 337 45876.2528 -17595.4776 338 42518.4372 45876.2528 339 32563.7714 42518.4372 340 3957.1427 32563.7714 341 45876.2528 3957.1427 342 45876.2528 45876.2528 343 45876.2528 45876.2528 344 8421.4511 45876.2528 345 45876.2528 8421.4511 346 27704.5864 45876.2528 347 45876.2528 27704.5864 348 45876.2528 45876.2528 349 2137.8717 45876.2528 350 38503.5036 2137.8717 351 34204.3680 38503.5036 352 30477.1743 34204.3680 353 40468.4630 30477.1743 354 45876.2528 40468.4630 355 45876.2528 45876.2528 356 22299.5254 45876.2528 357 45876.2528 22299.5254 358 45876.2528 45876.2528 359 45876.2528 45876.2528 360 40711.9569 45876.2528 361 45876.2528 40711.9569 362 45876.2528 45876.2528 363 -35785.4212 45876.2528 364 31340.9823 -35785.4212 365 45876.2528 31340.9823 366 45876.2528 45876.2528 367 45876.2528 45876.2528 368 45876.2528 45876.2528 369 32258.4971 45876.2528 370 45876.2528 32258.4971 371 35044.5333 45876.2528 372 53147.9855 35044.5333 373 30250.6819 53147.9855 374 45876.2528 30250.6819 375 45876.2528 45876.2528 376 9032.9946 45876.2528 377 30551.8278 9032.9946 378 -10839.0923 30551.8278 379 45876.2528 -10839.0923 380 45876.2528 45876.2528 381 12992.8254 45876.2528 382 -19634.9329 12992.8254 383 45876.2528 -19634.9329 384 30641.7681 45876.2528 385 -123138.4747 30641.7681 386 50836.5709 -123138.4747 387 11640.2014 50836.5709 388 39263.8138 11640.2014 389 15615.9244 39263.8138 390 31253.7532 15615.9244 391 5065.9793 31253.7532 392 -30129.8069 5065.9793 393 75257.3307 -30129.8069 394 -71988.3779 75257.3307 395 -8508.7294 -71988.3779 396 -166757.5998 -8508.7294 397 -170833.8386 -166757.5998 398 6175.9286 -170833.8386 399 -138197.8894 6175.9286 400 7711.3690 -138197.8894 401 36736.6607 7711.3690 402 6333.6046 36736.6607 403 -50925.8076 6333.6046 404 -14562.6528 -50925.8076 405 17503.4568 -14562.6528 406 19791.6817 17503.4568 407 -2024.5758 19791.6817 408 7417.6489 -2024.5758 409 -57365.0314 7417.6489 410 -62068.8376 -57365.0314 411 -171624.2565 -62068.8376 412 -44517.9813 -171624.2565 413 65182.9517 -44517.9813 414 -160213.4970 65182.9517 415 -88829.3916 -160213.4970 416 36889.0594 -88829.3916 417 -29612.3537 36889.0594 418 -52885.4894 -29612.3537 419 -38685.5501 -52885.4894 420 -39444.3383 -38685.5501 421 -186172.9561 -39444.3383 422 -692206.1533 -186172.9561 423 -479941.8205 -692206.1533 424 -72026.6461 -479941.8205 425 -271484.3057 -72026.6461 426 -341158.8368 -271484.3057 427 -1055853.5887 -341158.8368 428 -259264.3938 -1055853.5887 429 -340495.0209 -259264.3938 430 -760308.8934 -340495.0209 > 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/7l64k1292875566.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/8l64k1292875566.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/9eg3n1292875566.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/10eg3n1292875566.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/11hykb1292875566.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/122yih1292875566.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/13mtnw1292875566.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/14k9fe1292875566.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/15nadk1292875566.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/16j1ba1292875566.tab") + } > > try(system("convert tmp/17eoc1292875566.ps tmp/17eoc1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/20o6x1292875566.ps tmp/20o6x1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/30o6x1292875566.ps tmp/30o6x1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/40o6x1292875566.ps tmp/40o6x1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/5axnh1292875566.ps tmp/5axnh1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/6axnh1292875566.ps tmp/6axnh1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/7l64k1292875566.ps tmp/7l64k1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/8l64k1292875566.ps tmp/8l64k1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/9eg3n1292875566.ps tmp/9eg3n1292875566.png",intern=TRUE)) character(0) > try(system("convert tmp/10eg3n1292875566.ps tmp/10eg3n1292875566.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 11.253 2.312 23.541