R version 2.8.1 (2008-12-22) 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(464,675,703,887,1139,1077,1318,1260,1120,963,996,960,530,883,894,1045,1199,1287,1565,1577,1076,918,1008,1063,544,635,804,980,1018,1064,1404,1286,1104,999,996,1015,615,722,832,977,1270,1437,1520,1708,1151,934,1159,1209,699,830,996,1124,1458,1270,1753,2258,1208,1241,1265,1828,809,997,1164,1205,1538,1513,1378,2083,1357,1536,1526,1376) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 464 530 544 615 699 809 NA [2,] 675 883 635 722 830 997 NA [3,] 703 894 804 832 996 1164 NA [4,] 887 1045 980 977 1124 1205 NA [5,] 1139 1199 1018 1270 1458 1538 NA [6,] 1077 1287 1064 1437 1270 1513 NA [7,] 1318 1565 1404 1520 1753 1378 NA [8,] 1260 1577 1286 1708 2258 2083 NA [9,] 1120 1076 1104 1151 1208 1357 NA [10,] 963 918 999 934 1241 1536 NA [11,] 996 1008 996 1159 1265 1526 NA [12,] 960 1063 1015 1209 1828 1376 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 211 353 91 107 131 188 NA [2,] 28 11 169 110 166 167 NA [3,] 184 151 176 145 128 41 NA [4,] 252 154 38 293 334 333 NA [5,] -62 88 46 167 -188 -25 NA [6,] 241 278 340 83 483 -135 NA [7,] -58 12 -118 188 505 705 NA [8,] -140 -501 -182 -557 -1050 -726 NA [9,] -157 -158 -105 -217 33 179 NA [10,] 33 90 -3 225 24 -10 NA [11,] -36 55 19 50 563 -150 NA [12,] -430 -519 -400 -510 -1019 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/1wy5q1241095215.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/2n2wr1241095215.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/3bn9g1241095215.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/4rzrc1241095215.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] [,12] [1,] 464.0 635 703 887.0 1018.0 1064.0 1318 1260.0 1076.0 918 996.0 960 [2,] 530.0 675 804 977.0 1139.0 1077.0 1378 1286.0 1104.0 934 996.0 1015 [3,] 579.5 776 863 1012.5 1234.5 1278.5 1462 1642.5 1135.5 981 1083.5 1136 [4,] 699.0 883 996 1124.0 1458.0 1437.0 1565 2083.0 1208.0 1241 1265.0 1376 [5,] 809.0 997 1164 1205.0 1538.0 1513.0 1753 2258.0 1357.0 1536 1526.0 1828 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 470.4895 641.8333 739.1538 917.6803 1028.735 1046.288 1341.379 1128.409 [2,] 688.5105 910.1667 986.8462 1107.3197 1440.265 1510.712 1582.621 2156.591 [,9] [,10] [,11] [,12] [1,] 1068.417 782.9751 909.9863 903.1433 [2,] 1202.583 1179.0249 1257.0137 1368.8567 $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(464, 530, 579.5, 699, 809, 635, 675, 776, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/55bqg1241095215.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] [,11] [,12] [1,] 91.0 11 128 38.0 -188.0 -135.0 -118 -1050 -217 -10.0 -150.0 -519 [2,] 107.0 28 128 154.0 -62.0 83.0 -58 -726 -158 -3.0 -36.0 -519 [3,] 159.5 138 148 272.5 10.5 259.5 100 -529 -131 28.5 34.5 -510 [4,] 211.0 167 176 333.0 88.0 340.0 505 -182 33 90.0 55.0 -430 [5,] 353.0 169 184 334.0 167.0 483.0 705 -140 179 225.0 55.0 -400 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 92.41664 48.34051 117.0384 157.0392 -86.25484 93.7267 -263.1532 [2,] 226.58336 227.65949 178.9616 387.9608 107.25484 425.2733 463.1532 [,8] [,9] [,10] [,11] [,12] [1,] -879.8976 -254.201169 -31.48800 -24.19794 -572.8872 [2,] -178.1024 -7.798831 88.48800 93.19794 -447.1128 $out [1] 41 563 -1019 $group [1] 3 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(91, 107, 159.5, 211, 353, 11, 28, 138, 167, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6d1jj1241095215.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,] 464.0 530 635 615.0 699.0 809.0 NA [2,] 795.0 906 892 883.0 1060.0 1184.5 NA [3,] 979.5 1054 1007 1155.0 1253.0 1377.0 NA [4,] 1129.5 1243 1084 1353.5 1605.5 1531.0 NA [5,] 1318.0 1577 1286 1708.0 2258.0 1538.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 826.9323 900.292 919.4275 940.4018 1004.194 1218.959 NA [2,] 1132.0677 1207.708 1094.5725 1369.5982 1501.806 1535.041 NA $out [1] 544 1404 2083 $group [1] 3 3 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(464, 795, 979.5, 1129.5, 1318, 530, 906, 1054, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/79cwp1241095215.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,] 610.1667 579.5 605.750 [2,] 967.5833 922.0 965.875 [3,] 1163.8333 1109.5 1134.812 [4,] 1272.5000 1256.5 1272.438 [5,] 1695.3333 1642.5 1674.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1024.759 956.9323 994.9873 [2,] 1302.908 1262.0677 1274.6377 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(610.166666666667, 967.583333333333, 1163.83333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1wy5q1241095215.ps tmp/1wy5q1241095215.png") > system("convert tmp/2n2wr1241095215.ps tmp/2n2wr1241095215.png") > system("convert tmp/3bn9g1241095215.ps tmp/3bn9g1241095215.png") > system("convert tmp/4rzrc1241095215.ps tmp/4rzrc1241095215.png") > system("convert tmp/55bqg1241095215.ps tmp/55bqg1241095215.png") > system("convert tmp/6d1jj1241095215.ps tmp/6d1jj1241095215.png") > system("convert tmp/79cwp1241095215.ps tmp/79cwp1241095215.png") > > > proc.time() user system elapsed 1.185 0.929 1.476