R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. 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(76,75,74,72,70,69,70,72,73,73,74,76,74,67,66,58,55,58,64,68,66,76,75,88,85,83,77,66,65,65,63,62,57,68,69,79,74,76,82,75,75,76,78,77,67,74,68,87,76,88,95,96,96,105,108,113,101,107,102,116,105,121,134,140,131,141,131,128,123,129,125,144,135,141,156,159,146,154,145,133,126,127,122,148) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 76 74 85 74 76 105 135 NA [2,] 75 67 83 76 88 121 141 NA [3,] 74 66 77 82 95 134 156 NA [4,] 72 58 66 75 96 140 159 NA [5,] 70 55 65 75 96 131 146 NA [6,] 69 58 65 76 105 141 154 NA [7,] 70 64 63 78 108 131 145 NA [8,] 72 68 62 77 113 128 133 NA [9,] 73 66 57 67 101 123 126 NA [10,] 73 76 68 74 107 129 127 NA [11,] 74 75 69 68 102 125 122 NA [12,] 76 88 79 87 116 144 148 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 -7 -2 2 12 16 6 NA [2,] -1 -1 -6 6 7 13 15 NA [3,] -2 -8 -11 -7 1 6 3 NA [4,] -2 -3 -1 0 0 -9 -13 NA [5,] -1 3 0 1 9 10 8 NA [6,] 1 6 -2 2 3 -10 -9 NA [7,] 2 4 -1 -1 5 -3 -12 NA [8,] 1 -2 -5 -10 -12 -5 -7 NA [9,] 0 10 11 7 6 6 1 NA [10,] 1 -1 1 -6 -5 -4 -5 NA [11,] 2 13 10 19 14 19 26 NA [12,] -2 -3 -5 -11 -11 -9 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/freestat/rcomp/tmp/1iw6v1282169093.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/freestat/rcomp/tmp/2iw6v1282169093.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/freestat/rcomp/tmp/3iw6v1282169093.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/freestat/rcomp/tmp/4fqd41282169093.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] [,11] [,12] [1,] 74 67.0 66.0 58 55.0 58 63.0 62.0 57.0 68.0 68.0 76 [2,] 75 75.5 75.5 69 67.5 67 67.0 70.0 66.5 73.5 71.5 83 [3,] 76 83.0 82.0 75 75.0 76 78.0 77.0 73.0 76.0 75.0 88 [4,] 95 104.5 114.5 118 113.5 123 119.5 120.5 112.0 117.0 112.0 130 [5,] 105 141.0 156.0 159 146.0 154 145.0 133.0 126.0 129.0 125.0 148 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 64.05632 65.68167 58.70983 45.73799 47.52954 42.5577 46.64785 [2,] 87.94368 100.31833 105.29017 104.26201 102.47046 109.4423 109.35215 [,8] [,9] [,10] [,11] [,12] [1,] 46.84221 45.82813 50.0225 50.81405 59.93236 [2,] 107.15779 100.17187 101.9775 99.18595 116.06764 $out [1] 135 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(74, 75, 76, 95, 105, 67, 75.5, 83, 104.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5fqd41282169093.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,] -7.0 -6 -11.0 -13.0 -1.0 -10.0 -3 -12.0 0.0 -6 2.0 -11 [2,] -1.5 -1 -7.5 -6.0 0.5 -5.5 -2 -8.5 3.5 -5 11.5 -11 [3,] 2.0 6 -2.0 -2.0 3.0 1.0 -1 -5.0 6.0 -4 14.0 -7 [4,] 9.0 10 2.0 -0.5 8.5 2.5 3 -3.5 8.5 0 19.0 -3 [5,] 16.0 15 6.0 0.0 10.0 6.0 5 1.0 11.0 1 26.0 -2 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4.270431 -0.5690225 -7.673247 -5.284511 -1.777471 -3.777471 -3.985919 [2,] 8.270431 12.5690225 3.673247 1.284511 7.777471 5.777471 1.985919 [,8] [,9] [,10] [,11] [,12] [1,] -7.985919 3.014081 -6.985919 9.521121 -12.160258 [2,] -2.014081 8.985919 -1.014081 18.478879 -1.839742 $out [1] -12 $group [1] 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-7, -1.5, 2, 9, 16, -6, -1, 6, 10, 15, -11, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/63wmi1282169093.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] [1,] 69.0 55.0 57 74.0 88.0 105 122 NA [2,] 71.0 61.0 64 74.0 95.5 124 130 NA [3,] 73.0 66.5 67 75.5 101.5 130 143 NA [4,] 74.5 74.5 78 77.5 107.5 137 151 NA [5,] 76.0 88.0 85 82.0 116.0 144 159 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 71.40363 60.34256 60.6145 73.90363 96.02672 124.0706 133.4218 NA [2,] 74.59637 72.65744 73.3855 77.09637 106.97328 135.9294 152.5782 NA $out [1] 67 68 87 76 $group [1] 4 4 4 5 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(69, 71, 73, 74.5, 76, 55, 61, 66.5, 74.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/73wmi1282169093.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,] 87.57143 73 85.000 [2,] 90.92857 75 90.250 [3,] 93.35714 76 93.375 [4,] 95.28571 80 95.125 [5,] 97.71429 83 95.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 91.36982 73.71947 91.15148 [2,] 95.34446 78.28053 95.59852 $out [1] 105.4286 88.0000 106.5000 $group [1] 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(87.5714285714286, 90.9285714285714, 93.3571428571429, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1iw6v1282169093.ps tmp/1iw6v1282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/2iw6v1282169093.ps tmp/2iw6v1282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/3iw6v1282169093.ps tmp/3iw6v1282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/4fqd41282169093.ps tmp/4fqd41282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/5fqd41282169093.ps tmp/5fqd41282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/63wmi1282169093.ps tmp/63wmi1282169093.png",intern=TRUE)) character(0) > try(system("convert tmp/73wmi1282169093.ps tmp/73wmi1282169093.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.008 1.523 2.328