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(589,561,640,656,727,697,640,599,568,577,553,582,600,566,653,673,742,716,660,617,583,587,565,598,628,618,688,705,770,736,678,639,604,611,594,634,658,622,709,722,782,756,702,653,615,621,602,635,677,635,736,755,811,798,735,697,661,667,645,688,713,667,762,784,837,817,767,722,681,687,660,698,717,696,775,796,858,826,783,740,701,706,677,711,734,690,785,805,871,845,801,764,725,723,690,734,750,707,807,824,886,859,819,783,740,747,711,751,804,756,860,878,942,913,869,834,790,800,763,800,826,799,890,900,961,935,894,855,809,810,766,805,821,773,883,898,957,924,881,837,784,791,760,802,828,778,889,902,969,947,908,867,815,812,773,813,834,782,892,903,966,937,896,858,817,827,797,843) > 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] 168 > (np <- floor(n / par1)) [1] 14 > 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] 14 14 14 14 14 14 14 14 14 14 14 14 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 589 600 628 658 677 713 717 734 750 804 826 821 828 [2,] 561 566 618 622 635 667 696 690 707 756 799 773 778 [3,] 640 653 688 709 736 762 775 785 807 860 890 883 889 [4,] 656 673 705 722 755 784 796 805 824 878 900 898 902 [5,] 727 742 770 782 811 837 858 871 886 942 961 957 969 [6,] 697 716 736 756 798 817 826 845 859 913 935 924 947 [7,] 640 660 678 702 735 767 783 801 819 869 894 881 908 [8,] 599 617 639 653 697 722 740 764 783 834 855 837 867 [9,] 568 583 604 615 661 681 701 725 740 790 809 784 815 [10,] 577 587 611 621 667 687 706 723 747 800 810 791 812 [11,] 553 565 594 602 645 660 677 690 711 763 766 760 773 [12,] 582 598 634 635 688 698 711 734 751 800 805 802 813 [,14] [,15] [1,] 834 NA [2,] 782 NA [3,] 892 NA [4,] 903 NA [5,] 966 NA [6,] 937 NA [7,] 896 NA [8,] 858 NA [9,] 817 NA [10,] 827 NA [11,] 797 NA [12,] 843 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -28 -34 -10 -36 -42 -46 -21 -44 -43 -48 -27 -48 -50 [2,] 79 87 70 87 101 95 79 95 100 104 91 110 111 [3,] 16 20 17 13 19 22 21 20 17 18 10 15 13 [4,] 71 69 65 60 56 53 62 66 62 64 61 59 67 [5,] -30 -26 -34 -26 -13 -20 -32 -26 -27 -29 -26 -33 -22 [6,] -57 -56 -58 -54 -63 -50 -43 -44 -40 -44 -41 -43 -39 [7,] -41 -43 -39 -49 -38 -45 -43 -37 -36 -35 -39 -44 -41 [8,] -31 -34 -35 -38 -36 -41 -39 -39 -43 -44 -46 -53 -52 [9,] 9 4 7 6 6 6 5 -2 7 10 1 7 -3 [10,] -24 -22 -17 -19 -22 -27 -29 -33 -36 -37 -44 -31 -39 [11,] 29 33 40 33 43 38 34 44 40 37 39 42 40 [12,] 18 30 24 42 25 19 23 16 53 26 16 26 21 [,14] [,15] [1,] -52 NA [2,] 110 NA [3,] 11 NA [4,] 63 NA [5,] -29 NA [6,] -41 NA [7,] -38 NA [8,] -41 NA [9,] 10 NA [10,] -30 NA [11,] 46 NA [12,] 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/112sq1275129767.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/212sq1275129767.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/312sq1275129767.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/4ut9b1275129767.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,] 589.0 561 640 656.0 727.0 697.0 640 599 568 577.0 553.0 582.0 [2,] 658.0 622 709 722.0 782.0 756.0 702 653 615 621.0 602.0 635.0 [3,] 725.5 693 780 800.5 864.5 835.5 792 752 713 714.5 683.5 722.5 [4,] 821.0 773 883 898.0 957.0 924.0 881 837 790 800.0 763.0 802.0 [5,] 834.0 799 892 903.0 969.0 947.0 908 867 817 827.0 797.0 843.0 $n [1] 14 14 14 14 14 14 14 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 656.6695 629.2368 706.5245 726.18 790.6023 764.5582 716.4132 674.3018 [2,] 794.3305 756.7632 853.4755 874.82 938.3977 906.4418 867.5868 829.6982 [,9] [,10] [,11] [,12] [1,] 639.1023 638.9132 615.5141 651.9804 [2,] 786.8977 790.0868 751.4859 793.0196 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(589, 658, 725.5, 821, 834, 561, 622, 693, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5ut9b1275129767.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] [,12] [1,] -52.0 70 10 53.0 -34.0 -63 -49 -53 1 -44.0 29.0 16 [2,] -48.0 87 13 60.0 -30.0 -56 -43 -44 4 -36.0 34.0 19 [3,] -42.5 95 17 62.5 -26.5 -44 -40 -40 6 -29.5 39.5 24 [4,] -28.0 104 20 66.0 -26.0 -41 -38 -36 7 -22.0 42.0 26 [5,] -10.0 111 22 71.0 -20.0 -39 -35 -31 10 -17.0 46.0 30 $n [1] 14 14 14 14 14 14 14 14 14 14 14 13 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -50.94546 87.82136 14.04409 59.96636 -28.18909 -50.33409 -42.11136 [2,] -34.05454 102.17864 19.95591 65.03364 -24.81091 -37.66591 -37.88864 [,8] [,9] [,10] [,11] [,12] [1,] -43.37818 4.733182 -35.41182 36.12182 20.93251 [2,] -36.62182 7.266818 -23.58818 42.87818 27.06749 $out [1] -13 -2 -3 42 53 $group [1] 5 9 9 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-52, -48, -42.5, -28, -10, 70, 87, 95, 104, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6nkrw1275129767.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,] 553.0 565.0 594.0 602.0 635.0 660.0 677.0 690 707.0 756.0 766.0 760.0 [2,] 572.5 585.0 614.5 621.5 664.0 684.0 703.5 724 743.5 795.0 807.0 787.5 [3,] 594.0 608.5 636.5 655.5 692.5 717.5 728.5 749 767.0 819.0 840.5 829.0 [4,] 648.0 666.5 696.5 715.5 745.5 775.5 789.5 803 821.5 873.5 897.0 890.5 [5,] 727.0 742.0 770.0 782.0 811.0 837.0 858.0 871 886.0 942.0 961.0 957.0 [,13] [,14] [,15] [1,] 773.0 782.0 NA [2,] 812.5 822.0 NA [3,] 847.5 850.5 NA [4,] 905.0 899.5 NA [5,] 969.0 966.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 559.5639 571.3273 599.0992 612.626 655.3273 675.7662 689.2748 712.9676 [2,] 628.4361 645.6727 673.9008 698.374 729.6727 759.2338 767.7252 785.0324 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 731.4237 783.1956 799.4504 782.021 805.3101 815.1517 NA [2,] 802.5763 854.8044 881.5496 875.979 889.6899 885.8483 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" NA Warning message: In bxp(list(stats = c(553, 572.5, 594, 648, 727, 565, 585, 608.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7nkrw1275129767.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,] 682.5714 683.50 687.5000 [2,] 709.2500 713.75 711.3125 [3,] 737.2857 738.75 744.9375 [4,] 794.0714 796.25 804.0625 [5,] 862.7857 864.50 871.2500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 698.5981 701.1212 702.6336 [2,] 775.9733 776.3788 787.2414 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(682.571428571429, 709.25, 737.285714285714, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/112sq1275129767.ps tmp/112sq1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/212sq1275129767.ps tmp/212sq1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/312sq1275129767.ps tmp/312sq1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/4ut9b1275129767.ps tmp/4ut9b1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/5ut9b1275129767.ps tmp/5ut9b1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/6nkrw1275129767.ps tmp/6nkrw1275129767.png",intern=TRUE)) character(0) > try(system("convert tmp/7nkrw1275129767.ps tmp/7nkrw1275129767.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.846 1.462 2.080