R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-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. Natural language support but running in an English locale 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(30790,30800,31025,30835,31110,31270,31090,30755,31460,32135,32680,32700,32515,32275,32200,31835,31985,31875,31795,32260,33255,33160,32195,33130,33950,34210,33855,33735,34175,34265,33915,33660,33720,33810,33590,33545,33660,33165,33800,33880,33975,33930,33905,33890,33640,34395,34245,33940,34295,33745,33535,33715,33600,34120,34330,34130,33755,32910,32910,32850,32780,32565,31905,31975,31380,31355,31440,30310,31410,31300,31070,31075,31815) > 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] 73 > (np <- floor(n / par1)) [1] 6 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 7 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 30790 32515 33950 33660 34295 32780 31815 [2,] 30800 32275 34210 33165 33745 32565 NA [3,] 31025 32200 33855 33800 33535 31905 NA [4,] 30835 31835 33735 33880 33715 31975 NA [5,] 31110 31985 34175 33975 33600 31380 NA [6,] 31270 31875 34265 33930 34120 31355 NA [7,] 31090 31795 33915 33905 34330 31440 NA [8,] 30755 32260 33660 33890 34130 30310 NA [9,] 31460 33255 33720 33640 33755 31410 NA [10,] 32135 33160 33810 34395 32910 31300 NA [11,] 32680 32195 33590 34245 32910 31070 NA [12,] 32700 33130 33545 33940 32850 31075 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/www/html/freestat/rcomp/tmp/1udgs1304963919.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/2trkj1304963919.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/31tig1304963919.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4sfat1304963919.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] [1,] 30790 30800 31025.0 30835 31110.0 31270.0 31090 30310 31410.0 31300 31070 [2,] 32165 32275 31905.0 31835 31380.0 31355.0 31440 30755 31460.0 32135 32195 [3,] 32780 32865 32867.5 32845 32792.5 32902.5 32850 32960 33447.5 33035 32795 [4,] 33805 33745 33800.0 33735 33975.0 34120.0 33915 33890 33720.0 33810 33590 [5,] 34295 34210 33855.0 33880 34175.0 34265.0 34330 34130 33755.0 34395 34245 [,12] [1,] 32700 [2,] 32700 [3,] 32990 [4,] 33545 [5,] 33940 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 31800.62 31916.8 31645.16 31619.44 31118.64 31118.99 31253.55 30937.82 [2,] 33759.38 33813.2 34089.84 34070.56 34466.36 34686.01 34446.45 34982.18 [,9] [,10] [,11] [,12] [1,] 31989.73 31954.57 31895.18 32444.95 [2,] 34905.27 34115.43 33694.82 33535.05 $out [1] 31075 $group [1] 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(30790, 32165, 32780, 33805, 34295, 30800, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5fd491304963919.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,] 30755.0 31795.0 33545.0 33640.0 32850.0 30310.0 31815 [2,] 30817.5 31930.0 33690.0 33730.0 33222.5 31187.5 31815 [3,] 31100.0 32230.0 33832.5 33897.5 33730.0 31395.0 31815 [4,] 31797.5 32822.5 34062.5 33957.5 34125.0 31940.0 31815 [5,] 32700.0 33255.0 34265.0 34245.0 34330.0 32780.0 31815 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 30653.02 31822.92 33662.6 33793.74 33318.36 31051.78 31815 [2,] 31546.98 32637.08 34002.4 34001.26 34141.64 31738.22 31815 $out [1] 33165 34395 $group [1] 4 4 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(30755, 30817.5, 31100, 31797.5, 32700, 31795, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6oup01304963919.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 991.9056 2345 703.750 [2,] 1123.6678 2930 1285.625 [3,] 1231.2778 3080 1773.125 [4,] 1408.1483 3325 2366.875 [5,] 1663.3292 3820 2701.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1101.524 2899.838 1279.960 [2,] 1361.031 3260.162 2266.290 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(991.90557346285, 1123.66775790108, 1231.27779055085, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1udgs1304963919.ps tmp/1udgs1304963919.png",intern=TRUE)) character(0) > try(system("convert tmp/2trkj1304963919.ps tmp/2trkj1304963919.png",intern=TRUE)) character(0) > try(system("convert tmp/31tig1304963919.ps tmp/31tig1304963919.png",intern=TRUE)) character(0) > try(system("convert tmp/4sfat1304963919.ps tmp/4sfat1304963919.png",intern=TRUE)) character(0) > try(system("convert tmp/5fd491304963919.ps tmp/5fd491304963919.png",intern=TRUE)) character(0) > try(system("convert tmp/6oup01304963919.ps tmp/6oup01304963919.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.642 0.178 1.899