R version 2.9.0 (2009-04-17) 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(745.52,1962.64,2092.88,2034.73,568.28,1482.67,1282.26,1534.16,1621.42,1465.32,1373.15,1386.28,756.10,1359.02,1600.69,1299.97,161.31,622.74,552.60,570.36,806.43,1212.24,1806.54,1433.47,543.97,1453.00,1795.89,1608.99,484.29,1541.01,1094.14,1412.93,1612.35,1309.29,1626.17,874.45,794.25,1416.23,1575.25,1255.68,-787.86,1038.52,934.94,581.78,1069.98,1104.77,1152.28,880.48,169.49,953.93,1009.01,297.91,259.98,1164.90,924.60,840.77,1042.04,1026.08,637.53,1338.82) > 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] 60 > (np <- floor(n / par1)) [1] 5 > 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] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 745.52 756.10 543.97 794.25 169.49 NA [2,] 1962.64 1359.02 1453.00 1416.23 953.93 NA [3,] 2092.88 1600.69 1795.89 1575.25 1009.01 NA [4,] 2034.73 1299.97 1608.99 1255.68 297.91 NA [5,] 568.28 161.31 484.29 -787.86 259.98 NA [6,] 1482.67 622.74 1541.01 1038.52 1164.90 NA [7,] 1282.26 552.60 1094.14 934.94 924.60 NA [8,] 1534.16 570.36 1412.93 581.78 840.77 NA [9,] 1621.42 806.43 1612.35 1069.98 1042.04 NA [10,] 1465.32 1212.24 1309.29 1104.77 1026.08 NA [11,] 1373.15 1806.54 1626.17 1152.28 637.53 NA [12,] 1386.28 1433.47 874.45 880.48 1338.82 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1217.12 602.92 909.03 621.98 784.44 NA [2,] 130.24 241.67 342.89 159.02 55.08 NA [3,] -58.15 -300.72 -186.90 -319.57 -711.10 NA [4,] -1466.45 -1138.66 -1124.70 -2043.54 -37.93 NA [5,] 914.39 461.43 1056.72 1826.38 904.92 NA [6,] -200.41 -70.14 -446.87 -103.58 -240.30 NA [7,] 251.90 17.76 318.79 -353.16 -83.83 NA [8,] 87.26 236.07 199.42 488.20 201.27 NA [9,] -156.10 405.81 -303.06 34.79 -15.96 NA [10,] -92.17 594.30 316.88 47.51 -388.55 NA [11,] 13.13 -373.07 -751.72 -271.80 701.29 NA [12,] -630.18 -889.50 -80.20 -710.99 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/1xh0d1258294351.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/2zvav1258294351.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/3bo051258294351.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/4cdew1258294351.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] [1,] 543.97 1359.02 1575.25 1255.68 161.31 622.74 924.60 570.36 806.43 [2,] 543.97 1359.02 1575.25 1255.68 161.31 1038.52 924.60 581.78 1042.04 [3,] 745.52 1416.23 1600.69 1299.97 259.98 1164.90 934.94 840.77 1069.98 [4,] 756.10 1453.00 1795.89 1608.99 484.29 1482.67 1094.14 1412.93 1612.35 [5,] 794.25 1453.00 2092.88 2034.73 568.28 1541.01 1282.26 1534.16 1621.42 [,10] [,11] [,12] [1,] 1026.08 637.53 874.45 [2,] 1104.77 1152.28 880.48 [3,] 1212.24 1373.15 1338.82 [4,] 1309.29 1626.17 1386.28 [5,] 1465.32 1806.54 1433.47 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 595.6295 1349.824 1444.786 1050.322 31.76315 851.0647 815.1435 [2,] 895.4105 1482.636 1756.594 1549.618 488.19685 1478.7353 1054.7365 [,8] [,9] [,10] [,11] [,12] [1,] 253.4815 667.0004 1067.727 1038.301 981.423 [2,] 1428.0585 1472.9596 1356.753 1707.999 1696.217 $out [1] 169.49 1962.64 953.93 1009.01 297.91 -787.86 552.60 $group [1] 1 2 2 3 4 5 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(543.97, 543.97, 745.52, 756.1, 794.25, 1359.02, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/53sq01258294352.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] [1,] 602.92 55.08 -319.57 -1466.45 904.92 -240.30 -353.16 199.42 -303.06 [2,] 621.98 130.24 -319.57 -1466.45 904.92 -240.30 -83.83 199.42 -156.10 [3,] 784.44 159.02 -300.72 -1138.66 914.39 -200.41 17.76 201.27 -15.96 [4,] 909.03 241.67 -186.90 -1124.70 1056.72 -103.58 251.90 236.07 34.79 [5,] 1217.12 342.89 -58.15 -1124.70 1056.72 -70.14 318.79 236.07 34.79 [,10] [,11] [,12] [1,] -388.55 -751.72 -889.500 [2,] -92.17 -373.07 -800.245 [3,] 47.51 -271.80 -670.585 [4,] 316.88 13.13 -355.190 [5,] 594.30 13.13 -80.200 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 581.6112 80.28384 -394.4643 -1380.1397 807.1285 -297.016 -219.4660 [2,] 987.2688 237.75616 -206.9757 -897.1803 1021.6515 -103.804 254.9860 [,8] [,9] [,10] [,11] [,12] [1,] 175.3732 -150.8424 -241.5237 -544.687947 -1022.1785 [2,] 227.1668 118.9224 336.5437 1.087947 -318.9916 $out [1] -711.10 -2043.54 -37.93 461.43 1826.38 -446.87 87.26 488.20 [9] 405.81 701.29 $group [1] 3 4 4 5 5 6 8 8 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(602.92, 621.98, 784.44, 909.03, 1217.12, 55.08, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6f8fc1258294352.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] [1,] 745.520 161.310 484.290 581.780 169.490 NA [2,] 1327.705 596.550 984.295 837.365 467.720 NA [3,] 1473.995 1009.335 1432.965 1054.250 939.265 NA [4,] 1792.030 1396.245 1610.670 1203.980 1034.060 NA [5,] 2092.880 1806.540 1795.890 1575.250 1338.820 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1262.213 644.5887 1147.271 887.0344 680.9535 NA [2,] 1685.777 1374.0813 1718.659 1221.4656 1197.5765 NA $out [1] 568.28 -787.86 $group [1] 1 4 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(745.52, 1327.705, 1473.995, 1792.03, 2092.88, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7w9wl1258294352.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,] 601.866 259.980 650.035 [2,] 972.854 887.855 1003.363 [3,] 1203.120 1188.570 1233.812 [4,] 1309.295 1355.985 1397.617 [5,] 1614.744 1600.690 1685.570 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1049.667 975.0528 1053.990 [2,] 1356.573 1402.0872 1413.635 $out [1] 137.2 322.8 $group [1] 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(601.866, 972.854, 1203.12, 1309.295, 1614.744, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1xh0d1258294351.ps tmp/1xh0d1258294351.png") > system("convert tmp/2zvav1258294351.ps tmp/2zvav1258294351.png") > system("convert tmp/3bo051258294351.ps tmp/3bo051258294351.png") > system("convert tmp/4cdew1258294351.ps tmp/4cdew1258294351.png") > system("convert tmp/53sq01258294352.ps tmp/53sq01258294352.png") > system("convert tmp/6f8fc1258294352.ps tmp/6f8fc1258294352.png") > system("convert tmp/7w9wl1258294352.ps tmp/7w9wl1258294352.png") > > > proc.time() user system elapsed 1.188 0.931 1.531