R version 2.7.2 (2008-08-25) 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(2319,2383,2513,2553,2341,2540,2371,2122,2301,2512,3145,2741,2548,1987,2281,2016,2434,2637,1831,1851,1839,2609,2417,2394,2372,2717,2998,2538,3007,2475,2175,2465,2279,2323,2746,2601,2486,2718,2646,2551,2712,2606,2365,3533,3509,2912,3599,2719,2869,4085,2686,2545,3071,3388,2652,3190,2884,3295,3818,3226,3953,3810,2877,3515,3708,3450,3360,4110,4384,3729,4263,3505,3674,3911,2951,3317,3417,3498,2768,2899,3171,3004,3481,3016,2595,3509,2833,3125,2556,3628,2876,2575,2903,3438,2926,3068,3015) > 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] 97 > (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] 9 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2319 2548 2372 2486 2869 3953 3674 2595 3015 [2,] 2383 1987 2717 2718 4085 3810 3911 3509 NA [3,] 2513 2281 2998 2646 2686 2877 2951 2833 NA [4,] 2553 2016 2538 2551 2545 3515 3317 3125 NA [5,] 2341 2434 3007 2712 3071 3708 3417 2556 NA [6,] 2540 2637 2475 2606 3388 3450 3498 3628 NA [7,] 2371 1831 2175 2365 2652 3360 2768 2876 NA [8,] 2122 1851 2465 3533 3190 4110 2899 2575 NA [9,] 2301 1839 2279 3509 2884 4384 3171 2903 NA [10,] 2512 2609 2323 2912 3295 3729 3004 3438 NA [11,] 3145 2417 2746 3599 3818 4263 3481 2926 NA [12,] 2741 2394 2601 2719 3226 3505 3016 3068 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 64 -561 345 232 1216 -143 237 914 NA [2,] 130 294 281 -72 -1399 -933 -960 -676 NA [3,] 40 -265 -460 -95 -141 638 366 292 NA [4,] -212 418 469 161 526 193 100 -569 NA [5,] 199 203 -532 -106 317 -258 81 1072 NA [6,] -169 -806 -300 -241 -736 -90 -730 -752 NA [7,] -249 20 290 1168 538 750 131 -301 NA [8,] 179 -12 -186 -24 -306 274 272 328 NA [9,] 211 770 44 -597 411 -655 -167 535 NA [10,] 633 -192 423 687 523 534 477 -512 NA [11,] -404 -23 -145 -880 -592 -758 -465 142 NA [12,] -193 -22 -115 150 727 169 -421 -53 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/10qph1225745944.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/2z2401225745944.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/3hm8g1225745944.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/4esgu1225745944.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] [,11] [1,] 2319 1987.0 2281.0 2016.0 2341.0 2475.0 1831.0 1851.0 1839.0 2323.0 2417.0 [2,] 2486 2550.0 2579.5 2541.5 2495.0 2573.0 2270.0 2293.5 2290.0 2560.5 2836.0 [3,] 2595 3113.5 2759.5 2552.0 2859.5 3012.5 2511.5 2737.0 2893.5 2958.0 3313.0 [4,] 3015 3860.5 2914.0 3221.0 3244.0 3474.0 2822.0 3361.5 3340.0 3366.5 3708.5 [5,] 3674 4085.0 2998.0 3515.0 3708.0 3628.0 3360.0 4110.0 4384.0 3729.0 4263.0 [,12] [1,] 2394.0 [2,] 2660.0 [3,] 2878.5 [4,] 3147.0 [5,] 3505.0 $n [1] 9 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2316.393 2381.436 2572.643 2172.422 2441.098 2509.188 2203.145 2140.4 [2,] 2873.607 3845.564 2946.357 2931.578 3277.902 3515.812 2819.855 3333.6 [,9] [,10] [,11] [,12] [1,] 2306.955 2507.757 2825.609 2606.455 [2,] 3480.045 3408.243 3800.391 3150.545 $out [1] 3953 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2319, 2486, 2595, 3015, 3674, 1987, 2550, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/50qqz1225745944.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,] -561.0 -1399.0 -460.0 -569.0 -532 -806 -301.0 -306.0 -655.0 -512.0 -880.0 [2,] -39.5 -946.5 -203.0 -56.0 -182 -744 -114.5 -105.0 -382.0 115.5 -675.0 [3,] 234.5 -374.0 -27.5 177.0 140 -515 210.5 83.5 127.5 500.0 -434.5 [4,] 629.5 205.5 329.0 443.5 260 -205 644.0 273.0 473.0 583.5 -84.0 [5,] 1216.0 294.0 638.0 526.0 317 -90 1168.0 328.0 770.0 687.0 142.0 [,12] [1,] -421.0 [2,] -154.0 [3,] -37.5 [4,] 159.5 [5,] 169.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -139.213 -1017.5237 -324.6828 -102.0279 -106.9075 -816.0931 -213.209 [2,] 608.213 269.5237 269.6828 456.0279 386.9075 -213.9069 634.209 [,8] [,9] [,10] [,11] [,12] [1,] -127.6562 -350.1153 238.5685 -764.6411 -212.6256 [2,] 294.6562 605.1153 761.4315 -104.3589 137.6256 $out [1] 1072 727 $group [1] 5 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-561, -39.5, 234.5, 629.5, 1216, -1399, -946.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/62e3u1225745944.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,] 2122.0 1831.0 2175.0 2365.0 2545.0 2877.0 2768.0 2556.0 3015 [2,] 2330.0 1919.0 2347.5 2578.5 2777.5 3477.5 2977.5 2714.0 3015 [3,] 2447.5 2337.5 2506.5 2715.0 3130.5 3718.5 3244.0 2914.5 3015 [4,] 2546.5 2491.0 2731.5 3210.5 3341.5 4031.5 3489.5 3281.5 3015 [5,] 2741.0 2637.0 3007.0 3599.0 4085.0 4384.0 3911.0 3628.0 3015 $n [1] 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2348.753 2076.607 2331.355 2426.741 2873.256 3465.817 3010.473 2655.659 [2,] 2546.247 2598.393 2681.645 3003.259 3387.744 3971.183 3477.527 3173.341 [,9] [1,] 3015 [2,] 3015 $out [1] 3145 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(2122, 2330, 2447.5, 2546.5, 2741, 1831, 1919, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7lxvi1225745944.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,] 2549.750 2511.50 2556.250 [2,] 2806.562 2666.00 2764.812 [3,] 2907.250 2869.00 2849.812 [4,] 3002.750 2985.25 2991.750 [5,] 3140.000 3313.00 3267.375 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2817.768 2723.388 2746.305 [2,] 2996.732 3014.612 2953.320 $out [1] 3299.375 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(2549.75, 2806.5625, 2907.25, 3002.75, 3140, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/10qph1225745944.ps tmp/10qph1225745944.png") > system("convert tmp/2z2401225745944.ps tmp/2z2401225745944.png") > system("convert tmp/3hm8g1225745944.ps tmp/3hm8g1225745944.png") > system("convert tmp/4esgu1225745944.ps tmp/4esgu1225745944.png") > system("convert tmp/50qqz1225745944.ps tmp/50qqz1225745944.png") > system("convert tmp/62e3u1225745944.ps tmp/62e3u1225745944.png") > system("convert tmp/7lxvi1225745944.ps tmp/7lxvi1225745944.png") > > > proc.time() user system elapsed 2.182 1.486 2.435