R version 2.8.0 (2008-10-20) Copyright (C) 2008 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(-183.9235445 + ,-177.0726091 + ,-228.6351091 + ,-237.4476091 + ,-127.7601091 + ,-193.0101091 + ,-220.6351091 + ,-164.5101091 + ,-268.3226091 + ,-333.6976091 + ,-34.26010911 + ,-154.8851091 + ,-97.74528053 + ,101.1056549 + ,2.543154874 + ,-43.26934513 + ,-163.5818451 + ,-162.8318451 + ,46.54315487 + ,26.66815487 + ,-107.1443451 + ,42.48065487 + ,76.91815487 + ,196.2931549 + ,201.4329835 + ,12.28391886 + ,-0.278581137 + ,42.90891886 + ,87.59641886 + ,84.34641886 + ,57.72141886 + ,173.8464189 + ,-185.9660811 + ,47.65891886 + ,89.09641886 + ,-68.52858114 + ,272.6112475 + ,146.4621829 + ,162.8996829 + ,10.08718285 + ,279.7746829 + ,212.5246829 + ,248.8996829 + ,-41.97531715 + ,-5.787817149 + ,52.83718285 + ,274.2746829 + ,414.6496829 + ,310.7895114 + ,362.6404468 + ,26.07794684 + ,403.2654468 + ,327.9529468 + ,193.7029468 + ,317.0779468 + ,202.2029468 + ,321.3904468 + ,178.0154468 + ,16.45294684 + ,-68.17205316 + ,-157.0322246 + ,-76.18128917 + ,-81.74378917 + ,-134.5562892 + ,77.13121083 + ,199.8812108 + ,105.2562108 + ,198.3812108 + ,262.5687108 + ,196.1937108 + ,11.63121083 + ,-145.9937892 + ,-166.8539606 + ,-202.0030252 + ,43.43447482 + ,-113.3780252 + ,-113.6905252 + ,-155.9405252 + ,-210.5655252 + ,-124.4405252 + ,-64.25302518 + ,-298.6280252 + ,-154.1905252 + ,23.18447482 + ,-249.6756966 + ,118.1752388 + ,-180.3872612 + ,-79.19976119 + ,-81.51226119 + ,-246.7622612 + ,-105.3872612 + ,-319.2622612 + ,-72.07476119 + ,-90.44976119 + ,-80.01226119 + ,119.3627388 + ,-53.49743261 + ,-114.6464972 + ,-155.2089972 + ,-50.02149721 + ,-196.3339972 + ,-14.58399721 + ,-82.20899721 + ,17.91600279 + ,-162.8964972 + ,-132.2714972 + ,-16.83399721 + ,81.54100279 + ,275.6808314 + ,-32.46823322 + ,17.96926678 + ,27.15676678 + ,-123.1557332 + ,108.5942668 + ,67.96926678 + ,34.09426678 + ,-13.71823322 + ,-113.0932332 + ,54.34426678 + ,149.7192668 + ,153.8590954 + ,-28.28996923 + ,238.1475308 + ,50.33503077 + ,8.022530771 + ,-61.22746923 + ,-140.8524692 + ,-28.72746923 + ,9.460030771 + ,-121.9149692 + ,41.52253077 + ,115.8975308 + ,27.03735936 + ,-91.11170524 + ,3.325794759 + ,-29.48670524 + ,-73.79920524 + ,50.95079476 + ,-86.67420524 + ,-9.54920524 + ,-66.36170524 + ,73.26329476 + ,-216.2992052 + ,-128.9242052 + ,-142.7843767 + ,27.06655875 + ,60.50405875 + ,35.69155875 + ,16.37905875 + ,-64.87094125 + ,115.5040587 + ,-30.37094125 + ,87.81655875 + ,205.4415587 + ,-64.12094125 + ,-322.7459413 + ,-139.6061127 + ,35.24482274 + ,-4.317677263 + ,17.86982274 + ,2.557322737 + ,129.3073227 + ,-16.31767726 + ,164.8073227 + ,21.99482274 + ,138.6198227 + ,87.05732274 + ,51.43232274 + ,-80.42784867 + ,-105.1918797 + ,5.245620328 + ,68.43312033 + ,-0.879379672 + ,-105.1293797 + ,-82.75437967 + ,-132.6293797 + ,102.5581203 + ,23.18312033 + ,-180.3793797 + ,-267.0043797 + ,30.13544892 + ,23.98638432 + ,90.42388432 + ,31.61138432 + ,81.29888432 + ,25.04888432 + ,-13.57611568 + ,33.54888432 + ,140.7363843 + ,132.3613843 + ,94.79888432 + ,4.173884316) > 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] 192 > (np <- floor(n / par1)) [1] 16 > 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 16 16 16 16 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -183.92354 -97.745281 201.4329835 272.611247 310.78951 -157.03222 [2,] -177.07261 101.105655 12.2839189 146.462183 362.64045 -76.18129 [3,] -228.63511 2.543155 -0.2785811 162.899683 26.07795 -81.74379 [4,] -237.44761 -43.269345 42.9089189 10.087183 403.26545 -134.55629 [5,] -127.76011 -163.581845 87.5964189 279.774683 327.95295 77.13121 [6,] -193.01011 -162.831845 84.3464189 212.524683 193.70295 199.88121 [7,] -220.63511 46.543155 57.7214189 248.899683 317.07795 105.25621 [8,] -164.51011 26.668155 173.8464189 -41.975317 202.20295 198.38121 [9,] -268.32261 -107.144345 -185.9660811 -5.787817 321.39045 262.56871 [10,] -333.69761 42.480655 47.6589189 52.837183 178.01545 196.19371 [11,] -34.26011 76.918155 89.0964189 274.274683 16.45295 11.63121 [12,] -154.88511 196.293155 -68.5285811 414.649683 -68.17205 -145.99379 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -166.85396 -249.67570 -53.49743 275.68083 153.85910 27.037359 [2,] -202.00303 118.17524 -114.64650 -32.46823 -28.28997 -91.111705 [3,] 43.43447 -180.38726 -155.20900 17.96927 238.14753 3.325795 [4,] -113.37803 -79.19976 -50.02150 27.15677 50.33503 -29.486705 [5,] -113.69053 -81.51226 -196.33400 -123.15573 8.02253 -73.799205 [6,] -155.94053 -246.76226 -14.58400 108.59427 -61.22747 50.950795 [7,] -210.56553 -105.38726 -82.20900 67.96927 -140.85247 -86.674205 [8,] -124.44053 -319.26226 17.91600 34.09427 -28.72747 -9.549205 [9,] -64.25303 -72.07476 -162.89650 -13.71823 9.46003 -66.361705 [10,] -298.62803 -90.44976 -132.27150 -113.09323 -121.91497 73.263295 [11,] -154.19053 -80.01226 -16.83400 54.34427 41.52253 -216.299205 [12,] 23.18447 119.36274 81.54100 149.71927 115.89753 -128.924205 [,13] [,14] [,15] [,16] [,17] [1,] -142.78438 -139.606113 -80.4278487 30.135449 NA [2,] 27.06656 35.244823 -105.1918797 23.986384 NA [3,] 60.50406 -4.317677 5.2456203 90.423884 NA [4,] 35.69156 17.869823 68.4331203 31.611384 NA [5,] 16.37906 2.557323 -0.8793797 81.298884 NA [6,] -64.87094 129.307323 -105.1293797 25.048884 NA [7,] 115.50406 -16.317677 -82.7543797 -13.576116 NA [8,] -30.37094 164.807323 -132.6293797 33.548884 NA [9,] 87.81656 21.994823 102.5581203 140.736384 NA [10,] 205.44156 138.619823 23.1831203 132.361384 NA [11,] -64.12094 87.057323 -180.3793797 94.798884 NA [12,] -322.74594 51.432323 -267.0043797 4.173884 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 6.850935 198.850935 -189.1491 -126.1491 51.85094 80.85094 [2,] -51.562500 -98.562500 -12.5625 16.4375 -336.56250 -5.56250 [3,] -8.812500 -45.812500 43.1875 -152.8125 377.18750 -52.81250 [4,] 109.687500 -120.312500 44.6875 269.6875 -75.31250 211.68750 [5,] -65.250000 0.750000 -3.2500 -67.2500 -134.25000 122.75000 [6,] -27.625000 209.375000 -26.6250 36.3750 123.37500 -94.62500 [7,] 56.125000 -19.875000 116.1250 -290.8750 -114.87500 93.12500 [8,] -103.812500 -133.812500 -359.8125 36.1875 119.18750 64.18750 [9,] -65.375000 149.625000 233.6250 58.6250 -143.37500 -66.37500 [10,] 299.437500 34.437500 41.4375 221.4375 -161.56250 -184.56250 [11,] -120.625000 119.375000 -157.6250 140.3750 -84.62500 -157.62500 [12,] 57.139829 5.139829 341.1398 -103.8602 -88.86017 -20.86017 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -35.14906 367.8509 -61.14906 -308.149065 -182.14906 -118.14906 [2,] 245.43750 -298.5625 -40.56250 50.437500 266.43750 94.43750 [3,] -156.81250 101.1875 105.18750 9.187500 -187.81250 -32.81250 [4,] -0.31250 -2.3125 -146.31250 -150.312500 -42.31250 -44.31250 [5,] -42.25000 -165.2500 181.75000 231.750000 -69.25000 124.75000 [6,] -54.62500 141.3750 -67.62500 -40.625000 -79.62500 -137.62500 [7,] 86.12500 -213.8750 100.12500 -33.875000 112.12500 77.12500 [8,] 60.18750 247.1875 -180.81250 -47.812500 38.18750 -56.81250 [9,] -234.37500 -18.3750 30.62500 -99.375000 -131.37500 139.62500 [10,] 144.43750 10.4375 115.43750 167.437500 163.43750 -289.56250 [11,] 177.37500 199.3750 98.37500 95.375000 74.37500 87.37500 [12,] -272.86017 -172.8602 194.13983 4.139829 -88.86017 -13.86017 [,13] [,14] [,15] [,16] [,17] [1,] 169.8509 174.8509 -24.76403 -6.149065 NA [2,] 33.4375 -39.5625 110.43750 66.437500 NA [3,] -24.8125 22.1875 63.18750 -58.812500 NA [4,] -19.3125 -15.3125 -69.31250 49.687500 NA [5,] -81.2500 126.7500 -104.25000 -56.250000 NA [6,] 180.3750 -145.6250 22.37500 -38.625000 NA [7,] -145.8750 181.1250 -49.87500 47.125000 NA [8,] 118.1875 -142.8125 235.18750 107.187500 NA [9,] 117.6250 116.6250 -79.37500 -8.375000 NA [10,] -269.5625 -51.5625 -203.56250 -37.562500 NA [11,] -258.6250 -35.6250 -86.62500 -90.625000 NA [12,] 183.1398 -131.8602 297.13983 NA 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/1qbea1227790328.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/2tk721227790328.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/3ym7l1227790328.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/4hfhq1227790328.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,] -249.67570 -202.003025 -180.387261 -134.55629 -196.3339972 -246.762261 [2,] -149.90830 -98.151792 -43.030733 -64.61063 -118.4231292 -130.534952 [3,] -66.96264 -8.003025 4.285708 13.97850 0.8389715 5.232444 [4,] 177.64604 68.175239 51.969267 39.30024 79.2150476 118.950795 [5,] 310.78951 146.462183 162.899683 68.43312 327.9529468 212.524683 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -220.63511 -319.262261 -268.322609 -333.69761 -216.29921 -322.74594 [2,] -96.03073 -83.207921 -89.609553 -117.50410 -72.06660 -137.45900 [3,] -14.94690 4.183399 -9.753025 45.06979 14.04208 13.67918 [4,] 86.61274 99.450795 95.187340 135.49060 81.98774 117.63013 [5,] 317.07795 202.202947 321.390447 205.44156 274.27468 414.64968 $n [1] 16 16 16 16 16 16 16 16 16 16 16 16 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -196.34660 -73.70220 -33.23929 -27.06629 -77.22811 -93.31443 -87.09107 [2,] 62.42132 57.69615 41.81071 55.02330 78.90605 103.77931 57.19727 [,8] [,9] [,10] [,11] [,12] [1,] -67.96679 -82.74780 -54.86312 -46.80939 -87.08103 [2,] 76.33359 63.24175 145.00270 74.89354 114.43939 $out [1] 362.6404 -228.6351 238.1475 -237.4476 403.2654 $group [1] 2 3 3 4 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-249.6756966, -149.90830065, -66.96264064, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5nwtc1227790328.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] [,7] [1,] -308.14906 -98.5625 -187.8125 -150.3125 -165.25 -145.625 -290.875 [2,] -122.14906 -46.0625 -55.8125 -72.3125 -75.25 -73.625 -82.375 [3,] -15.45655 5.4375 -16.8125 -17.3125 -49.25 -33.125 51.625 [4,] 125.35094 80.4375 53.1875 47.1875 123.75 79.875 96.625 [5,] 367.85094 266.4375 105.1875 211.6875 231.75 209.375 181.125 [,8] [,9] [,10] [,11] [,12] [1,] -359.8125 -234.375 -289.5625 -258.625 -272.86017 [2,] -118.8125 -89.375 -173.0625 -105.625 -96.36017 [3,] 37.1875 -13.375 22.4375 19.375 -13.86017 [4,] 112.6875 117.125 153.9375 108.875 120.13983 [5,] 247.1875 233.625 299.4375 199.375 341.13983 $n [1] 16 16 16 16 16 16 16 16 16 16 16 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -113.21905 -44.530 -59.8675 -64.515 -127.855 -93.7575 -19.08 -54.255 [2,] 82.30595 55.405 26.2425 29.890 29.355 27.5075 122.33 128.630 [,9] [,10] [,11] [,12] [1,] -94.9425 -106.7275 -65.3525 -102.18227 [2,] 68.1925 151.6025 104.1025 74.46192 $out [1] -336.5625 -298.5625 377.1875 269.6875 $group [1] 2 2 3 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-308.14906462, -122.1490646, -15.456547815, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ws4a1227790328.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] [,7] [1,] -333.6976 -163.58185 -68.528581 -41.97532 -68.17205 -157.03222 -298.62803 [2,] -233.0414 -102.44481 6.002669 31.46218 102.04670 -108.15004 -184.42849 [3,] -188.4668 14.60565 52.690169 187.71218 256.49623 44.38121 -139.31553 [4,] -159.6976 61.73065 88.346419 273.44297 324.67170 197.28746 -88.81553 [5,] -127.7601 196.29315 201.432984 414.64968 403.26545 262.56871 43.43447 [,8] [,9] [,10] [,11] [,12] [,13] [1,] -319.26226 -196.33400 -123.15573 -140.85247 -216.29921 -142.78438 [2,] -213.57476 -143.74025 -23.09323 -44.97747 -88.89296 -64.49594 [3,] -85.98101 -67.85321 30.62552 8.74128 -47.92421 21.72281 [4,] -75.63726 -15.70900 88.28177 83.11628 15.18158 74.16031 [5,] 119.36274 81.54100 149.71927 238.14753 73.26329 205.44156 [,14] [,15] [,16] [,17] [1,] -139.6061127 -267.00438 -13.57612 NA [2,] -0.8801773 -118.91063 24.51763 NA [3,] 28.6198227 -81.59111 32.58013 NA [4,] 108.1823227 14.21437 92.61138 NA [5,] 164.8073227 102.55812 140.73638 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -221.9194 -60.27588 15.13263 77.34312 154.9555 -94.93088 -182.9252 [2,] -155.0143 89.48719 90.24771 298.08124 358.0370 183.69330 -95.7058 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -148.89523 -126.249127 -20.17337 -49.68314 -95.3932980 -41.51924 -21.12432 [2,] -23.06679 -9.457302 81.42440 67.16570 -0.4551124 84.96486 78.36396 [,15] [,16] [,17] [1,] -142.31032 1.522118 NA [2,] -20.87191 63.638151 NA $out [1] -34.26011 -185.96608 275.68083 -322.74594 $group [1] 1 3 10 13 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" NA Warning message: In bxp(list(stats = c(-333.6976091, -233.0413591, -188.4668268, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7urla1227790328.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,] -7.500005e-09 -14.946896 -21.46085904 [2,] -3.750001e-09 -8.878025 -9.06198321 [3,] -1.500001e-09 4.234553 0.03031907 [4,] 3.625000e-09 13.828841 7.49662558 [5,] 6.312499e-09 45.069787 12.01383278 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -4.863789e-09 -6.122201 -7.522174 [2,] 1.863786e-09 14.591307 7.582812 $out [1] -66.96264 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-7.50000461735567e-09, -3.75000114294366e-09, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1qbea1227790328.ps tmp/1qbea1227790328.png") > system("convert tmp/2tk721227790328.ps tmp/2tk721227790328.png") > system("convert tmp/3ym7l1227790328.ps tmp/3ym7l1227790328.png") > system("convert tmp/4hfhq1227790328.ps tmp/4hfhq1227790328.png") > system("convert tmp/5nwtc1227790328.ps tmp/5nwtc1227790328.png") > system("convert tmp/6ws4a1227790328.ps tmp/6ws4a1227790328.png") > system("convert tmp/7urla1227790328.ps tmp/7urla1227790328.png") > > > proc.time() user system elapsed 1.185 0.920 1.838