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(1664.81,2397.53,2840.71,3547.29,3752.96,3714.74,4349.61,3566.34,5021.82,6423.48,7600.60,19756.21,2499.81,5198.24,7225.14,4806.03,5900.88,4951.34,6179.12,4752.15,5496.43,5835.10,12600.08,28541.72,4717.02,5702.63,9957.58,5304.78,6492.43,6630.80,7349.62,8176.62,8573.17,9690.50,15151.84,34061.01,5921.10,5814.58,12421.25,6369.77,7609.12,7224.75,8121.22,7979.25,8093.06,8476.70,17914.66,30114.41,4826.64,6470.23,9638.77,8821.17,8722.37,10209.48,11276.55,12552.22,11637.39,13606.89,21822.11,45060.69,7615.03,9849.69,14558.40,11587.33,9332.56,13082.09,16732.78,19888.61,23933.38,25391.35,36024.80,80721.71,10243.24,11266.88,21826.84,17357.33,15997.79,18601.53,26155.15,28586.52,30505.41,30821.33,46634.38,104660.67) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1664.81 2499.81 4717.02 5921.10 4826.64 7615.03 10243.24 NA [2,] 2397.53 5198.24 5702.63 5814.58 6470.23 9849.69 11266.88 NA [3,] 2840.71 7225.14 9957.58 12421.25 9638.77 14558.40 21826.84 NA [4,] 3547.29 4806.03 5304.78 6369.77 8821.17 11587.33 17357.33 NA [5,] 3752.96 5900.88 6492.43 7609.12 8722.37 9332.56 15997.79 NA [6,] 3714.74 4951.34 6630.80 7224.75 10209.48 13082.09 18601.53 NA [7,] 4349.61 6179.12 7349.62 8121.22 11276.55 16732.78 26155.15 NA [8,] 3566.34 4752.15 8176.62 7979.25 12552.22 19888.61 28586.52 NA [9,] 5021.82 5496.43 8573.17 8093.06 11637.39 23933.38 30505.41 NA [10,] 6423.48 5835.10 9690.50 8476.70 13606.89 25391.35 30821.33 NA [11,] 7600.60 12600.08 15151.84 17914.66 21822.11 36024.80 46634.38 NA [12,] 19756.21 28541.72 34061.01 30114.41 45060.69 80721.71 104660.67 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 732.72 2698.43 985.61 -106.52 1643.59 2234.66 1023.64 NA [2,] 443.18 2026.90 4254.95 6606.67 3168.54 4708.71 10559.96 NA [3,] 706.58 -2419.11 -4652.80 -6051.48 -817.60 -2971.07 -4469.51 NA [4,] 205.67 1094.85 1187.65 1239.35 -98.80 -2254.77 -1359.54 NA [5,] -38.22 -949.54 138.37 -384.37 1487.11 3749.53 2603.74 NA [6,] 634.87 1227.78 718.82 896.47 1067.07 3650.69 7553.62 NA [7,] -783.27 -1426.97 827.00 -141.97 1275.67 3155.83 2431.37 NA [8,] 1455.48 744.28 396.55 113.81 -914.83 4044.77 1918.89 NA [9,] 1401.66 338.67 1117.33 383.64 1969.50 1457.97 315.92 NA [10,] 1177.12 6764.98 5461.34 9437.96 8215.22 10633.45 15813.05 NA [11,] 12155.61 15941.64 18909.17 12199.75 23238.58 44696.91 58026.29 NA [12,] -17256.40 -23824.70 -28139.91 -25287.77 -37445.66 -70478.47 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/13i3h1244298299.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/2bk441244298299.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/3w9p41244298299.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/4n5yp1244298299.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] [1,] 1664.810 2397.530 2840.710 3547.290 3752.960 3714.74 4349.61 [2,] 3608.415 5450.435 8431.955 5055.405 6196.655 5791.07 6764.37 [3,] 4826.640 5814.580 9957.580 6369.770 7609.120 7224.75 8121.22 [4,] 6768.065 8159.960 13489.825 10204.250 9027.465 11645.78 14004.66 [5,] 10243.240 11266.880 14558.400 17357.330 9332.560 18601.53 16732.78 [,8] [,9] [,10] [,11] [,12] [1,] 3566.34 5021.820 5835.10 7600.60 19756.21 [2,] 6365.70 6794.745 7450.09 13875.96 29328.07 [3,] 8176.62 8573.170 9690.50 17914.66 34061.01 [4,] 16220.42 17785.385 19499.12 28923.46 62891.20 [5,] 28586.52 30505.410 30821.33 46634.38 104660.67 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2939.748 4196.495 6937.102 3294.963 5918.606 3728.409 3797.433 [2,] 6713.532 7432.665 12978.058 9444.577 9299.634 10721.091 12445.007 [,8] [,9] [,10] [,11] [,12] [1,] 2291.543 2009.737 2495.014 8928.539 14017.65 [2,] 14061.697 15136.603 16885.986 26900.781 54104.37 $out [1] 21826.84 15997.79 26155.15 $group [1] 3 5 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1664.81, 3608.415, 4826.64, 6768.065, 10243.24, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5x6nb1244298299.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] [1,] -106.520 443.18 -6051.480 -2254.77 -949.540 634.870 -1426.97 -914.830 [2,] 859.165 2597.72 -4561.155 -729.17 -211.295 807.645 -462.62 255.180 [3,] 1023.640 4254.95 -2971.070 205.67 138.370 1067.070 827.00 744.280 [4,] 1939.125 5657.69 -1618.355 1141.25 2045.425 2439.235 1853.52 1687.185 [5,] 2698.430 6606.67 706.580 1239.35 3749.530 3650.690 3155.83 1918.890 [,9] [,10] [,11] [,12] [1,] 315.920 1177.12 12155.61 -37445.66 [2,] 361.155 6113.16 14070.69 -37445.66 [3,] 1117.330 8215.22 18909.17 -26713.84 [4,] 1429.815 10035.71 33967.75 -23824.70 [5,] 1969.500 15813.05 58026.29 -17256.40 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 378.7053 2427.585 -4728.463 -911.3146 -1209.307 92.71077 -556.1614 [2,] 1668.5747 6082.315 -1213.677 1322.6546 1486.047 2041.42923 2210.1614 [,8] [,9] [,10] [,11] [,12] [1,] -110.8903 479.1435 5872.739 7026.973 -35499.80 [2,] 1599.4503 1755.5165 10557.701 30791.367 -17927.88 $out [1] 10559.96 7553.62 4044.77 -70478.47 $group [1] 2 6 8 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-106.520000000000, 859.165, 1023.64, 1939.125, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/60zd51244298299.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] [1,] 1664.81 2499.810 4717.02 5814.580 4826.64 7615.03 10243.24 NA [2,] 3194.00 4878.685 6097.53 6797.260 8771.77 10718.51 16677.56 NA [3,] 3733.85 5665.765 7763.12 8036.155 10743.01 15645.59 23991.00 NA [4,] 5722.65 6702.130 9824.04 10448.975 13079.56 24662.36 30663.37 NA [5,] 7600.60 7225.140 15151.84 12421.250 13606.89 36024.80 46634.38 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2580.516 4834.079 6063.434 6370.583 8778.205 9285.704 17611.97 NA [2,] 4887.184 6497.451 9462.806 9701.727 12707.825 22005.476 30370.02 NA $out [1] 19756.21 12600.08 28541.72 34061.01 17914.66 30114.41 21822.11 [8] 45060.69 80721.71 104660.67 $group [1] 1 2 2 3 4 4 5 5 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(1664.81, 3194, 3733.85, 5722.65, 7600.6, 2499.81, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/74tta1244298299.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,] 5355.379 4826.64 5188.240 [2,] 8257.272 6797.26 7620.944 [3,] 11330.910 8148.92 10672.704 [4,] 13821.858 9824.04 12882.335 [5,] 14320.764 9957.58 13474.605 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8792.865 6768.385 8272.948 [2,] 13868.955 9529.455 13072.460 $out [1] 22535.50 48988.06 17914.66 34061.01 21399.71 46109.63 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(5355.37857142857, 8257.27214285714, 11330.91, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/13i3h1244298299.ps tmp/13i3h1244298299.png") > system("convert tmp/2bk441244298299.ps tmp/2bk441244298299.png") > system("convert tmp/3w9p41244298299.ps tmp/3w9p41244298299.png") > system("convert tmp/4n5yp1244298299.ps tmp/4n5yp1244298299.png") > system("convert tmp/5x6nb1244298299.ps tmp/5x6nb1244298299.png") > system("convert tmp/60zd51244298299.ps tmp/60zd51244298299.png") > system("convert tmp/74tta1244298299.ps tmp/74tta1244298299.png") > > > proc.time() user system elapsed 1.832 1.418 2.110