R version 2.6.2 (2008-02-08) 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(10.893,10.756,10.940,10.997,10.827,10.166,10.186,10.457,10.368,10.244,10.511,10.812,10.738,10.171,9.721,9.897,9.828,9.924,10.371,10.846,10.413,10.709,10.662,10.570,10.297,10.635,10.872,10.296,10.383,10.431,10.574,10.653,10.805,10.872,10.625,10.407,10.463,10.556,10.646,10.702,11.353,11.346,11.451,11.964,12.574,13.031,13.812,14.544,14.931,14.886,16.005,17.064,15.168,16.050,15.839,15.137,14.954,15.648,15.305,15.579,16.348,15.928,16.171,15.937,15.713,15.594,15.683,16.438,17.032,17.696,17.745,19.394) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 10.893 10.738 10.297 10.463 14.931 16.348 NA [2,] 10.756 10.171 10.635 10.556 14.886 15.928 NA [3,] 10.940 9.721 10.872 10.646 16.005 16.171 NA [4,] 10.997 9.897 10.296 10.702 17.064 15.937 NA [5,] 10.827 9.828 10.383 11.353 15.168 15.713 NA [6,] 10.166 9.924 10.431 11.346 16.050 15.594 NA [7,] 10.186 10.371 10.574 11.451 15.839 15.683 NA [8,] 10.457 10.846 10.653 11.964 15.137 16.438 NA [9,] 10.368 10.413 10.805 12.574 14.954 17.032 NA [10,] 10.244 10.709 10.872 13.031 15.648 17.696 NA [11,] 10.511 10.662 10.625 13.812 15.305 17.745 NA [12,] 10.812 10.570 10.407 14.544 15.579 19.394 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.137 -0.567 0.338 0.093 -0.045 -0.420 NA [2,] 0.184 -0.450 0.237 0.090 1.119 0.243 NA [3,] 0.057 0.176 -0.576 0.056 1.059 -0.234 NA [4,] -0.170 -0.069 0.087 0.651 -1.896 -0.224 NA [5,] -0.661 0.096 0.048 -0.007 0.882 -0.119 NA [6,] 0.020 0.447 0.143 0.105 -0.211 0.089 NA [7,] 0.271 0.475 0.079 0.513 -0.702 0.755 NA [8,] -0.089 -0.433 0.152 0.610 -0.183 0.594 NA [9,] -0.124 0.296 0.067 0.457 0.694 0.664 NA [10,] 0.267 -0.047 -0.247 0.781 -0.343 0.049 NA [11,] 0.301 -0.092 -0.218 0.732 0.274 1.649 NA [12,] -0.074 -0.273 0.056 0.387 0.769 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/1dyig1209312455.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/2t49e1209312455.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/3hpc21209312455.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/4e4cm1209312455.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] [1,] 10.2970 10.1710 9.721 9.8970 9.828 9.9240 10.1860 10.457 10.3680 [2,] 10.4630 10.5560 10.646 10.2960 10.383 10.1660 10.3710 10.653 10.4130 [3,] 10.8155 10.6955 10.906 10.8495 11.090 10.8885 11.0125 11.405 11.6895 [4,] 14.9310 14.8860 16.005 15.9370 15.168 15.5940 15.6830 15.137 14.9540 [5,] 16.3480 15.9280 16.171 17.0640 15.713 16.0500 15.8390 16.438 17.0320 [,10] [,11] [,12] [1,] 10.2440 10.511 10.407 [2,] 10.7090 10.625 10.570 [3,] 11.9515 12.237 12.678 [4,] 15.6480 15.305 15.579 [5,] 17.6960 17.745 19.394 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 7.933496 7.90251 7.449272 7.210873 8.00352 7.387265 7.586088 [2,] 13.697504 13.48849 14.362728 14.488127 14.17648 14.389735 14.438912 [,8] [,9] [,10] [,11] [,12] [1,] 8.512675 8.760408 8.765685 9.218249 9.447033 [2,] 14.297325 14.618592 15.137315 15.255751 15.908967 $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(10.297, 10.463, 10.8155, 14.931, 16.348, 10.171, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/539wd1209312455.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,] -0.567 0.0900 -0.5760 -0.2240 -0.1190 0.020 0.079 -0.4330 -0.1240 -0.343 [2,] -0.420 0.0900 -0.2340 -0.2240 -0.1190 0.020 0.079 -0.1830 0.0670 -0.247 [3,] -0.091 0.2105 0.0565 -0.1195 0.0205 0.097 0.373 0.0315 0.3765 0.001 [4,] 0.093 0.2430 0.1760 0.0870 0.0960 0.143 0.513 0.5940 0.6640 0.267 [5,] 0.338 0.2430 0.1760 0.0870 0.0960 0.143 0.755 0.6100 0.6940 0.781 [,11] [,12] [1,] -0.2180 -0.273 [2,] -0.0920 -0.074 [3,] 0.2875 0.056 [4,] 0.7320 0.387 [5,] 1.6490 0.769 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.4219016 0.1118101 -0.2079632 -0.32010504 -0.1181819 0.01766103 [2,] 0.2399016 0.3091899 0.3209632 0.08110504 0.1591819 0.17633897 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.09305598 -0.4696901 -0.008584282 -0.3305466 -0.2440066 -0.2697414 [2,] 0.65294402 0.5326901 0.761584282 0.3325466 0.8190066 0.3817414 $out [1] -0.450 1.119 1.059 0.651 -1.896 -0.661 0.882 0.447 -0.211 -0.702 $group [1] 2 2 3 4 4 5 5 6 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.567, -0.419999999999998, -0.0910000000000002, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/69nj31209312455.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] [1,] 10.1660 9.7210 10.2960 10.4630 14.8860 15.5940 NA [2,] 10.3060 9.9105 10.3950 10.6740 15.0455 15.8205 NA [3,] 10.6335 10.3920 10.5995 11.4020 15.4420 16.2595 NA [4,] 10.8600 10.6855 10.7290 12.8025 15.9220 17.3640 NA [5,] 10.9970 10.8460 10.8720 14.5440 17.0640 19.3940 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 10.38082 10.03852 10.44716 10.43118 15.04222 15.5555 NA [2,] 10.88618 10.74548 10.75184 12.37282 15.84178 16.9635 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(10.166, 10.306, 10.6335, 10.86, 10.997, 9.721, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/785bz1209312455.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,] 12.15533 10.69550 12.21462 [2,] 12.26508 10.86900 12.36812 [3,] 12.43733 11.05125 12.52294 [4,] 12.86217 11.82050 12.75181 [5,] 13.55100 12.67800 12.97537 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 12.16500 10.61726 12.34794 [2,] 12.70967 11.48524 12.69794 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(12.1553333333333, 12.2650833333333, 12.4373333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1dyig1209312455.ps tmp/1dyig1209312455.png") > system("convert tmp/2t49e1209312455.ps tmp/2t49e1209312455.png") > system("convert tmp/3hpc21209312455.ps tmp/3hpc21209312455.png") > system("convert tmp/4e4cm1209312455.ps tmp/4e4cm1209312455.png") > system("convert tmp/539wd1209312455.ps tmp/539wd1209312455.png") > system("convert tmp/69nj31209312455.ps tmp/69nj31209312455.png") > system("convert tmp/785bz1209312455.ps tmp/785bz1209312455.png") > > > proc.time() user system elapsed 2.135 1.410 2.290