R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(112,118,129,99,116,168,118,129,205,147,150,267,126,129,124,97,102,127,222,214,118,141,154,226,89,77,82,97,127,121,117,117,106,112,134,169,75,108,115,85,101,108,109,124,105,95,135,164,88,85,112,87,91,87,87,142,95,108,139,159,61,82,124,93,108,75,87,103,90,108,123,129,57,65,67,71,76,67,110,118,99,85,107,141,58,65,70,86,93,74,87,73,101,100,96,157,63,115,70,66,67,83,79,77,102,116,100,135,71,60,89,74,73,91,86,74,87,87,109,137,43,69,73,77,69,76,78,70,83,65,110,132,54,55,66,65,60,65,96,55,71,63,74,106,34,47,56,53,53,55,67,52,46,51,58,91,33,40,46,45,41,55,57,54,46,52,48,77,30,35,42,48,44,45,46) > 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] 175 > (np <- floor(n / par1)) [1] 14 > 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] 15 15 15 15 15 15 15 14 14 14 14 14 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 112 126 89 75 88 61 57 58 63 71 43 54 34 [2,] 118 129 77 108 85 82 65 65 115 60 69 55 47 [3,] 129 124 82 115 112 124 67 70 70 89 73 66 56 [4,] 99 97 97 85 87 93 71 86 66 74 77 65 53 [5,] 116 102 127 101 91 108 76 93 67 73 69 60 53 [6,] 168 127 121 108 87 75 67 74 83 91 76 65 55 [7,] 118 222 117 109 87 87 110 87 79 86 78 96 67 [8,] 129 214 117 124 142 103 118 73 77 74 70 55 52 [9,] 205 118 106 105 95 90 99 101 102 87 83 71 46 [10,] 147 141 112 95 108 108 85 100 116 87 65 63 51 [11,] 150 154 134 135 139 123 107 96 100 109 110 74 58 [12,] 267 226 169 164 159 129 141 157 135 137 132 106 91 [,14] [,15] [1,] 33 30 [2,] 40 35 [3,] 46 42 [4,] 45 48 [5,] 41 44 [6,] 55 45 [7,] 57 46 [8,] 54 NA [9,] 46 NA [10,] 52 NA [11,] 48 NA [12,] 77 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 6 3 -12 33 -3 21 8 7 52 -11 26 1 13 [2,] 11 -5 5 7 27 42 2 5 -45 29 4 11 9 [3,] -30 -27 15 -30 -25 -31 4 16 -4 -15 4 -1 -3 [4,] 17 5 30 16 4 15 5 7 1 -1 -8 -5 0 [5,] 52 25 -6 7 -4 -33 -9 -19 16 18 7 5 2 [6,] -50 95 -4 1 0 12 43 13 -4 -5 2 31 12 [7,] 11 -8 0 15 55 16 8 -14 -2 -12 -8 -41 -15 [8,] 76 -96 -11 -19 -47 -13 -19 28 25 13 13 16 -6 [9,] -58 23 6 -10 13 18 -14 -1 14 0 -18 -8 5 [10,] 3 13 22 40 31 15 22 -4 -16 22 45 11 7 [11,] 117 72 35 29 20 6 34 61 35 28 22 32 33 [12,] -141 -137 -94 -76 -98 -72 -83 -94 -64 -94 -78 -72 -58 [,14] [,15] [1,] 7 5 [2,] 6 7 [3,] -1 6 [4,] -4 -4 [5,] 14 1 [6,] 2 1 [7,] -3 NA [8,] -8 NA [9,] 6 NA [10,] -4 NA [11,] 29 NA [12,] -47 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/rcomp/tmp/1jd4p1272370587.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/rcomp/tmp/2u4ls1272370587.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/rcomp/tmp/3u4ls1272370587.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/rcomp/tmp/45e3d1272370587.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,] 30.0 35.0 42.0 45.0 41.0 45.0 46.0 52 71 51.0 48.0 77 [2,] 48.5 57.5 66.5 65.5 63.5 66.0 78.5 70 83 65.0 96.0 129 [3,] 61.0 69.0 73.0 77.0 76.0 76.0 87.0 90 97 97.5 109.5 139 [4,] 81.5 96.5 113.5 90.0 101.5 99.5 109.5 124 105 112.0 135.0 164 [5,] 126.0 129.0 129.0 99.0 127.0 127.0 118.0 142 118 147.0 154.0 169 $n [1] 15 15 15 15 15 15 15 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 47.53751 53.08978 53.82615 67.00512 60.49774 62.33353 74.35342 67.19727 [2,] 74.46249 84.91022 92.17385 86.99488 91.50226 89.66647 99.64658 112.80273 [,9] [,10] [,11] [,12] [1,] 87.71 77.65318 93.03136 124.2205 [2,] 106.29 117.34682 125.96864 153.7795 $out [1] 168 222 214 205 46 46 267 226 $group [1] 6 7 8 9 9 9 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(30, 48.5, 61, 81.5, 126, 35, 57.5, 69, 96.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/55e3d1272370587.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,] -12 -5.0 -31 -8.0 -33 -5.0 -41.0 -47 -18.0 -16 20.0 -98.0 [2,] 2 4.5 -26 -2.5 -5 -2.0 -12.0 -19 -10.0 3 28.0 -94.0 [3,] 7 7.0 -3 4.0 5 2.0 -2.5 -7 2.5 14 32.5 -80.5 [4,] 17 11.0 4 11.0 15 12.5 11.0 16 13.0 22 35.0 -72.0 [5,] 33 11.0 16 30.0 25 31.0 16.0 28 23.0 45 35.0 -47.0 $n [1] 15 15 15 15 15 15 14 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.8806863 4.348297 -15.238627 -1.507382 -3.159085 -3.915337 -12.212274 [2,] 13.1193137 9.651703 9.238627 9.507382 13.159085 7.915337 7.212274 [,8] [,9] [,10] [,11] [,12] [1,] -21.779547 -7.212274 5.976818 29.54409 -89.79 [2,] 7.779547 12.212274 22.023182 35.45591 -71.21 $out [1] 52 27 42 -45 29 52 -50 95 43 55 76 -96 -58 117 72 [16] 6 61 -141 -137 $group [1] 1 2 2 2 2 5 6 6 6 7 8 8 9 11 11 11 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12, 2, 7, 17, 33, -5, 4.5, 7, 11, 11, -31, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/65e3d1272370587.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] [,12] [1,] 99 97 77.0 75.0 85.0 61.0 57.0 58.0 63.0 60.0 65.0 54.0 [2,] 117 121 93.0 98.0 87.0 84.5 67.0 71.5 68.5 73.5 69.0 57.5 [3,] 129 128 114.5 108.0 93.0 98.0 80.5 86.5 81.0 86.5 74.5 65.0 [4,] 159 184 124.0 119.5 125.5 115.5 108.5 98.0 108.5 90.0 80.5 72.5 [5,] 205 226 169.0 135.0 159.0 129.0 141.0 101.0 135.0 109.0 83.0 74.0 [,13] [,14] [,15] [1,] 46 33.0 30.0 [2,] 49 43.0 38.5 [3,] 53 47.0 44.0 [4,] 57 54.5 45.5 [5,] 67 57.0 48.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 109.8435 99.26528 100.3607 98.1937 75.43989 83.8607 61.57157 74.41317 [2,] 148.1565 156.73472 128.6393 117.8063 110.56011 112.1393 99.42843 98.58683 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 62.75573 78.97424 69.25477 58.1584 49.35115 41.75477 39.81971 [2,] 99.24427 94.02576 79.74523 71.8416 56.64885 52.24523 48.18029 $out [1] 267 164 157 137 43 110 132 96 106 34 91 77 $group [1] 1 4 8 10 11 11 11 12 12 13 13 14 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" NA Warning message: In bxp(list(stats = c(99, 117, 129, 159, 205, 97, 121, 128, 184, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7xn2y1272370587.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,] 66.26667 61.00 65.000 [2,] 79.03333 74.50 80.125 [3,] 90.73333 82.00 90.250 [4,] 98.42857 97.25 95.375 [5,] 109.78571 109.50 115.875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 81.88704 71.62357 83.29437 [2,] 99.57963 92.37643 97.20563 $out [1] 149.2857 139.0000 146.2500 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(66.2666666666667, 79.0333333333333, 90.7333333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1jd4p1272370587.ps tmp/1jd4p1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/2u4ls1272370587.ps tmp/2u4ls1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/3u4ls1272370587.ps tmp/3u4ls1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/45e3d1272370587.ps tmp/45e3d1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/55e3d1272370587.ps tmp/55e3d1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/65e3d1272370587.ps tmp/65e3d1272370587.png",intern=TRUE)) character(0) > try(system("convert tmp/7xn2y1272370587.ps tmp/7xn2y1272370587.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.217 0.946 1.524