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(105.46,104.66,103.52,103.71,103.78,103.67,103.66,102.76,102,101.5,101.5,99.22,98.97,98.9,99.78,104.4,106.21,105.46,108.33,111.72,111.88,112.86,113.09,116.9,114.62,118.86,124.71,122.53,127.89,136.16,134.12,130.26,135.35,131.43,129.61,123.96,121.1,125.38,123.1,129.92,136.68,131.17,124.82,122.47,126.15,118.74,116.8,116.64,116.53,117.68,119.46,126.19,124.39,121.9,122.53,122.93,124.66,124.41,120.93,120.18,123.44,126.1,125.82,122.18,117.27,117.86,119.09,123.08,125.42,121.81,121.66,121.27,120.92,122.16,124.17,127.26,134.16,134.09,135.57,136.13,136.23,140.6,136.5,130.59,129.5,135.25,138.06,146.28,145.04,147.96,156.71,160.97,168.17,163.91,153.05,151.76) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 105.46 98.97 114.62 121.10 116.53 123.44 120.92 129.50 NA [2,] 104.66 98.90 118.86 125.38 117.68 126.10 122.16 135.25 NA [3,] 103.52 99.78 124.71 123.10 119.46 125.82 124.17 138.06 NA [4,] 103.71 104.40 122.53 129.92 126.19 122.18 127.26 146.28 NA [5,] 103.78 106.21 127.89 136.68 124.39 117.27 134.16 145.04 NA [6,] 103.67 105.46 136.16 131.17 121.90 117.86 134.09 147.96 NA [7,] 103.66 108.33 134.12 124.82 122.53 119.09 135.57 156.71 NA [8,] 102.76 111.72 130.26 122.47 122.93 123.08 136.13 160.97 NA [9,] 102.00 111.88 135.35 126.15 124.66 125.42 136.23 168.17 NA [10,] 101.50 112.86 131.43 118.74 124.41 121.81 140.60 163.91 NA [11,] 101.50 113.09 129.61 116.80 120.93 121.66 136.50 153.05 NA [12,] 99.22 116.90 123.96 116.64 120.18 121.27 130.59 151.76 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -0.80 -0.07 4.24 4.28 1.15 2.66 1.24 5.75 NA [2,] -1.14 0.88 5.85 -2.28 1.78 -0.28 2.01 2.81 NA [3,] 0.19 4.62 -2.18 6.82 6.73 -3.64 3.09 8.22 NA [4,] 0.07 1.81 5.36 6.76 -1.80 -4.91 6.90 -1.24 NA [5,] -0.11 -0.75 8.27 -5.51 -2.49 0.59 -0.07 2.92 NA [6,] -0.01 2.87 -2.04 -6.35 0.63 1.23 1.48 8.75 NA [7,] -0.90 3.39 -3.86 -2.35 0.40 3.99 0.56 4.26 NA [8,] -0.76 0.16 5.09 3.68 1.73 2.34 0.10 7.20 NA [9,] -0.50 0.98 -3.92 -7.41 -0.25 -3.61 4.37 -4.26 NA [10,] 0.00 0.23 -1.82 -1.94 -3.48 -0.15 -4.10 -10.86 NA [11,] -2.28 3.81 -5.65 -0.16 -0.75 -0.39 -5.91 -1.29 NA [12,] -0.25 -2.28 -2.86 -0.11 3.26 -0.35 -1.09 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/rcomp/tmp/1mzca1250436401.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/2ce8l1250436401.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/3mado1250436401.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/4eyu91250436401.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,] 98.970 98.90 99.780 103.71 103.78 103.670 103.660 102.760 102.000 [2,] 110.040 111.17 111.490 113.29 111.74 111.660 113.710 117.095 118.270 [3,] 118.725 120.51 123.635 124.36 126.14 126.535 123.675 123.005 125.785 [4,] 122.270 125.74 125.265 128.59 135.42 135.125 134.845 133.195 135.790 [5,] 129.500 135.25 138.060 146.28 145.04 147.960 156.710 136.130 136.230 [,10] [,11] [,12] [1,] 101.500 101.500 116.640 [2,] 115.800 114.945 116.770 [3,] 123.110 121.295 120.725 [4,] 136.015 133.055 127.275 [5,] 163.910 153.050 130.590 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 111.8931 112.371 115.9401 115.8132 112.9120 113.4271 111.8687 114.0113 [2,] 125.5569 128.649 131.3299 132.9068 139.3680 139.6429 135.4813 131.9987 [,9] [,10] [,11] [,12] [1,] 115.9981 111.8176 111.1785 114.8568 [2,] 135.5719 134.4024 131.4115 126.5932 $out [1] 160.97 168.17 99.22 151.76 $group [1] 8 9 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(98.97, 110.04, 118.725, 122.27, 129.5, 98.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5c4kp1250436401.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] [1,] -0.80 -2.28 -3.640 -4.91 -5.510 -2.040 -3.860 -0.760 -7.410 -4.100 -5.910 [2,] 0.54 -0.71 -0.995 -1.52 -1.620 -1.025 -1.625 0.130 -4.090 -3.790 -3.965 [3,] 1.95 1.33 3.855 0.94 -0.090 0.930 0.480 2.035 -2.055 -1.880 -1.020 [4,] 4.26 2.41 6.775 6.06 1.755 2.175 3.690 4.385 0.365 -0.075 -0.275 [5,] 5.75 5.85 8.220 6.90 2.920 2.870 4.260 7.200 4.370 0.230 3.810 [,12] [1,] -2.860 [2,] -1.685 [3,] -0.350 [4,] -0.180 [5,] -0.110 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.1280454 -0.4128768 -0.4854336 -3.294297 -1.975323 -0.857566 -2.489035 [2,] 4.0280454 3.0728768 8.1954336 5.174297 1.795323 2.717566 3.449035 [,8] [,9] [,10] [,11] [,12] [1,] -0.3419041 -4.543627 -3.9552523 -3.081287 -1.2487617 [2,] 4.4119041 0.433627 0.1952523 1.041287 0.5487617 $out [1] 8.27 -6.35 8.75 -10.86 3.26 $group [1] 5 6 6 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.799999999999997, 0.540000000000006, 1.95000000000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6887j1250436401.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] [1,] 99.220 98.90 114.620 116.640 116.530 117.270 120.920 129.50 NA [2,] 101.750 102.09 123.245 119.920 119.820 120.180 125.715 141.55 NA [3,] 103.590 107.27 128.750 123.960 122.215 121.995 134.125 149.86 NA [4,] 103.745 112.37 132.775 128.035 124.400 124.430 136.180 158.84 NA [5,] 105.460 116.90 136.160 136.680 126.190 126.100 140.600 168.17 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 102.6801 102.5812 124.4033 120.2587 120.1260 120.0565 129.3518 141.9739 [2,] 104.4999 111.9588 133.0967 127.6613 124.3040 123.9335 138.8982 157.7461 [,9] [1,] NA [2,] NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(99.22, 101.75, 103.59, 103.745, 105.46, 98.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7fziz1250436401.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,] 116.3175 118.7250 117.0075 [2,] 121.1962 121.0100 120.7287 [3,] 124.2850 123.3725 124.1250 [4,] 125.9469 125.0725 125.4688 [5,] 128.7325 126.5350 128.5175 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 122.1182 121.5196 121.9631 [2,] 126.4518 125.2254 126.2869 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(116.3175, 121.19625, 124.285, 125.946875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1mzca1250436401.ps tmp/1mzca1250436401.png") > system("convert tmp/2ce8l1250436401.ps tmp/2ce8l1250436401.png") > system("convert tmp/3mado1250436401.ps tmp/3mado1250436401.png") > system("convert tmp/4eyu91250436401.ps tmp/4eyu91250436401.png") > system("convert tmp/5c4kp1250436401.ps tmp/5c4kp1250436401.png") > system("convert tmp/6887j1250436401.ps tmp/6887j1250436401.png") > system("convert tmp/7fziz1250436401.ps tmp/7fziz1250436401.png") > > > proc.time() user system elapsed 1.187 0.924 1.769