R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(611,639,630,586,695,552,619,681,421,307,754,690,644,643,608,651,691,627,634,731,475,337,803,722,590,724,627,696,825,677,656,785,412,352,839,729,696,641,695,638,762,635,721,854,418,367,824,687,601,676,740,691,683,594,729,731,386,331,706,715,657,653,642,643,718,654,632,731,392,344,792,852,649,629,685,617,715,715,629,916,531,357,917,828,708,858,775,785,1006,789,734,906,532,387,991,841,892,782,811,792,978,773,796,946,594,438,1023,868,791,760,779,852,1001,734,996,869,599,426,1138,1091) > 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] [,11] [1,] 611 644 590 696 601 657 649 708 892 791 NA [2,] 639 643 724 641 676 653 629 858 782 760 NA [3,] 630 608 627 695 740 642 685 775 811 779 NA [4,] 586 651 696 638 691 643 617 785 792 852 NA [5,] 695 691 825 762 683 718 715 1006 978 1001 NA [6,] 552 627 677 635 594 654 715 789 773 734 NA [7,] 619 634 656 721 729 632 629 734 796 996 NA [8,] 681 731 785 854 731 731 916 906 946 869 NA [9,] 421 475 412 418 386 392 531 532 594 599 NA [10,] 307 337 352 367 331 344 357 387 438 426 NA [11,] 754 803 839 824 706 792 917 991 1023 1138 NA [12,] 690 722 729 687 715 852 828 841 868 1091 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 28 -1 134 -55 75 -4 -20 150 -110 -31 NA [2,] -9 -35 -97 54 64 -11 56 -83 29 19 NA [3,] -44 43 69 -57 -49 1 -68 10 -19 73 NA [4,] 109 40 129 124 -8 75 98 221 186 149 NA [5,] -143 -64 -148 -127 -89 -64 0 -217 -205 -267 NA [6,] 67 7 -21 86 135 -22 -86 -55 23 262 NA [7,] 62 97 129 133 2 99 287 172 150 -127 NA [8,] -260 -256 -373 -436 -345 -339 -385 -374 -352 -270 NA [9,] -114 -138 -60 -51 -55 -48 -174 -145 -156 -173 NA [10,] 447 466 487 457 375 448 560 604 585 712 NA [11,] -64 -81 -110 -137 9 60 -89 -150 -155 -47 NA [12,] -46 -132 -33 -86 -58 -203 -120 51 -77 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/rcomp/tmp/16ryx1324577689.ps",horizontal=F,onefile=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/rcomp/tmp/2mc5m1324577689.ps",horizontal=F,onefile=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/rcomp/tmp/3rbuk1324577689.ps",horizontal=F,onefile=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/rcomp/tmp/47bsl1324577689.ps",horizontal=F,onefile=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] [,12] [1,] 590 629.0 608 586 683 552.0 619.0 681.0 386 307.0 706.0 687.0 [2,] 611 641.0 630 638 695 627.0 632.0 731.0 412 337.0 792.0 715.0 [3,] 653 664.5 690 671 740 665.5 688.5 819.5 448 354.5 831.5 778.5 [4,] 708 760.0 775 785 978 734.0 734.0 906.0 532 387.0 991.0 852.0 [5,] 791 858.0 811 852 1006 789.0 796.0 946.0 599 438.0 1138.0 868.0 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 604.5349 605.0429 617.5522 597.5529 598.6019 612.0385 637.5367 732.063 [2,] 701.4651 723.9571 762.4478 744.4471 881.3981 718.9615 739.4633 906.937 [,9] [,10] [,11] [,12] [1,] 388.0432 329.518 732.0717 710.0493 [2,] 507.9568 379.482 930.9283 846.9507 $out [1] 892 996 1091 $group [1] 1 7 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(590, 611, 653, 708, 791, 629, 641, 664.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5hhmk1324577689.ps",horizontal=F,onefile=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,] -110.0 -97 -68 -8.0 -267 -86 2 -436.0 -174 375.0 -155 -203 [2,] -31.0 -35 -49 75.0 -205 -22 62 -374.0 -156 448.0 -137 -120 [3,] -2.5 5 -9 116.5 -135 15 114 -348.5 -126 476.5 -85 -77 [4,] 75.0 54 43 149.0 -64 86 150 -270.0 -55 585.0 -47 -46 [5,] 150.0 64 73 221.0 0 135 172 -256.0 -48 712.0 60 51 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -55.46183 -39.46795 -54.96687 79.52665 -205.44922 -38.96111 70.03169 [2,] 50.46183 49.46795 36.96687 153.47335 -64.55078 68.96111 157.96831 [,8] [,9] [,10] [,11] [,12] [1,] -400.4625 -176.46363 408.0493 -129.96759 -115.97333 [2,] -296.5375 -75.53637 544.9507 -40.03241 -38.02667 $out [1] 262 287 -127 $group [1] 6 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-110, -31, -2.5, 75, 150, -97, -35, 5, 54, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6km201324577689.ps",horizontal=F,onefile=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] [,10] [,11] [1,] 421.0 608.0 412.0 635.0 594.0 632.0 531.0 532 594.0 426.0 NA [2,] 569.0 617.5 608.5 636.5 597.5 637.0 623.0 721 777.5 747.0 NA [3,] 624.5 643.5 686.5 691.0 687.0 653.5 667.0 787 803.5 821.5 NA [4,] 685.5 706.5 757.0 741.5 722.0 724.5 771.5 882 919.0 998.5 NA [5,] 754.0 803.0 839.0 854.0 740.0 852.0 917.0 1006 1023.0 1138.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 571.3636 602.9065 618.7682 643.1088 630.2147 613.5907 599.2682 713.5668 [2,] 677.6364 684.0935 754.2318 738.8912 743.7853 693.4093 734.7318 860.4332 [,9] [,10] [,11] [1,] 738.9609 706.7892 NA [2,] 868.0391 936.2108 NA $out [1] 307 475 337 352 418 367 386 331 392 344 357 387 438 $group [1] 1 2 2 3 4 4 5 5 6 6 7 8 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(421, 569, 624.5, 685.5, 754, 608, 617.5, 643.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/762731324577689.ps",horizontal=F,onefile=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,] 675.00 653.00 662.1250 [2,] 679.45 658.75 670.6250 [3,] 699.85 679.75 697.9375 [4,] 804.85 759.25 798.4375 [5,] 878.70 831.50 883.6250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 642.6542 633.9113 639.6414 [2,] 757.0458 725.5887 756.2336 $out [1] 476.000 364.600 448.000 354.500 472.625 360.375 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(675, 679.45, 699.85, 804.85, 878.7, 653, 658.75, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/16ryx1324577689.ps tmp/16ryx1324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/2mc5m1324577689.ps tmp/2mc5m1324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/3rbuk1324577689.ps tmp/3rbuk1324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/47bsl1324577689.ps tmp/47bsl1324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/5hhmk1324577689.ps tmp/5hhmk1324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/6km201324577689.ps tmp/6km201324577689.png",intern=TRUE)) character(0) > try(system("convert tmp/762731324577689.ps tmp/762731324577689.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.640 0.260 2.897