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(6802.96,7132.68,7073.29,7264.5,7105.33,7218.71,7225.72,7354.25,7745.46,8070.26,8366.33,8667.51,8854.34,9218.1,9332.9,9358.31,9248.66,9401.2,9652.04,9957.38,10110.63,10169.26,10343.78,10750.21,11337.5,11786.96,12083.04,12007.74,11745.93,11051.51,11445.9,11924.88,12247.63,12690.91,12910.7,13202.12,13654.67,13862.82,13523.93,14211.17,14510.35,14289.23,14111.82,13086.59,13351.54,13747.69,12855.61,12926.93,12121.95,11731.65,11639.51,12163.78,12029.53,11234.18,9852.13,9709.04,9332.75,7108.6,6691.49,6143.05) > 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,] 6802.96 8854.34 11337.50 13654.67 12121.95 NA [2,] 7132.68 9218.10 11786.96 13862.82 11731.65 NA [3,] 7073.29 9332.90 12083.04 13523.93 11639.51 NA [4,] 7264.50 9358.31 12007.74 14211.17 12163.78 NA [5,] 7105.33 9248.66 11745.93 14510.35 12029.53 NA [6,] 7218.71 9401.20 11051.51 14289.23 11234.18 NA [7,] 7225.72 9652.04 11445.90 14111.82 9852.13 NA [8,] 7354.25 9957.38 11924.88 13086.59 9709.04 NA [9,] 7745.46 10110.63 12247.63 13351.54 9332.75 NA [10,] 8070.26 10169.26 12690.91 13747.69 7108.60 NA [11,] 8366.33 10343.78 12910.70 12855.61 6691.49 NA [12,] 8667.51 10750.21 13202.12 12926.93 6143.05 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 329.72 363.76 449.46 208.15 -390.30 NA [2,] -59.39 114.80 296.08 -338.89 -92.14 NA [3,] 191.21 25.41 -75.30 687.24 524.27 NA [4,] -159.17 -109.65 -261.81 299.18 -134.25 NA [5,] 113.38 152.54 -694.42 -221.12 -795.35 NA [6,] 7.01 250.84 394.39 -177.41 -1382.05 NA [7,] 128.53 305.34 478.98 -1025.23 -143.09 NA [8,] 391.21 153.25 322.75 264.95 -376.29 NA [9,] 324.80 58.63 443.28 396.15 -2224.15 NA [10,] 296.07 174.52 219.79 -892.08 -417.11 NA [11,] 301.18 406.43 291.42 71.32 -548.44 NA [12,] 186.83 587.29 452.55 -804.98 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/1w0qi1257850671.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/28tjk1257850671.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/34hxg1257850671.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/4s4f81257850671.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,] 6802.96 7132.68 7073.29 7264.50 7105.33 7218.71 7225.72 7354.25 [2,] 8854.34 9218.10 9332.90 9358.31 9248.66 9401.20 9652.04 9709.04 [3,] 11337.50 11731.65 11639.51 12007.74 11745.93 11051.51 9852.13 9957.38 [4,] 12121.95 11786.96 12083.04 12163.78 12029.53 11234.18 11445.90 11924.88 [5,] 13654.67 13862.82 13523.93 14211.17 14510.35 11234.18 14111.82 13086.59 [,9] [,10] [,11] [,12] [1,] 7745.46 7108.60 6691.49 6143.05 [2,] 9332.75 8070.26 8366.33 8667.51 [3,] 10110.63 10169.26 10343.78 10750.21 [4,] 12247.63 12690.91 12855.61 12926.93 [5,] 13351.54 13747.69 12910.70 13202.12 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 9028.615 9916.5 9696.268 10025.40 9780.974 9756.33 8584.593 [2,] 13646.385 13546.8 13582.752 13990.08 13710.886 12346.69 11119.667 [,8] [,9] [,10] [,11] [,12] [1,] 8391.673 8050.983 6904.32 7171.666 7740.515 [2,] 11523.087 12170.277 13434.20 13515.894 13759.905 $out [1] 14289.23 $group [1] 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(6802.96, 8854.34, 11337.5, 12121.95, 13654.67, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5kjv81257850671.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,] 208.15 -338.89 -75.30 -159.17 -795.35 -177.41 -143.09 153.25 58.63 [2,] 208.15 -92.14 25.41 -159.17 -694.42 -177.41 -143.09 153.25 58.63 [3,] 329.72 -59.39 191.21 -134.25 -221.12 7.01 128.53 264.95 324.80 [4,] 363.76 114.80 524.27 -109.65 113.38 250.84 305.34 322.75 396.15 [5,] 449.46 296.08 687.24 -109.65 152.54 394.39 478.98 391.21 443.28 [,10] [,11] [,12] [1,] -892.08 71.32 -804.980 [2,] -417.11 71.32 -309.075 [3,] 174.52 291.42 319.690 [4,] 219.79 301.18 519.920 [5,] 296.07 406.43 587.290 $n [1] 5 5 5 5 5 5 5 5 5 5 5 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 219.7664 -205.61328 -161.2832 -169.2407 -791.9094 -295.5904 -188.3295 [2,] 439.6736 86.83328 543.7032 -99.2593 349.6694 309.6104 445.3895 [,8] [,9] [,10] [,11] [,12] [1,] 145.1817 86.30922 -275.5119 129.0015 -335.2161 [2,] 384.7183 563.29078 624.5519 453.8385 974.5960 $out [1] -390.30 -261.81 299.18 -1382.05 -1025.23 -376.29 -2224.15 -548.44 $group [1] 1 4 4 6 7 8 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(208.150000000000, 208.150000000000, 329.72, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6u1ta1257850671.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,] 6802.960 8854.34 11051.51 12855.61 6143.050 NA [2,] 7119.005 9290.78 11595.92 13219.07 8220.675 NA [3,] 7245.110 9526.62 11966.31 13701.18 10543.155 NA [4,] 7907.860 10139.94 12469.27 14161.49 11880.590 NA [5,] 8667.510 10750.21 13202.12 14510.35 12163.780 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 6885.308 9139.31 11567.97 13271.33 8873.843 NA [2,] 7604.912 9913.93 12364.65 14131.03 12212.467 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(6802.96, 7119.005, 7245.11, 7907.86, 8667.51, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/734qv1257850671.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,] 10233.58 9852.13 10317.69 [2,] 10381.89 10139.94 10495.34 [3,] 10555.94 10900.86 10625.03 [4,] 10738.49 11685.58 10775.62 [5,] 11001.10 12007.74 10816.96 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 10393.29 10195.89 10497.19 [2,] 10718.59 11605.83 10752.87 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/1w0qi1257850671.ps tmp/1w0qi1257850671.png") > system("convert tmp/28tjk1257850671.ps tmp/28tjk1257850671.png") > system("convert tmp/34hxg1257850671.ps tmp/34hxg1257850671.png") > system("convert tmp/4s4f81257850671.ps tmp/4s4f81257850671.png") > system("convert tmp/5kjv81257850671.ps tmp/5kjv81257850671.png") > system("convert tmp/6u1ta1257850671.ps tmp/6u1ta1257850671.png") > system("convert tmp/734qv1257850671.ps tmp/734qv1257850671.png") > > > proc.time() user system elapsed 1.233 0.980 3.178