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(84.9,81.9,95.9,81,89.2,102.5,89.8,88.8,83.2,90.2,100.4,187.1,87.6,85.4,86.1,86.7,89.1,103.7,86.9,85.2,80.8,91.2,102.8,182.5,80.9,83.1,88.3,86.6,93,105.3,93.8,86.4,87,96.7,100.5,196.7,86.8,88.2,93.8,85,90.4,115.9,94.9,87.7,91.7,95.9,106.8,204.5,90.2,90.5,93.2,97.8,99.4,120,108.2,98.5,104.3,102.9,111.1,188.1,93.8,94.5,112.4,102.5,115.8,136.5,122.1,110.6,116.4,112.6,121.5,199.3) > 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,] 84.9 87.6 80.9 86.8 90.2 93.8 NA [2,] 81.9 85.4 83.1 88.2 90.5 94.5 NA [3,] 95.9 86.1 88.3 93.8 93.2 112.4 NA [4,] 81.0 86.7 86.6 85.0 97.8 102.5 NA [5,] 89.2 89.1 93.0 90.4 99.4 115.8 NA [6,] 102.5 103.7 105.3 115.9 120.0 136.5 NA [7,] 89.8 86.9 93.8 94.9 108.2 122.1 NA [8,] 88.8 85.2 86.4 87.7 98.5 110.6 NA [9,] 83.2 80.8 87.0 91.7 104.3 116.4 NA [10,] 90.2 91.2 96.7 95.9 102.9 112.6 NA [11,] 100.4 102.8 100.5 106.8 111.1 121.5 NA [12,] 187.1 182.5 196.7 204.5 188.1 199.3 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3.0 -2.2 2.2 1.4 0.3 0.7 NA [2,] 14.0 0.7 5.2 5.6 2.7 17.9 NA [3,] -14.9 0.6 -1.7 -8.8 4.6 -9.9 NA [4,] 8.2 2.4 6.4 5.4 1.6 13.3 NA [5,] 13.3 14.6 12.3 25.5 20.6 20.7 NA [6,] -12.7 -16.8 -11.5 -21.0 -11.8 -14.4 NA [7,] -1.0 -1.7 -7.4 -7.2 -9.7 -11.5 NA [8,] -5.6 -4.4 0.6 4.0 5.8 5.8 NA [9,] 7.0 10.4 9.7 4.2 -1.4 -3.8 NA [10,] 10.2 11.6 3.8 10.9 8.2 8.9 NA [11,] 86.7 79.7 96.2 97.7 77.0 77.8 NA [12,] -99.5 -101.6 -109.9 -114.3 -94.3 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/1onse1290006779.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/2onse1290006779.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/3ze9z1290006779.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/4ze9z1290006779.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,] 80.9 81.9 86.1 81.00 89.1 102.5 86.90 85.20 80.80 90.2 100.4 182.5 [2,] 84.9 83.1 88.3 85.00 89.2 103.7 89.80 86.40 83.20 91.2 100.5 187.1 [3,] 87.2 86.8 93.5 86.65 91.7 110.6 94.35 88.25 89.35 96.3 104.8 192.4 [4,] 90.2 90.5 95.9 97.80 99.4 120.0 108.20 98.50 104.30 102.9 111.1 199.3 [5,] 93.8 94.5 95.9 102.50 99.4 136.5 122.10 110.60 116.40 112.6 121.5 204.5 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 83.78133 82.02676 88.59775 78.39359 85.12067 100.0860 82.4814 80.44511 [2,] 90.61867 91.57324 98.40225 94.90641 98.27933 121.1140 106.2186 96.05489 [,9] [,10] [,11] [,12] [1,] 75.73982 88.75312 97.96266 184.5306 [2,] 102.96018 103.84688 111.63734 200.2694 $out [1] 112.4 115.8 $group [1] 3 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(80.9, 84.9, 87.2, 90.2, 93.8, 81.9, 83.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5sn821290006779.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,] -3.0 0.7 -14.90 1.6 12.3 -21.00 -11.5 -5.6 -3.8 8.20 77.0 -114.3 [2,] -2.2 2.7 -9.90 2.4 13.3 -16.80 -9.7 -4.4 -1.4 8.20 77.8 -109.9 [3,] 0.5 5.4 -5.25 5.9 17.6 -13.55 -7.3 2.3 5.6 9.55 83.2 -101.6 [4,] 1.4 14.0 0.60 8.2 20.7 -11.80 -1.7 5.8 9.7 10.90 96.2 -99.5 [5,] 2.2 17.9 4.60 13.3 25.5 -11.50 -1.0 5.8 10.4 11.60 97.7 -94.3 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.822116 -1.888865 -12.022839 2.158813 12.82676 -16.77516 -12.460258 [2,] 2.822116 12.688865 1.522839 9.641187 22.37324 -10.32484 -2.139742 [,8] [,9] [,10] [,11] [,12] [1,] -4.279329 -1.559859 7.808413 71.3314 -108.94861 [2,] 8.879329 12.759859 11.291587 95.0686 -94.25139 $out [1] 3.8 $group [1] 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3, -2.19999999999999, 0.5, 1.40000000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6sn821290006779.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,] 81.00 80.80 80.90 85.00 90.20 93.80 NA [2,] 84.05 85.75 86.50 87.95 95.50 106.55 NA [3,] 89.50 87.25 90.65 92.75 101.15 114.20 NA [4,] 98.15 97.00 98.60 101.35 109.65 121.80 NA [5,] 102.50 103.70 105.30 115.90 120.00 136.50 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 83.0689 82.1188 85.13111 86.63817 94.69609 107.2444 NA [2,] 95.9311 92.3812 96.16889 98.86183 107.60391 121.1556 NA $out [1] 187.1 182.5 196.7 204.5 188.1 199.3 $group [1] 1 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(81, 84.05, 89.5, 98.15, 102.5, 80.8, 85.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7sn821290006779.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,] 87.26667 86.650 86.80000 [2,] 91.40000 87.725 90.80625 [3,] 95.55000 92.600 93.15000 [4,] 103.23333 100.550 101.69375 [5,] 113.98333 110.600 111.53750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 90.15274 86.75043 88.18414 [2,] 100.94726 98.44957 98.11586 $out [1] 193.0333 192.4000 193.0000 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(87.2666666666667, 91.4, 95.55, 103.233333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1onse1290006779.ps tmp/1onse1290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/2onse1290006779.ps tmp/2onse1290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/3ze9z1290006779.ps tmp/3ze9z1290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/4ze9z1290006779.ps tmp/4ze9z1290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/5sn821290006779.ps tmp/5sn821290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/6sn821290006779.ps tmp/6sn821290006779.png",intern=TRUE)) character(0) > try(system("convert tmp/7sn821290006779.ps tmp/7sn821290006779.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.944 1.478 2.179