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(112,118,132,129,121,135,148,148,136,119,104,118,115,126,141,135,125,149,170,170,158,133,114,140,145,150,178,163,172,178,199,199,184,162,146,166,171,180,193,181,183,218,230,242,209,191,172,194,196,196,236,235,229,243,264,272,237,211,180,201,204,188,235,227,234,264,302,293,259,229,203,229,242,233,267,269,270,315,364,347,312,274,237,278,284,277,317,313,318,374,413,405,355,306,271,306,315,301,356,348,355,422,465,467,404,347,305,336,340,318,362,348,363,435,491,505,404,359,310,337,360,342,406,396,420,472,548,559,463,407,362,405,417,391,419,461,472,535,622,606,508,461,390,432) > 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] 144 > (np <- floor(n / par1)) [1] 12 > 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] 12 12 12 12 12 12 12 12 12 12 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 112 115 145 171 196 204 242 284 315 340 360 417 NA [2,] 118 126 150 180 196 188 233 277 301 318 342 391 NA [3,] 132 141 178 193 236 235 267 317 356 362 406 419 NA [4,] 129 135 163 181 235 227 269 313 348 348 396 461 NA [5,] 121 125 172 183 229 234 270 318 355 363 420 472 NA [6,] 135 149 178 218 243 264 315 374 422 435 472 535 NA [7,] 148 170 199 230 264 302 364 413 465 491 548 622 NA [8,] 148 170 199 242 272 293 347 405 467 505 559 606 NA [9,] 136 158 184 209 237 259 312 355 404 404 463 508 NA [10,] 119 133 162 191 211 229 274 306 347 359 407 461 NA [11,] 104 114 146 172 180 203 237 271 305 310 362 390 NA [12,] 118 140 166 194 201 229 278 306 336 337 405 432 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 6 11 5 9 0 -16 -9 -7 -14 -22 -18 -26 NA [2,] 14 15 28 13 40 47 34 40 55 44 64 28 NA [3,] -3 -6 -15 -12 -1 -8 2 -4 -8 -14 -10 42 NA [4,] -8 -10 9 2 -6 7 1 5 7 15 24 11 NA [5,] 14 24 6 35 14 30 45 56 67 72 52 63 NA [6,] 13 21 21 12 21 38 49 39 43 56 76 87 NA [7,] 0 0 0 12 8 -9 -17 -8 2 14 11 -16 NA [8,] -12 -12 -15 -33 -35 -34 -35 -50 -63 -101 -96 -98 NA [9,] -17 -25 -22 -18 -26 -30 -38 -49 -57 -45 -56 -47 NA [10,] -15 -19 -16 -19 -31 -26 -37 -35 -42 -49 -45 -71 NA [11,] 14 26 20 22 21 26 41 35 31 27 43 42 NA [12,] -3 5 5 2 3 13 6 9 4 23 12 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/1uew31289575915.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/2uew31289575915.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/3mnvn1289575915.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/4mnvn1289575915.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,] 112.0 118.0 132.0 129 121.0 135.0 148.0 148.0 136.0 119.0 104.0 118.0 [2,] 158.0 165.0 185.5 172 177.5 198.0 214.5 220.5 196.5 176.5 159.0 180.0 [3,] 223.0 214.5 251.5 252 252.0 289.5 333.0 320.0 285.5 251.5 220.0 253.5 [4,] 327.5 309.5 359.0 348 359.0 428.5 478.0 486.0 404.0 353.0 307.5 336.5 [5,] 417.0 391.0 419.0 461 472.0 535.0 622.0 606.0 508.0 461.0 390.0 432.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 145.6899 148.5926 172.3655 171.7252 169.2166 184.3674 212.8159 198.9037 [2,] 300.3101 280.4074 330.6345 332.2748 334.7834 394.6326 453.1841 441.0963 [,9] [,10] [,11] [,12] [1,] 190.8579 170.9972 152.2682 182.1193 [2,] 380.1421 332.0028 287.7318 324.8807 $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(112, 158, 223, 327.5, 417, 118, 165, 214.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5mnvn1289575915.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,] -26.0 13.0 -15 -10.0 6.0 12.0 -17.0 -101.0 -57.0 -71.0 14.0 -3.0 [2,] -17.0 21.5 -11 -2.5 19.0 21.0 -8.5 -79.5 -48.0 -43.5 21.5 3.5 [3,] -8.0 37.0 -7 6.0 40.0 38.5 0.0 -35.0 -34.0 -33.0 26.5 5.0 [4,] 5.5 45.5 -2 10.0 59.5 52.5 9.5 -24.0 -23.5 -19.0 38.0 10.5 [5,] 11.0 64.0 2 24.0 72.0 87.0 14.0 -12.0 -17.0 -15.0 43.0 13.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -18.262401 26.05344 -11.104960 0.2986661 21.52768 24.13264 -8.209921 [2,] 2.262401 47.94656 -2.895040 11.7013339 58.47232 52.86736 8.209921 [,8] [,9] [,10] [,11] [,12] [1,] -60.313923 -45.17461 -44.17461 18.97424 1.665285 [2,] -9.686077 -22.82539 -21.82539 34.02576 8.334715 $out [1] 42 23 $group [1] 3 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-26, -17, -8, 5.5, 11, 13, 21.5, 37, 45.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6xxdr1289575915.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] [,12] [1,] 104.0 114.0 145 171.0 180.0 188.0 233.0 271.0 301.0 310.0 342.0 390.0 [2,] 118.0 125.5 156 180.5 198.5 215.5 254.5 295.0 325.5 338.5 379.0 418.0 [3,] 125.0 137.5 169 192.0 232.0 231.5 272.0 315.0 351.5 360.5 406.5 461.0 [4,] 135.5 153.5 181 213.5 240.0 261.5 313.5 364.5 413.0 419.5 467.5 521.5 [5,] 148.0 170.0 199 242.0 272.0 302.0 364.0 413.0 467.0 505.0 559.0 622.0 [,13] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 117.0181 124.729 157.5973 176.9485 213.0716 210.5191 245.0897 283.3006 [2,] 132.9819 150.271 180.4027 207.0515 250.9284 252.4809 298.9103 346.6994 [,9] [,10] [,11] [,12] [,13] [1,] 311.5907 323.5554 366.1346 413.7930 NA [2,] 391.4093 397.4446 446.8654 508.2070 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(104, 118, 125, 135.5, 148, 114, 125.5, 137.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7xxdr1289575915.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,] 232.8333 214.50 235.875 [2,] 251.7917 237.25 252.250 [3,] 268.6250 252.00 267.750 [4,] 307.0417 287.50 310.000 [5,] 351.3333 333.00 353.875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 243.4251 229.0806 241.4098 [2,] 293.8249 274.9194 294.0902 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(232.833333333333, 251.791666666667, 268.625, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1uew31289575915.ps tmp/1uew31289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/2uew31289575915.ps tmp/2uew31289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/3mnvn1289575915.ps tmp/3mnvn1289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/4mnvn1289575915.ps tmp/4mnvn1289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/5mnvn1289575915.ps tmp/5mnvn1289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/6xxdr1289575915.ps tmp/6xxdr1289575915.png",intern=TRUE)) character(0) > try(system("convert tmp/7xxdr1289575915.ps tmp/7xxdr1289575915.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.963 1.494 2.267