R version 2.13.0 (2011-04-13) Copyright (C) 2011 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(616.38,614.62,618,621.38,642.62,658.12,639,635.38,634.62,624.38,626.12,629.25,647.62,646.12,642.62,621.38,635.38,642.38,633.75,630.62,630.88,617,629.62,624.88,621.38,623.88,621.75,623.88,657.5,614.62,612.12,608.88,615.12,626.38,611.12,615.88,612.38,606,658.5,632.75,652.38,654.5,646.12,642.62,637.62,662.12,673.12,676.62,688.5,693.38,697.75,704,703.88,731.62,728,729.12,738.12,753.5,763.62,765.62,775.88,765.12,788.12,785.88,789.75,770.62,762,754.38,746.62,755,737.62,732.12,731.75,706.12,725.5,725.62,710.38,711.12,703.62,703.12,684.5,682.88,651.88,673.62,677.75,673,704.88,715.38,681,674.88,672.88,697.38,699.88,692.12,687.62,683.62,693.38,680.75,687.38,700.62,693.62,681.38,695.38,707.88,716,667.5,638.12,647.38,650,635.25,624.12,635.25,625.62,611.62) > 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] 114 > (np <- floor(n / par1)) [1] 9 > 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 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 616.38 647.62 621.38 612.38 688.50 775.88 731.75 677.75 693.38 650.00 [2,] 614.62 646.12 623.88 606.00 693.38 765.12 706.12 673.00 680.75 635.25 [3,] 618.00 642.62 621.75 658.50 697.75 788.12 725.50 704.88 687.38 624.12 [4,] 621.38 621.38 623.88 632.75 704.00 785.88 725.62 715.38 700.62 635.25 [5,] 642.62 635.38 657.50 652.38 703.88 789.75 710.38 681.00 693.62 625.62 [6,] 658.12 642.38 614.62 654.50 731.62 770.62 711.12 674.88 681.38 611.62 [7,] 639.00 633.75 612.12 646.12 728.00 762.00 703.62 672.88 695.38 NA [8,] 635.38 630.62 608.88 642.62 729.12 754.38 703.12 697.38 707.88 NA [9,] 634.62 630.88 615.12 637.62 738.12 746.62 684.50 699.88 716.00 NA [10,] 624.38 617.00 626.38 662.12 753.50 755.00 682.88 692.12 667.50 NA [11,] 626.12 629.62 611.12 673.12 763.62 737.62 651.88 687.62 638.12 NA [12,] 629.25 624.88 615.88 676.62 765.62 732.12 673.62 683.62 647.38 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -1.76 -1.50 2.50 -6.38 4.88 -10.76 -25.63 -4.75 -12.63 -14.75 [2,] 3.38 -3.50 -2.13 52.50 4.37 23.00 19.38 31.88 6.63 -11.13 [3,] 3.38 -21.24 2.13 -25.75 6.25 -2.24 0.12 10.50 13.24 11.13 [4,] 21.24 14.00 33.62 19.63 -0.12 3.87 -15.24 -34.38 -7.00 -9.63 [5,] 15.50 7.00 -42.88 2.12 27.74 -19.13 0.74 -6.12 -12.24 -14.00 [6,] -19.12 -8.63 -2.50 -8.38 -3.62 -8.62 -7.50 -2.00 14.00 NA [7,] -3.62 -3.13 -3.24 -3.50 1.12 -7.62 -0.50 24.50 12.50 NA [8,] -0.76 0.26 6.24 -5.00 9.00 -7.76 -18.62 2.50 8.12 NA [9,] -10.24 -13.88 11.26 24.50 15.38 8.38 -1.62 -7.76 -48.50 NA [10,] 1.74 12.62 -15.26 11.00 10.12 -17.38 -31.00 -4.50 -29.38 NA [11,] 3.13 -4.74 4.76 3.50 2.00 -5.50 21.74 -4.00 9.26 NA [12,] 18.37 -3.50 -3.50 11.88 10.26 -0.37 4.13 9.76 2.62 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/wessaorg/rcomp/tmp/1nqhn1321368568.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/wessaorg/rcomp/tmp/2qy2w1321368568.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/wessaorg/rcomp/tmp/3b5xx1321368568.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/wessaorg/rcomp/tmp/4ytnv1321368568.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] [1,] 612.380 606.00 618.00 621.380 625.62 611.62 612.12 608.88 615.12 617.00 [2,] 621.380 623.88 624.12 623.880 642.62 642.38 639.00 635.38 634.62 626.38 [3,] 663.875 659.56 672.94 667.935 669.25 666.50 672.88 697.38 684.50 667.50 [4,] 693.380 693.38 704.88 715.380 703.88 711.12 703.62 707.88 716.00 692.12 [5,] 775.880 765.12 788.12 785.880 789.75 770.62 762.00 754.38 746.62 755.00 [,11] [,12] [1,] 611.12 615.88 [2,] 629.62 629.25 [3,] 651.88 673.62 [4,] 687.62 683.62 [5,] 763.62 732.12 $n [1] 10 10 10 10 10 10 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 627.9009 624.835 632.5891 622.218 638.6421 632.1548 638.8468 659.1967 [2,] 699.8491 694.285 713.2909 713.652 699.8579 700.8452 706.9132 735.5633 [,9] [,10] [,11] [,12] [1,] 641.6399 632.8769 621.3333 644.9851 [2,] 727.3601 702.1231 682.4267 702.2549 $out [1] 765.62 $group [1] 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(612.38, 621.38, 663.875, 693.38, 775.88, 606, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/54bg41321368568.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] [1,] -25.630 -11.13 -21.240 -34.380 -42.88 -8.63 -7.62 -18.62 -13.88 -31.00 [2,] -12.630 -2.13 -2.240 -9.630 -14.00 -8.62 -3.50 -5.00 -10.24 -17.38 [3,] -5.565 5.50 2.755 1.875 -2.69 -7.50 -3.13 0.26 -1.62 -4.50 [4,] -1.500 23.00 10.500 19.630 7.00 -2.50 1.12 6.24 11.26 10.12 [5,] 4.880 52.50 13.240 33.620 27.74 -2.00 1.12 9.00 24.50 12.62 [,11] [,12] [1,] -5.50 -3.50 [2,] -4.00 -0.37 [3,] 3.13 4.13 [4,] 4.76 10.26 [5,] 9.26 18.37 $n [1] 10 10 10 10 10 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -11.125991757 -7.05595 -3.610412 -12.74446 -13.182437 -10.7232 -5.5632 [2,] -0.004008243 18.05595 9.120412 16.49446 7.802437 -4.2768 -0.6968 [,8] [,9] [,10] [,11] [,12] [1,] -5.659733 -12.943333 -18.983333 -1.4836 -1.468467 [2,] 6.179733 9.703333 9.983333 7.7436 9.728467 $out [1] -25.75 -19.12 14.00 24.50 12.50 -48.50 21.74 $group [1] 3 6 6 7 7 9 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-25.63, -12.63, -5.565, -1.5, 4.88, -11.13, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6o6pf1321368568.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] [1,] 614.620 617.000 608.88 606.000 688.500 732.12 651.88 672.880 647.380 [2,] 619.690 627.250 613.37 635.185 700.815 750.50 683.69 676.315 674.125 [3,] 627.685 632.315 618.63 649.250 728.560 763.56 704.87 685.620 690.380 [4,] 637.190 642.500 623.88 660.310 745.810 780.88 718.31 698.630 698.000 [5,] 658.120 647.620 626.38 676.620 765.620 789.75 731.75 715.380 716.000 [,10] [1,] 611.620 [2,] 624.120 [3,] 630.435 [4,] 635.250 [5,] 650.000 $n [1] 12 12 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 619.7031 625.3594 613.8363 637.7903 708.0375 749.7035 689.0796 675.442 [2,] 635.6669 639.2706 623.4237 660.7097 749.0825 777.4165 720.6604 695.798 [,9] [,10] [1,] 679.4905 623.2558 [2,] 701.2695 637.6142 $out [1] 657.50 638.12 $group [1] 3 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(614.62, 619.69, 627.685, 637.19, 658.12, 617, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/741lq1321368568.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,] 664.4240 659.5600 656.4350 [2,] 671.8060 665.1875 658.9350 [3,] 676.1337 668.5925 667.6187 [4,] 677.5683 673.2800 672.4088 [5,] 679.2130 684.5000 675.3100 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 673.5054 664.9015 661.4733 [2,] 678.7619 672.2835 673.7642 $out [1] 697.38 651.88 $group [1] 2 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(664.424, 671.806, 676.133666666667, 677.568333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1nqhn1321368568.ps tmp/1nqhn1321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/2qy2w1321368568.ps tmp/2qy2w1321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/3b5xx1321368568.ps tmp/3b5xx1321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/4ytnv1321368568.ps tmp/4ytnv1321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/54bg41321368568.ps tmp/54bg41321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/6o6pf1321368568.ps tmp/6o6pf1321368568.png",intern=TRUE)) character(0) > try(system("convert tmp/741lq1321368568.ps tmp/741lq1321368568.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.542 0.315 1.871