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(335,334,333,331,329,328,329,331,332,332,333,335,335,333,325,322,322,315,321,324,329,332,322,324,324,323,309,306,305,300,301,302,308,311,301,301,308,302,290,286,286,275,284,289,292,293,285,280,281,280,265,260,254,238,247,246,247,237,222,216,212,209,185,186,178,158,166,162,164,147,132,124,117,120,89,81,71,52,63,62,74,67,53,42) > 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,] 335 335 324 308 281 212 117 NA [2,] 334 333 323 302 280 209 120 NA [3,] 333 325 309 290 265 185 89 NA [4,] 331 322 306 286 260 186 81 NA [5,] 329 322 305 286 254 178 71 NA [6,] 328 315 300 275 238 158 52 NA [7,] 329 321 301 284 247 166 63 NA [8,] 331 324 302 289 246 162 62 NA [9,] 332 329 308 292 247 164 74 NA [10,] 332 332 311 293 237 147 67 NA [11,] 333 322 301 285 222 132 53 NA [12,] 335 324 301 280 216 124 42 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -1 -2 -1 -6 -1 -3 3 NA [2,] -1 -8 -14 -12 -15 -24 -31 NA [3,] -2 -3 -3 -4 -5 1 -8 NA [4,] -2 0 -1 0 -6 -8 -10 NA [5,] -1 -7 -5 -11 -16 -20 -19 NA [6,] 1 6 1 9 9 8 11 NA [7,] 2 3 1 5 -1 -4 -1 NA [8,] 1 5 6 3 1 2 12 NA [9,] 0 3 3 1 -10 -17 -7 NA [10,] 1 -10 -10 -8 -15 -15 -14 NA [11,] 2 2 0 -5 -6 -8 -11 NA [12,] 0 0 7 1 -4 -7 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/1oru41281523830.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/2oru41281523830.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/3zitp1281523830.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/4zitp1281523830.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,] 212.0 120.0 89 186 71.0 52.0 63.0 62 74.0 67.0 53.0 42.0 [2,] 246.5 244.5 225 223 216.0 198.0 206.5 204 205.5 192.0 177.0 170.0 [3,] 308.0 302.0 290 286 286.0 275.0 284.0 289 292.0 293.0 285.0 280.0 [4,] 329.5 328.0 317 314 313.5 307.5 311.0 313 318.5 321.5 311.5 312.5 [5,] 335.0 334.0 333 331 329.0 328.0 329.0 331 332.0 332.0 333.0 335.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,] 258.4337 252.1351 235.0591 231.6563 227.7746 209.6084 221.5943 223.9070 [2,] 357.5663 351.8649 344.9409 340.3437 344.2254 340.3916 346.4057 354.0930 [,9] [,10] [,11] [,12] [1,] 224.5182 215.6647 204.6788 194.9013 [2,] 359.4818 370.3353 365.3212 365.0987 $out [1] 117 81 $group [1] 1 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(212, 246.5, 308, 329.5, 335, 120, 244.5, 302, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5zitp1281523830.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,] -3.0 -31.0 -5.0 -10.0 -20.0 1.0 -4.0 1.0 -17.0 -15.0 -11 -7 [2,] -2.5 -19.5 -4.5 -7.0 -17.5 3.5 -1.0 1.5 -8.5 -14.5 -7 -4 [3,] -1.0 -14.0 -3.0 -2.0 -11.0 8.0 1.0 3.0 0.0 -10.0 -5 0 [4,] -1.0 -10.0 -2.5 -0.5 -6.0 9.0 2.5 5.5 2.0 -9.0 1 1 [5,] -1.0 -1.0 -2.0 0.0 -1.0 11.0 5.0 6.0 3.0 -8.0 2 7 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.8957758 -19.673247 -4.194368 -5.881695 -17.867614 4.715489 -1.090144 [2,] -0.1042242 -8.326753 -1.805632 1.881695 -4.132386 11.284511 3.090144 [,8] [,9] [,10] [,11] [,12] [1,] 0.6112645 -6.27043 -13.284511 -9.7774709 -3.225161 [2,] 5.3887355 6.27043 -6.715489 -0.2225291 3.225161 $out [1] -6 3 1 -8 12 1 $group [1] 1 1 3 3 8 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-3, -2.5, -1, -1, -1, -31, -19.5, -14, -10, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/69ras1281523830.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,] 328.0 315.0 300.0 275.0 216.0 124.0 42.0 NA [2,] 330.0 322.0 301.0 284.5 237.5 152.5 57.5 NA [3,] 332.0 324.0 305.5 287.5 247.0 165.0 69.0 NA [4,] 333.5 330.5 310.0 292.5 262.5 185.5 85.0 NA [5,] 335.0 335.0 323.0 302.0 281.0 212.0 120.0 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 330.4036 320.1231 301.3950 283.8511 235.5973 149.9485 56.45707 NA [2,] 333.5964 327.8769 309.6050 291.1489 258.4027 180.0515 81.54293 NA $out [1] 324 308 $group [1] 3 4 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(328, 330, 332, 333.5, 335, 315, 322, 324, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/79ras1281523830.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,] 231.7143 275.0 241.250 [2,] 241.2143 284.5 254.750 [3,] 247.4286 287.5 260.375 [4,] 254.8571 292.5 269.750 [5,] 273.1429 302.0 288.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 241.2060 283.8511 253.5334 [2,] 253.6512 291.1489 267.2166 $out [1] 308 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(231.714285714286, 241.214285714286, 247.428571428571, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1oru41281523830.ps tmp/1oru41281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/2oru41281523830.ps tmp/2oru41281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/3zitp1281523830.ps tmp/3zitp1281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/4zitp1281523830.ps tmp/4zitp1281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/5zitp1281523830.ps tmp/5zitp1281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/69ras1281523830.ps tmp/69ras1281523830.png",intern=TRUE)) character(0) > try(system("convert tmp/79ras1281523830.ps tmp/79ras1281523830.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.91 0.69 2.59