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(158,157,156,154,152,151,152,154,155,155,156,158,156,152,145,141,140,145,143,141,144,139,141,142,141,132,122,122,127,128,122,123,128,128,128,129,124,121,109,110,107,107,104,110,114,118,117,122,113,106,102,111,106,110,105,104,106,110,107,111,101,105,108,124,122,128,124,121,125,134,126,126,111,117,118,128,127,129,124,113,120,127,114,107) > 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,] 158 156 141 124 113 101 111 NA [2,] 157 152 132 121 106 105 117 NA [3,] 156 145 122 109 102 108 118 NA [4,] 154 141 122 110 111 124 128 NA [5,] 152 140 127 107 106 122 127 NA [6,] 151 145 128 107 110 128 129 NA [7,] 152 143 122 104 105 124 124 NA [8,] 154 141 123 110 104 121 113 NA [9,] 155 144 128 114 106 125 120 NA [10,] 155 139 128 118 110 134 127 NA [11,] 156 141 128 117 107 126 114 NA [12,] 158 142 129 122 111 126 107 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 -4 -9 -3 -7 4 6 NA [2,] -1 -7 -10 -12 -4 3 1 NA [3,] -2 -4 0 1 9 16 10 NA [4,] -2 -1 5 -3 -5 -2 -1 NA [5,] -1 5 1 0 4 6 2 NA [6,] 1 -2 -6 -3 -5 -4 -5 NA [7,] 2 -2 1 6 -1 -3 -11 NA [8,] 1 3 5 4 2 4 7 NA [9,] 0 -5 0 4 4 9 7 NA [10,] 1 2 0 -1 -3 -8 -13 NA [11,] 2 1 1 5 4 0 -7 NA [12,] -2 -1 -5 -9 -10 -15 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/1a4yc1280407705.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/2a4yc1280407705.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/32eff1280407705.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/42eff1280407705.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,] 101.0 105.0 102.0 110.0 106.0 107 104.0 104.0 106 110.0 107.0 107.0 [2,] 112.0 111.5 108.5 116.5 114.5 119 113.5 111.5 117 122.5 115.5 116.5 [3,] 124.0 121.0 118.0 124.0 127.0 128 124.0 121.0 125 128.0 126.0 126.0 [4,] 148.5 142.0 133.5 134.5 133.5 137 133.5 132.0 136 136.5 134.5 135.5 [5,] 158.0 157.0 156.0 154.0 152.0 151 152.0 154.0 155 155.0 156.0 158.0 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 102.2028 102.7859 103.0704 113.2507 115.6535 117.2507 112.0563 108.7577 [2,] 145.7972 139.2141 132.9296 134.7493 138.3465 138.7493 135.9437 133.2423 [,9] [,10] [,11] [,12] [1,] 113.6535 119.6394 114.6535 114.6535 [2,] 136.3465 136.3606 137.3465 137.3465 $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(101, 112, 124, 148.5, 158, 105, 111.5, 121, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/52eff1280407705.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.0 -12.0 -4.0 -3.0 -1.0 -6.0 -3.0 1.0 -5.0 -13.0 0.0 -15 [2,] -5.5 -8.5 -1.0 -2.5 0.5 -5.0 -2.5 2.5 0.0 -5.5 0.5 -10 [3,] -3.0 -4.0 1.0 -2.0 2.0 -4.0 -1.0 4.0 4.0 -1.0 1.0 -7 [4,] 1.5 0.0 9.5 -1.0 4.5 -2.5 1.5 4.5 5.5 0.5 3.0 -2 [5,] 6.0 3.0 16.0 -1.0 6.0 1.0 6.0 7.0 9.0 2.0 5.0 -1 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -7.180287 -9.076063 -5.270431 -2.895776 -0.3887355 -5.492960 -3.388735 [2,] 1.180287 1.076063 7.270431 -1.104224 4.3887355 -2.507040 1.388735 [,8] [,9] [,10] [,11] [,12] [1,] 2.805632 0.7154887 -4.583103 -0.4929597 -12.160258 [2,] 5.194368 7.2845113 2.583103 2.4929597 -1.839742 $out [1] 5 -5 -11 -7 $group [1] 4 4 7 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-9, -5.5, -3, 1.5, 6, -12, -8.5, -4, 0, 3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6v5e01280407705.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,] 151.0 139.0 122.0 104.0 102.0 101.0 107.0 NA [2,] 153.0 141.0 122.5 108.0 105.5 114.5 113.5 NA [3,] 155.0 142.5 128.0 112.0 106.5 124.0 119.0 NA [4,] 156.5 145.0 128.5 119.5 110.5 126.0 127.0 NA [5,] 158.0 145.0 132.0 124.0 113.0 134.0 129.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 153.4036 140.6756 125.2634 106.7548 104.2195 118.7548 112.8426 NA [2,] 156.5964 144.3244 130.7366 117.2452 108.7805 129.2452 125.1574 NA $out [1] 156 152 141 $group [1] 2 2 3 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(151, 153, 155, 156.5, 158, 139, 141, 142.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7v5e01280407705.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,] 122.8571 118.0 121.000 [2,] 125.3571 122.5 123.750 [3,] 127.1429 124.5 125.750 [4,] 128.0714 126.5 127.375 [5,] 130.1429 128.0 130.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 125.9049 122.6756 124.0966 [2,] 128.3809 126.3244 127.4034 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(122.857142857143, 125.357142857143, 127.142857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1a4yc1280407705.ps tmp/1a4yc1280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/2a4yc1280407705.ps tmp/2a4yc1280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/32eff1280407705.ps tmp/32eff1280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/42eff1280407705.ps tmp/42eff1280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/52eff1280407705.ps tmp/52eff1280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/6v5e01280407705.ps tmp/6v5e01280407705.png",intern=TRUE)) character(0) > try(system("convert tmp/7v5e01280407705.ps tmp/7v5e01280407705.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.995 1.532 2.233