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(25000,25284,12434.5,33955,14980.5,50831,4198.5,34566,35000,11055.5,20807,21887.29,16977.5,19613.5,14570,24416.5,16825.5,13980,21450.5,27239.5,19078.5,20459.1,20373.5,19306.5,16723.16,11638,20917,17903.5,28218.5,15268,21555,23143,16691,17932.5,30512,41931.5,10853.5,25939.5,14900,25127.76,22063.5,25306.5,31217.5,23201.5,38148,26264,16359,27945.5,16218.5,36003.5,20323.5,20100.5,18741,24426.75,19174.5,13766,18999,21745,34469,13248,16218.5,36003.5,20323.5,20100.5,18741,24426.75,19174.5,13766,18999,21745,34469,13248) > 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,] 25000.00 16977.5 16723.16 10853.50 16218.50 16218.50 NA [2,] 25284.00 19613.5 11638.00 25939.50 36003.50 36003.50 NA [3,] 12434.50 14570.0 20917.00 14900.00 20323.50 20323.50 NA [4,] 33955.00 24416.5 17903.50 25127.76 20100.50 20100.50 NA [5,] 14980.50 16825.5 28218.50 22063.50 18741.00 18741.00 NA [6,] 50831.00 13980.0 15268.00 25306.50 24426.75 24426.75 NA [7,] 4198.50 21450.5 21555.00 31217.50 19174.50 19174.50 NA [8,] 34566.00 27239.5 23143.00 23201.50 13766.00 13766.00 NA [9,] 35000.00 19078.5 16691.00 38148.00 18999.00 18999.00 NA [10,] 11055.50 20459.1 17932.50 26264.00 21745.00 21745.00 NA [11,] 20807.00 20373.5 30512.00 16359.00 34469.00 34469.00 NA [12,] 21887.29 19306.5 41931.50 27945.50 13248.00 13248.00 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 284.00 2636.00 -5085.16 15086.00 19785.00 19785.00 NA [2,] -12849.50 -5043.50 9279.00 -11039.50 -15680.00 -15680.00 NA [3,] 21520.50 9846.50 -3013.50 10227.76 -223.00 -223.00 NA [4,] -18974.50 -7591.00 10315.00 -3064.26 -1359.50 -1359.50 NA [5,] 35850.50 -2845.50 -12950.50 3243.00 5685.75 5685.75 NA [6,] -46632.50 7470.50 6287.00 5911.00 -5252.25 -5252.25 NA [7,] 30367.50 5789.00 1588.00 -8016.00 -5408.50 -5408.50 NA [8,] 434.00 -8161.00 -6452.00 14946.50 5233.00 5233.00 NA [9,] -23944.50 1380.60 1241.50 -11884.00 2746.00 2746.00 NA [10,] 9751.50 -85.60 12579.50 -9905.00 12724.00 12724.00 NA [11,] 1080.29 -1067.00 11419.50 11586.50 -21221.00 -21221.00 NA [12,] -4909.79 -2583.34 -31078.00 -11727.00 2970.50 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/131201275657646.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/231201275657646.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/331201275657646.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/4eb131275657646.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,] 16218.50 11638.00 12434.50 17903.50 14980.5 13980.00 19174.5 13766.00 [2,] 16218.50 19613.50 14570.00 20100.50 16825.5 15268.00 19174.5 13766.00 [3,] 16470.83 25611.75 17611.75 22258.50 18741.0 24426.75 20312.5 23172.25 [4,] 16977.50 36003.50 20323.50 25127.76 22063.5 25306.50 21555.0 27239.50 [5,] 16977.50 36003.50 20917.00 25127.76 28218.5 25306.50 21555.0 34566.00 [,9] [,10] [,11] [,12] [1,] 16691.00 17932.50 16359.0 13248.00 [2,] 18999.00 17932.50 20373.5 13248.00 [3,] 19038.75 21102.05 25659.5 20596.90 [4,] 35000.00 21745.00 34469.0 27945.50 [5,] 38148.00 26264.00 34469.0 41931.50 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15981.25 15039.67 13900.56 19015.75 15362.32 17951.59 18777 14481.41 [2,] 16960.41 36183.83 21322.94 25501.25 22119.68 30901.91 21848 31863.09 [,9] [,10] [,11] [,12] [1,] 8717.588 18642.86 16567.45 11116.53 [2,] 29359.912 23561.24 34751.55 30077.26 $out [1] 25000.0 10853.5 33955.0 50831.0 4198.5 31217.5 11055.5 $group [1] 1 1 4 6 7 7 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(16218.5, 16218.5, 16470.83, 16977.5, 16977.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5eb131275657646.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,] -5085.16 -15680.0 -3013.50 -7591.00 -12950.500 -5252.250 -8016.00 -8161.0 [2,] 284.00 -15680.0 -223.00 -7591.00 -2845.500 -5252.250 -5408.50 -6452.0 [3,] 8861.00 -11944.5 4811.75 -2211.88 4464.375 329.375 -1910.25 2833.5 [4,] 19785.00 -5043.5 10227.76 -1359.50 5685.750 6287.000 5789.00 5233.0 [5,] 19785.00 9279.0 21520.50 -1359.50 5685.750 7470.500 5789.00 14946.5 [,9] [,10] [,11] [,12] [1,] -23944.50 -9905.0 -21221.000 -11727.00 [2,] -11884.00 -85.6 -21221.000 -11727.00 [3,] 1311.05 11165.5 6.645 -4909.79 [4,] 2746.00 12724.0 11419.500 -2583.34 [5,] 2746.00 12724.0 11586.500 2970.50 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3717.775 -18805.386 -1929.328 -6231.399 -1038.557 -7113.814 -9132.999 [2,] 21439.775 -5083.614 11552.828 1807.639 9967.307 7772.564 5312.499 [,8] [,9] [,10] [,11] [,12] [1,] -4703.702 -8125.773 2902.894 -21047.53 -11370.677 [2,] 10370.702 10747.873 19428.106 21060.82 1551.097 $out [1] -18974.5 10315.0 35850.5 -46632.5 30367.5 -31078.0 $group [1] 4 4 5 6 7 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-5085.16, 284, 8861, 19785, 19785, -15680, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6okjo1275657646.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,] 4198.50 13980.0 11638.00 10853.50 13248.00 13248.00 NA [2,] 13707.50 16901.5 16707.08 19211.25 17479.75 17479.75 NA [3,] 23443.65 19460.0 19424.75 25217.13 19637.50 19637.50 NA [4,] 34260.50 20954.8 25680.75 27104.75 23085.88 23085.88 NA [5,] 50831.00 24416.5 30512.00 38148.00 24426.75 24426.75 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 14069.28 17611.26 15331.8 21616.85 17080.51 17080.51 NA [2,] 32818.01 21308.74 23517.7 28817.41 22194.49 22194.49 NA $out [1] 27239.5 41931.5 36003.5 34469.0 36003.5 34469.0 $group [1] 2 3 5 5 6 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(4198.5, 13707.5, 23443.645, 34260.5, 50831, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7okjo1275657646.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,] 16998.53 16470.83 16566.21 [2,] 19664.30 18889.88 19711.60 [3,] 22770.73 20849.47 20883.46 [4,] 25096.21 23799.50 23767.27 [5,] 26164.92 25659.50 27259.31 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 20293.20 18610.16 19033.64 [2,] 25248.26 23088.79 22733.27 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(16998.5266666667, 19664.3, 22770.7325, 25096.2083333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/131201275657646.ps tmp/131201275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/231201275657646.ps tmp/231201275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/331201275657646.ps tmp/331201275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/4eb131275657646.ps tmp/4eb131275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/5eb131275657646.ps tmp/5eb131275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/6okjo1275657646.ps tmp/6okjo1275657646.png",intern=TRUE)) character(0) > try(system("convert tmp/7okjo1275657646.ps tmp/7okjo1275657646.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.197 0.947 12.749