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(100.34,115.78,114.6,114.2,115.88,125.22,161.71,165.01,135.78,153.67,125.52,135.29,103.05,120.79,120.17,119.62,121.17,129.86,167.8,167.14,140.55,158.44,131.07,140.55,106.15,123.65,122.8,122.25,123.88,132.96,171.82,173.69,149.5,164.44,133.37,143.77,69.49,84.5,82.3,78.8,79.47,88.93,138.13,139.69,114.43,128.65,95.92,98.22,56.65,69.6,66.91,63.76,64,35.24,45.3,43.02,43.08,43.17,46.38,70.85,72.81,59.51,67.54,56.51,53.82,112.55,127.65,126.51,126.08,127.34) > 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] 70 > (np <- floor(n / par1)) [1] 5 > 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] 6 6 6 6 6 6 6 6 6 6 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.34 103.05 106.15 69.49 56.65 72.81 [2,] 115.78 120.79 123.65 84.50 69.60 59.51 [3,] 114.60 120.17 122.80 82.30 66.91 67.54 [4,] 114.20 119.62 122.25 78.80 63.76 56.51 [5,] 115.88 121.17 123.88 79.47 64.00 53.82 [6,] 125.22 129.86 132.96 88.93 35.24 112.55 [7,] 161.71 167.80 171.82 138.13 45.30 127.65 [8,] 165.01 167.14 173.69 139.69 43.02 126.51 [9,] 135.78 140.55 149.50 114.43 43.08 126.08 [10,] 153.67 158.44 164.44 128.65 43.17 127.34 [11,] 125.52 131.07 133.37 95.92 46.38 NA [12,] 135.29 140.55 143.77 98.22 70.85 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/rcomp/tmp/1rof51262534858.ps",horizontal=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/rcomp/tmp/28hzv1262534858.ps",horizontal=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/rcomp/tmp/3lbv61262534858.ps",horizontal=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/rcomp/tmp/4k8wf1262534858.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] [,10] [1,] 56.650 59.51 66.91 56.51 53.820 35.240 127.65 126.51 114.43 127.34 [2,] 69.490 69.60 67.54 63.76 64.000 88.930 127.65 126.51 114.43 127.34 [3,] 86.575 100.14 98.45 96.50 97.675 118.885 149.92 152.35 130.93 141.16 [4,] 103.050 120.79 120.17 119.62 121.170 129.860 167.80 167.14 140.55 158.44 [5,] 106.150 123.65 122.80 122.25 123.880 132.960 171.82 173.69 149.50 164.44 [,11] [,12] [1,] 46.38 70.85 [2,] 95.92 98.22 [3,] 125.52 135.29 [4,] 131.07 140.55 [5,] 133.37 143.77 $n [1] 6 6 6 6 6 6 6 6 6 6 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 64.92772 67.1208 64.50195 60.4685 60.7985 92.48383 124.0220 126.1423 [2,] 108.22228 133.1592 132.39805 132.5315 134.5515 145.28617 175.8180 178.5577 [,9] [,10] [,11] [,12] [1,] 114.0818 121.0995 100.6831 105.3797 [2,] 147.7782 161.2205 150.3569 165.2003 $out [1] 45.30 43.02 43.08 43.17 $group [1] 7 8 9 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(56.65, 69.49, 86.575, 103.05, 106.15, 59.51, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5tjj11262534858.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] [1,] 100.340 103.050 106.150 69.490 35.240 53.82 [2,] 115.190 120.480 123.225 80.885 43.125 59.51 [3,] 125.370 130.465 133.165 92.425 51.515 92.68 [4,] 144.725 149.495 156.970 121.540 65.455 126.51 [5,] 165.010 167.800 173.690 139.690 70.850 127.65 $n [1] 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 111.8989 117.2311 117.7737 73.88198 41.33014 59.20413 [2,] 138.8411 143.6989 148.5563 110.96802 61.69986 126.15587 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(100.34, 115.19, 125.37, 144.725, 165.01, 103.05, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6djav1262534859.ps",horizontal=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,] 20.98065 49.500 22.01500 [2,] 28.85138 64.940 32.95875 [3,] 34.37420 79.955 36.40500 [4,] 41.84746 113.845 46.88000 [5,] 48.95023 130.670 51.98000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 28.44660 57.6491 30.05542 [2,] 40.30180 102.2609 42.75458 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(20.9806543431483, 28.8513833671048, 34.3742010263974, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1rof51262534858.ps tmp/1rof51262534858.png",intern=TRUE)) character(0) > try(system("convert tmp/28hzv1262534858.ps tmp/28hzv1262534858.png",intern=TRUE)) character(0) > try(system("convert tmp/3lbv61262534858.ps tmp/3lbv61262534858.png",intern=TRUE)) character(0) > try(system("convert tmp/4k8wf1262534858.ps tmp/4k8wf1262534858.png",intern=TRUE)) character(0) > try(system("convert tmp/5tjj11262534858.ps tmp/5tjj11262534858.png",intern=TRUE)) character(0) > try(system("convert tmp/6djav1262534859.ps tmp/6djav1262534859.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.078 0.832 1.352