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(80,79,78,76,74,73,74,76,77,77,78,80,90,90,89,82,78,76,74,78,81,82,88,99,117,113,106,100,97,96,100,104,104,111,117,118,140,147,134,126,116,114,120,122,117,119,132,134,154,152,132,130,123,129,124,128,128,129,141,138,155,160,142,133,131,140,134,134,134,136,145,137,152,168,160,157,147,161,159,164,163,158,175,163) > 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,] 80 90 117 140 154 155 152 NA [2,] 79 90 113 147 152 160 168 NA [3,] 78 89 106 134 132 142 160 NA [4,] 76 82 100 126 130 133 157 NA [5,] 74 78 97 116 123 131 147 NA [6,] 73 76 96 114 129 140 161 NA [7,] 74 74 100 120 124 134 159 NA [8,] 76 78 104 122 128 134 164 NA [9,] 77 81 104 117 128 134 163 NA [10,] 77 82 111 119 129 136 158 NA [11,] 78 88 117 132 141 145 175 NA [12,] 80 99 118 134 138 137 163 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 0 -4 7 -2 5 16 NA [2,] -1 -1 -7 -13 -20 -18 -8 NA [3,] -2 -7 -6 -8 -2 -9 -3 NA [4,] -2 -4 -3 -10 -7 -2 -10 NA [5,] -1 -2 -1 -2 6 9 14 NA [6,] 1 -2 4 6 -5 -6 -2 NA [7,] 2 4 4 2 4 0 5 NA [8,] 1 3 0 -5 0 0 -1 NA [9,] 0 1 7 2 1 2 -5 NA [10,] 1 6 6 13 12 9 17 NA [11,] 2 11 1 2 -3 -8 -12 NA [12,] 10 18 22 20 17 15 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/167bd1282041806.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/267bd1282041806.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/367bd1282041806.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/4hhtg1282041806.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,] 80.0 79.0 78.0 76.0 74.0 73.0 74 76 77.0 77.0 78.0 80.0 [2,] 103.5 101.5 97.5 91.0 87.5 86.0 87 91 92.5 96.5 102.5 108.5 [3,] 140.0 147.0 132.0 126.0 116.0 114.0 120 122 117.0 119.0 132.0 134.0 [4,] 153.0 156.0 138.0 131.5 127.0 134.5 129 131 131.0 132.5 143.0 137.5 [5,] 155.0 168.0 160.0 157.0 147.0 161.0 159 164 163.0 158.0 175.0 163.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 110.4394 114.4535 107.8141 101.8141 92.41124 85.03658 94.91828 [2,] 169.5606 179.5465 156.1859 150.1859 139.58876 142.96342 145.08172 [,8] [,9] [,10] [,11] [,12] [1,] 98.11265 94.00842 97.50138 107.8141 116.6817 [2,] 145.88735 139.99158 140.49862 156.1859 151.3183 $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(80, 103.5, 140, 153, 155, 79, 101.5, 147, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5hhtg1282041806.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,] -4.0 -20.0 -9.0 -10.0 -2.0 -6.0 0 -1.0 0.0 1.0 -12.0 10.0 [2,] -1.5 -15.5 -7.5 -8.5 -1.5 -3.5 2 -0.5 0.5 6.0 -5.5 15.0 [3,] 0.0 -8.0 -6.0 -4.0 -1.0 -2.0 4 0.0 1.0 9.0 1.0 17.5 [4,] 6.0 -4.0 -2.5 -2.5 7.5 2.5 4 0.5 2.0 12.5 2.0 20.0 [5,] 16.0 -1.0 -2.0 -2.0 14.0 6.0 5 1.0 2.0 17.0 11.0 22.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4.478879 -14.867614 -8.985919 -7.5831032 -6.374655 -5.583103 2.805632 [2,] 4.478879 -1.132386 -3.014081 -0.4168968 4.374655 1.583103 5.194368 [,8] [,9] [,10] [,11] [,12] [1,] -0.5971839 0.1042242 5.118305 -3.478879 14.27484 [2,] 0.5971839 1.8957758 12.881695 5.478879 20.72516 $out [1] 3 -5 7 -5 $group [1] 8 8 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4, -1.5, 0, 6, 16, -20, -15.5, -8, -4, -1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/69qs11282041806.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,] 73.0 74.0 96 114 123.0 131.0 152.0 NA [2,] 75.0 78.0 100 118 128.0 134.0 157.5 NA [3,] 77.0 82.0 105 124 129.5 136.5 160.5 NA [4,] 78.5 89.5 115 134 139.5 143.5 163.5 NA [5,] 80.0 99.0 118 147 154.0 155.0 168.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 75.40363 76.75477 98.1584 116.7023 124.2548 132.1670 157.7634 NA [2,] 78.59637 87.24523 111.8416 131.2977 134.7452 140.8330 163.2366 NA $out [1] 160 147 175 $group [1] 6 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(73, 75, 77, 78.5, 80, 74, 78, 82, 89.5, 99, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/79qs11282041806.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,] 109.4286 114 107.250 [2,] 113.7857 118 110.625 [3,] 115.5714 124 113.125 [4,] 124.6429 133 122.875 [5,] 129.8571 147 128.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 110.6194 117.1584 107.5377 [2,] 120.5234 130.8416 118.7123 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(109.428571428571, 113.785714285714, 115.571428571429, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/167bd1282041806.ps tmp/167bd1282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/267bd1282041806.ps tmp/267bd1282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/367bd1282041806.ps tmp/367bd1282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/4hhtg1282041806.ps tmp/4hhtg1282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/5hhtg1282041806.ps tmp/5hhtg1282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/69qs11282041806.ps tmp/69qs11282041806.png",intern=TRUE)) character(0) > try(system("convert tmp/79qs11282041806.ps tmp/79qs11282041806.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.993 1.500 2.349