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(76,75,74,72,70,69,70,72,73,73,74,76,74,67,66,58,55,58,64,68,66,76,75,88,85,83,77,66,65,65,63,62,57,68,69,79,74,76,82,75,75,76,78,77,67,74,68,87,76,88,95,96,96,105,108,113,101,107,102,116,105,121,134,140,131,141,131,128,123,129,125,144,135,141,156,159,146,154,145,133,126,127,122,148) > 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,] 76 74 85 74 76 105 135 NA [2,] 75 67 83 76 88 121 141 NA [3,] 74 66 77 82 95 134 156 NA [4,] 72 58 66 75 96 140 159 NA [5,] 70 55 65 75 96 131 146 NA [6,] 69 58 65 76 105 141 154 NA [7,] 70 64 63 78 108 131 145 NA [8,] 72 68 62 77 113 128 133 NA [9,] 73 66 57 67 101 123 126 NA [10,] 73 76 68 74 107 129 127 NA [11,] 74 75 69 68 102 125 122 NA [12,] 76 88 79 87 116 144 148 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 -7 -2 2 12 16 6 NA [2,] -1 -1 -6 6 7 13 15 NA [3,] -2 -8 -11 -7 1 6 3 NA [4,] -2 -3 -1 0 0 -9 -13 NA [5,] -1 3 0 1 9 10 8 NA [6,] 1 6 -2 2 3 -10 -9 NA [7,] 2 4 -1 -1 5 -3 -12 NA [8,] 1 -2 -5 -10 -12 -5 -7 NA [9,] 0 10 11 7 6 6 1 NA [10,] 1 -1 1 -6 -5 -4 -5 NA [11,] 2 13 10 19 14 19 26 NA [12,] -2 -3 -5 -11 -11 -9 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/162721282053318.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/262721282053318.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/3hto51282053318.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/4hto51282053318.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,] 74 67.0 66.0 58 55.0 58 63.0 62.0 57.0 68.0 68.0 76 [2,] 75 75.5 75.5 69 67.5 67 67.0 70.0 66.5 73.5 71.5 83 [3,] 76 83.0 82.0 75 75.0 76 78.0 77.0 73.0 76.0 75.0 88 [4,] 95 104.5 114.5 118 113.5 123 119.5 120.5 112.0 117.0 112.0 130 [5,] 105 141.0 156.0 159 146.0 154 145.0 133.0 126.0 129.0 125.0 148 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 64.05632 65.68167 58.70983 45.73799 47.52954 42.5577 46.64785 [2,] 87.94368 100.31833 105.29017 104.26201 102.47046 109.4423 109.35215 [,8] [,9] [,10] [,11] [,12] [1,] 46.84221 45.82813 50.0225 50.81405 59.93236 [2,] 107.15779 100.17187 101.9775 99.18595 116.06764 $out [1] 135 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(74, 75, 76, 95, 105, 67, 75.5, 83, 104.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5sk5q1282053318.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,] -7.0 -6 -11.0 -13.0 -1.0 -10.0 -3 -12.0 0.0 -6 2.0 -11 [2,] -1.5 -1 -7.5 -6.0 0.5 -5.5 -2 -8.5 3.5 -5 11.5 -11 [3,] 2.0 6 -2.0 -2.0 3.0 1.0 -1 -5.0 6.0 -4 14.0 -7 [4,] 9.0 10 2.0 -0.5 8.5 2.5 3 -3.5 8.5 0 19.0 -3 [5,] 16.0 15 6.0 0.0 10.0 6.0 5 1.0 11.0 1 26.0 -2 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4.270431 -0.5690225 -7.673247 -5.284511 -1.777471 -3.777471 -3.985919 [2,] 8.270431 12.5690225 3.673247 1.284511 7.777471 5.777471 1.985919 [,8] [,9] [,10] [,11] [,12] [1,] -7.985919 3.014081 -6.985919 9.521121 -12.160258 [2,] -2.014081 8.985919 -1.014081 18.478879 -1.839742 $out [1] -12 $group [1] 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-7, -1.5, 2, 9, 16, -6, -1, 6, 10, 15, -11, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6sk5q1282053318.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,] 69.0 55.0 57 74.0 88.0 105 122 NA [2,] 71.0 61.0 64 74.0 95.5 124 130 NA [3,] 73.0 66.5 67 75.5 101.5 130 143 NA [4,] 74.5 74.5 78 77.5 107.5 137 151 NA [5,] 76.0 88.0 85 82.0 116.0 144 159 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 71.40363 60.34256 60.6145 73.90363 96.02672 124.0706 133.4218 NA [2,] 74.59637 72.65744 73.3855 77.09637 106.97328 135.9294 152.5782 NA $out [1] 67 68 87 76 $group [1] 4 4 4 5 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(69, 71, 73, 74.5, 76, 55, 61, 66.5, 74.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7sk5q1282053318.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,] 87.57143 73 85.000 [2,] 90.92857 75 90.250 [3,] 93.35714 76 93.375 [4,] 95.28571 80 95.125 [5,] 97.71429 83 95.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 91.36982 73.71947 91.15148 [2,] 95.34446 78.28053 95.59852 $out [1] 105.4286 88.0000 106.5000 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(87.5714285714286, 90.9285714285714, 93.3571428571429, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/162721282053318.ps tmp/162721282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/262721282053318.ps tmp/262721282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/3hto51282053318.ps tmp/3hto51282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/4hto51282053318.ps tmp/4hto51282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/5sk5q1282053318.ps tmp/5sk5q1282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/6sk5q1282053318.ps tmp/6sk5q1282053318.png",intern=TRUE)) character(0) > try(system("convert tmp/7sk5q1282053318.ps tmp/7sk5q1282053318.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.860 1.430 3.254