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(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,707,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) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 627 695 740 642 685 775 NA [2,] 696 638 691 643 617 785 NA [3,] 825 762 683 718 715 1006 NA [4,] 677 635 594 654 715 789 NA [5,] 656 721 729 632 629 734 NA [6,] 785 854 731 731 916 906 NA [7,] 412 418 386 392 531 532 NA [8,] 352 367 331 344 357 387 NA [9,] 839 824 707 792 917 991 NA [10,] 729 687 715 852 828 841 NA [11,] 696 601 657 649 708 892 NA [12,] 641 676 653 629 858 782 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 69 -57 -49 1 -68 10 NA [2,] 129 124 -8 75 98 221 NA [3,] -148 -127 -89 -64 0 -217 NA [4,] -21 86 135 -22 -86 -55 NA [5,] 129 133 2 99 287 172 NA [6,] -373 -436 -345 -339 -385 -374 NA [7,] -60 -51 -55 -48 -174 -145 NA [8,] 487 457 376 448 560 604 NA [9,] -110 -137 8 60 -89 -150 NA [10,] -33 -86 -58 -203 -120 51 NA [11,] -55 75 -4 -20 150 -110 NA [12,] 54 64 -11 56 -83 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/1zb931291982708.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/2zb931291982708.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/3zb931291982708.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/4skq61291982708.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,] 627 617 683 594.0 629.0 731.0 386 331.0 707.0 687.0 601.0 629.0 [2,] 642 638 715 635.0 632.0 731.0 392 344.0 792.0 715.0 649.0 641.0 [3,] 690 667 740 665.5 688.5 819.5 415 354.5 831.5 778.5 676.5 664.5 [4,] 740 696 825 715.0 729.0 906.0 531 367.0 917.0 841.0 708.0 782.0 [5,] 775 696 825 789.0 734.0 916.0 532 387.0 991.0 852.0 708.0 858.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 626.7868 629.5881 669.0464 613.8974 625.9319 706.6193 325.3405 339.6643 [2,] 753.2132 704.4119 810.9536 717.1026 751.0681 932.3807 504.6595 369.3357 [,9] [,10] [,11] [,12] [1,] 750.871 697.226 638.4431 573.5504 [2,] 912.129 859.774 714.5569 755.4496 $out [1] 785 1006 892 $group [1] 2 3 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(627, 642, 690, 740, 775, 617, 638, 667, 696, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5skq61291982708.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,] -68 75 -217 -86.0 2 -436.0 -174.0 376 -150.0 -203 -110 -83 [2,] -57 75 -148 -55.0 99 -385.0 -145.0 448 -137.0 -120 -55 -11 [3,] -24 111 -108 -21.5 131 -373.5 -57.5 472 -99.5 -72 -12 54 [4,] 10 129 -64 86.0 172 -345.0 -51.0 560 8.0 -33 75 56 [5,] 69 129 0 135.0 172 -339.0 -48.0 604 60.0 51 150 64 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -67.21716 76.16826 -162.18271 -112.44955 83.91264 -399.3013 -118.133036 [2,] 19.21716 145.83174 -53.81729 69.44955 178.08736 -347.6987 3.133036 [,8] [,9] [,10] [,11] [,12] [1,] 399.7564 -193.029683 -128.11781 -95.8542 6.657969 [2,] 544.2436 -5.970317 -15.88219 71.8542 101.342031 $out [1] -8 221 287 $group [1] 2 2 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-68, -57, -24, 10, 69, 75, 75, 111, 129, 129, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6skq61291982708.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] [1,] 627.0 601.0 594.0 629.0 357.0 734.0 NA [2,] 634.0 618.0 623.5 630.5 623.0 754.5 NA [3,] 686.5 681.5 687.0 646.0 711.5 787.0 NA [4,] 757.0 741.5 722.0 724.5 843.0 899.0 NA [5,] 839.0 854.0 740.0 852.0 917.0 1006.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 630.3989 625.1708 642.0735 603.126 611.1565 721.0926 NA [2,] 742.6011 737.8292 731.9265 688.874 811.8435 852.9074 NA $out [1] 412 352 418 367 386 331 392 344 532 387 $group [1] 1 1 2 2 3 3 4 4 6 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(627, 634, 686.5, 757, 839, 601, 618, 681.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/72bq91291982708.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,] 677.3333 664.50 667.0000 [2,] 677.8333 665.00 669.8125 [3,] 697.2500 682.50 686.6250 [4,] 780.0833 759.25 770.3125 [5,] 845.0000 831.50 848.7500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 650.6131 639.512 640.7863 [2,] 743.8869 725.488 732.4637 $out [1] 445.1667 356.3333 415.0000 354.5000 449.8750 355.2500 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(677.333333333333, 677.833333333333, 697.25, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1zb931291982708.ps tmp/1zb931291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/2zb931291982708.ps tmp/2zb931291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/3zb931291982708.ps tmp/3zb931291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/4skq61291982708.ps tmp/4skq61291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/5skq61291982708.ps tmp/5skq61291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/6skq61291982708.ps tmp/6skq61291982708.png",intern=TRUE)) character(0) > try(system("convert tmp/72bq91291982708.ps tmp/72bq91291982708.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.540 1.150 2.712