R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) 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(3004,3080,3017,3114,3057,3032,3127,3050,2910,2671,2638,2672,2654,2568,2467,2419,2363,2291,2560,2527,2370,2310,2231,2367,2346,2286,2249,2226,2108,2131,2387,2358,2284,2312,2293,2576,2665,2749,2926,2886,2893,2944,3060,3045,2894,2955,2954,3243,3120,3074,3034,2981,2876,2835,2978,2881,2768,2722,2630,2753,2771,2652,2584,2501,2449,2445,2620,2579,2460,2434,2392,1037,1212,1232,1174,1158,1140,1118,1212,1207,1186,608,627,626,649,619,612,643,623,649,699,693,659,669,668,693) > par1 = '12' > 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,] 3004 2654 2346 2665 3120 2771 1212 649 NA [2,] 3080 2568 2286 2749 3074 2652 1232 619 NA [3,] 3017 2467 2249 2926 3034 2584 1174 612 NA [4,] 3114 2419 2226 2886 2981 2501 1158 643 NA [5,] 3057 2363 2108 2893 2876 2449 1140 623 NA [6,] 3032 2291 2131 2944 2835 2445 1118 649 NA [7,] 3127 2560 2387 3060 2978 2620 1212 699 NA [8,] 3050 2527 2358 3045 2881 2579 1207 693 NA [9,] 2910 2370 2284 2894 2768 2460 1186 659 NA [10,] 2671 2310 2312 2955 2722 2434 608 669 NA [11,] 2638 2231 2293 2954 2630 2392 627 668 NA [12,] 2672 2367 2576 3243 2753 1037 626 693 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 76 -86 -60 84 -46 -119 20 -30 NA [2,] -63 -101 -37 177 -40 -68 -58 -7 NA [3,] 97 -48 -23 -40 -53 -83 -16 31 NA [4,] -57 -56 -118 7 -105 -52 -18 -20 NA [5,] -25 -72 23 51 -41 -4 -22 26 NA [6,] 95 269 256 116 143 175 94 50 NA [7,] -77 -33 -29 -15 -97 -41 -5 -6 NA [8,] -140 -157 -74 -151 -113 -119 -21 -34 NA [9,] -239 -60 28 61 -46 -26 -578 10 NA [10,] -33 -79 -19 -1 -92 -42 19 -1 NA [11,] 34 136 283 289 123 -1355 -1 25 NA [12,] -18 -21 89 -123 18 175 23 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/wessaorg/rcomp/tmp/1ana91419845393.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2junc1419845393.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/32ryl1419845393.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4hpim1419845393.ps",horizontal=F,onefile=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,] 649.0 619.0 612.0 643.0 623.0 649.0 699.0 693.0 659 608.0 627.0 [2,] 1779.0 1759.0 1711.5 1692.0 1624.0 1624.5 1799.5 1782.5 1735 1489.5 1449.5 [3,] 2659.5 2610.0 2525.5 2460.0 2406.0 2368.0 2590.0 2553.0 2415 2373.0 2342.5 [4,] 2887.5 2911.5 2971.5 2933.5 2884.5 2889.5 3019.0 2963.0 2831 2696.5 2634.0 [5,] 3120.0 3080.0 3034.0 3114.0 3057.0 3032.0 3127.0 3050.0 2910 2955.0 2954.0 [,12] [1,] 626.0 [2,] 865.0 [3,] 2471.5 [4,] 2712.5 [5,] 3243.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,] 2040.276 1966.197 1821.646 1766.48 1701.867 1661.353 1908.77 1893.556 [2,] 3278.724 3253.803 3229.354 3153.52 3110.133 3074.647 3271.23 3212.444 [,9] [,10] [,11] [,12] [1,] 1802.759 1698.752 1680.821 1439.46 [2,] 3027.241 3047.248 3004.179 3503.54 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(649, 1779, 2659.5, 2887.5, 3120, 619, 1759, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5il941419845393.ps",horizontal=F,onefile=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] [,11] [,12] [1,] -119 -101.0 -83.0 -118 -72.0 50.0 -97.0 -157.0 -239.0 -92.0 -1.0 -123.0 [2,] -73 -65.5 -50.5 -81 -33.0 94.5 -59.0 -145.5 -149.5 -60.5 12.0 -19.5 [3,] -38 -49.0 -31.5 -54 -13.0 129.5 -31.0 -116.0 -36.0 -26.0 78.5 18.0 [4,] 48 -22.0 7.5 -19 24.5 215.5 -10.5 -54.0 19.0 -1.0 209.5 56.0 [5,] 84 -7.0 31.0 7 51.0 269.0 -5.0 -21.0 61.0 19.0 289.0 89.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -105.59234 -73.29972 -63.8996327 -88.63409 -45.12033 61.90766 -58.092796 [2,] 29.59234 -24.70028 0.8996327 -19.36591 19.12033 197.09234 -3.907204 [,8] [,9] [,10] [,11] [,12] [1,] -167.11321 -130.12652 -59.237554 -31.82634 -27.08738 [2,] -64.88679 58.12652 7.237554 188.82634 63.08738 $out [1] 177 97 -578 -1355 175 $group [1] 2 3 9 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-119, -73, -38, 48, 84, -101, -65.5, -49, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6utce1419845393.ps",horizontal=F,onefile=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,] 2638.0 2231.0 2108.0 2749.0 2630.0 2392.0 608.0 612 NA [2,] 2791.0 2336.5 2237.5 2889.5 2760.5 2439.5 872.5 633 NA [3,] 3024.5 2394.5 2289.5 2935.0 2878.5 2480.5 1166.0 654 NA [4,] 3068.5 2543.5 2352.0 3000.0 3007.5 2602.0 1209.5 681 NA [5,] 3127.0 2654.0 2387.0 3060.0 3120.0 2771.0 1232.0 699 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2897.93 2300.086 2237.276 2884.6 2765.842 2406.383 1012.292 632.1069 NA [2,] 3151.07 2488.914 2341.724 2985.4 2991.158 2554.617 1319.708 675.8931 NA $out [1] 2576 2665 3243 1037 $group [1] 3 4 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(2638, 2791, 3024.5, 3068.5, 3127, 2231, 2336.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7bj5e1419845393.ps",horizontal=F,onefile=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,] 1995.875 2342.50 2236.125 [2,] 2132.875 2389.50 2330.875 [3,] 2216.188 2465.75 2415.438 [4,] 2287.500 2571.50 2455.188 [5,] 2330.375 2659.50 2545.875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2145.662 2382.739 2358.738 [2,] 2286.713 2548.761 2472.137 $out [1] 1821.625 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1995.875, 2132.875, 2216.1875, 2287.5, 2330.375, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ana91419845393.ps tmp/1ana91419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/2junc1419845393.ps tmp/2junc1419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/32ryl1419845393.ps tmp/32ryl1419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/4hpim1419845393.ps tmp/4hpim1419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/5il941419845393.ps tmp/5il941419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/6utce1419845393.ps tmp/6utce1419845393.png",intern=TRUE)) character(0) > try(system("convert tmp/7bj5e1419845393.ps tmp/7bj5e1419845393.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.685 0.512 3.230