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(81.28,69.39,67.63,51.25,103.97,133.83,162.37,172.91,163.01,151.50,111.73,88.58,74.29,63.98,61.18,76.48,107.98,124.97,145.57,140.20,143.84,138.80,104.06,74.70,60.18,55.16,35.62,56.18,85.44,114.08,133.64,67.14,95.58,89.37,75.24,69.18,54.49,57.50,62.16,76.67,110.04,127.38,156.47,167.56,153.54,124.08,100.97,79.17,68.13,61.77,54.31,60.30,84.18,104.05,114.66,105.55,96.61,70.94,63.91,58.61,44.53,49.58,57.39,76.76,104.57,125.41,143.11,136.35,135.15,131.70,96.87,70.63,66.29,63.49,62.97,66.43,101.49,127.69,133.21,158.72,148.61,134.31,100.99,75.16,59.74,52.87,52.07,57.38,79.43,101.40,120.19,134.38,135.97,113.83,84.38,70.28,65.96,56.36,49.57,68.33,90.32,117.06,134.69,131.67,129.25,118.77,88.44,76.79,75.28,73.89,76.24,88.58,105.83,115.84,127.76,131.75,119.63,93.38,75.55,51.79) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 81.28 74.29 60.18 54.49 68.13 44.53 66.29 59.74 65.96 75.28 [2,] 69.39 63.98 55.16 57.50 61.77 49.58 63.49 52.87 56.36 73.89 [3,] 67.63 61.18 35.62 62.16 54.31 57.39 62.97 52.07 49.57 76.24 [4,] 51.25 76.48 56.18 76.67 60.30 76.76 66.43 57.38 68.33 88.58 [5,] 103.97 107.98 85.44 110.04 84.18 104.57 101.49 79.43 90.32 105.83 [6,] 133.83 124.97 114.08 127.38 104.05 125.41 127.69 101.40 117.06 115.84 [7,] 162.37 145.57 133.64 156.47 114.66 143.11 133.21 120.19 134.69 127.76 [8,] 172.91 140.20 67.14 167.56 105.55 136.35 158.72 134.38 131.67 131.75 [9,] 163.01 143.84 95.58 153.54 96.61 135.15 148.61 135.97 129.25 119.63 [10,] 151.50 138.80 89.37 124.08 70.94 131.70 134.31 113.83 118.77 93.38 [11,] 111.73 104.06 75.24 100.97 63.91 96.87 100.99 84.38 88.44 75.55 [12,] 88.58 74.70 69.18 79.17 58.61 70.63 75.16 70.28 76.79 51.79 [,11] [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] [1,] -11.89 -10.31 -5.02 3.01 -6.36 5.05 -2.80 -6.87 -9.60 -1.39 [2,] -1.76 -2.80 -19.54 4.66 -7.46 7.81 -0.52 -0.80 -6.79 2.35 [3,] -16.38 15.30 20.56 14.51 5.99 19.37 3.46 5.31 18.76 12.34 [4,] 52.72 31.50 29.26 33.37 23.88 27.81 35.06 22.05 21.99 17.25 [5,] 29.86 16.99 28.64 17.34 19.87 20.84 26.20 21.97 26.74 10.01 [6,] 28.54 20.60 19.56 29.09 10.61 17.70 5.52 18.79 17.63 11.92 [7,] 10.54 -5.37 -66.50 11.09 -9.11 -6.76 25.51 14.19 -3.02 3.99 [8,] -9.90 3.64 28.44 -14.02 -8.94 -1.20 -10.11 1.59 -2.42 -12.12 [9,] -11.51 -5.04 -6.21 -29.46 -25.67 -3.45 -14.30 -22.14 -10.48 -26.25 [10,] -39.77 -34.74 -14.13 -23.11 -7.03 -34.83 -33.32 -29.45 -30.33 -17.83 [11,] -23.15 -29.36 -6.06 -21.80 -5.30 -26.24 -25.83 -14.10 -11.65 -23.76 [12,] -14.29 -14.52 -14.69 -11.04 -14.08 -4.34 -15.42 -4.32 -1.51 NA [,11] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] 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/1c0mj1272301850.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/2c0mj1272301850.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/3c0mj1272301850.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/4nr341272301850.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,] 44.530 49.580 49.570 51.25 79.43 101.400 114.660 105.550 95.58 70.940 [2,] 59.740 55.160 52.070 57.38 85.44 114.080 127.760 131.670 119.63 93.380 [3,] 66.125 59.635 59.285 67.38 102.73 121.015 134.165 135.365 135.56 121.425 [4,] 74.290 63.980 62.970 76.67 105.83 127.380 145.570 158.720 148.61 134.310 [5,] 81.280 73.890 76.240 88.58 110.04 133.830 162.370 172.910 163.01 151.500 [,11] [,12] [1,] 63.910 58.610 [2,] 75.550 69.180 [3,] 92.655 72.665 [4,] 100.990 76.790 [5,] 111.730 79.170 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 58.85524 55.22818 53.83893 57.74195 92.54234 114.3698 125.2664 121.8497 [2,] 73.39476 64.04182 64.73107 77.01805 112.91766 127.6602 143.0636 148.8803 [,9] [,10] [,11] [,12] [1,] 121.0804 100.9747 79.94416 68.86274 [2,] 150.0396 141.8753 105.36584 76.46726 $out [1] 35.62 67.14 88.58 51.79 $group [1] 3 8 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(44.53, 59.74, 66.125, 74.29, 81.28, 49.58, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5nr341272301850.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] [1,] -11.89 -19.54 3.460 17.250 10.010 5.520 -9.110 -14.02 -29.460 -39.77 [2,] -9.60 -6.79 5.310 22.050 17.340 11.920 -6.760 -10.11 -25.670 -34.74 [3,] -5.69 -1.28 13.425 28.535 21.405 18.245 0.485 -5.68 -12.905 -29.89 [4,] -1.39 2.35 18.760 33.370 26.740 20.600 11.090 1.59 -6.210 -17.83 [5,] 5.05 7.81 20.560 35.060 29.860 29.090 25.510 3.64 -3.450 -7.03 [,11] [,12] [1,] -29.360 -15.42 [2,] -25.830 -14.52 [3,] -22.475 -14.08 [4,] -11.650 -4.34 [5,] -5.300 -1.51 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -9.792043 -5.846708 6.704844 22.87908 16.70839 13.90813 -8.433572 [2,] -1.587957 3.286708 20.145156 34.19092 26.10161 22.58187 9.403572 [,8] [,9] [,10] [,11] [,12] [1,] -11.5257865 -22.627992 -38.33891 -29.55989 -19.441467 [2,] 0.1657865 -3.182008 -21.44109 -15.39011 -8.718533 $out [1] -16.38 52.72 -66.50 28.44 $group [1] 3 4 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-11.89, -9.6, -5.69, -1.39, 5.05, -19.54, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6nr341272301850.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,] 51.250 61.180 35.620 54.490 54.310 44.530 62.97 52.070 49.570 [2,] 75.335 74.495 58.180 69.415 61.035 64.010 66.36 58.560 67.145 [3,] 107.850 106.020 72.210 105.505 69.535 100.720 101.24 81.905 89.380 [4,] 156.935 139.500 92.475 140.460 100.330 133.425 133.76 117.010 124.010 [5,] 172.910 145.570 133.640 167.560 114.660 143.110 158.72 135.970 134.690 [,10] [,11] [1,] 51.790 NA [2,] 75.415 NA [3,] 90.980 NA [4,] 117.735 NA [5,] 131.750 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 70.63169 76.37078 56.56782 73.1009 51.61229 69.05935 70.49841 [2,] 145.06831 135.66922 87.85218 137.9091 87.45771 132.38065 131.98159 [,8] [,9] [,10] [,11] [1,] 55.24556 63.44349 71.67756 NA [2,] 108.56444 115.31651 110.28244 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(51.25, 75.335, 107.85, 156.935, 172.91, 61.18, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7g0271272301850.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,] 57.9140 59.2850 57.69875 [2,] 66.4265 66.7525 66.83312 [3,] 93.7695 97.6925 92.72938 [4,] 125.6450 127.7950 127.71500 [5,] 137.1670 135.5600 142.89000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 66.75954 69.8506 64.96074 [2,] 120.77946 125.5344 120.49801 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(57.914, 66.4265, 93.7695, 125.645, 137.167, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1c0mj1272301850.ps tmp/1c0mj1272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/2c0mj1272301850.ps tmp/2c0mj1272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/3c0mj1272301850.ps tmp/3c0mj1272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/4nr341272301850.ps tmp/4nr341272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/5nr341272301850.ps tmp/5nr341272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/6nr341272301850.ps tmp/6nr341272301850.png",intern=TRUE)) character(0) > try(system("convert tmp/7g0271272301850.ps tmp/7g0271272301850.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.890 1.480 2.215