R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-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(679,687,638,628,604,713,712,693,697,555,486,470,465,426,384,379,381,380,351,346,339,336,333,324,324,321,304,343,407,389,361,353,361,387,692,704,742,721,843,847,945,946,946,945,1082,1075,820,832,851,1090,1203,1239,1535,1527,1480,1452,1383,1381,1429,1376,1602,1597,2003,1958,1997,1986,2129,2115,2297,2250,2309,2648,2627,2711,2732,2825,2932,2910,2969,2999,2965,2846,2847,2751) > par1 = '52' > par1 <- as.numeric(par1) > (n <- length(x)) [1] 84 > (np <- floor(n / par1)) [1] 1 > 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] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 > arr [,1] [,2] [1,] 679 1535 [2,] 687 1527 [3,] 638 1480 [4,] 628 1452 [5,] 604 1383 [6,] 713 1381 [7,] 712 1429 [8,] 693 1376 [9,] 697 1602 [10,] 555 1597 [11,] 486 2003 [12,] 470 1958 [13,] 465 1997 [14,] 426 1986 [15,] 384 2129 [16,] 379 2115 [17,] 381 2297 [18,] 380 2250 [19,] 351 2309 [20,] 346 2648 [21,] 339 2627 [22,] 336 2711 [23,] 333 2732 [24,] 324 2825 [25,] 324 2932 [26,] 321 2910 [27,] 304 2969 [28,] 343 2999 [29,] 407 2965 [30,] 389 2846 [31,] 361 2847 [32,] 353 2751 [33,] 361 NA [34,] 387 NA [35,] 692 NA [36,] 704 NA [37,] 742 NA [38,] 721 NA [39,] 843 NA [40,] 847 NA [41,] 945 NA [42,] 946 NA [43,] 946 NA [44,] 945 NA [45,] 1082 NA [46,] 1075 NA [47,] 820 NA [48,] 832 NA [49,] 851 NA [50,] 1090 NA [51,] 1203 NA [52,] 1239 NA > darr [,1] [,2] [1,] 8 -8 [2,] -49 -47 [3,] -10 -28 [4,] -24 -69 [5,] 109 -2 [6,] -1 48 [7,] -19 -53 [8,] 4 226 [9,] -142 -5 [10,] -69 406 [11,] -16 -45 [12,] -5 39 [13,] -39 -11 [14,] -42 143 [15,] -5 -14 [16,] 2 182 [17,] -1 -47 [18,] -29 59 [19,] -5 339 [20,] -7 -21 [21,] -3 84 [22,] -3 21 [23,] -9 93 [24,] 0 107 [25,] -3 -22 [26,] -17 59 [27,] 39 30 [28,] 64 -34 [29,] -18 -119 [30,] -28 1 [31,] -8 -96 [32,] 8 NA [33,] 26 NA [34,] 305 NA [35,] 12 NA [36,] 38 NA [37,] -21 NA [38,] 122 NA [39,] 4 NA [40,] 98 NA [41,] 1 NA [42,] 0 NA [43,] -1 NA [44,] 137 NA [45,] -7 NA [46,] -255 NA [47,] 12 NA [48,] 19 NA [49,] 239 NA [50,] 113 NA [51,] 36 NA [52,] 296 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/fisher/rcomp/tmp/1o12m1382176621.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/fisher/rcomp/tmp/200ua1382176621.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/fisher/rcomp/tmp/36w1l1382176621.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/fisher/rcomp/tmp/4fx161382176621.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] [,12] [1,] 679 687 638 628 604.0 713 712.0 693.0 697.0 555 486.0 470 [2,] 679 687 638 628 604.0 713 712.0 693.0 697.0 555 486.0 470 [3,] 1107 1107 1059 1040 993.5 1047 1070.5 1034.5 1149.5 1076 1244.5 1214 [4,] 1535 1527 1480 1452 1383.0 1381 1429.0 1376.0 1602.0 1597 2003.0 1958 [5,] 1535 1527 1480 1452 1383.0 1381 1429.0 1376.0 1602.0 1597 2003.0 1958 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 465 426 384.0 379 381 380 351 346 339 336.0 333.0 [2,] 465 426 384.0 379 381 380 351 346 339 336.0 333.0 [3,] 1231 1206 1256.5 1247 1339 1315 1330 1497 1483 1523.5 1532.5 [4,] 1997 1986 2129.0 2115 2297 2250 2309 2648 2627 2711.0 2732.0 [5,] 1997 1986 2129.0 2115 2297 2250 2309 2648 2627 2711.0 2732.0 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [1,] 324.0 324 321.0 304.0 343 407 389.0 361 353 361 387 [2,] 324.0 324 321.0 304.0 343 407 389.0 361 353 361 387 [3,] 1574.5 1628 1615.5 1636.5 1671 1686 1617.5 1604 1552 361 387 [4,] 2825.0 2932 2910.0 2969.0 2999 2965 2846.0 2847 2751 361 387 [5,] 2825.0 2932 2910.0 2969.0 2999 2965 2846.0 2847 2751 361 387 [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] 692 704 742 721 843 847 945 946 946 945 1082 1075 [2,] 692 704 742 721 843 847 945 946 946 945 1082 1075 [3,] 692 704 742 721 843 847 945 946 946 945 1082 1075 [4,] 692 704 742 721 843 847 945 946 946 945 1082 1075 [5,] 692 704 742 721 843 847 945 946 946 945 1082 1075 [,47] [,48] [,49] [,50] [,51] [,52] [1,] 820 832 851 1090 1203 1239 [2,] 820 832 851 1090 1203 1239 [3,] 820 832 851 1090 1203 1239 [4,] 820 832 851 1090 1203 1239 [5,] 820 832 851 1090 1203 1239 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 150.6522 168.5279 118.2934 119.4035 123.1788 300.6912 269.447 [2,] 2063.3478 2045.4721 1999.7066 1960.5965 1863.8212 1793.3088 1871.553 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 271.4328 138.408 -88.15232 -450.336 -448.4363 -480.5944 -536.8768 [2,] 1797.5672 2160.592 2240.15232 2939.336 2876.4363 2942.5944 2948.8768 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -693.0641 -692.509 -801.6102 -774.2177 -857.5338 -1074.861 -1073.219 [2,] 3206.0641 3186.509 3479.6102 3404.2177 3517.5338 4068.861 4039.219 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -1129.918 -1147.732 -1219.689 -1285.732 -1277.005 -1340.915 -1296.359 [2,] 4176.918 4212.732 4368.689 4541.732 4508.005 4613.915 4638.359 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] -1171.871 -1127.531 -1173.431 -1127.114 361 387 692 704 742 [2,] 4543.871 4362.531 4381.431 4231.114 361 387 692 704 742 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 721 843 847 945 946 946 945 1082 1075 820 832 851 [2,] 721 843 847 945 946 946 945 1082 1075 820 832 851 [,50] [,51] [,52] [1,] 1090 1203 1239 [2,] 1090 1203 1239 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" Warning message: In bxp(list(stats = c(679, 679, 1107, 1535, 1535, 687, 687, 1107, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/5aary1382176621.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] [,13] [1,] -8 -49 -28 -69.0 -2.0 -1.0 -53 4 -142.0 -69.0 -45.0 -5 -39 [2,] -8 -49 -28 -69.0 -2.0 -1.0 -53 4 -142.0 -69.0 -45.0 -5 -39 [3,] 0 -48 -19 -46.5 53.5 23.5 -36 115 -73.5 168.5 -30.5 17 -25 [4,] 8 -47 -10 -24.0 109.0 48.0 -19 226 -5.0 406.0 -16.0 39 -11 [5,] 8 -47 -10 -24.0 109.0 48.0 -19 226 -5.0 406.0 -16.0 39 -11 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -42.0 -14.0 2 -47 -29 -5 -21 -3.0 -3 -9 0.0 -22.0 [2,] -42.0 -14.0 2 -47 -29 -5 -21 -3.0 -3 -9 0.0 -22.0 [3,] 50.5 -9.5 92 -24 15 167 -14 40.5 9 42 53.5 -12.5 [4,] 143.0 -5.0 182 -1 59 339 -7 84.0 21 93 107.0 -3.0 [5,] 143.0 -5.0 182 -1 59 339 -7 84.0 21 93 107.0 -3.0 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] -17 30.0 -34 -119.0 -28.0 -96 8 26 305 12 38 -21 [2,] -17 30.0 -34 -119.0 -28.0 -96 8 26 305 12 38 -21 [3,] 21 34.5 15 -68.5 -13.5 -52 8 26 305 12 38 -21 [4,] 59 39.0 64 -18.0 1.0 -8 8 26 305 12 38 -21 [5,] 59 39.0 64 -18.0 1.0 -8 8 26 305 12 38 -21 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 122 4 98 1 0 -1 137 -7 -255 12 19 239 [2,] 122 4 98 1 0 -1 137 -7 -255 12 19 239 [3,] 122 4 98 1 0 -1 137 -7 -255 12 19 239 [4,] 122 4 98 1 0 -1 137 -7 -255 12 19 239 [5,] 122 4 98 1 0 -1 137 -7 -255 12 19 239 [,50] [,51] [,52] [1,] 113 36 296 [2,] 113 36 296 [3,] 113 36 296 [4,] 113 36 296 [5,] 113 36 296 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -17.87566 -50.23446 -39.110117 -96.775292 -70.51239 -31.24421 -73.985776 [2,] 17.87566 -45.76554 1.110117 3.775292 177.51239 78.24421 1.985776 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -133.0248 -226.56033 -362.1836 -62.899633 -32.15806 -56.282404 -156.1873 [2,] 363.0248 79.56033 699.1836 1.899633 66.15806 6.282404 257.1873 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -19.5550584 -109.1012 -75.39252 -83.31613 -217.3267 -29.641202 -56.6989 [2,] 0.5550584 293.1012 27.39252 113.31613 551.3267 1.641202 137.6989 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -17.81349 -71.95733 -66.04347 -33.727346 -63.90938 24.44494 -94.48841 [2,] 35.81349 155.95733 173.04347 8.727346 105.90938 44.55506 124.48841 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -181.3401 -45.89963 -150.31613 8 26 305 12 38 -21 122 [2,] 44.3401 18.89963 46.31613 8 26 305 12 38 -21 122 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] 4 98 1 0 -1 137 -7 -255 12 19 239 113 [2,] 4 98 1 0 -1 137 -7 -255 12 19 239 113 [,51] [,52] [1,] 36 296 [2,] 36 296 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" Warning message: In bxp(list(stats = c(-8, -8, 0, 8, 8, -49, -49, -48, -47, -47, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/6zr3n1382176621.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] [1,] 304 1376.0 [2,] 370 1566.0 [3,] 616 2189.5 [4,] 826 2788.0 [5,] 1239 2999.0 $n [1] 52 32 $conf [,1] [,2] [1,] 516.0874 1848.187 [2,] 715.9126 2530.813 $out numeric(0) $group numeric(0) $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/7zj371382176621.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,] 361.0 361.0 361.0 [2,] 945.5 945.5 945.5 [3,] 1107.0 1107.0 1107.0 [4,] 1411.0 1411.0 1411.0 [5,] 1686.0 1686.0 1686.0 $n [1] 52 52 52 $conf [,1] [,2] [,3] [1,] 1005.006 1005.006 1005.006 [2,] 1208.994 1208.994 1208.994 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/1o12m1382176621.ps tmp/1o12m1382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/200ua1382176621.ps tmp/200ua1382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/36w1l1382176621.ps tmp/36w1l1382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/4fx161382176621.ps tmp/4fx161382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/5aary1382176621.ps tmp/5aary1382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/6zr3n1382176621.ps tmp/6zr3n1382176621.png",intern=TRUE)) character(0) > try(system("convert tmp/7zj371382176621.ps tmp/7zj371382176621.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.175 0.538 3.706