R version 2.10.1 (2009-12-14) 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(2851,2672,2755,2721,2946,3036,2282,2212,2922,4301,5764,7132,2541,2475,3031,3266,3776,3230,3028,1759,3595,4474,6838,8357,3113,3006,4047,3523,3937,3986,3260,1573,3528,5211,7614,9254,5375,3088,3718,4514,4520,4539,3663,1643,4739,5428,8314,10651,3633,4292,4154,4121,4647,4753,3965,1723,5048,6922,9858,11331,4016,3975,4510,4276,4968,4677,3523,1821,5222,6873,10803,13916,2639,2899,3370,3740,2927,3986,4217,1738,5221,6424,9842,13076,3934,3162,4286,4676,5010,4874,4633,1659,5951,6981,9851,12670) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2851 2541 3113 5375 3633 4016 2639 3934 NA [2,] 2672 2475 3006 3088 4292 3975 2899 3162 NA [3,] 2755 3031 4047 3718 4154 4510 3370 4286 NA [4,] 2721 3266 3523 4514 4121 4276 3740 4676 NA [5,] 2946 3776 3937 4520 4647 4968 2927 5010 NA [6,] 3036 3230 3986 4539 4753 4677 3986 4874 NA [7,] 2282 3028 3260 3663 3965 3523 4217 4633 NA [8,] 2212 1759 1573 1643 1723 1821 1738 1659 NA [9,] 2922 3595 3528 4739 5048 5222 5221 5951 NA [10,] 4301 4474 5211 5428 6922 6873 6424 6981 NA [11,] 5764 6838 7614 8314 9858 10803 9842 9851 NA [12,] 7132 8357 9254 10651 11331 13916 13076 12670 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -179 -66 -107 -2287 659 -41 260 -772 NA [2,] 83 556 1041 630 -138 535 471 1124 NA [3,] -34 235 -524 796 -33 -234 370 390 NA [4,] 225 510 414 6 526 692 -813 334 NA [5,] 90 -546 49 19 106 -291 1059 -136 NA [6,] -754 -202 -726 -876 -788 -1154 231 -241 NA [7,] -70 -1269 -1687 -2020 -2242 -1702 -2479 -2974 NA [8,] 710 1836 1955 3096 3325 3401 3483 4292 NA [9,] 1379 879 1683 689 1874 1651 1203 1030 NA [10,] 1463 2364 2403 2886 2936 3930 3418 2870 NA [11,] 1368 1519 1640 2337 1473 3113 3234 2819 NA [12,] -4591 -5244 -3879 -7018 -7315 -11277 -9142 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/rcomp/tmp/18fq01274005322.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/rcomp/tmp/28fq01274005322.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/rcomp/tmp/38fq01274005322.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/rcomp/tmp/4j6p31274005322.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] [,7] [,8] [,9] [,10] [,11] [1,] 2541 2475.0 2755.0 2721.0 2927.0 3036.0 2282 1573.0 2922.0 4301.0 5764.0 [2,] 2745 2785.5 3200.5 3394.5 3361.0 3608.0 3144 1651.0 3561.5 4842.5 7226.0 [3,] 3373 3047.0 3882.5 3930.5 4228.5 4262.5 3593 1730.5 4893.5 5926.0 9078.0 [4,] 3975 3568.5 4220.0 4395.0 4807.5 4715.0 4091 1790.0 5221.5 6897.5 9854.5 [5,] 5375 4292.0 4510.0 4676.0 5010.0 4874.0 4633 1821.0 5951.0 6981.0 10803.0 [,12] [1,] 7132.0 [2,] 8805.5 [3,] 10991.0 [4,] 12873.0 [5,] 13916.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2685.904 2609.605 3312.993 3371.606 3420.464 3644.114 3063.992 1652.853 [2,] 4060.096 3484.395 4452.007 4489.394 5036.536 4880.886 4122.008 1808.147 [,9] [,10] [,11] [,12] [1,] 3966.2 4778.047 7609.682 8718.836 [2,] 5820.8 7073.953 10546.318 13263.164 $out [1] 2212 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2541, 2745, 3373, 3975, 5375, 2475, 2785.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5j6p31274005322.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] [,8] [,9] [,10] [1,] -772.0 -138.0 -524 6.0 -546.0 -1154.0 -2974.0 710.0 689.0 1463.0 [2,] -475.5 277.0 -134 115.5 -213.5 -832.0 -2360.5 1895.5 954.5 2383.5 [3,] -86.5 545.5 101 374.0 34.0 -740.0 -1861.0 3210.5 1291.0 2878.0 [4,] 109.5 835.5 380 518.0 98.0 -221.5 -1478.0 3442.0 1667.0 3177.0 [5,] 659.0 1124.0 796 692.0 106.0 231.0 -1269.0 4292.0 1874.0 3930.0 [,11] [,12] [1,] 1368.0 -11277.0 [2,] 1496.0 -8228.5 [3,] 1988.5 -7018.0 [4,] 2966.0 -4917.5 [5,] 3234.0 -3879.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -413.2894 233.5139 -186.1278 149.1577 -140.0084 -1081.0341 -2353.977 [2,] 240.2894 857.4861 388.1278 598.8423 208.0084 -398.9659 -1368.023 [,8] [,9] [,10] [,11] [,12] [1,] 2346.603 892.9873 2434.740 1167.337 -8995.276 [2,] 4074.397 1689.0127 3321.260 2809.663 -5040.724 $out [1] -2287 -813 1059 -70 $group [1] 1 4 5 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-772, -475.5, -86.5, 109.5, 659, -138, 277, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6bxoo1274005322.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] [,8] [,9] [1,] 2212.0 1759.0 1573.0 1643.0 1723.0 1821.0 1738.0 1659 NA [2,] 2696.5 2784.5 3186.5 3690.5 4043.0 3995.5 2913.0 4110 NA [3,] 2886.5 3248.0 3732.5 4529.5 4469.5 4593.5 3863.0 4775 NA [4,] 3668.5 4125.0 4629.0 5401.5 5985.0 6047.5 5822.5 6466 NA [5,] 4301.0 4474.0 5211.0 5428.0 6922.0 6873.0 9842.0 9851 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2443.164 2636.589 3074.566 3749.101 3583.741 3657.569 2535.958 3700.413 [2,] 3329.836 3859.411 4390.434 5309.899 5355.259 5529.431 5190.042 5849.587 [,9] [1,] NA [2,] NA $out [1] 5764 7132 6838 8357 7614 9254 8314 10651 9858 11331 10803 13916 [13] 13076 12670 $group [1] 1 1 2 2 3 3 4 4 5 5 6 6 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(2212, 2696.5, 2886.5, 3668.5, 4301, 1759, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7bxoo1274005322.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,] 1766.000 1730.50 1714.750 [2,] 3542.062 3483.00 3495.625 [3,] 3973.000 4079.50 4022.500 [4,] 5177.500 5409.75 5177.875 [5,] 5826.750 5926.00 5956.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 3227.066 3200.696 3255.214 [2,] 4718.934 4958.304 4789.786 $out [1] 8610.500 10798.375 9078.000 10991.000 8636.375 10900.625 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1766, 3542.0625, 3973, 5177.5, 5826.75, 1730.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/18fq01274005322.ps tmp/18fq01274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/28fq01274005322.ps tmp/28fq01274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/38fq01274005322.ps tmp/38fq01274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/4j6p31274005322.ps tmp/4j6p31274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/5j6p31274005322.ps tmp/5j6p31274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/6bxoo1274005322.ps tmp/6bxoo1274005322.png",intern=TRUE)) character(0) > try(system("convert tmp/7bxoo1274005322.ps tmp/7bxoo1274005322.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.620 1.550 2.199