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(104.65,109.2,103.93,101.93,101.58,102.63,106.83,105.57,105,104.16,106.93,106.5,106.47,112.33,106.81,103.49,104.13,103.2,105.55,102.8,106.68,105.43,109.01,112.24,112.87,110.98,112.85,112.08,110.72,109.69,112.53,112.99,111.74,111.15,114.82,117.38,117.81,122.85,116.96,119.16,117.74,118.84,123.81,120.33,119.2,117.32,128.58,129.2,126.19,132.1,128.12,122.28,122.36,123.13,125.97,126.14,122.7,122.67,129.19,133.01,123.96,128.96,127.32,131.38,125.25,127.91,130.42,128.44,125.86,125.71,130.63,131.78,125.61,131.84,122.14,127.13,124.49,125.48,129.86,126.32,125.56,125.64,128.26,125.47,134.4,134.5,131.22,121.62,124.16,127.5,132.86,127.87,124.07,124.25,131.16,129.24,129.24,135.51,128.97,126.89,127.52,130.31,132.39,132.69,128.73,129.66,127.72,132.63,129.74,138.46,134.31,128.8,129.95,134.15,136.01,135.1,132.27,132.49,130.06,136.11,131.47,140.61,141.65,126.75,133.9,138.75,141.86,141.13,138.76,138.65,138.59,147.09,140.17,152.38,144.02,139.55,141.1,145.85,147.88,145.17) > 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] 140 > (np <- floor(n / par1)) [1] 11 > 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 11 11 11 11 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 104.65 106.47 112.87 117.81 126.19 123.96 125.61 134.40 129.24 129.74 [2,] 109.20 112.33 110.98 122.85 132.10 128.96 131.84 134.50 135.51 138.46 [3,] 103.93 106.81 112.85 116.96 128.12 127.32 122.14 131.22 128.97 134.31 [4,] 101.93 103.49 112.08 119.16 122.28 131.38 127.13 121.62 126.89 128.80 [5,] 101.58 104.13 110.72 117.74 122.36 125.25 124.49 124.16 127.52 129.95 [6,] 102.63 103.20 109.69 118.84 123.13 127.91 125.48 127.50 130.31 134.15 [7,] 106.83 105.55 112.53 123.81 125.97 130.42 129.86 132.86 132.39 136.01 [8,] 105.57 102.80 112.99 120.33 126.14 128.44 126.32 127.87 132.69 135.10 [9,] 105.00 106.68 111.74 119.20 122.70 125.86 125.56 124.07 128.73 132.27 [10,] 104.16 105.43 111.15 117.32 122.67 125.71 125.64 124.25 129.66 132.49 [11,] 106.93 109.01 114.82 128.58 129.19 130.63 128.26 131.16 127.72 130.06 [12,] 106.50 112.24 117.38 129.20 133.01 131.78 125.47 129.24 132.63 136.11 [,11] [,12] [1,] 131.47 140.17 [2,] 140.61 152.38 [3,] 141.65 144.02 [4,] 126.75 139.55 [5,] 133.90 141.10 [6,] 138.75 145.85 [7,] 141.86 147.88 [8,] 141.13 145.17 [9,] 138.76 NA [10,] 138.65 NA [11,] 138.59 NA [12,] 147.09 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 4.55 5.86 -1.89 5.04 5.91 5.00 6.23 0.10 6.27 8.72 9.14 12.21 [2,] -5.27 -5.52 1.87 -5.89 -3.98 -1.64 -9.70 -3.28 -6.54 -4.15 1.04 -8.36 [3,] -2.00 -3.32 -0.77 2.20 -5.84 4.06 4.99 -9.60 -2.08 -5.51 -14.90 -4.47 [4,] -0.35 0.64 -1.36 -1.42 0.08 -6.13 -2.64 2.54 0.63 1.15 7.15 1.55 [5,] 1.05 -0.93 -1.03 1.10 0.77 2.66 0.99 3.34 2.79 4.20 4.85 4.75 [6,] 4.20 2.35 2.84 4.97 2.84 2.51 4.38 5.36 2.08 1.86 3.11 2.03 [7,] -1.26 -2.75 0.46 -3.48 0.17 -1.98 -3.54 -4.99 0.30 -0.91 -0.73 -2.71 [8,] -0.57 3.88 -1.25 -1.13 -3.44 -2.58 -0.76 -3.80 -3.96 -2.83 -2.37 NA [9,] -0.84 -1.25 -0.59 -1.88 -0.03 -0.15 0.08 0.18 0.93 0.22 -0.11 NA [10,] 2.77 3.58 3.67 11.26 6.52 4.92 2.62 6.91 -1.94 -2.43 -0.06 NA [11,] -0.43 3.23 2.56 0.62 3.82 1.15 -2.79 -1.92 4.91 6.05 8.50 NA [12,] -0.03 0.63 0.43 -3.01 -9.05 -6.17 8.93 0.00 -2.89 -4.64 -6.92 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/1lfp61259492551.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/2b2gd1259492551.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/3joir1259492551.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/4jmpy1259492551.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,] 104.650 109.200 103.930 101.930 101.580 102.630 105.550 102.800 105.000 [2,] 115.340 117.590 114.905 115.620 114.230 114.265 118.170 116.660 115.470 [3,] 125.900 131.970 127.720 124.515 124.325 126.490 130.140 127.095 124.070 [4,] 130.605 136.985 132.765 127.965 128.735 132.230 134.435 133.895 127.295 [5,] 140.170 152.380 144.020 139.550 141.100 145.850 147.880 145.170 138.760 [,10] [,11] [,12] [1,] 104.160 109.010 106.500 [2,] 114.235 121.270 121.425 [3,] 124.250 128.580 129.240 [4,] 127.685 130.345 132.820 [5,] 138.650 138.590 147.090 $n [1] 12 12 12 12 12 12 12 12 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 118.9375 123.1238 119.5739 118.8844 117.7092 118.2960 122.7214 119.234 [2,] 132.8625 140.8162 135.8661 130.1456 130.9408 134.6840 137.5586 134.956 [,9] [,10] [,11] [,12] [1,] 118.4367 117.8426 124.2568 123.8116 [2,] 129.7033 130.6574 132.9032 134.6684 $out [1] 106.93 $group [1] 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(104.65, 115.34, 125.9, 130.605, 140.17, 109.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5rtdx1259492551.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,] 4.550 -9.700 -14.900 -2.640 -1.03 1.860 -4.990 -3.960 -1.880 -2.43 -2.790 [2,] 4.775 -6.215 -5.675 -1.390 0.88 2.215 -3.115 -3.135 -0.715 1.28 0.095 [3,] 5.885 -4.710 -2.700 0.355 1.88 2.840 -1.620 -2.370 -0.110 3.58 2.560 [4,] 7.495 -2.460 0.715 1.350 3.77 4.290 -0.280 -0.945 0.130 5.72 4.365 [5,] 9.140 1.870 4.990 2.540 4.85 5.360 0.460 -0.570 0.930 11.26 8.500 [,12] [1,] -9.050 [2,] -5.405 [3,] -2.890 [4,] 0.215 [5,] 0.630 $n [1] 12 12 12 12 12 12 12 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 4.64439 -6.422681 -5.6145219 -0.8947324 0.5618516 1.893579 -2.9130625 [2,] 7.12561 -2.997319 0.2145219 1.6047324 3.1981484 3.786421 -0.3269375 [,8] [,9] [,10] [,11] [,12] [1,] -3.413290 -0.5125478 1.464838 0.5258236 -5.5673001 [2,] -1.326710 0.2925478 5.695162 4.5941764 -0.2126999 $out [1] -1.89 0.10 12.21 -6.13 7.15 3.88 8.93 $group [1] 1 1 1 4 4 8 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(4.55, 4.77500000000001, 5.885, 7.49499999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ukul1259492551.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,] 101.580 102.80 109.690 116.960 122.280 123.960 124.490 121.620 126.890 [2,] 103.280 103.81 111.065 117.775 122.685 125.785 125.475 124.205 128.225 [3,] 104.825 106.01 112.305 119.180 126.055 128.175 125.625 128.555 129.450 [4,] 106.665 107.91 112.930 123.330 128.655 130.525 127.695 132.040 132.510 [5,] 109.200 112.33 114.820 129.200 133.010 131.780 129.860 134.500 135.510 [,10] [,11] [,12] [1,] 128.800 131.470 139.550 [2,] 130.005 136.245 140.635 [3,] 133.320 138.755 144.595 [4,] 135.555 141.390 146.865 [5,] 138.460 147.090 152.380 $n [1] 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 103.2811 104.1400 111.4544 116.6463 123.3320 126.0131 124.6124 124.9814 [2,] 106.3689 107.8800 113.1556 121.7137 128.7780 130.3369 126.6376 132.1286 [,9] [,10] [,11] [,12] [1,] 127.4956 130.7886 136.4083 141.1148 [2,] 131.4044 135.8514 141.1017 148.0752 $out [1] 117.38 131.84 122.14 126.75 $group [1] 3 7 7 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(101.58, 103.28, 104.825, 106.665, 109.2, 102.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7pbde1259492551.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,] 121.5573 124.0700 120.9600 [2,] 121.8892 124.4200 122.2625 [3,] 124.4058 126.7925 123.9369 [4,] 126.2717 128.9100 126.5081 [5,] 129.1433 131.9700 128.2337 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 122.4069 124.7446 122.0004 [2,] 126.4047 128.8404 125.8733 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(121.557272727273, 121.889166666667, 124.405833333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1lfp61259492551.ps tmp/1lfp61259492551.png") > system("convert tmp/2b2gd1259492551.ps tmp/2b2gd1259492551.png") > system("convert tmp/3joir1259492551.ps tmp/3joir1259492551.png") > system("convert tmp/4jmpy1259492551.ps tmp/4jmpy1259492551.png") > system("convert tmp/5rtdx1259492551.ps tmp/5rtdx1259492551.png") > system("convert tmp/6ukul1259492551.ps tmp/6ukul1259492551.png") > system("convert tmp/7pbde1259492551.ps tmp/7pbde1259492551.png") > > > proc.time() user system elapsed 1.186 0.918 3.126