R version 2.11.1 (2010-05-31) Copyright (C) 2010 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(297,296,295,293,291,290,291,293,294,294,295,297,302,297,301,298,295,287,290,288,288,287,274,282,296,292,298,296,292,296,293,295,294,291,279,284,299,296,299,299,291,298,288,284,277,270,251,257,269,271,268,268,258,261,255,251,239,229,210,218,226,227,222,215,203,205,194,190,182,179,158,163,165,169,163,154,142,146,133,131,128,120,88,95) > 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,] 297 302 296 299 269 226 165 NA [2,] 296 297 292 296 271 227 169 NA [3,] 295 301 298 299 268 222 163 NA [4,] 293 298 296 299 268 215 154 NA [5,] 291 295 292 291 258 203 142 NA [6,] 290 287 296 298 261 205 146 NA [7,] 291 290 293 288 255 194 133 NA [8,] 293 288 295 284 251 190 131 NA [9,] 294 288 294 277 239 182 128 NA [10,] 294 287 291 270 229 179 120 NA [11,] 295 274 279 251 210 158 88 NA [12,] 297 282 284 257 218 163 95 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 -5 -4 -3 2 1 4 NA [2,] -1 4 6 3 -3 -5 -6 NA [3,] -2 -3 -2 0 0 -7 -9 NA [4,] -2 -3 -4 -8 -10 -12 -12 NA [5,] -1 -8 4 7 3 2 4 NA [6,] 1 3 -3 -10 -6 -11 -13 NA [7,] 2 -2 2 -4 -4 -4 -2 NA [8,] 1 0 -1 -7 -12 -8 -3 NA [9,] 0 -1 -3 -7 -10 -3 -8 NA [10,] 1 -13 -12 -19 -19 -21 -32 NA [11,] 2 8 5 6 8 5 7 NA [12,] 5 14 15 12 8 2 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/rcomp/tmp/13dld1279960560.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/rcomp/tmp/23dld1279960560.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/rcomp/tmp/3w42y1279960560.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/rcomp/tmp/4w42y1279960560.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,] 226.0 227 222.0 215.0 142.0 146 133.0 131.0 128.0 120 88.0 95.0 [2,] 247.5 249 245.0 241.5 230.5 233 224.5 220.5 210.5 204 184.0 190.5 [3,] 296.0 292 295.0 293.0 291.0 287 288.0 284.0 277.0 270 251.0 257.0 [4,] 298.0 296 298.5 297.0 291.5 293 290.5 290.5 291.0 289 276.5 283.0 [5,] 302.0 297 301.0 299.0 295.0 298 293.0 295.0 294.0 294 295.0 297.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,] 265.8422 263.9324 263.0507 259.8563 254.5718 251.1690 248.5859 242.1971 [2,] 326.1578 320.0676 326.9493 326.1437 327.4282 322.8310 327.4141 325.8029 [,9] [,10] [,11] [,12] [1,] 228.9267 219.2394 195.7605 201.7605 [2,] 325.0733 320.7606 306.2395 312.2395 $out [1] 165 169 163 154 $group [1] 1 2 3 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(226, 247.5, 296, 298, 302, 227, 249, 292, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5w42y1279960560.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,] -5.0 -6.0 -9 -12.0 -1.0 -13.0 -4 -12.0 -10.0 -21.0 2.0 2 [2,] -3.5 -4.0 -5 -11.0 0.5 -10.5 -4 -7.5 -7.5 -20.0 5.0 5 [3,] -1.0 -1.0 -2 -8.0 3.0 -6.0 -2 -3.0 -3.0 -19.0 6.0 10 [4,] 1.5 3.5 -1 -3.5 4.0 -1.0 0 -0.5 -2.0 -12.5 7.5 14 [5,] 4.0 6.0 0 -2.0 7.0 3.0 2 1.0 0.0 -12.0 8.0 15 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3.985919 -5.478879 -4.3887355 -12.478879 0.9098565 -11.6732467 -4.3887355 [2,] 1.985919 3.478879 0.3887355 -3.521121 5.0901435 -0.3267533 0.3887355 [,8] [,9] [,10] [,11] [,12] [1,] -7.180287 -6.2845113 -23.47888 4.50704 4.194709 [2,] 1.180287 0.2845113 -14.52112 7.49296 15.805291 $out [1] -8 1 -32 $group [1] 5 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-5, -3.5, -1, 1.5, 4, -6, -4, -1, 3.5, 6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/66d1j1279960560.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,] 290.0 274.0 291.0 251.0 210.0 158.0 88.0 NA [2,] 292.0 287.0 291.5 273.5 234.0 180.5 124.0 NA [3,] 294.0 289.0 293.5 289.5 256.5 198.5 137.5 NA [4,] 295.5 297.5 296.0 298.5 268.0 218.5 158.5 NA [5,] 297.0 302.0 298.0 299.0 271.0 227.0 169.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 292.4036 284.2109 291.4475 278.0973 240.9924 181.1679 121.7643 NA [2,] 295.5964 293.7891 295.5525 300.9027 272.0076 215.8321 153.2357 NA $out [1] 279 284 $group [1] 3 3 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(290, 292, 294, 295.5, 297, 274, 287, 289, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/76d1j1279960560.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,] 222.1429 251.0 230.250 [2,] 240.8571 273.5 248.625 [3,] 251.1429 287.5 259.250 [4,] 262.0714 292.5 270.500 [5,] 264.8571 296.0 272.750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 241.4669 278.834 249.2727 [2,] 260.8188 296.166 269.2273 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(222.142857142857, 240.857142857143, 251.142857142857, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/13dld1279960560.ps tmp/13dld1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/23dld1279960560.ps tmp/23dld1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/3w42y1279960560.ps tmp/3w42y1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/4w42y1279960560.ps tmp/4w42y1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/5w42y1279960560.ps tmp/5w42y1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/66d1j1279960560.ps tmp/66d1j1279960560.png",intern=TRUE)) character(0) > try(system("convert tmp/76d1j1279960560.ps tmp/76d1j1279960560.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.430 0.120 1.559