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(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/html/freestat/rcomp/tmp/13cku1281785526.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/23cku1281785526.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/3wl1x1281785526.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/4wl1x1281785526.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/html/freestat/rcomp/tmp/5wl1x1281785526.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/html/freestat/rcomp/tmp/67v111281785526.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/html/freestat/rcomp/tmp/77v111281785526.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/13cku1281785526.ps tmp/13cku1281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/23cku1281785526.ps tmp/23cku1281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/3wl1x1281785526.ps tmp/3wl1x1281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/4wl1x1281785526.ps tmp/4wl1x1281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/5wl1x1281785526.ps tmp/5wl1x1281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/67v111281785526.ps tmp/67v111281785526.png",intern=TRUE)) character(0) > try(system("convert tmp/77v111281785526.ps tmp/77v111281785526.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.921 1.444 2.173