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. 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(9.733,9.259,9.864,9.215,10.103,9.380,9.896,10.117,9.451,9.700,9.081,9.084,9.743,8.587,9.731,9.563,9.998,9.437,10.038,9.918,9.252,9.737,9.035,9.133,9.487,8.700,9.627,8.947,9.283,8.829,9.947,9.628,9.318,9.605,8.640,9.214,9.567,8.547,9.185,9.470,9.123,9.278,10.170,9.434,9.655,9.429,8.739,9.552,9.687,9.019,9.672,9.206,9.069,9.788,10.312,10.105,9.863,9.656,9.295,9.946,9.701,9.049,10.190,9.706,9.765,9.893,9.994,10.433,10.073,10.112,9.266,9.820,10.097,9.115,10.411,9.678,10.408,10.153,10.368,10.581,10.597,10.680,9.738,9.556) > 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,] 9.733 9.743 9.487 9.567 9.687 9.701 10.097 NA [2,] 9.259 8.587 8.700 8.547 9.019 9.049 9.115 NA [3,] 9.864 9.731 9.627 9.185 9.672 10.190 10.411 NA [4,] 9.215 9.563 8.947 9.470 9.206 9.706 9.678 NA [5,] 10.103 9.998 9.283 9.123 9.069 9.765 10.408 NA [6,] 9.380 9.437 8.829 9.278 9.788 9.893 10.153 NA [7,] 9.896 10.038 9.947 10.170 10.312 9.994 10.368 NA [8,] 10.117 9.918 9.628 9.434 10.105 10.433 10.581 NA [9,] 9.451 9.252 9.318 9.655 9.863 10.073 10.597 NA [10,] 9.700 9.737 9.605 9.429 9.656 10.112 10.680 NA [11,] 9.081 9.035 8.640 8.739 9.295 9.266 9.738 NA [12,] 9.084 9.133 9.214 9.552 9.946 9.820 9.556 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -0.474 -1.156 -0.787 -1.020 -0.668 -0.652 -0.982 NA [2,] 0.605 1.144 0.927 0.638 0.653 1.141 1.296 NA [3,] -0.649 -0.168 -0.680 0.285 -0.466 -0.484 -0.733 NA [4,] 0.888 0.435 0.336 -0.347 -0.137 0.059 0.730 NA [5,] -0.723 -0.561 -0.454 0.155 0.719 0.128 -0.255 NA [6,] 0.516 0.601 1.118 0.892 0.524 0.101 0.215 NA [7,] 0.221 -0.120 -0.319 -0.736 -0.207 0.439 0.213 NA [8,] -0.666 -0.666 -0.310 0.221 -0.242 -0.360 0.016 NA [9,] 0.249 0.485 0.287 -0.226 -0.207 0.039 0.083 NA [10,] -0.619 -0.702 -0.965 -0.690 -0.361 -0.846 -0.942 NA [11,] 0.003 0.098 0.574 0.813 0.651 0.554 -0.182 NA [12,] 0.659 0.354 0.353 0.135 -0.245 0.277 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/1hj6e1227894450.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/2pbc91227894450.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/34uxa1227894450.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/4lhqj1227894450.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] [1,] 9.487 8.5470 9.1850 8.9470 9.0690 8.8290 9.8960 9.434 9.2520 9.4290 [2,] 9.627 8.6435 9.6495 9.2105 9.2030 9.3290 9.9705 9.773 9.3845 9.6305 [3,] 9.701 9.0190 9.7310 9.4700 9.7650 9.4370 10.0380 10.105 9.6550 9.7000 [4,] 9.738 9.0820 10.0270 9.6205 10.0505 9.8405 10.2410 10.275 9.9680 9.9245 [5,] 9.743 9.2590 10.4110 9.7060 10.4080 10.1530 10.3680 10.581 10.5970 10.1120 [,11] [,12] [1,] 8.6400 9.0840 [2,] 8.8870 9.1735 [3,] 9.0810 9.5520 [4,] 9.2805 9.6880 [5,] 9.7380 9.9460 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9.634713 8.757135 9.505563 9.225155 9.258887 9.13154 9.876462 9.805214 [2,] 9.767287 9.280865 9.956437 9.714845 10.271113 9.74246 10.199538 10.404786 [,9] [,10] [,11] [,12] [1,] 9.306543 9.524428 8.846008 9.244749 [2,] 10.003457 9.875572 9.315992 9.859251 $out [1] 10.097 10.680 $group [1] 1 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9.487, 9.627, 9.701, 9.738, 9.743, 8.547, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5aggi1227894450.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] [1,] -1.156 0.6050 -0.7330 -0.3470 -0.7230 0.1010 -0.736 -0.666 -0.226 -0.9650 [2,] -1.001 0.6455 -0.6645 -0.0390 -0.5075 0.3655 -0.263 -0.513 -0.084 -0.8940 [3,] -0.787 0.9270 -0.4840 0.3360 -0.2550 0.5240 -0.120 -0.310 0.083 -0.7020 [4,] -0.660 1.1425 -0.3170 0.5825 0.1415 0.7465 0.217 -0.113 0.268 -0.6545 [5,] -0.474 1.2960 -0.1680 0.8880 0.7190 1.1180 0.439 0.221 0.485 -0.3610 [,11] [,12] [1,] -0.1820 0.135 [2,] 0.0505 0.135 [3,] 0.5540 0.315 [4,] 0.6125 0.354 [5,] 0.8130 0.659 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.9906397 0.6301996 -0.6915214 -0.03514977 -0.6425723 0.2964729 [2,] -0.5833603 1.2238004 -0.2764786 0.70714977 0.1325723 0.7515271 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.4066483 -0.54887355 -0.1272087 -0.8450255 0.2183827 0.1737379 [2,] 0.1666483 -0.07112645 0.2932087 -0.5589745 0.8896173 0.4562621 $out [1] 0.285 -0.245 $group [1] 3 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1.156, -1.001, -0.787000000000001, -0.66, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/65hlh1227894450.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,] 9.0810 8.5870 8.6400 8.5470 9.0190 9.2660 9.1150 NA [2,] 9.2370 9.1925 8.8880 9.1540 9.2505 9.7035 9.7080 NA [3,] 9.5755 9.6470 9.3005 9.4315 9.6795 9.8565 10.2605 NA [4,] 9.8800 9.8305 9.6160 9.5595 9.9045 10.0925 10.4960 NA [5,] 10.1170 10.0380 9.9470 9.6550 10.3120 10.4330 10.6800 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9.282223 9.356004 8.968454 9.246549 9.381206 9.679074 9.901088 NA [2,] 9.868777 9.937996 9.632546 9.616451 9.977794 10.033926 10.619912 NA $out [1] 10.170 9.049 $group [1] 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(9.081, 9.237, 9.5755, 9.88, 10.117, 8.587, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7oy5l1227894450.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,] 8.896571 9.0190 8.862750 [2,] 9.435000 9.4535 9.423125 [3,] 9.697429 9.6775 9.651500 [4,] 9.828500 9.7480 9.807875 [5,] 10.103571 10.1050 10.105750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 9.517951 9.543177 9.476013 [2,] 9.876907 9.811823 9.826987 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(8.89657142857143, 9.435, 9.69742857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1hj6e1227894450.ps tmp/1hj6e1227894450.png") > system("convert tmp/2pbc91227894450.ps tmp/2pbc91227894450.png") > system("convert tmp/34uxa1227894450.ps tmp/34uxa1227894450.png") > system("convert tmp/4lhqj1227894450.ps tmp/4lhqj1227894450.png") > system("convert tmp/5aggi1227894450.ps tmp/5aggi1227894450.png") > system("convert tmp/65hlh1227894450.ps tmp/65hlh1227894450.png") > system("convert tmp/7oy5l1227894450.ps tmp/7oy5l1227894450.png") > > > proc.time() user system elapsed 1.883 1.483 2.094