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(727,817,918,786,803,756,725,523,538,587,505,521,498,550,637,622,668,669,670,499,539,593,429,622,533,655,835,686,706,869,777,739,637,597,629,940,444,496,801,659,767,876,601,697,745,655,572,628,650,677,900,780,896,1092,823,735,770,915,645,566,707,785,762,712,714,823,609,620,619,638,483,535,617,698,804,824,878,1019,974,773,734,827,804,721,659,732,839,994,828,1039,1072,803,1035,922,834,1739,359,513,699,741,793,877,750,752,675,682,583,632,606,645,980,847,941,1066,936,880,808,741,780,675,782,795,873,727,998,768,714,782,578,664,560,516,752,597,716,691,752,718,737,621,472,719,497,536,653,605,637,743,719,653,675,590,527,534,463,542,568,501,678,774,665,742,715,638,656,606,498,587,677,547,871,731,752,862,619,700,667,667,650,547,637,655,703,886,896,831,741,833,750,779,655,739,845,795,1021,726,1045,915,852,772,729,755,691,729,702,702,894,765,753,876,781,776,606,775,663,649,821,771,635,1070,693,779) > 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] 222 > (np <- floor(n / par1)) [1] 18 > 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] 19 19 19 19 19 19 18 18 18 18 18 18 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 727 498 533 444 650 707 617 659 359 606 782 752 653 [2,] 817 550 655 496 677 785 698 732 513 645 795 597 605 [3,] 918 637 835 801 900 762 804 839 699 980 873 716 637 [4,] 786 622 686 659 780 712 824 994 741 847 727 691 743 [5,] 803 668 706 767 896 714 878 828 793 941 998 752 719 [6,] 756 669 869 876 1092 823 1019 1039 877 1066 768 718 653 [7,] 725 670 777 601 823 609 974 1072 750 936 714 737 675 [8,] 523 499 739 697 735 620 773 803 752 880 782 621 590 [9,] 538 539 637 745 770 619 734 1035 675 808 578 472 527 [10,] 587 593 597 655 915 638 827 922 682 741 664 719 534 [11,] 505 429 629 572 645 483 804 834 583 780 560 497 463 [12,] 521 622 940 628 566 535 721 1739 632 675 516 536 542 [,14] [,15] [,16] [,17] [,18] [,19] [1,] 568 677 637 845 702 821 [2,] 501 547 655 795 702 771 [3,] 678 871 703 1021 894 635 [4,] 774 731 886 726 765 1070 [5,] 665 752 896 1045 753 693 [6,] 742 862 831 915 876 779 [7,] 715 619 741 852 781 NA [8,] 638 700 833 772 776 NA [9,] 656 667 750 729 606 NA [10,] 606 667 779 755 775 NA [11,] 498 650 655 691 663 NA [12,] 587 547 739 729 649 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 90 52 122 52 27 78 81 73 154 39 13 -155 -48 [2,] 101 87 180 305 223 -23 106 107 186 335 78 119 32 [3,] -132 -15 -149 -142 -120 -50 20 155 42 -133 -146 -25 106 [4,] 17 46 20 108 116 2 54 -166 52 94 271 61 -24 [5,] -47 1 163 109 196 109 141 211 84 125 -230 -34 -66 [6,] -31 1 -92 -275 -269 -214 -45 33 -127 -130 -54 19 22 [7,] -202 -171 -38 96 -88 11 -201 -269 2 -56 68 -116 -85 [8,] 15 40 -102 48 35 -1 -39 232 -77 -72 -204 -149 -63 [9,] 49 54 -40 -90 145 19 93 -113 7 -67 86 247 7 [10,] -82 -164 32 -83 -270 -155 -23 -88 -99 39 -104 -222 -71 [11,] 16 193 311 56 -79 52 -83 905 49 -105 -44 39 79 [12,] -23 -89 -496 22 141 82 -62 -1380 -26 107 236 117 26 [,14] [,15] [,16] [,17] [,18] [,19] [1,] -67 -130 18 -50 0 -50 [2,] 177 324 48 226 192 -136 [3,] 96 -140 183 -295 -129 435 [4,] -109 21 10 319 -12 -377 [5,] 77 110 -65 -130 123 86 [6,] -27 -243 -90 -63 -95 NA [7,] -77 81 92 -80 -5 NA [8,] 18 -33 -83 -43 -170 NA [9,] -50 0 29 26 169 NA [10,] -108 -17 -124 -64 -112 NA [11,] 89 -103 84 38 -14 NA [12,] 90 90 106 -27 172 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/rcomp/tmp/1dy7i1229340498.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/rcomp/tmp/24t5k1229340498.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/rcomp/tmp/32fw21229340498.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/rcomp/tmp/452xu1229340498.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] [,12] [1,] 444 496.0 635.0 622 665.0 653 601 499 472.0 534.0 429 516 [2,] 587 573.5 701.0 719 716.5 762 675 621 578.0 606.0 498 542 [3,] 653 655.0 804.0 743 767.0 862 739 737 661.5 674.5 606 625 [4,] 717 751.5 883.5 805 887.0 896 823 776 745.0 775.0 663 721 [5,] 845 817.0 1021.0 886 1045.0 1092 974 880 808.0 922.0 834 940 $n [1] 19 19 19 19 19 19 18 18 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 605.878 590.4791 737.848 711.827 705.1977 813.4281 683.8834 679.2765 [2,] 700.122 719.5209 870.152 774.173 828.8023 910.5719 794.1166 794.7235 [,9] [,10] [,11] [,12] [1,] 599.3076 611.5628 544.5524 558.3387 [2,] 723.6924 737.4372 667.4476 691.6613 $out [1] 359 994 1070 1072 1035 1739 $group [1] 1 4 4 7 9 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(444, 587, 653, 717, 845, 496, 573.5, 655, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5jn1u1229340498.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,] -155.0 -23.0 -295.0 -109.0 -230.0 -275.0 -269.0 -204 -113.0 -164.0 -105 [2,] -49.0 82.5 -136.5 -5.0 -40.5 -130.0 -116.0 -83 -40.0 -124.0 -44 [3,] 27.0 119.0 -50.0 21.0 86.0 -76.5 -66.5 -41 22.5 -93.5 44 [4,] 75.5 207.5 69.0 77.5 124.0 -27.0 11.0 18 86.0 -64.0 84 [5,] 154.0 335.0 183.0 116.0 211.0 33.0 96.0 48 247.0 -17.0 193 [,12] [1,] -89 [2,] -27 [3,] 54 [4,] 107 [5,] 236 $n [1] 19 19 19 19 19 18 18 18 18 18 18 18 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -18.12837 73.69039 -124.48899 -8.904341 26.37256 -114.85819 -113.79602 [2,] 72.12837 164.30961 24.48899 50.904341 145.62744 -38.14181 -19.20398 [,8] [,9] [,10] [,11] [,12] [1,] -78.613367 -24.42361 -115.84457 -3.668425 4.097117 [2,] -3.386633 69.42361 -71.15543 91.668425 103.902883 $out [1] -136 435 -166 271 319 -377 232 32 -270 39 -222 311 [13] 905 -496 -1380 $group [1] 2 3 4 4 4 4 8 10 10 10 10 11 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-155, -49, 27, 75.5, 154, -23, 82.5, 119, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/67xi71229340498.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] [,10] [,11] [,12] [1,] 505.0 429.0 533 444.0 566.0 483.0 617.0 659.0 513.0 606.0 516.0 472.0 [2,] 530.5 519.0 633 586.5 663.5 614.0 727.5 815.5 607.5 708.0 621.0 566.5 [3,] 726.0 607.5 696 657.0 775.0 672.5 804.0 880.5 690.5 827.5 747.5 703.5 [4,] 794.5 652.5 806 756.0 898.0 738.0 852.5 1037.0 751.0 938.5 788.5 728.0 [5,] 918.0 670.0 940 876.0 1092.0 823.0 1019.0 1072.0 877.0 1066.0 998.0 752.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 463 498.0 547.0 637.0 691.0 606.0 635 [2,] 538 577.5 634.5 679.0 729.0 682.5 693 [3,] 621 647.0 672.0 745.5 783.5 759.0 775 [4,] 664 696.5 741.5 832.0 883.5 778.5 821 [5,] 743 774.0 871.0 896.0 1045.0 894.0 821 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 605.5878 546.6098 617.0935 579.6899 668.043 615.9428 746.9867 779.4724 [2,] 846.4122 668.3902 774.9065 734.3101 881.957 729.0572 861.0133 981.5276 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 625.0487 722.3674 671.1021 629.8388 563.5306 592.7233 623.1966 675.7157 [2,] 755.9513 932.6326 823.8979 777.1612 678.4694 701.2767 720.8034 815.2843 [,17] [,18] [,19] [1,] 713.0315 715.2138 692.4359 [2,] 853.9685 802.7862 857.5641 $out [1] 1739 359 1070 $group [1] 8 9 19 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" NA Warning message: In bxp(list(stats = c(505, 530.5, 726, 794.5, 918, 429, 519, 607.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/74bl31229340498.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,] 607.8333 606.00 580.375 [2,] 665.5892 654.00 657.250 [3,] 705.2500 705.75 696.125 [4,] 788.6053 755.00 777.125 [5,] 854.2105 862.00 829.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 649.1415 659.6832 641.4492 [2,] 761.3585 751.8168 750.8008 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(607.833333333333, 665.58918128655, 705.25, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1dy7i1229340498.ps tmp/1dy7i1229340498.png") > system("convert tmp/24t5k1229340498.ps tmp/24t5k1229340498.png") > system("convert tmp/32fw21229340498.ps tmp/32fw21229340498.png") > system("convert tmp/452xu1229340498.ps tmp/452xu1229340498.png") > system("convert tmp/5jn1u1229340498.ps tmp/5jn1u1229340498.png") > system("convert tmp/67xi71229340498.ps tmp/67xi71229340498.png") > system("convert tmp/74bl31229340498.ps tmp/74bl31229340498.png") > > > proc.time() user system elapsed 1.258 0.982 1.912