R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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(190,189,188,186,184,183,184,186,187,187,188,190,190,190,197,187,185,182,182,191,183,192,178,181,179,175,183,179,178,175,170,179,169,178,161,168,167,165,181,181,184,181,177,183,162,166,151,162,159,152,164,158,160,161,151,149,131,138,130,147,151,140,149,143,145,139,136,133,118,130,121,142,148,131,137,128,130,119,107,113,93,106,98,118) > 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,] 190 190 179 167 159 151 148 NA [2,] 189 190 175 165 152 140 131 NA [3,] 188 197 183 181 164 149 137 NA [4,] 186 187 179 181 158 143 128 NA [5,] 184 185 178 184 160 145 130 NA [6,] 183 182 175 181 161 139 119 NA [7,] 184 182 170 177 151 136 107 NA [8,] 186 191 179 183 149 133 113 NA [9,] 187 183 169 162 131 118 93 NA [10,] 187 192 178 166 138 130 106 NA [11,] 188 178 161 151 130 121 98 NA [12,] 190 181 168 162 147 142 118 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 0 -4 -2 -7 -11 -17 NA [2,] -1 7 8 16 12 9 6 NA [3,] -2 -10 -4 0 -6 -6 -9 NA [4,] -2 -2 -1 3 2 2 2 NA [5,] -1 -3 -3 -3 1 -6 -11 NA [6,] 1 0 -5 -4 -10 -3 -12 NA [7,] 2 9 9 6 -2 -3 6 NA [8,] 1 -8 -10 -21 -18 -15 -20 NA [9,] 0 9 9 4 7 12 13 NA [10,] 1 -14 -17 -15 -8 -9 -8 NA [11,] 2 3 7 11 17 21 20 NA [12,] 0 -2 -1 -3 4 6 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/freestat/rcomp/tmp/1pgpd1281547553.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/freestat/rcomp/tmp/2pgpd1281547553.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/freestat/rcomp/tmp/3pgpd1281547553.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/freestat/rcomp/tmp/4ip7y1281547553.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,] 148.0 131 137.0 128.0 130.0 119.0 107.0 113.0 93.0 106.0 98.0 118.0 [2,] 155.0 146 156.5 150.5 152.5 150.0 143.5 141.0 124.5 134.0 125.5 144.5 [3,] 167.0 165 181.0 179.0 178.0 175.0 170.0 179.0 162.0 166.0 151.0 162.0 [4,] 184.5 182 185.5 183.5 184.0 181.5 179.5 184.5 176.0 182.5 169.5 174.5 [5,] 190.0 190 197.0 187.0 185.0 183.0 184.0 191.0 187.0 192.0 188.0 190.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,] 149.3831 143.5014 163.6817 159.2929 159.1887 156.1887 148.5014 153.0225 [2,] 184.6169 186.4986 198.3183 198.7071 196.8113 193.8113 191.4986 204.9775 [,9] [,10] [,11] [,12] [1,] 131.2450 137.0366 124.7239 144.0845 [2,] 192.7550 194.9634 177.2761 179.9155 $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(148, 155, 167, 184.5, 190, 131, 146, 165, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5ip7y1281547553.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,] -17.0 6.0 -10.0 -2.0 -6.0 -12.0 -3.0 -21 0.0 -17.0 2.0 -3.0 [2,] -9.0 6.5 -7.5 -1.5 -4.5 -7.5 0.0 -19 5.5 -14.5 5.0 -2.0 [3,] -4.0 8.0 -6.0 2.0 -3.0 -4.0 6.0 -15 9.0 -9.0 11.0 -0.5 [4,] -1.5 10.5 -3.0 2.0 -2.0 -1.5 7.5 -9 10.5 -8.0 18.5 4.0 [5,] 0.0 16.0 0.0 3.0 1.0 1.0 9.0 1 13.0 1.0 21.0 6.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -8.478879 5.611265 -8.687327 -0.09014354 -4.492960 -7.5831032 1.521121 [2,] 0.478879 10.388735 -3.312673 4.09014354 -1.507040 -0.4168968 10.478879 [,8] [,9] [,10] [,11] [,12] [1,] -20.971839 6.014081 -12.881695 2.938018 -4.370194 [2,] -9.028161 11.985919 -5.118305 19.061982 3.370194 $out [1] -1 -11 $group [1] 2 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-17, -9, -4, -1.5, 0, 6, 6.5, 8, 10.5, 16, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6ah611281547553.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,] 183.0 178.0 161.0 151.0 130.0 118.0 93.0 NA [2,] 185.0 182.0 169.5 163.5 142.5 131.5 106.5 NA [3,] 187.0 186.0 176.5 172.0 151.5 139.5 118.5 NA [4,] 188.5 190.5 179.0 181.0 159.5 144.0 130.5 NA [5,] 190.0 197.0 183.0 184.0 164.0 151.0 148.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 185.4036 182.1231 172.167 164.0181 143.7462 133.7987 107.5534 NA [2,] 188.5964 189.8769 180.833 179.9819 159.2538 145.2013 129.4466 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(183, 185, 187, 188.5, 190, 178, 182, 186, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7ah611281547553.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,] 146.7143 151.0 147.500 [2,] 157.4286 163.5 158.875 [3,] 162.4286 168.5 163.375 [4,] 166.2857 178.5 167.625 [5,] 171.2857 181.0 171.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 158.3888 161.6584 159.3841 [2,] 166.4684 175.3416 167.3659 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(146.714285714286, 157.428571428571, 162.428571428571, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1pgpd1281547553.ps tmp/1pgpd1281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/2pgpd1281547553.ps tmp/2pgpd1281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/3pgpd1281547553.ps tmp/3pgpd1281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/4ip7y1281547553.ps tmp/4ip7y1281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/5ip7y1281547553.ps tmp/5ip7y1281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/6ah611281547553.ps tmp/6ah611281547553.png",intern=TRUE)) character(0) > try(system("convert tmp/7ah611281547553.ps tmp/7ah611281547553.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.984 1.497 2.472