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(75,74,73,71,91,90,75,65,66,66,67,69,75,79,75,77,100,100,94,83,83,83,84,88,89,98,94,84,111,98,98,83,79,78,80,94,98,104,94,90,115,104,114,99,96,98,104,111,117,125,117,118,151,145,155,133,124,125,131,133,136,141,130,137,177,183,191,166,156,153,164,164,168,173,164,165,205,207,215,190,169,175,188,188,196,201,194,197,237,236,244,222,195,199,207,204,212,222,214,217,258,256,251,223,198,206,214,212,227,238,228,235,275,278,278,251,225,232,238,239) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 75 75 89 98 117 136 168 196 212 227 NA [2,] 74 79 98 104 125 141 173 201 222 238 NA [3,] 73 75 94 94 117 130 164 194 214 228 NA [4,] 71 77 84 90 118 137 165 197 217 235 NA [5,] 91 100 111 115 151 177 205 237 258 275 NA [6,] 90 100 98 104 145 183 207 236 256 278 NA [7,] 75 94 98 114 155 191 215 244 251 278 NA [8,] 65 83 83 99 133 166 190 222 223 251 NA [9,] 66 83 79 96 124 156 169 195 198 225 NA [10,] 66 83 78 98 125 153 175 199 206 232 NA [11,] 67 84 80 104 131 164 188 207 214 238 NA [12,] 69 88 94 111 133 164 188 204 212 239 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1 4 9 6 8 5 5 5 10 11 NA [2,] -1 -4 -4 -10 -8 -11 -9 -7 -8 -10 NA [3,] -2 2 -10 -4 1 7 1 3 3 7 NA [4,] 20 23 27 25 33 40 40 40 41 40 NA [5,] -1 0 -13 -11 -6 6 2 -1 -2 3 NA [6,] -15 -6 0 10 10 8 8 8 -5 0 NA [7,] -10 -11 -15 -15 -22 -25 -25 -22 -28 -27 NA [8,] 1 0 -4 -3 -9 -10 -21 -27 -25 -26 NA [9,] 0 0 -1 2 1 -3 6 4 8 7 NA [10,] 1 1 2 6 6 11 13 8 8 6 NA [11,] 2 4 14 7 2 0 0 -3 -2 1 NA [12,] 6 1 4 6 3 4 8 8 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/12pyp1281872740.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/22pyp1281872740.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/3vyxs1281872740.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/4vyxs1281872740.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,] 75.0 74 73.0 71.0 91 90 75 65.0 66 66 67.0 69.0 [2,] 89.0 98 94.0 84.0 111 100 98 83.0 83 83 84.0 94.0 [3,] 126.5 133 123.5 127.5 164 164 173 149.5 140 139 147.5 148.5 [4,] 196.0 201 194.0 197.0 237 236 244 222.0 195 199 207.0 204.0 [5,] 227.0 238 228.0 235.0 275 278 278 251.0 225 232 238.0 239.0 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 73.03853 81.5371 73.53601 71.0407 101.0454 96.04898 100.0526 [2,] 179.96147 184.4629 173.46399 183.9593 226.9546 231.95102 245.9474 [,8] [,9] [,10] [,11] [,12] [1,] 80.05006 84.04033 81.04178 86.0443 93.53961 [2,] 218.94994 195.95967 196.95822 208.9557 203.46039 $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(75, 89, 126.5, 196, 227, 74, 98, 133, 201, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5vyxs1281872740.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,] -1.0 -11 -4.0 20.0 -13 -15 -28 -27.0 -3.0 1 -3.0 1 [2,] 5.0 -10 -2.0 25.0 -6 -5 -25 -25.0 0.0 2 0.0 4 [3,] 5.5 -8 1.5 36.5 -1 4 -22 -9.5 1.5 6 1.5 6 [4,] 9.0 -4 3.0 40.0 2 8 -15 -3.0 6.0 8 4.0 8 [5,] 11.0 -1 7.0 41.0 6 10 -10 1.0 8.0 13 7.0 8 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3.501441 -10.997839 -0.9981994 29.0054 -4.997119 -2.495318 -26.9964 [2,] 7.498559 -5.002161 3.9981994 43.9946 2.997119 10.495318 -17.0036 [,8] [,9] [,10] [,11] [,12] [1,] -20.492077 -1.497839 3.002161 -0.4985595 3.893333 [2,] 1.492077 4.497839 8.997839 3.4985595 8.106667 $out [1] -10 14 15 $group [1] 3 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1, 5, 5.5, 9, 11, -11, -10, -8, -4, -1, -4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/66qfd1281872740.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] [,9] [,10] [,11] [1,] 65.0 75 78.0 90.0 117 130.0 164.0 194.0 198.0 225 NA [2,] 66.5 78 81.5 97.0 121 139.0 168.5 196.5 212.0 230 NA [3,] 72.0 83 91.5 101.5 128 160.0 181.5 202.5 215.5 238 NA [4,] 75.0 91 98.0 107.5 139 171.5 197.5 229.0 237.0 263 NA [5,] 75.0 100 111.0 115.0 155 191.0 215.0 244.0 258.0 278 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 68.12309 77.07061 83.97424 96.71088 119.7901 145.1765 168.2729 187.6765 [2,] 75.87691 88.92939 99.02576 106.28912 136.2099 174.8235 194.7271 217.3235 [,9] [,10] [,11] [1,] 204.0973 222.9485 NA [2,] 226.9027 253.0515 NA $out [1] 91 90 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(65, 66.5, 72, 75, 75, 75, 78, 83, 91, 100, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/76qfd1281872740.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,] 138.3 123.50 137.2500 [2,] 139.2 130.25 140.0000 [3,] 146.6 143.75 146.1875 [4,] 160.6 156.75 157.6875 [5,] 172.0 173.00 170.5000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 136.8393 131.6632 138.1201 [2,] 156.3607 155.8368 154.2549 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(138.3, 139.2, 146.6, 160.6, 172, 123.5, 130.25, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/12pyp1281872740.ps tmp/12pyp1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/22pyp1281872740.ps tmp/22pyp1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/3vyxs1281872740.ps tmp/3vyxs1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/4vyxs1281872740.ps tmp/4vyxs1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/5vyxs1281872740.ps tmp/5vyxs1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/66qfd1281872740.ps tmp/66qfd1281872740.png",intern=TRUE)) character(0) > try(system("convert tmp/76qfd1281872740.ps tmp/76qfd1281872740.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.009 1.525 2.220