R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) 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(116,111,104,100,93,91,119,139,134,124,113,109,109,106,101,98,93,91,122,139,140,132,117,114,113,110,107,103,98,98,137,148,147,139,130,128,127,123,118,114,108,111,151,159,158,148,138,137,136,133,126,120,114,116,153,162,161,149,139,135,130,127,122,117,112,113,149,157,157,147,137,132,125,123,117,114,111,112,144,150,149,134,123,116,117,111,105,102,95,93,124,130,124,115,106,105,105,101,95,93,84,87,116,120,117,109,105,107,109,109,108,107,99,103,131,137,135,124,118,121,121,118,113,107,100,102,130,136,133,120,112,109,110) > 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] 133 > (np <- floor(n / par1)) [1] 11 > 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] 12 11 11 11 11 11 11 11 11 11 11 11 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 116 109 113 127 136 130 125 117 105 109 121 110 [2,] 111 106 110 123 133 127 123 111 101 109 118 NA [3,] 104 101 107 118 126 122 117 105 95 108 113 NA [4,] 100 98 103 114 120 117 114 102 93 107 107 NA [5,] 93 93 98 108 114 112 111 95 84 99 100 NA [6,] 91 91 98 111 116 113 112 93 87 103 102 NA [7,] 119 122 137 151 153 149 144 124 116 131 130 NA [8,] 139 139 148 159 162 157 150 130 120 137 136 NA [9,] 134 140 147 158 161 157 149 124 117 135 133 NA [10,] 124 132 139 148 149 147 134 115 109 124 120 NA [11,] 113 117 130 138 139 137 123 106 105 118 112 NA [12,] 109 114 128 137 135 132 116 105 107 121 109 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -5 -3 -3 -4 -3 -3 -2 -6 -4 0 -3 NA [2,] -7 -5 -3 -5 -7 -5 -6 -6 -6 -1 -5 NA [3,] -4 -3 -4 -4 -6 -5 -3 -3 -2 -1 -6 NA [4,] -7 -5 -5 -6 -6 -5 -3 -7 -9 -8 -7 NA [5,] -2 -2 0 3 2 1 1 -2 3 4 2 NA [6,] 28 31 39 40 37 36 32 31 29 28 28 NA [7,] 20 17 11 8 9 8 6 6 4 6 6 NA [8,] -5 1 -1 -1 -1 0 -1 -6 -3 -2 -3 NA [9,] -10 -8 -8 -10 -12 -10 -15 -9 -8 -11 -13 NA [10,] -11 -15 -9 -10 -10 -10 -11 -9 -4 -6 -8 NA [11,] -4 -3 -2 -1 -4 -5 -7 -1 2 3 -3 NA [12,] 0 -1 -1 -1 -5 -7 1 0 2 0 1 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/wessaorg/rcomp/tmp/146sf1305815090.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/27v8f1305815090.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3uqrg1305815090.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4fs791305815090.ps",horizontal=F,onefile=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,] 105.0 101.0 95.0 93 84.0 87.0 116.0 120.0 117.0 109 105.0 105 [2,] 109.5 109.5 104.5 101 94.0 92.0 123.0 136.5 133.5 122 112.5 109 [3,] 116.5 111.0 108.0 107 99.0 102.0 131.0 139.0 140.0 132 118.0 116 [4,] 126.0 123.0 117.5 114 109.5 111.5 146.5 153.5 153.0 143 133.5 130 [5,] 136.0 133.0 126.0 120 114.0 116.0 153.0 162.0 161.0 149 139.0 137 $n [1] 12 11 11 11 11 11 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 108.9742 104.5688 101.8070 100.8070 91.61599 92.71044 119.8049 130.9014 [2,] 124.0258 117.4312 114.1930 113.1930 106.38401 111.28956 142.1951 147.0986 [,9] [,10] [,11] [,12] [1,] 130.7104 121.9959 107.9959 105.9959 [2,] 149.2896 142.0041 128.0041 126.0041 $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(105, 109.5, 116.5, 126, 136, 101, 109.5, 111, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/5nfic1305815090.ps",horizontal=F,onefile=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,] -5 -7 -6.0 -9 -2.0 28.0 4 -6 -15.0 -11.0 -7 -1.0 [2,] -4 -6 -4.5 -7 -1.0 28.5 6 -3 -11.5 -10.5 -4 -1.0 [3,] -3 -5 -4.0 -6 1.0 31.0 8 -1 -10.0 -10.0 -3 0.0 [4,] -3 -5 -3.0 -5 2.5 36.5 10 -1 -8.5 -8.5 -1 0.5 [5,] -2 -5 -1.0 -3 4.0 40.0 11 1 -8.0 -6.0 3 2.0 $n [1] 11 11 11 11 11 11 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3.476388 -5.476388 -4.714582 -6.952776 -0.6673577 27.18890 6.094448 [2,] -2.523612 -4.523612 -3.285418 -5.047224 2.6673577 34.81110 9.905552 [,8] [,9] [,10] [,11] [,12] [1,] -1.95277585 -11.429164 -10.952776 -4.429164 -0.7145819 [2,] -0.04722415 -8.570836 -9.047224 -1.570836 0.7145819 $out [1] -6 0 -3 -1 20 17 -15 -4 -5 -7 $group [1] 1 1 2 2 7 7 10 10 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-5, -4, -3, -3, -2, -7, -6, -5, -5, -5, -6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/64t2f1305815090.ps",horizontal=F,onefile=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] [,9] [,10] [,11] [,12] [1,] 91.0 91.0 98.0 108.0 114.0 112.0 111 93.0 84.0 99.0 100.0 110 [2,] 102.0 99.5 105.0 116.0 123.0 119.5 115 103.5 94.0 107.5 108.0 110 [3,] 112.0 111.5 120.5 132.0 135.5 131.0 123 108.5 105.0 113.5 115.5 110 [4,] 121.5 127.0 138.0 149.5 151.0 148.0 139 120.5 112.5 127.5 125.5 110 [5,] 139.0 140.0 148.0 159.0 162.0 157.0 150 130.0 120.0 137.0 136.0 110 $n [1] 12 12 12 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 103.1059 98.95707 105.4485 116.7204 122.729 118.0010 112.0534 100.7462 [2,] 120.8941 124.04293 135.5515 147.2796 148.271 143.9990 133.9466 116.2538 [,9] [,10] [,11] [,12] [1,] 96.56203 104.3779 107.5181 110 [2,] 113.43797 122.6221 123.4819 110 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(91, 102, 112, 121.5, 139, 91, 99.5, 111.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/7q3qy1305815090.ps",horizontal=F,onefile=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,] 100.6364 99.00 101.7500 [2,] 108.6818 107.50 109.2500 [3,] 118.7652 116.25 118.5625 [4,] 132.5909 131.50 133.6250 [5,] 143.3636 140.00 145.0000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 107.8601 105.3034 107.4449 [2,] 129.6702 127.1966 129.6801 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(100.636363636364, 108.681818181818, 118.765151515152, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/146sf1305815090.ps tmp/146sf1305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/27v8f1305815090.ps tmp/27v8f1305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/3uqrg1305815090.ps tmp/3uqrg1305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/4fs791305815090.ps tmp/4fs791305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/5nfic1305815090.ps tmp/5nfic1305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/64t2f1305815090.ps tmp/64t2f1305815090.png",intern=TRUE)) character(0) > try(system("convert tmp/7q3qy1305815090.ps tmp/7q3qy1305815090.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.490 0.530 5.154