R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(613.20,614.70,618.40,628.20,629.00,629.70,630.40,630.40,639.30,639.40,640.90,640.80,642.10,645.30,647.60,648.40,648.80,648.90,648.90,648.90,650.30,650.30,650.00,650.00,650.50,658.40,666.00,675.50,680.70,690.60,690.60,691.10,692.90,693.80,692.80,697.50,699.00,702.10,704.80,715.50,721.80,726.40,727.70,727.40,731.30,734.40,733.40,733.40,738.10,742.60,747.20,751.10,752.60,758.90,759.10,764.30,765.60,767.60,767.60,765.60,768.20,770.90,775.10,777.60,778.60,778.90,779.40,779.90,781.70,789.10,788.70,788.80,790.80,794.10,795.10,797.30,803.80,805.60,804.60,804.50,805.80,806.80,805.20,814.90,816.60,819.50,823.00,824.00,831.40,831.70,831.10,832.10,833.30,838.80,838.00,837.30,994.20,994.20,994.20,994.20,994.20,1092.60,1100.00,1100.00,1092.60,1000.70,1000.70,1000.50,1000.50,1000.50,1000.50,1000.50,1000.50,1087.70,1113.20,1116.00,1085.20,1031.30,1028.70,1027.50,1027.50,1027.50,1027.50,1027.50,1027.50,1152.20,1155.30,1154.00,1119.90,1079.30,1074.30,1069.80) > 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] 132 > (np <- floor(n / par1)) [1] 11 > 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] 11 11 11 11 11 11 11 11 11 11 11 11 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 613.2 642.1 650.5 699.0 738.1 768.2 790.8 816.6 994.2 1000.5 1027.5 [2,] 614.7 645.3 658.4 702.1 742.6 770.9 794.1 819.5 994.2 1000.5 1027.5 [3,] 618.4 647.6 666.0 704.8 747.2 775.1 795.1 823.0 994.2 1000.5 1027.5 [4,] 628.2 648.4 675.5 715.5 751.1 777.6 797.3 824.0 994.2 1000.5 1027.5 [5,] 629.0 648.8 680.7 721.8 752.6 778.6 803.8 831.4 994.2 1000.5 1027.5 [6,] 629.7 648.9 690.6 726.4 758.9 778.9 805.6 831.7 1092.6 1087.7 1152.2 [7,] 630.4 648.9 690.6 727.7 759.1 779.4 804.6 831.1 1100.0 1113.2 1155.3 [8,] 630.4 648.9 691.1 727.4 764.3 779.9 804.5 832.1 1100.0 1116.0 1154.0 [9,] 639.3 650.3 692.9 731.3 765.6 781.7 805.8 833.3 1092.6 1085.2 1119.9 [10,] 639.4 650.3 693.8 734.4 767.6 789.1 806.8 838.8 1000.7 1031.3 1079.3 [11,] 640.9 650.0 692.8 733.4 767.6 788.7 805.2 838.0 1000.7 1028.7 1074.3 [12,] 640.8 650.0 697.5 733.4 765.6 788.8 814.9 837.3 1000.5 1027.5 1069.8 [,12] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 1.5 3.2 7.9 3.1 4.5 2.7 3.3 2.9 0.0 0.0 0.0 NA [2,] 3.7 2.3 7.6 2.7 4.6 4.2 1.0 3.5 0.0 0.0 0.0 NA [3,] 9.8 0.8 9.5 10.7 3.9 2.5 2.2 1.0 0.0 0.0 0.0 NA [4,] 0.8 0.4 5.2 6.3 1.5 1.0 6.5 7.4 0.0 0.0 0.0 NA [5,] 0.7 0.1 9.9 4.6 6.3 0.3 1.8 0.3 98.4 87.2 124.7 NA [6,] 0.7 0.0 0.0 1.3 0.2 0.5 -1.0 -0.6 7.4 25.5 3.1 NA [7,] 0.0 0.0 0.5 -0.3 5.2 0.5 -0.1 1.0 0.0 2.8 -1.3 NA [8,] 8.9 1.4 1.8 3.9 1.3 1.8 1.3 1.2 -7.4 -30.8 -34.1 NA [9,] 0.1 0.0 0.9 3.1 2.0 7.4 1.0 5.5 -91.9 -53.9 -40.6 NA [10,] 1.5 -0.3 -1.0 -1.0 0.0 -0.4 -1.6 -0.8 0.0 -2.6 -5.0 NA [11,] -0.1 0.0 4.7 0.0 -2.0 0.1 9.7 -0.7 -0.2 -1.2 -4.5 NA [12,] 1.3 0.5 1.5 4.7 2.6 2.0 1.7 156.9 0.0 0.0 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/rcomp/tmp/19srn1250084328.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/2g72o1250084328.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/3c46j1250084328.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/494mt1250084328.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] [1,] 613.20 614.70 618.4 628.2 629.00 629.7 630.40 630.40 639.30 [2,] 674.75 680.25 685.4 695.5 701.25 708.5 709.15 709.25 712.10 [3,] 768.20 770.90 775.1 777.6 778.60 778.9 779.40 779.90 781.70 [4,] 905.40 906.85 908.6 909.1 912.80 959.7 965.55 966.05 959.25 [5,] 1027.50 1027.50 1027.5 1027.5 1027.50 1152.2 1155.30 1154.00 1119.90 [,10] [,11] [,12] [1,] 639.40 640.90 640.80 [2,] 714.10 713.10 715.45 [3,] 789.10 788.70 788.80 [4,] 919.75 919.35 918.90 [5,] 1079.30 1074.30 1069.80 $n [1] 11 11 11 11 11 11 11 11 11 11 11 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 658.3211 662.9505 668.7702 675.8435 677.8201 659.2314 657.2541 657.5636 [2,] 878.0789 878.8495 881.4298 879.3565 879.3799 898.5686 901.5459 902.2364 [,9] [,10] [,11] [,12] [1,] 663.9607 691.1308 690.445 691.8789 [2,] 899.4393 887.0692 886.955 885.7211 $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(613.2, 674.75, 768.2, 905.4, 1027.5, 614.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5uyum1250084328.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,] 0.00 0.00 0.0 0.00 0.10 -1.0 -0.30 -7.4 -53.90 -2.60 -2.00 0.0 [2,] 0.75 0.50 0.4 0.20 0.50 0.0 -0.05 -3.1 -20.30 -1.30 -0.95 0.5 [3,] 2.90 2.70 2.2 1.00 4.60 0.5 0.00 1.3 0.90 -0.80 -0.10 1.6 [4,] 3.25 3.95 6.7 5.75 48.55 2.2 0.75 1.8 2.55 -0.15 0.05 2.6 [5,] 4.50 7.60 10.7 7.40 98.40 3.1 1.00 8.9 7.40 1.50 0.10 4.7 $n [1] 11 11 11 11 11 11 11 11 11 11 11 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.709030 1.056462 -0.801244 -1.643953 -18.29044 -0.5480534 -0.3811103 [2,] 4.090970 4.343538 5.201244 3.643953 27.49044 1.5480534 0.3811103 [,8] [,9] [,10] [,11] [,12] [1,] -1.034301 -9.985464 -1.3478461 -0.5763879 0.5507563 [2,] 3.634301 11.785464 -0.2521539 0.3763879 2.6492437 $out [1] 7.9 124.7 7.4 25.5 5.2 2.8 -1.3 -30.8 -34.1 -91.9 -5.0 4.7 [13] 9.7 -4.5 156.9 $group [1] 1 5 6 6 7 7 7 8 8 9 10 11 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0, 0.75, 2.89999999999998, 3.25, 4.5, 0, 0.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6getn1250084328.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] [1,] 613.20 645.3 650.50 699.00 738.10 768.20 790.80 816.60 994.2 1000.50 [2,] 623.30 648.0 670.75 710.15 749.15 776.35 796.20 823.50 994.2 1000.50 [3,] 630.05 648.9 690.60 726.90 759.00 779.15 804.55 831.55 1000.6 1028.10 [4,] 639.35 650.0 692.85 732.35 765.60 785.20 805.70 835.30 1092.6 1086.45 [5,] 640.90 650.3 697.50 734.40 767.60 789.10 814.90 838.80 1100.0 1116.00 [,11] [,12] [1,] 1027.50 NA [2,] 1027.50 NA [3,] 1072.05 NA [4,] 1136.05 NA [5,] 1155.30 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 622.7295 647.9878 680.52 716.7744 751.497 775.1135 800.217 826.168 [2,] 637.3705 649.8122 700.68 737.0256 766.503 783.1865 808.883 836.932 [,9] [,10] [,11] [,12] [1,] 955.7191 988.8976 1022.540 NA [2,] 1045.4809 1067.3024 1121.560 NA $out [1] 642.1 $group [1] 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" NA Warning message: In bxp(list(stats = c(613.2, 623.3, 630.05, 639.35, 640.9, 645.3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/70x5o1250084328.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,] 794.6091 768.20 790.0750 [2,] 801.7818 776.35 799.6500 [3,] 820.2909 779.15 816.5750 [4,] 836.4136 785.20 834.8875 [5,] 840.7818 789.10 837.6500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 804.4951 775.1135 800.503 [2,] 836.0867 783.1865 832.647 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(794.609090909091, 801.781818181818, 820.290909090909, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/19srn1250084328.ps tmp/19srn1250084328.png") > system("convert tmp/2g72o1250084328.ps tmp/2g72o1250084328.png") > system("convert tmp/3c46j1250084328.ps tmp/3c46j1250084328.png") > system("convert tmp/494mt1250084328.ps tmp/494mt1250084328.png") > system("convert tmp/5uyum1250084328.ps tmp/5uyum1250084328.png") > system("convert tmp/6getn1250084328.ps tmp/6getn1250084328.png") > system("convert tmp/70x5o1250084328.ps tmp/70x5o1250084328.png") > > > proc.time() user system elapsed 1.183 0.922 2.071