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. 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(113.9000,112.0000,113.8500,113.0800,111.7200,110.6900,113.5300,113.9900,112.7400,112.1500,115.8200,118.3800,118.8100,123.8500,117.9600,120.1600,118.7400,119.8400,124.8100,121.3300,120.2000,118.3200,129.5800,130.2000,127.1900,133.1000,129.1200,123.2800,123.3600,124.1300,126.9700,127.1400,123.7000,123.6700,130.1900,134.0100,124.9600,129.9600,128.3200,132.3800,126.2500,128.9100,131.4200,129.4400,126.8600,126.7100,131.6300,132.7800,126.6100,132.8400,123.1400,128.1300,125.4900,126.4800,130.8600,127.3200,126.5600,126.6400,129.2600,126.4700,135.4000,135.5000,132.2200,122.6200,125.1600,128.5000,133.8600,128.8700,125.0700,125.2500,132.1600,130.2400) > 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,] 113.90 118.81 127.19 124.96 126.61 135.40 NA [2,] 112.00 123.85 133.10 129.96 132.84 135.50 NA [3,] 113.85 117.96 129.12 128.32 123.14 132.22 NA [4,] 113.08 120.16 123.28 132.38 128.13 122.62 NA [5,] 111.72 118.74 123.36 126.25 125.49 125.16 NA [6,] 110.69 119.84 124.13 128.91 126.48 128.50 NA [7,] 113.53 124.81 126.97 131.42 130.86 133.86 NA [8,] 113.99 121.33 127.14 129.44 127.32 128.87 NA [9,] 112.74 120.20 123.70 126.86 126.56 125.07 NA [10,] 112.15 118.32 123.67 126.71 126.64 125.25 NA [11,] 115.82 129.58 130.19 131.63 129.26 132.16 NA [12,] 118.38 130.20 134.01 132.78 126.47 130.24 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.90 5.04 5.91 5.00 6.23 0.10 NA [2,] 1.85 -5.89 -3.98 -1.64 -9.70 -3.28 NA [3,] -0.77 2.20 -5.84 4.06 4.99 -9.60 NA [4,] -1.36 -1.42 0.08 -6.13 -2.64 2.54 NA [5,] -1.03 1.10 0.77 2.66 0.99 3.34 NA [6,] 2.84 4.97 2.84 2.51 4.38 5.36 NA [7,] 0.46 -3.48 0.17 -1.98 -3.54 -4.99 NA [8,] -1.25 -1.13 -3.44 -2.58 -0.76 -3.80 NA [9,] -0.59 -1.88 -0.03 -0.15 0.08 0.18 NA [10,] 3.67 11.26 6.52 4.92 2.62 6.91 NA [11,] 2.56 0.62 3.82 1.15 -2.79 -1.92 NA [12,] 0.43 -3.01 -9.05 -6.17 8.93 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/1bofb1227962917.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/26jlc1227962917.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/3f3as1227962917.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/4ktx71227962917.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,] 113.900 112.00 113.85 113.08 111.72 110.690 124.810 113.99 112.740 112.15 [2,] 118.810 123.85 117.96 120.16 118.74 119.840 124.810 121.33 120.200 118.32 [3,] 125.785 131.40 125.73 122.95 124.26 125.305 128.915 127.23 124.385 124.46 [4,] 127.190 133.10 129.12 128.13 125.49 128.500 131.420 128.87 126.560 126.64 [5,] 135.400 135.50 132.22 132.38 126.25 128.910 133.860 129.44 126.860 126.71 [,11] [,12] [1,] 129.260 118.38 [2,] 129.260 126.47 [3,] 129.885 130.22 [4,] 131.630 132.78 [5,] 132.160 134.01 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 120.3796 125.4335 118.5314 117.8091 119.9060 119.7190 124.6513 122.3665 [2,] 131.1904 137.3665 132.9286 128.0909 128.6140 130.8910 133.1787 132.0935 [,9] [,10] [,11] [,12] [1,] 120.2826 119.0933 128.3563 126.1498 [2,] 128.4874 129.8267 131.4137 134.2902 $out [1] 113.53 115.82 $group [1] 7 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(113.9, 118.81, 125.785, 127.19, 135.4, 112, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/524y31227962917.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,] -1.90 -9.70 -9.600 -6.13 -1.030 2.51 -4.99 -3.800 -0.59 2.62 -2.790 -9.05 [2,] 0.10 -5.89 -5.840 -2.64 0.770 2.84 -3.54 -3.440 -0.59 3.67 -1.920 -6.17 [3,] 5.02 -3.63 0.715 -1.39 1.045 3.61 -2.73 -1.915 -0.09 5.72 0.885 -3.01 [4,] 5.91 -1.64 4.060 0.08 2.660 4.97 0.17 -1.130 0.08 6.91 2.560 0.43 [5,] 6.23 1.85 4.990 2.54 3.340 5.36 0.46 -0.760 0.18 11.26 3.820 8.93 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.272362 -6.3713873 -5.67082 -3.1444879 -0.1741110 2.236081 -5.1230698 [2,] 8.767638 -0.8886127 7.10082 0.3644879 2.2641110 4.983919 -0.3369302 [,8] [,9] [,10] [,11] [,12] [1,] -3.4050246 -0.5221716 3.630095 -2.004745 -7.673543 [2,] -0.4249754 0.3421716 7.809905 3.774745 1.653543 $out [1] -1.88 $group [1] 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1.90000000000001, 0.0999999999999943, 5.02, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6ya3x1227962917.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,] 110.690 117.960 123.280 124.960 125.490 122.620 NA [2,] 112.075 118.775 123.685 126.785 126.475 125.205 NA [3,] 113.305 120.180 127.055 129.175 126.625 129.555 NA [4,] 113.945 124.330 129.655 131.525 128.695 133.040 NA [5,] 115.820 130.200 134.010 132.780 130.860 135.500 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 112.4521 117.6463 124.3320 127.0131 125.6124 125.9814 NA [2,] 114.1579 122.7137 129.7780 131.3369 127.6376 133.1286 NA $out [1] 118.38 132.84 123.14 $group [1] 1 5 5 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(110.69, 112.075, 113.305, 113.945, 115.82, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7yglf1227962917.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.7867 122.9500 122.6513 [2,] 122.8067 124.4225 123.6637 [3,] 124.2900 125.7575 124.2706 [4,] 127.3917 129.4000 128.7606 [5,] 128.6800 131.4000 130.3050 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 122.1988 123.4872 121.9459 [2,] 126.3812 128.0278 126.5953 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(121.786666666667, 122.806666666667, 124.29, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1bofb1227962917.ps tmp/1bofb1227962917.png") > system("convert tmp/26jlc1227962917.ps tmp/26jlc1227962917.png") > system("convert tmp/3f3as1227962917.ps tmp/3f3as1227962917.png") > system("convert tmp/4ktx71227962917.ps tmp/4ktx71227962917.png") > system("convert tmp/524y31227962917.ps tmp/524y31227962917.png") > system("convert tmp/6ya3x1227962917.ps tmp/6ya3x1227962917.png") > system("convert tmp/7yglf1227962917.ps tmp/7yglf1227962917.png") > > > proc.time() user system elapsed 1.911 1.502 2.157