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(213,212,211,209,229,228,213,203,204,204,205,207,205,208,201,201,220,218,203,185,179,182,182,185,183,192,177,172,188,182,162,150,141,135,139,148,142,156,143,134,146,142,117,106,104,99,105,106,96,104,96,85,91,98,73,70,62,60,70,82,72,73,68,53,61,73,46,50,52,45,58,73,58,49,44,35,46,61,29,33,37,31,44,57,42,34,27,22,30,47,12,13,18,11,26,41,21,24,30,34,48,64,35,44,55,53,73,94,73,78,87,87,91,104,73,84,103,111,131,155) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 213 205 183 142 96 72 58 42 21 73 NA [2,] 212 208 192 156 104 73 49 34 24 78 NA [3,] 211 201 177 143 96 68 44 27 30 87 NA [4,] 209 201 172 134 85 53 35 22 34 87 NA [5,] 229 220 188 146 91 61 46 30 48 91 NA [6,] 228 218 182 142 98 73 61 47 64 104 NA [7,] 213 203 162 117 73 46 29 12 35 73 NA [8,] 203 185 150 106 70 50 33 13 44 84 NA [9,] 204 179 141 104 62 52 37 18 55 103 NA [10,] 204 182 135 99 60 45 31 11 53 111 NA [11,] 205 182 139 105 70 58 44 26 73 131 NA [12,] 207 185 148 106 82 73 57 41 94 155 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1 3 9 14 8 1 -9 -8 3 5 NA [2,] -1 -7 -15 -13 -8 -5 -5 -7 6 9 NA [3,] -2 0 -5 -9 -11 -15 -9 -5 4 0 NA [4,] 20 19 16 12 6 8 11 8 14 4 NA [5,] -1 -2 -6 -4 7 12 15 17 16 13 NA [6,] -15 -15 -20 -25 -25 -27 -32 -35 -29 -31 NA [7,] -10 -18 -12 -11 -3 4 4 1 9 11 NA [8,] 1 -6 -9 -2 -8 2 4 5 11 19 NA [9,] 0 3 -6 -5 -2 -7 -6 -7 -2 8 NA [10,] 1 0 4 6 10 13 13 15 20 20 NA [11,] 2 3 9 1 12 15 13 15 21 24 NA [12,] -2 -2 -6 -10 -10 -15 -15 -20 -21 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/14nfi1278856756.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/2exe31278856756.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/3exe31278856756.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/476vo1278856756.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,] 21.0 24 27.0 22 30 47 12 13 18.0 11.0 26 41 [2,] 58.0 49 44.0 35 48 64 35 44 52.0 45.0 58 73 [3,] 84.5 91 91.5 86 91 101 73 77 82.5 79.5 89 100 [4,] 183.0 192 177.0 172 188 182 162 150 141.0 135.0 139 155 [5,] 213.0 212 211.0 209 229 228 213 203 204.0 204.0 205 207 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 22.04502 19.5515 25.04790 17.54934 21.05042 42.04250 9.545736 [2,] 146.95498 162.4485 157.95210 154.45066 160.94958 159.95750 136.454264 [,8] [,9] [,10] [,11] [,12] [1,] 24.03817 38.03205 34.53241 48.52917 59.02953 [2,] 129.96183 126.96795 124.46759 129.47083 140.97047 $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(21, 58, 84.5, 183, 213, 24, 49, 91, 192, 212, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/576vo1278856756.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,] -9 -15 -15 4.0 -6.0 -35 -18 -9.0 -7.0 0.0 1.0 -21 [2,] -1 -8 -9 8.0 -2.0 -31 -11 -6.0 -6.0 4.0 3.0 -15 [3,] 3 -6 -5 11.5 9.5 -26 -1 1.5 -3.5 11.5 12.5 -10 [4,] 8 -1 0 16.0 15.0 -20 4 5.0 0.0 15.0 15.0 -6 [5,] 14 9 4 20.0 17.0 -15 11 19.0 8.0 20.0 24.0 -2 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.496759 -9.497479 -9.4967588 7.502881 1.006122 -31.49604 -8.494598 [2,] 7.496759 -2.502521 -0.5032412 15.497119 17.993878 -20.50396 6.494598 [,8] [,9] [,10] [,11] [,12] [1,] -3.996039 -6.4978392 6.003961 6.504322 -14.74 [2,] 6.996039 -0.5021608 16.996039 18.495678 -5.26 $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(-9, -1, 3, 8, 14, -15, -8, -6, -1, 9, -15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/60xcr1278856756.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] [,11] [1,] 203.0 179.0 135.0 99.0 60.0 45.0 29 11.0 21.0 73.0 NA [2,] 204.5 183.5 144.5 105.5 70.0 51.0 34 15.5 32.0 81.0 NA [3,] 210.0 201.0 167.0 125.5 83.5 59.5 44 26.5 46.0 89.0 NA [4,] 213.0 206.5 182.5 142.5 96.0 72.5 53 37.5 59.5 107.5 NA [5,] 213.0 220.0 192.0 156.0 104.0 73.0 61 47.0 94.0 131.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 206.1231 190.5095 149.6679 108.6241 71.64123 49.69371 35.33397 16.46565 [2,] 213.8769 211.4905 184.3321 142.3759 95.35877 69.30629 52.66603 36.53435 [,9] [,10] [,11] [1,] 33.45707 76.91317 NA [2,] 58.54293 101.08683 NA $out [1] 229 228 155 $group [1] 1 1 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(203, 204.5, 210, 213, 213, 179, 183.5, 201, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/70xcr1278856756.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,] 93.10 73.00 88.000 [2,] 95.90 81.00 93.250 [3,] 105.85 87.50 105.125 [4,] 113.90 91.25 115.750 [5,] 121.70 101.00 119.125 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 97.64008 82.8249 94.8626 [2,] 114.05992 92.1751 115.3874 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(93.1, 95.9, 105.85, 113.9, 121.7, 73, 81, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/14nfi1278856756.ps tmp/14nfi1278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/2exe31278856756.ps tmp/2exe31278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/3exe31278856756.ps tmp/3exe31278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/476vo1278856756.ps tmp/476vo1278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/576vo1278856756.ps tmp/576vo1278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/60xcr1278856756.ps tmp/60xcr1278856756.png",intern=TRUE)) character(0) > try(system("convert tmp/70xcr1278856756.ps tmp/70xcr1278856756.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.012 1.540 2.240