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 <- c(7.03783512969226 + ,-98.3312467372175 + ,-117.325401144406 + ,-42.7948033030729 + ,-165.254264496858 + ,-196.555340415592 + ,145.254550435451 + ,-79.0414259157364 + ,-34.0832910828553 + ,171.087420498627 + ,-28.1542367670882 + ,-68.9435615590894 + ,-156.001393418334 + ,-104.456934875882 + ,130.549461497950 + ,34.8922283561497 + ,-123.098349508062 + ,74.5576548712037 + ,67.0810176316715 + ,152.706476785095 + ,91.2303867678834 + ,-110.882335556812 + ,-63.1515441793628 + ,11.9161441760431 + ,56.111817635295 + ,34.6434739833669 + ,0.459168959230808 + ,124.384105899383 + ,-277.557054985246 + ,77.3190779022267 + ,86.6189780737299 + ,-118.689668627819 + ,267.580038935061 + ,42.1259919911072 + ,45.9789517089369 + ,-99.395391658839 + ,227.000383506021 + ,84.4157532050709 + ,105.331312172398 + ,-192.382641815791 + ,-56.2564209143736 + ,39.5829396160272 + ,241.089876240609 + ,294.048754836687 + ,80.3314860927041 + ,154.577463913597 + ,-204.181032235406 + ,292.712866704525 + ,137.895615077045 + ,-22.3900588496519 + ,164.690640469680 + ,14.2498538101749 + ,164.397885385949 + ,-5.26345425644624 + ,-131.777413552913 + ,-128.88093775439 + ,-169.431718103336 + ,-12.2335813289512 + ,-32.2123203622983 + ,-91.8937474036534 + ,149.134033642112 + ,195.332136919032 + ,18.3981744397026 + ,121.462520305033 + ,153.777651802018 + ,45.3498183649984 + ,-124.422532349626 + ,-204.260726112894 + ,-145.205342416386 + ,-119.804957070208 + ,148.122321310936 + ,-93.2058969118374 + ,-65.2388449926761 + ,-88.3620370284215 + ,-129.532157838408 + ,-9.2469890086225 + ,28.3388095679010 + ,-240.907523502917 + ,-16.6782243489416 + ,139.666336757527 + ,-234.487641945176 + ,228.394801661146 + ,-184.933176456224 + ,-28.666192801664 + ,-14.8453193979652 + ,-199.483852991607 + ,13.6465274997104 + ,-230.29179207403 + ,81.8579620197308 + ,4.86081765309311 + ,-17.5590342194350 + ,179.618783151966 + ,-94.3631378156688 + ,-106.097639570226 + ,-107.534979176420 + ,21.9000068305806 + ,-142.522113796213 + ,81.2097507416726 + ,-34.3238154153577 + ,62.2956437618072 + ,-147.206044248034 + ,-67.4179711660492 + ,66.1124152079459 + ,117.156704601882 + ,244.820327089704 + ,-145.019708301844 + ,-22.6714338685113 + ,8.34530597612882 + ,-142.513254377861 + ,154.846929269182 + ,47.7414082166392 + ,-5.25010190300413 + ,-35.8451116699208 + ,-117.492487104458 + ,99.5934034310243 + ,149.795051334379 + ,79.603783152614 + ,-106.925290146154 + ,213.396122633205 + ,-51.5978154550942 + ,-52.7767772308038 + ,-76.7178470217067 + ,-124.740443525246 + ,38.5615824089955 + ,50.7924078971276 + ,-109.513841122060 + ,97.3231116440577 + ,125.340786977475 + ,-30.8742409164725 + ,-112.250837824287 + ,28.0870498207617 + ,-21.8358842027445 + ,-54.7201336644692 + ,93.2607812073961 + ,-87.500795519974 + ,24.4728882163302 + ,-41.541757987332 + ,107.428710097945 + ,-224.631733186458 + ,-49.2094660310937 + ,-63.6220721404223 + ,118.202235380549 + ,87.9365313771192 + ,16.8857286629086 + ,3.40429185321661 + ,-73.7608136638874 + ,141.323268192647 + ,-59.2442489060646 + ,88.2489237339183 + ,181.402622666532 + ,-155.569657072819 + ,-328.526166234309 + ,-19.9557930759967 + ,149.911989238589 + ,22.4025271335545 + ,31.1270716215838 + ,8.7700149451982 + ,132.375377611128 + ,-59.4000352570887 + ,156.352868592634 + ,-35.91415873356 + ,105.827992809706 + ,32.1302880061727 + ,-4.83085047993892 + ,-126.535939617738 + ,-225.420029592616 + ,1.76168052830480 + ,70.0224321866176 + ,-21.7339751364453 + ,-105.449404941814 + ,-36.659878508723 + ,-76.9976316993702 + ,179.970481462103 + ,20.4389668756409 + ,-191.810150339689 + ,-192.074837155084 + ,154.993155076085 + ,74.334476707483 + ,95.1535385324919 + ,2.62893383453497 + ,60.9978202261871 + ,-7.13200682666083 + ,-32.6574342387546 + ,39.5365173594474 + ,135.289040674639 + ,80.3820167101494 + ,27.2451043992194 + ,-52.8188161337465) > par1 = '12' > #'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!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 188 > (np <- floor(n / par1)) [1] 15 > arr <- array(NA,dim=c(par1,np+1)) > darr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > dx <- diff(x) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + darr[j,ari[j]] <- dx[i] + if (j == par1) j = 0 + } > ari [1] 16 16 16 16 16 16 16 16 15 15 15 15 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 7.037835 -156.00139 56.111818 227.00038 137.895615 149.13403 [2,] -98.331247 -104.45693 34.643474 84.41575 -22.390059 195.33214 [3,] -117.325401 130.54946 0.459169 105.33131 164.690640 18.39817 [4,] -42.794803 34.89223 124.384106 -192.38264 14.249854 121.46252 [5,] -165.254264 -123.09835 -277.557055 -56.25642 164.397885 153.77765 [6,] -196.555340 74.55765 77.319078 39.58294 -5.263454 45.34982 [7,] 145.254550 67.08102 86.618978 241.08988 -131.777414 -124.42253 [8,] -79.041426 152.70648 -118.689669 294.04875 -128.880938 -204.26073 [9,] -34.083291 91.23039 267.580039 80.33149 -169.431718 -145.20534 [10,] 171.087420 -110.88234 42.125992 154.57746 -12.233581 -119.80496 [11,] -28.154237 -63.15154 45.978952 -204.18103 -32.212320 148.12232 [12,] -68.943562 11.91614 -99.395392 292.71287 -91.893747 -93.20590 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -65.238845 -14.845319 -142.522114 -142.513254 -52.77678 -54.72013 [2,] -88.362037 -199.483853 81.209751 154.846929 -76.71785 93.26078 [3,] -129.532158 13.646527 -34.323815 47.741408 -124.74044 -87.50080 [4,] -9.246989 -230.291792 62.295644 -5.250102 38.56158 24.47289 [5,] 28.338810 81.857962 -147.206044 -35.845112 50.79241 -41.54176 [6,] -240.907524 4.860818 -67.417971 -117.492487 -109.51384 107.42871 [7,] -16.678224 -17.559034 66.112415 99.593403 97.32311 -224.63173 [8,] 139.666337 179.618783 117.156705 149.795051 125.34079 -49.20947 [9,] -234.487642 -94.363138 244.820327 79.603783 -30.87424 -63.62207 [10,] 228.394802 -106.097640 -145.019708 -106.925290 -112.25084 118.20224 [11,] -184.933176 -107.534979 -22.671434 213.396123 28.08705 87.93653 [12,] -28.666193 21.900007 8.345306 -51.597815 -21.83588 16.88573 [,13] [,14] [,15] [,16] [1,] 3.404292 8.770015 -21.733975 60.997820 [2,] -73.760814 132.375378 -105.449405 -7.132007 [3,] 141.323268 -59.400035 -36.659879 -32.657434 [4,] -59.244249 156.352869 -76.997632 39.536517 [5,] 88.248924 -35.914159 179.970481 135.289041 [6,] 181.402623 105.827993 20.438967 80.382017 [7,] -155.569657 32.130288 -191.810150 27.245104 [8,] -328.526166 -4.830850 -192.074837 -52.818816 [9,] -19.955793 -126.535940 154.993155 NA [10,] 149.911989 -225.420030 74.334477 NA [11,] 22.402527 1.761681 95.153539 NA [12,] 31.127072 70.022432 2.628934 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -105.36908 51.544459 -21.46834 -142.58463 -160.285674 46.19810 [2,] -18.99415 235.006396 -34.18431 20.91556 187.080699 -176.93396 [3,] 74.53060 -95.657233 123.92494 -297.71395 -150.440787 103.06435 [4,] -122.45946 -157.990578 -401.94116 136.12622 150.148032 32.31513 [5,] -31.30108 197.656004 354.87613 95.83936 -169.661340 -108.42783 [6,] 341.80989 -7.476637 9.29990 201.50694 -126.513959 -169.77235 [7,] -224.29598 85.625459 -205.30865 52.95888 2.896476 -79.83819 [8,] 44.95813 -61.476090 386.26971 -213.71727 -40.550780 59.05538 [9,] 205.17071 -202.112722 -225.45405 74.24598 157.198137 25.40039 [10,] -199.24166 47.730791 3.85296 -358.75850 -19.978739 267.92728 [11,] -40.78932 75.067688 -145.37434 496.89390 -59.681427 -241.32822 [12,] -87.05783 44.195673 326.39578 -154.81725 241.027781 27.96705 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -23.12319 -184.638534 223.73186 297.360184 -23.94107 147.98091 [2,] -41.17012 213.130380 -115.53357 -107.105521 -48.02260 -180.76158 [3,] 120.28517 -243.938320 96.61946 -52.991510 163.30203 111.97368 [4,] 37.58580 312.149754 -209.50169 -30.595010 12.23083 -66.01465 [5,] -269.24633 -76.997144 79.78807 -81.647375 -160.30625 148.97047 [6,] 224.22930 -22.419852 133.53039 217.085891 206.83695 -332.06044 [7,] 156.34456 197.177817 51.04429 50.201648 28.01768 175.42227 [8,] -374.15398 -273.981921 127.66362 -70.191268 -156.21503 -14.41261 [9,] 462.88244 -11.734502 -389.84004 -186.529073 -81.37660 181.82431 [10,] -413.32798 -1.437340 122.34827 320.321413 140.33789 -30.26570 [11,] 156.26698 129.434986 31.01674 -264.993938 -49.92293 -71.05080 [12,] 13.82087 -164.422121 -150.85856 -1.178962 -32.88425 -13.48144 [,13] [,14] [,15] [,16] [1,] -77.165106 123.60536 -83.7154298 -68.12983 [2,] 215.084082 -191.77541 68.7895264 -25.52543 [3,] -200.567517 215.75290 -40.3377532 72.19395 [4,] 147.493173 -192.26703 256.9681132 95.75252 [5,] 93.153699 141.74215 -159.5315146 -54.90702 [6,] -336.972280 -73.69770 -212.2491172 -53.13691 [7,] -172.956509 -36.96114 -0.2646868 -80.06392 [8,] 308.570373 -121.70509 347.0679922 NA [9,] 169.867782 -98.88409 -80.6586784 NA [10,] -127.509462 227.18171 20.8190618 NA [11,] 8.724544 68.26075 -92.5246047 NA [12,] -22.357057 -91.75641 58.3688864 NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/1jf721261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2lylv1261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/37bvx1261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4ke9u1261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -156.001393 -199.48385 -129.53216 -192.38264 -277.557055 -240.90752 [2,] -59.979489 -93.34664 -73.45042 -51.01953 -89.677385 -88.46591 [3,] -5.720514 -14.76103 -16.09913 19.36137 -3.753151 30.01095 [4,] 58.554819 88.83827 76.53636 50.91608 111.768982 78.85055 [5,] 227.000384 195.33214 164.69064 156.35287 179.970481 181.40262 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -224.63173 -328.52617 -234.48764 -225.42003 -204.181032 -99.395392 [2,] -128.09997 -123.78530 -110.44954 -111.56659 -47.681932 -60.270689 [3,] 29.68770 -27.02016 -30.87424 -12.23358 1.761681 2.628934 [4,] 91.97104 144.73069 85.78094 134.05711 66.957742 19.392868 [5,] 241.08988 294.04875 267.58004 228.39480 213.396123 70.022432 $n [1] 16 16 16 16 16 16 16 16 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -52.54157 -86.72407 -75.34391 -20.90319 -83.32447 -36.07905 -57.24036 [2,] 41.10054 57.20201 43.14564 59.62594 75.81816 96.10095 116.61575 [,8] [,9] [,10] [,11] [,12] [1,] -133.08398 -110.92730 -112.43681 -45.00606 -29.87015 [2,] 79.04366 49.17881 87.96965 48.52942 35.12802 $out [1] -230.2918 292.7129 $group [1] 4 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-156.001393418334, -59.9794893285727, -5.7205137723743, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5wilk1261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -184.63853 -191.77541 -297.71395 -401.94116 -269.24633 -336.97228 [2,] -94.54226 -111.31954 -123.04901 -140.22502 -133.97967 -148.14316 [3,] -23.53213 -29.85487 73.36227 22.27298 -43.10405 -14.94824 [4,] 87.57491 127.93511 116.12943 141.80970 118.79076 204.17194 [5,] 297.36018 235.00640 215.75290 312.14975 354.87613 341.80989 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -224.29598 -374.15398 -389.84004 -358.75850 -264.99394 -164.42212 [2,] -79.95106 -138.96006 -142.70658 -78.88758 -81.78770 -89.40712 [3,] 15.45708 -40.55078 -11.73450 3.85296 -40.78932 -13.48144 [4,] 69.29217 93.35950 163.53296 131.34308 71.66422 36.08136 [5,] 197.17782 386.26971 462.88244 320.32141 156.26698 58.36889 $n [1] 16 16 16 16 16 16 16 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -95.46841 -124.36046 -21.11321 -89.13073 -142.94837 -154.1127 -43.49400 [2,] 48.40415 64.65072 167.83776 133.67669 56.74027 124.2162 74.40815 [,8] [,9] [,10] [,11] [,12] [1,] -135.32653 -136.6662 -81.91153 -103.39069 -64.67500 [2,] 54.22497 113.1972 89.61745 21.81204 37.71212 $out [1] -413.3280 496.8939 326.3958 241.0278 $group [1] 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-184.638533593642, -94.5422558361392, -23.5321309133241, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/61ryl1261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -196.555340 -156.00139 -118.68967 -204.18103 -169.43172 -204.26073 [2,] -107.828324 -107.66964 -49.46811 -8.33674 -110.38734 -122.11374 [3,] -55.869182 23.40419 44.05247 94.87353 -17.31182 31.87400 [4,] -10.558201 82.89402 81.96903 234.04513 76.07273 148.62818 [5,] 7.037835 152.70648 267.58004 294.04875 164.69064 195.33214 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -240.90752 -230.29179 -147.206044 -142.51325 -124.74044 -224.63173 [2,] -157.23267 -106.81631 -104.970042 -79.26155 -93.11584 -59.17110 [3,] -46.95252 -16.20218 -7.163064 21.24565 -26.35506 -12.32801 [4,] 9.54591 17.77327 73.661083 124.69423 44.67700 90.59866 [5,] 228.39480 179.61878 244.820327 213.39612 125.34079 118.20224 [,13] [,14] [,15] [,16] [1,] -328.52617 -225.420030 -192.07484 -52.81882 [2,] -66.50253 -47.657097 -91.22352 -19.89472 [3,] 12.90341 5.265848 -9.55252 33.39081 [4,] 114.78610 87.925212 84.74401 70.68992 [5,] 181.40262 156.352869 179.97048 135.28904 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -100.23474 -63.51318 -15.89689 -15.67847 -102.35751 -91.61321 -123.02135 [2,] -11.50363 110.32155 104.00183 205.42553 67.73387 155.36120 29.11631 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -73.02832 -88.63792 -71.77995 -89.20330 -80.63901 -69.78355 -56.57415 [2,] 40.62397 74.31179 114.27125 36.49318 55.98298 95.59037 67.10585 [,15] [,16] [1,] -89.81249 -17.21107 [2,] 70.70745 83.99269 $out [1] 145.2546 171.0874 -277.5571 $group [1] 1 1 3 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] NA Warning message: In bxp(list(stats = c(-196.555340415592, -107.828323940812, -55.8691824310811, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7j5t81261344666.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] -4.079700e-14 -30.874241 -20.4389104 [2,] -2.808864e-14 -15.430083 -7.2387359 [3,] -9.850454e-15 -4.736832 -0.4243397 [4,] 1.211863e-15 10.995152 10.0592097 [5,] 4.091564e-14 30.010953 13.5210250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -2.321461e-14 -16.789560 -8.314049 [2,] 3.513704e-15 7.315895 7.465369 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-4.07969954115591e-14, -2.80886425230165e-14, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > try(system("convert tmp/1jf721261344666.ps tmp/1jf721261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/2lylv1261344666.ps tmp/2lylv1261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/37bvx1261344666.ps tmp/37bvx1261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/4ke9u1261344666.ps tmp/4ke9u1261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/5wilk1261344666.ps tmp/5wilk1261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/61ryl1261344666.ps tmp/61ryl1261344666.png",intern=TRUE)) character(0) > try(system("convert tmp/7j5t81261344666.ps tmp/7j5t81261344666.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.210 0.943 1.537