R version 2.12.1 (2010-12-16) 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(505.7,55.7,735.7,575.9,545.8,905.8,765.8,945.7,15.7,645.7,155.9,416,825.8,725.9,925.9,556,116.1,876.3,336.2,186.1,286.1,26,915.8,405.7,965.7,395.6,425.8,545.6,65.6,445.6,895.5,175.4,715.4,865.5,57.4,145.4,315.3,635.4,5.2,515.2,515.1,955,955,634.9,205,275,425,84.9,534.7,4.8,704.7,684.7,884.6,994.6,294.7,524.7,914.5,564.4,984.5,934.4,514.6,474.5,784.4,504.5,824.4,414.6,964.7,64.6,244.7,344.7,34.7,685,425,484.8,785.1,704.9,245.4,285.6,218.8,706.1,856.2,456.6,606.8,527.3,657.8,948.2,486.6,238.9,289.4,969.5,589.5,189.7,639.8,9710.1,969.9,939.9,859.7,679.9,879.9,329.8,349.6,39.5,849.5,449.6,749.6,249.7,649.8,619.4,939,778.9) > 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] 110 > (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 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 505.7 825.8 965.7 315.3 534.7 514.6 425.0 657.8 859.7 939.0 [2,] 55.7 725.9 395.6 635.4 4.8 474.5 484.8 948.2 679.9 778.9 [3,] 735.7 925.9 425.8 5.2 704.7 784.4 785.1 486.6 879.9 NA [4,] 575.9 556.0 545.6 515.2 684.7 504.5 704.9 238.9 329.8 NA [5,] 545.8 116.1 65.6 515.1 884.6 824.4 245.4 289.4 349.6 NA [6,] 905.8 876.3 445.6 955.0 994.6 414.6 285.6 969.5 39.5 NA [7,] 765.8 336.2 895.5 955.0 294.7 964.7 218.8 589.5 849.5 NA [8,] 945.7 186.1 175.4 634.9 524.7 64.6 706.1 189.7 449.6 NA [9,] 15.7 286.1 715.4 205.0 914.5 244.7 856.2 639.8 749.6 NA [10,] 645.7 26.0 865.5 275.0 564.4 344.7 456.6 9710.1 249.7 NA [11,] 155.9 915.8 57.4 425.0 984.5 34.7 606.8 969.9 649.8 NA [12,] 416.0 405.7 145.4 84.9 934.4 685.0 527.3 939.9 619.4 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -450.0 -99.9 -570.1 320.1 -529.9 -40.1 59.8 290.4 -179.8 -160.1 [2,] 680.0 200.0 30.2 -630.2 699.9 309.9 300.3 -461.6 200.0 NA [3,] -159.8 -369.9 119.8 510.0 -20.0 -279.9 -80.2 -247.7 -550.1 NA [4,] -30.1 -439.9 -480.0 -0.1 199.9 319.9 -459.5 50.5 19.8 NA [5,] 360.0 760.2 380.0 439.9 110.0 -409.8 40.2 680.1 -310.1 NA [6,] -140.0 -540.1 449.9 0.0 -699.9 550.1 -66.8 -380.0 810.0 NA [7,] 179.9 -150.1 -720.1 -320.1 230.0 -900.1 487.3 -399.8 -399.9 NA [8,] -930.0 100.0 540.0 -429.9 389.8 180.1 150.1 450.1 300.0 NA [9,] 630.0 -260.1 150.1 70.0 -350.1 100.0 -399.6 9070.3 -499.9 NA [10,] -489.8 889.8 -808.1 150.0 420.1 -310.0 150.2 -8740.2 400.1 NA [11,] 260.1 -510.1 88.0 -340.1 -50.1 650.3 -79.5 -30.0 -30.4 NA [12,] 409.8 560.0 169.9 449.8 -419.8 -260.0 130.5 -80.2 319.6 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/1vfhh1306483679.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/20mzq1306483679.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/3x8ox1306483679.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/43njw1306483679.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,] 315.30 4.8 425.8 504.5 65.6 39.5 218.8 64.6 15.7 26.0 34.7 84.9 [2,] 505.70 395.6 486.6 504.5 245.4 414.6 336.2 186.1 244.7 275.0 155.9 405.7 [3,] 596.25 560.1 735.7 545.6 349.6 876.3 765.8 449.6 639.8 456.6 606.8 527.3 [4,] 859.70 725.9 785.1 575.9 545.8 955.0 895.5 634.9 749.6 645.7 915.8 685.0 [5,] 965.70 948.2 925.9 575.9 884.6 994.6 964.7 945.7 914.5 865.5 984.5 939.9 $n [1] 10 10 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 419.3775 395.0690 578.49 507.996 191.3893 591.6893 471.2353 213.232 [2,] 773.1225 725.1310 892.91 583.204 507.8107 1160.9107 1060.3647 685.968 [,9] [,10] [,11] [,12] [1,] 373.886 261.3647 206.586 380.202 [2,] 905.714 651.8353 1007.014 674.398 $out [1] 5.2 684.7 704.9 238.9 329.8 9710.1 $group [1] 3 4 4 4 4 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(315.3, 505.7, 596.25, 859.7, 965.7, 4.8, 395.6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/56u2o1306483679.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] [1,] -570.1 30.2 -550.1 -480.0 -409.8 -699.9 -900.1 100.0 -499.9 -808.1 -79.5 [2,] -450.0 30.2 -279.9 -439.9 40.2 -380.0 -399.9 100.0 -350.1 -489.8 -79.5 [3,] -130.0 200.0 -159.8 -0.1 360.0 -66.8 -320.1 180.1 70.0 150.0 -30.4 [4,] 59.8 309.9 -20.0 50.5 439.9 449.9 179.9 389.8 150.1 400.1 88.0 [5,] 320.1 699.9 119.8 319.9 760.2 810.0 487.3 540.0 630.0 889.8 260.1 [,12] [1,] -419.8 [2,] -80.2 [3,] 169.9 [4,] 409.8 [5,] 560.0 $n [1] 10 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -384.7164 52.69133 -296.68067 -258.3773 149.4913 -503.8807 -625.46133 [2,] 124.7164 347.30867 -22.91933 258.1773 570.5087 370.2807 -14.73867 [,8] [,9] [,10] [,11] [,12] [1,] 27.472 -193.4387 -318.6807 -118.61667 -88.16667 [2,] 332.728 333.4387 618.6807 57.81667 427.96667 $out [1] -630.2 -461.6 510.0 -930.0 -429.9 9070.3 -8740.2 -510.1 -340.1 [10] 650.3 $group [1] 2 2 3 8 8 9 10 11 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-570.1, -450, -130, 59.8, 320.1, 30.2, 30.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6d86v1306483679.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] [,10] [1,] 15.70 26.00 57.40 5.20 4.80 34.7 218.80 189.70 39.50 778.90 [2,] 285.95 236.10 160.40 240.00 529.70 294.7 355.30 388.00 339.70 778.90 [3,] 560.85 480.85 435.70 470.05 694.70 489.5 506.05 648.80 634.60 858.95 [4,] 750.75 851.05 790.45 635.15 924.45 734.7 705.50 958.85 799.55 939.00 [5,] 945.70 925.90 965.70 955.00 994.60 964.7 856.20 969.90 879.90 939.00 $n [1] 12 12 12 12 12 12 12 12 12 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 348.8516 200.3672 148.3300 289.8194 514.6519 288.8130 346.3214 388.4315 [2,] 772.8484 761.3328 723.0700 650.2806 874.7481 690.1870 665.7786 909.1685 [,9] [,10] [1,] 424.8593 680.0817 [2,] 844.3407 1037.8183 $out [1] 9710.1 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(15.7, 285.95, 560.85, 750.75, 945.7, 26, 236.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7rz7r1306483679.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,] 426.2222 349.600 395.600 [2,] 515.6944 491.950 478.750 [3,] 530.9889 578.175 542.775 [4,] 653.1222 687.750 625.850 [5,] 654.3300 876.300 684.800 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 468.3072 488.8693 475.6817 [2,] 593.6706 667.4807 609.8683 $out [1] 1459.744 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(426.222222222222, 515.694444444444, 530.988888888889, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1vfhh1306483679.ps tmp/1vfhh1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/20mzq1306483679.ps tmp/20mzq1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/3x8ox1306483679.ps tmp/3x8ox1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/43njw1306483679.ps tmp/43njw1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/56u2o1306483679.ps tmp/56u2o1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/6d86v1306483679.ps tmp/6d86v1306483679.png",intern=TRUE)) character(0) > try(system("convert tmp/7rz7r1306483679.ps tmp/7rz7r1306483679.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.936 0.388 2.309