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(96.1,96.5,96.9,97.8,98.9,100.2,101.2,101,101.6,102.4,103.7,103.7,104.6,104.5,104.5,105.6,106.1,107.6,107.7,108.3,108.1,108.1,108,108.2,108.9,109.8,109.9,109.8,110.9,111.1,112.2,112.7,114.6,114.2,114.7,114.7,116,116.3,116.4,116.6,118.1,117.2,108.3,109.5,110.5,110.6,111.2,111.1,111,112.4,112.5,112.4,111.8,111.6,112.9,112.8,113.7,113.8,114,113.8,113.9,114.4,114.4,114.5,113.8,114.3,115,115.4,115.3,114.9,114.3,114.5,115.5,115.8,115.8,116,114.9,114.1,114.1,113.5,115,114.7,115.4,116.1,116.6,117.2,118.2,118,117.7,118.5,117.5,118,117.7,116.3,115,115.7,113.6,114.8,114.9,117.3,117.3,117.7,120,119.6,119.2,117.3,117.5,119,112.5,118.9,118.4,119.4,120.6,118.6,122,122.6,120.6,117.4,116.4,122.2,121,122.4,124.9,126.1,124.5,123.2,126.4,123.9,116,126.6,125.9,126.6,116.7,126.4,129,128.7,128.4,129.2,133.3,128.9,132.7,127.7,131.8,133.9) > 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] 144 > (np <- floor(n / par1)) [1] 12 > 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] 12 12 12 12 12 12 12 12 12 12 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 96.1 104.6 108.9 116.0 111.0 113.9 115.5 116.6 113.6 112.5 121.0 116.7 [2,] 96.5 104.5 109.8 116.3 112.4 114.4 115.8 117.2 114.8 118.9 122.4 126.4 [3,] 96.9 104.5 109.9 116.4 112.5 114.4 115.8 118.2 114.9 118.4 124.9 129.0 [4,] 97.8 105.6 109.8 116.6 112.4 114.5 116.0 118.0 117.3 119.4 126.1 128.7 [5,] 98.9 106.1 110.9 118.1 111.8 113.8 114.9 117.7 117.3 120.6 124.5 128.4 [6,] 100.2 107.6 111.1 117.2 111.6 114.3 114.1 118.5 117.7 118.6 123.2 129.2 [7,] 101.2 107.7 112.2 108.3 112.9 115.0 114.1 117.5 120.0 122.0 126.4 133.3 [8,] 101.0 108.3 112.7 109.5 112.8 115.4 113.5 118.0 119.6 122.6 123.9 128.9 [9,] 101.6 108.1 114.6 110.5 113.7 115.3 115.0 117.7 119.2 120.6 116.0 132.7 [10,] 102.4 108.1 114.2 110.6 113.8 114.9 114.7 116.3 117.3 117.4 126.6 127.7 [11,] 103.7 108.0 114.7 111.2 114.0 114.3 115.4 115.0 117.5 116.4 125.9 131.8 [12,] 103.7 108.2 114.7 111.1 113.8 114.5 116.1 115.7 119.0 122.2 126.6 133.9 [,13] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.4 -0.1 0.9 0.3 1.4 0.5 0.3 0.6 1.2 6.4 1.4 9.7 NA [2,] 0.4 0.0 0.1 0.1 0.1 0.0 0.0 1.0 0.1 -0.5 2.5 2.6 NA [3,] 0.9 1.1 -0.1 0.2 -0.1 0.1 0.2 -0.2 2.4 1.0 1.2 -0.3 NA [4,] 1.1 0.5 1.1 1.5 -0.6 -0.7 -1.1 -0.3 0.0 1.2 -1.6 -0.3 NA [5,] 1.3 1.5 0.2 -0.9 -0.2 0.5 -0.8 0.8 0.4 -2.0 -1.3 0.8 NA [6,] 1.0 0.1 1.1 -8.9 1.3 0.7 0.0 -1.0 2.3 3.4 3.2 4.1 NA [7,] -0.2 0.6 0.5 1.2 -0.1 0.4 -0.6 0.5 -0.4 0.6 -2.5 -4.4 NA [8,] 0.6 -0.2 1.9 1.0 0.9 -0.1 1.5 -0.3 -0.4 -2.0 -7.9 3.8 NA [9,] 0.8 0.0 -0.4 0.1 0.1 -0.4 -0.3 -1.4 -1.9 -3.2 10.6 -5.0 NA [10,] 1.3 -0.1 0.5 0.6 0.2 -0.6 0.7 -1.3 0.2 -1.0 -0.7 4.1 NA [11,] 0.0 0.2 0.0 -0.1 -0.2 0.2 0.7 0.7 1.5 5.8 0.7 2.1 NA [12,] 0.9 0.7 1.3 -0.1 0.1 1.0 0.5 -2.1 -6.5 -1.2 -9.9 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/19w1v1289817033.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/29w1v1289817033.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/32oig1289817033.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/42oig1289817033.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,] 104.60 104.50 104.50 105.6 106.10 107.60 101.20 101.00 108.10 108.10 [2,] 109.95 111.10 111.20 111.1 111.35 111.35 110.25 111.10 112.10 112.20 [3,] 113.75 115.30 115.35 116.3 116.10 115.75 114.55 114.45 115.15 114.80 [4,] 116.30 118.05 118.30 118.7 119.35 118.55 121.00 121.10 118.45 117.35 [5,] 121.00 126.40 124.90 128.7 128.40 129.20 133.30 128.90 120.60 117.40 [,11] [,12] [1,] 108.00 103.70 [2,] 112.60 112.45 [3,] 114.85 115.20 [4,] 116.95 120.60 [5,] 117.50 126.60 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 110.8537 112.1301 112.1116 112.8336 112.4511 112.4660 109.6469 109.8889 [2,] 116.6463 118.4699 118.5884 119.7664 119.7489 119.0340 119.4531 119.0111 [,9] [,10] [,11] [,12] [1,] 112.2537 112.4511 112.8659 111.4827 [2,] 118.0463 117.1489 116.8341 118.9173 $out [1] 96.1 96.5 96.9 129.0 97.8 98.9 100.2 101.6 132.7 102.4 126.6 127.7 [13] 103.7 125.9 131.8 133.9 $group [1] 1 2 3 3 4 5 6 9 9 10 10 10 11 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(104.6, 109.95, 113.75, 116.3, 121, 104.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/52oig1289817033.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,] -0.10 -0.5 -0.30 -1.60 -2.00 -1.00 -0.60 -2.00 -3.20 -1.30 -0.20 -2.10 [2,] 0.35 0.0 -0.10 -0.65 -0.85 0.05 -0.50 -0.35 -1.65 -0.65 0.00 -1.65 [3,] 0.75 0.1 0.20 -0.15 0.30 1.05 0.15 0.25 -0.35 0.20 0.45 0.10 [4,] 1.40 0.7 1.05 1.10 0.80 2.75 0.55 1.25 0.10 0.65 1.10 0.80 [5,] 1.40 1.0 2.40 1.50 1.50 4.10 1.20 1.90 0.80 1.30 2.10 1.30 $n [1] 12 12 12 12 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.2710880 -0.2192747 -0.3245227 -0.9481867 -0.4525761 -0.1814881 [2,] 1.2289120 0.4192747 0.7245227 0.6481867 1.0525761 2.2814881 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.3289120 -0.4797707 -1.1481867 -0.3929387 -0.05171738 -1.067150 [2,] 0.6289120 0.9797707 0.4481867 0.7929387 0.95171738 1.267150 $out [1] 6.4 9.7 2.5 2.6 -8.9 -2.5 -4.4 -7.9 3.8 10.6 -5.0 4.1 5.8 -6.5 -9.9 $group [1] 1 1 2 2 6 7 7 8 8 9 9 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0999999999999943, 0.350000000000001, 0.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6vfzj1289817033.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,] 96.10 104.50 108.90 108.30 111.00 113.80 113.5 115.00 113.6 116.40 121.00 [2,] 97.35 105.10 109.85 110.55 112.10 114.30 114.4 116.45 116.1 117.90 122.80 [3,] 100.60 107.65 111.65 113.60 112.65 114.45 115.2 117.60 117.4 119.15 124.70 [4,] 102.00 108.10 114.40 116.50 113.75 114.95 115.8 118.00 119.1 121.30 126.25 [5,] 103.70 108.30 114.70 118.10 114.00 115.40 116.1 118.50 120.0 122.60 126.60 [,12] [,13] [1,] 126.40 NA [2,] 128.05 NA [3,] 128.95 NA [4,] 132.25 NA [5,] 133.90 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 98.4791 106.2817 109.5747 110.8862 111.8974 114.1535 114.5615 116.8930 [2,] 102.7209 109.0183 113.7253 116.3138 113.4026 114.7465 115.8385 118.3070 [,9] [,10] [,11] [,12] [,13] [1,] 116.0317 117.5992 123.1264 127.0344 NA [2,] 118.7683 120.7008 126.2736 130.8656 NA $out [1] 112.5 116.0 116.7 $group [1] 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(96.1, 97.35, 100.6, 102, 103.7, 104.5, 105.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7vfzj1289817033.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,] 114.1167 113.750 114.6875 [2,] 114.9167 114.675 114.9938 [3,] 115.3042 115.175 115.1000 [4,] 115.5875 115.550 115.6750 [5,] 115.8833 116.300 116.4625 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 114.9982 114.7759 114.7893 [2,] 115.6101 115.5741 115.4107 $out [1] 112.2000 116.6250 113.3125 $group [1] 1 1 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(114.116666666667, 114.916666666667, 115.304166666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/19w1v1289817033.ps tmp/19w1v1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/29w1v1289817033.ps tmp/29w1v1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/32oig1289817033.ps tmp/32oig1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/42oig1289817033.ps tmp/42oig1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/52oig1289817033.ps tmp/52oig1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/6vfzj1289817033.ps tmp/6vfzj1289817033.png",intern=TRUE)) character(0) > try(system("convert tmp/7vfzj1289817033.ps tmp/7vfzj1289817033.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.915 1.431 2.276