R version 2.8.0 (2008-10-20) 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(14532.2,15167,16071.1,14827.5,15082,14772.7,16083,14272.5,15223.3,14897.3,13062.6,12603.8,13629.8,14421.1,13978.3,12927.9,13429.9,13470.1,14785.8,14292,14308.8,14013,13240.9,12153.4,14289.7,15669.2,14169.5,14569.8,14469.1,14264.9,15320.9,14433.5,13691.5,14194.1,13519.2,11857.9,14616,15643.4,14077.2,14887.5,14159.9,14643,17192.5,15386.1,14287.1,17526.6,14497,14398.3,16629.6,16670.7,16614.8,16869.2,15663.9,16359.9,18447.7,16889,16505,18320.9,15052.1,15699.8,18135.3,16768.7,18883,19021,18101.9,17776.1,21489.9,17065.3,18690,18953.1,16398.9,16895.7,18553,19270,19422.1,17579.4,18637.3,18076.7,20438.6,18075.2,19563,19899.2,19227.5,17789.6,19220.8,22058.6,21230.8,19504.4,23913.1,23165.7,23574.3,25002,22603.9,23408.6) > 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] 94 > (np <- floor(n / par1)) [1] 7 > 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 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 14532.2 13629.8 14289.7 14616.0 16629.6 18135.3 18553.0 19220.8 [2,] 15167.0 14421.1 15669.2 15643.4 16670.7 16768.7 19270.0 22058.6 [3,] 16071.1 13978.3 14169.5 14077.2 16614.8 18883.0 19422.1 21230.8 [4,] 14827.5 12927.9 14569.8 14887.5 16869.2 19021.0 17579.4 19504.4 [5,] 15082.0 13429.9 14469.1 14159.9 15663.9 18101.9 18637.3 23913.1 [6,] 14772.7 13470.1 14264.9 14643.0 16359.9 17776.1 18076.7 23165.7 [7,] 16083.0 14785.8 15320.9 17192.5 18447.7 21489.9 20438.6 23574.3 [8,] 14272.5 14292.0 14433.5 15386.1 16889.0 17065.3 18075.2 25002.0 [9,] 15223.3 14308.8 13691.5 14287.1 16505.0 18690.0 19563.0 22603.9 [10,] 14897.3 14013.0 14194.1 17526.6 18320.9 18953.1 19899.2 23408.6 [11,] 13062.6 13240.9 13519.2 14497.0 15052.1 16398.9 19227.5 NA [12,] 12603.8 12153.4 11857.9 14398.3 15699.8 16895.7 17789.6 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 634.8 791.3 1379.5 1027.4 41.1 -1366.6 717.0 2837.8 [2,] 904.1 -442.8 -1499.7 -1566.2 -55.9 2114.3 152.1 -827.8 [3,] -1243.6 -1050.4 400.3 810.3 254.4 138.0 -1842.7 -1726.4 [4,] 254.5 502.0 -100.7 -727.6 -1205.3 -919.1 1057.9 4408.7 [5,] -309.3 40.2 -204.2 483.1 696.0 -325.8 -560.6 -747.4 [6,] 1310.3 1315.7 1056.0 2549.5 2087.8 3713.8 2361.9 408.6 [7,] -1810.5 -493.8 -887.4 -1806.4 -1558.7 -4424.6 -2363.4 1427.7 [8,] 950.8 16.8 -742.0 -1099.0 -384.0 1624.7 1487.8 -2398.1 [9,] -326.0 -295.8 502.6 3239.5 1815.9 263.1 336.2 804.7 [10,] -1834.7 -772.1 -674.9 -3029.6 -3268.8 -2554.2 -671.7 NA [11,] -458.8 -1087.5 -1661.3 -98.7 647.7 496.8 -1437.9 NA [12,] 1026.0 2136.3 2758.1 2231.3 2435.5 1657.3 1431.2 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/1h3261229890742.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/2uwcx1229890742.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/37cqy1229890742.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/45o2w1229890742.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] [1,] 13629.80 14421.10 13978.30 12927.90 13429.90 13470.10 14785.80 14272.50 [2,] 14410.95 15405.20 14123.35 14698.65 14314.50 14453.95 15701.95 14362.75 [3,] 15622.80 16169.95 16342.95 15878.35 15372.95 15566.30 17820.10 16137.55 [4,] 18344.15 18019.35 19152.55 18300.20 18369.60 17926.40 20964.25 17570.25 [5,] 19220.80 19270.00 21230.80 19504.40 23913.10 18076.70 23574.30 18075.20 [,9] [,10] [,11] [,12] [1,] 13691.50 14013.00 13062.60 11857.90 [2,] 14297.95 14545.70 13380.05 12378.60 [3,] 15864.15 17923.75 14497.00 14398.30 [4,] 19126.50 19426.15 15725.50 16297.75 [5,] 22603.90 23408.60 19227.50 17789.60 $n [1] 8 8 8 8 8 8 8 8 8 8 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 13425.66 14709.65 13533.57 13866.47 13107.71 13626.54 14880.50 14345.79 [2,] 17819.94 17630.25 19152.33 17890.23 17638.19 17506.06 20759.70 17929.31 [,9] [,10] [,11] [,12] [1,] 13166.85 15197.46 13096.34 12057.85 [2,] 18561.45 20650.04 15897.66 16738.75 $out [1] 22058.6 23165.7 25002.0 $group [1] 2 6 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(13629.8, 14410.95, 15622.8, 18344.15, 19220.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5k2yw1229890742.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] [1,] 41.10 -1566.20 -1842.70 -1205.30 -747.40 408.60 -2363.40 -2398.1 [2,] 337.95 -1163.75 -1485.00 -823.35 -443.20 1183.15 -2086.95 -920.5 [3,] 754.15 -249.35 -456.20 76.90 -256.75 1701.75 -1682.55 -183.6 [4,] 1203.45 528.10 327.35 779.95 261.65 2455.70 -690.60 1219.3 [5,] 1379.50 2114.30 810.30 1057.90 696.00 3713.80 -493.80 1624.7 [,9] [,10] [,11] [,12] [1,] -326.00 -3268.8 -1661.30 1026.00 [2,] -16.35 -2791.9 -1262.70 1544.25 [3,] 419.40 -1834.7 -458.80 2136.30 [4,] 1310.30 -723.5 199.05 2333.40 [5,] 3239.50 -671.7 647.70 2758.10 $n [1] 8 8 8 8 8 8 8 8 8 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 270.6693 -1194.4417 -1468.6047 -818.7264 -650.4893 990.8853 -2462.5712 [2,] 1237.6307 695.7417 556.2047 972.5264 136.9893 2412.6147 -902.5288 [,8] [,9] [,10] [,11] [,12] [1,] -1378.923 -321.6857 -3069.9151 -1331.7335 1665.032 [2,] 1011.723 1160.4857 -599.4849 414.1335 2607.568 $out [1] -1366.6 2837.8 4408.7 -4424.6 1427.7 $group [1] 1 1 4 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(41.1000000000022, 337.950000000001, 754.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/60id31229890742.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] [1,] 14272.50 12153.40 13519.20 14077.20 15052.10 16398.90 17579.40 19220.8 [2,] 14402.35 13335.40 13930.50 14342.70 16029.85 16980.50 18075.95 21230.8 [3,] 14862.40 13804.05 14277.30 14629.50 16622.20 18118.60 18932.40 22884.8 [4,] 15195.15 14300.40 14519.45 15514.75 16879.10 18918.05 19492.55 23574.3 [5,] 16083.00 14785.80 15320.90 17192.50 16889.00 21489.90 20438.60 25002.0 $n [1] 12 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 14500.8 13363.91 14008.68 14094.92 16234.85 17234.87 18286.28 21713.89 [2,] 15224.0 14244.19 14545.92 15164.08 17009.55 19002.33 19578.52 24055.71 $out [1] 13062.6 12603.8 15669.2 11857.9 17526.6 18447.7 18320.9 $group [1] 1 1 3 3 4 5 5 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(14272.5, 14402.35, 14862.4, 15195.15, 16083, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/767sf1229890743.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,] 16200.80 14398.30 15857.95 [2,] 16237.07 15469.62 16028.91 [3,] 16743.99 15871.25 16353.54 [4,] 16942.77 16256.45 16593.96 [5,] 17651.60 16342.95 16955.56 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 16422.12 15512.37 16095.82 [2,] 17065.87 16230.13 16611.27 $out [1] 18416.59 14999.74 14485.50 17820.10 17923.75 18296.95 14552.77 14338.17 $group [1] 1 1 1 2 2 3 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(16200.8, 16237.06875, 16743.99375, 16942.76875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1h3261229890742.ps tmp/1h3261229890742.png") > system("convert tmp/2uwcx1229890742.ps tmp/2uwcx1229890742.png") > system("convert tmp/37cqy1229890742.ps tmp/37cqy1229890742.png") > system("convert tmp/45o2w1229890742.ps tmp/45o2w1229890742.png") > system("convert tmp/5k2yw1229890742.ps tmp/5k2yw1229890742.png") > system("convert tmp/60id31229890742.ps tmp/60id31229890742.png") > system("convert tmp/767sf1229890743.ps tmp/767sf1229890743.png") > > > proc.time() user system elapsed 1.228 0.941 1.476