R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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(14724,14404,14058,13427,19946,19631,14724,11462,11778,11778,12093,12760,13742,13427,11462,11778,20929,22893,17666,14724,15387,15707,17351,18964,19315,16022,16369,12093,24222,27800,19631,17004,18649,20613,23555,27164,27164,24853,23871,17982,27800,32391,28462,24222,24853,27164,30426,34355,31724,30111,30111,24853,32391,37297,33373,29129,30426,35653,37964,41222,38595,34355,33373,25520,30742,36315,30111,26502,30111,33689,35653,40906,38280,31724,32391,26186,31409,36000,30742,27164,30426,34355,33689,41542,40244,35017,35333,28462,32706,39262,34355,31409,36315,39262,36982,47431,44835,38946,37297,29764,34035,37964,33053,33053,38595,41542,39924,51355,48413,42871,40560,32391,35333,40560,36631,35653,40244,44168,39924,50057) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 14724 13742 19315 27164 31724 38595 38280 40244 44835 48413 NA [2,] 14404 13427 16022 24853 30111 34355 31724 35017 38946 42871 NA [3,] 14058 11462 16369 23871 30111 33373 32391 35333 37297 40560 NA [4,] 13427 11778 12093 17982 24853 25520 26186 28462 29764 32391 NA [5,] 19946 20929 24222 27800 32391 30742 31409 32706 34035 35333 NA [6,] 19631 22893 27800 32391 37297 36315 36000 39262 37964 40560 NA [7,] 14724 17666 19631 28462 33373 30111 30742 34355 33053 36631 NA [8,] 11462 14724 17004 24222 29129 26502 27164 31409 33053 35653 NA [9,] 11778 15387 18649 24853 30426 30111 30426 36315 38595 40244 NA [10,] 11778 15707 20613 27164 35653 33689 34355 39262 41542 44168 NA [11,] 12093 17351 23555 30426 37964 35653 33689 36982 39924 39924 NA [12,] 12760 18964 27164 34355 41222 40906 41542 47431 51355 50057 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -320 -315 -3293 -2311 -1613 -4240 -6556 -5227 -5889 -5542 NA [2,] -346 -1965 347 -982 0 -982 667 316 -1649 -2311 NA [3,] -631 316 -4276 -5889 -5258 -7853 -6205 -6871 -7533 -8169 NA [4,] 6519 9151 12129 9818 7538 5222 5223 4244 4271 2942 NA [5,] -315 1964 3578 4591 4906 5573 4591 6556 3929 5227 NA [6,] -4907 -5227 -8169 -3929 -3924 -6204 -5258 -4907 -4911 -3929 NA [7,] -3262 -2942 -2627 -4240 -4244 -3609 -3578 -2946 0 -978 NA [8,] 316 663 1645 631 1297 3609 3262 4906 5542 4591 NA [9,] 0 320 1964 2311 5227 3578 3929 2947 2947 3924 NA [10,] 315 1644 2942 3262 2311 1964 -666 -2280 -1618 -4244 NA [11,] 667 1613 3609 3929 3258 5253 7853 10449 11431 10133 NA [12,] 982 351 0 -2631 -2627 -2626 -1298 -2596 -2942 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/wessaorg/rcomp/tmp/1r8rk1313682903.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/wessaorg/rcomp/tmp/2y2131313682903.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/wessaorg/rcomp/tmp/35kqh1313682903.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/wessaorg/rcomp/tmp/4zrf71313682903.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] [1,] 13742 13427.0 11462 11778.0 19946.0 19631.0 14724.0 11462 11778.0 11778 [2,] 19315 16022.0 16369 13427.0 24222.0 27800.0 19631.0 17004 18649.0 20613 [3,] 35002 30917.5 31251 25186.5 31075.5 36157.5 30426.5 26833 30268.5 34022 [4,] 40244 35017.0 35333 28462.0 32706.0 37964.0 33373.0 31409 36315.0 39262 [5,] 48413 42871.0 40560 32391.0 35333.0 40560.0 36631.0 35653 40244.0 44168 [,11] [,12] [1,] 12093 12760 [2,] 23555 27164 [3,] 34671 41064 [4,] 37964 47431 [5,] 39924 51355 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 24545.04 21426.84 21775.83 17674.41 26836.56 31079.16 23560.45 19635.69 [2,] 45458.96 40408.16 40726.17 32698.59 35314.44 41235.84 37292.55 34030.31 [,9] [,10] [,11] [,12] [1,] 21441.86 24704.22 27471.69 30937.8 [2,] 39095.14 43339.78 41870.31 51190.2 $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(13742, 19315, 35002, 40244, 48413, 13427, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5st2g1313682903.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] [1,] -6556.0 -2311 -8169 2942 1964 -6204 -4244 316.0 0 -4244.0 667 -2942 [2,] -5542.0 -1649 -7533 4271 3578 -5258 -3609 663.0 1964 -1618.0 3258 -2627 [3,] -3766.5 -664 -6047 5871 4591 -4909 -3104 2453.5 2947 979.5 4591 -2596 [4,] -1613.0 316 -4276 9151 5227 -3929 -2627 4591.0 3924 2311.0 10133 0 [5,] -315.0 667 316 12129 6556 -3924 -2627 5542.0 5227 3262.0 11431 982 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -5729.585 -1645.7923 -7674.327 3432.757 3767.094 -5573.021 -3594.646 [2,] -1803.415 317.7923 -4419.673 8309.243 5414.906 -4244.979 -2613.354 [,8] [,9] [,10] [,11] [,12] [1,] 490.9146 1967.706 -983.5851 1155.976 -3979.553 [2,] 4416.0854 3926.294 2942.5851 8026.024 -1212.447 $out [1] -315 -8169 0 -978 $group [1] 5 6 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-6556, -5542, -3766.5, -1613, -315, -2311, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/60m641313682903.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] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 11462.0 11462.0 12093.0 17982.0 24853.0 25520 26186.0 28462.0 29764.0 [2,] 11935.5 13584.5 16686.5 24537.5 30111.0 30111 30584.0 33530.5 33544.0 [3,] 13742.5 15547.0 19473.0 27164.0 32057.5 33531 32057.5 35824.0 38279.5 [4,] 14724.0 18315.0 23888.5 29444.0 36475.0 35984 35177.5 39262.0 40733.0 [5,] 14724.0 22893.0 27800.0 34355.0 41222.0 40906 41542.0 47431.0 51355.0 [,10] [,11] [1,] 32391.0 NA [2,] 36142.0 NA [3,] 40402.0 NA [4,] 43519.5 NA [5,] 50057.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12470.65 13389.39 16188.12 24926.11 29154.84 30852.29 29962.37 33209.82 [2,] 15014.35 17704.61 22757.88 29401.89 34960.16 36209.71 34152.63 38438.18 [,9] [,10] [,11] [1,] 35000.55 37037.07 NA [2,] 41558.45 43766.93 NA $out [1] 19946 19631 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(11462, 11935.5, 13742.5, 14724, 14724, 11462, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7zq961313682903.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,] 22245.60 25186.50 21229.38 [2,] 27580.45 30347.50 26542.19 [3,] 28562.15 31163.25 28218.88 [4,] 31229.85 34836.50 31025.06 [5,] 36575.60 41064.00 37460.25 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 26897.63 29115.79 26174.21 [2,] 30226.67 33210.71 30263.54 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(22245.6, 27580.45, 28562.15, 31229.85, 36575.6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1r8rk1313682903.ps tmp/1r8rk1313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/2y2131313682903.ps tmp/2y2131313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/35kqh1313682903.ps tmp/35kqh1313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/4zrf71313682903.ps tmp/4zrf71313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/5st2g1313682903.ps tmp/5st2g1313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/60m641313682903.ps tmp/60m641313682903.png",intern=TRUE)) character(0) > try(system("convert tmp/7zq961313682903.ps tmp/7zq961313682903.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.662 0.363 2.029