R version 2.11.1 (2010-05-31) Copyright (C) 2010 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(152,151,150,148,146,145,146,148,149,149,150,152,154,164,167,162,164,172,169,174,181,181,172,181,183,200,199,190,197,194,190,195,204,197,185,193,192,211,210,197,191,182,170,166,175,163,153,171,165,184,179,163,163,148,132,127,130,118,113,137,133,155,151,132,134,118,102,98,91,77,74,102,98,113,114,96,102,90,72,67,62,48,43,75) > 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] 84 > (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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 152 154 183 192 165 133 98 NA [2,] 151 164 200 211 184 155 113 NA [3,] 150 167 199 210 179 151 114 NA [4,] 148 162 190 197 163 132 96 NA [5,] 146 164 197 191 163 134 102 NA [6,] 145 172 194 182 148 118 90 NA [7,] 146 169 190 170 132 102 72 NA [8,] 148 174 195 166 127 98 67 NA [9,] 149 181 204 175 130 91 62 NA [10,] 149 181 197 163 118 77 48 NA [11,] 150 172 185 153 113 74 43 NA [12,] 152 181 193 171 137 102 75 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 10 17 19 19 22 15 NA [2,] -1 3 -1 -1 -5 -4 1 NA [3,] -2 -5 -9 -13 -16 -19 -18 NA [4,] -2 2 7 -6 0 2 6 NA [5,] -1 8 -3 -9 -15 -16 -12 NA [6,] 1 -3 -4 -12 -16 -16 -18 NA [7,] 2 5 5 -4 -5 -4 -5 NA [8,] 1 7 9 9 3 -7 -5 NA [9,] 0 0 -7 -12 -12 -14 -14 NA [10,] 1 -9 -12 -10 -5 -3 -5 NA [11,] 2 9 8 18 24 28 32 NA [12,] 2 2 -1 -6 -4 -4 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/rcomp/tmp/1wkye1281958578.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/rcomp/tmp/2wkye1281958578.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/rcomp/tmp/3wkye1281958578.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/rcomp/tmp/47cgh1281958578.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] [,10] [,11] [,12] [1,] 98.0 113 114.0 96.0 102.0 90.0 72.0 67.0 62.0 48.0 43.0 75.0 [2,] 142.5 153 150.5 140.0 140.0 131.5 117.0 112.5 110.5 97.5 93.5 119.5 [3,] 154.0 164 167.0 162.0 163.0 148.0 146.0 148.0 149.0 149.0 150.0 152.0 [4,] 174.0 192 189.0 176.5 177.5 177.0 169.5 170.0 178.0 172.0 162.5 176.0 [5,] 192.0 211 210.0 197.0 197.0 194.0 190.0 195.0 204.0 197.0 185.0 193.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 135.1887 140.7098 144.0084 140.2028 140.6056 120.8281 114.6478 113.6619 [2,] 172.8113 187.2902 189.9916 183.7972 185.3944 175.1719 177.3522 182.3381 [,9] [,10] [,11] [,12] [1,] 108.6901 104.5098 108.7943 118.2591 [2,] 189.3099 193.4902 191.2057 185.7409 $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(98, 142.5, 154, 174, 192, 113, 153, 164, 192, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/57cgh1281958578.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] [,10] [,11] [,12] [1,] 10.0 -5.0 -19 -6 -16.0 -18.0 -5.0 -7 -14.0 -12.0 2.0 -6.0 [2,] 12.5 -2.5 -17 -1 -13.5 -16.0 -4.5 -2 -13.0 -9.5 8.5 -4.0 [3,] 17.0 -1.0 -13 2 -9.0 -12.0 -4.0 3 -12.0 -5.0 18.0 -2.5 [4,] 19.0 0.0 -7 4 -2.0 -3.5 3.5 8 -3.5 -4.0 26.0 2.0 [5,] 22.0 3.0 -2 7 8.0 1.0 5.0 9 0.0 1.0 32.0 2.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 13.11830 -2.4929597 -18.971839 -0.9859193 -15.867614 -19.464798 -8.777471 [2,] 20.88170 0.4929597 -7.028161 4.9859193 -2.132386 -4.535202 0.777471 [,8] [,9] [,10] [,11] [,12] [1,] -2.971839 -17.673247 -8.284511 7.549282 -6.370194 [2,] 8.971839 -6.326753 -1.715489 28.450718 1.370194 $out [1] -1 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(10, 12.5, 17, 19, 22, -5, -2.5, -1, 0, 3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/603fk1281958578.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,] 145.0 154.0 183.0 153.0 113.0 74.0 43.0 NA [2,] 147.0 164.0 190.0 168.0 128.5 94.5 64.5 NA [3,] 149.0 170.5 194.5 178.5 142.5 110.0 82.5 NA [4,] 150.5 177.5 198.0 194.5 164.0 133.5 100.0 NA [5,] 152.0 181.0 204.0 211.0 184.0 155.0 114.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 147.4036 164.3426 190.8511 166.4132 126.3082 92.21184 66.30821 NA [2,] 150.5964 176.6574 198.1489 190.5868 158.6918 127.78816 98.69179 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(145, 147, 149, 150.5, 152, 154, 164, 170.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/703fk1281958578.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,] 127.1429 146.0 128.00 [2,] 139.7143 148.5 142.25 [3,] 147.1429 151.0 151.00 [4,] 156.0714 162.5 158.50 [5,] 168.2857 167.0 172.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 139.6823 144.6145 143.5883 [2,] 154.6035 157.3855 158.4117 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(127.142857142857, 139.714285714286, 147.142857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1wkye1281958578.ps tmp/1wkye1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/2wkye1281958578.ps tmp/2wkye1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/3wkye1281958578.ps tmp/3wkye1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/47cgh1281958578.ps tmp/47cgh1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/57cgh1281958578.ps tmp/57cgh1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/603fk1281958578.ps tmp/603fk1281958578.png",intern=TRUE)) character(0) > try(system("convert tmp/703fk1281958578.ps tmp/703fk1281958578.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.960 1.320 3.249