R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(442,441,440,438,458,457,442,432,433,433,434,436,439,439,441,436,460,453,435,421,412,408,402,409,410,410,416,410,437,431,411,398,394,395,389,404,397,401,402,383,406,400,377,372,362,365,361,372,355,365,367,341,370,366,333,320,298,306,293,313,293,304,304,286,320,313,283,272,251,262,247,268,251,257,261,242,274,272,243,234,217,231,209,226,208,214,222,194,230,226,197,188,175,190,165,176,159,169,170,141,170,164,132,123,113,125,101,99,87,90,89,66,102,97,65,54,33,49,30,34) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 442 439 410 397 355 293 251 208 159 87 NA [2,] 441 439 410 401 365 304 257 214 169 90 NA [3,] 440 441 416 402 367 304 261 222 170 89 NA [4,] 438 436 410 383 341 286 242 194 141 66 NA [5,] 458 460 437 406 370 320 274 230 170 102 NA [6,] 457 453 431 400 366 313 272 226 164 97 NA [7,] 442 435 411 377 333 283 243 197 132 65 NA [8,] 432 421 398 372 320 272 234 188 123 54 NA [9,] 433 412 394 362 298 251 217 175 113 33 NA [10,] 433 408 395 365 306 262 231 190 125 49 NA [11,] 434 402 389 361 293 247 209 165 101 30 NA [12,] 436 409 404 372 313 268 226 176 99 34 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1 0 0 4 10 11 6 6 10 3 NA [2,] -1 2 6 1 2 0 4 8 1 -1 NA [3,] -2 -5 -6 -19 -26 -18 -19 -28 -29 -23 NA [4,] 20 24 27 23 29 34 32 36 29 36 NA [5,] -1 -7 -6 -6 -4 -7 -2 -4 -6 -5 NA [6,] -15 -18 -20 -23 -33 -30 -29 -29 -32 -32 NA [7,] -10 -14 -13 -5 -13 -11 -9 -9 -9 -11 NA [8,] 1 -9 -4 -10 -22 -21 -17 -13 -10 -21 NA [9,] 0 -4 1 3 8 11 14 15 12 16 NA [10,] 1 -6 -6 -4 -13 -15 -22 -25 -24 -19 NA [11,] 2 7 15 11 20 21 17 11 -2 4 NA [12,] 3 1 -7 -17 -20 -17 -18 -17 -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/rcomp/tmp/1t0ov1280748511.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/rcomp/tmp/2t0ov1280748511.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/rcomp/tmp/3t0ov1280748511.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/rcomp/tmp/4395g1280748511.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,] 87 90.0 89.0 66.0 102 97.0 65 54 33.0 49 30 34.0 [2,] 208 214.0 222.0 194.0 230 226.0 197 188 175.0 190 165 176.0 [3,] 324 334.5 335.5 313.5 345 339.5 308 296 274.5 284 270 290.5 [4,] 410 410.0 416.0 410.0 437 431.0 411 398 394.0 395 389 404.0 [5,] 442 441.0 441.0 438.0 460 457.0 442 432 433.0 433 434 436.0 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 223.0727 236.5706 238.5699 205.5778 241.5745 237.0738 201.0771 191.0756 [2,] 424.9273 432.4294 432.4301 421.4222 448.4255 441.9262 414.9229 400.9244 [,9] [,10] [,11] [,12] [1,] 165.0789 181.5738 158.0807 176.5821 [2,] 383.9211 386.4262 381.9193 404.4179 $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(87, 208, 324, 410, 442, 90, 214, 334.5, 410, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5395g1280748511.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,] -1 -1.0 -29 20 -7.0 -33 -14.0 -22.0 -4.0 -25 -2 -20 [2,] 0 0.0 -26 24 -6.0 -32 -13.0 -21.0 1.0 -22 4 -17 [3,] 5 1.5 -19 29 -5.5 -29 -10.5 -11.5 9.5 -14 11 -17 [4,] 10 4.0 -6 34 -4.0 -20 -9.0 -9.0 14.0 -6 17 -7 [5,] 11 8.0 -2 36 -1.0 -15 -5.0 1.0 16.0 1 21 3 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.003601297 -0.4985595 -28.992797 24.0036 -6.49928 -34.99568 -12.49856 [2,] 9.996398703 3.4985595 -9.007203 33.9964 -4.50072 -23.00432 -8.50144 [,8] [,9] [,10] [,11] [,12] [1,] -17.495678 3.004682 -21.994238 4.504682 -22.26667 [2,] -5.504322 15.995318 -6.005762 17.495318 -11.73333 $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(-1, 0, 5, 10, 11, -1, 0, 1.5, 4, 8, -29, -26, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6395g1280748511.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] [1,] 432.0 402.0 389.0 361.0 293.0 247.0 209.0 165.0 99.0 30.0 NA [2,] 433.5 410.5 396.5 368.5 309.5 265.0 228.5 182.0 118.0 41.5 NA [3,] 439.0 435.5 410.0 380.0 337.0 284.5 242.5 195.5 136.5 65.5 NA [4,] 442.0 440.0 413.5 400.5 365.5 304.0 259.0 218.0 166.5 89.5 NA [5,] 442.0 460.0 437.0 406.0 370.0 320.0 274.0 230.0 170.0 102.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 435.1231 422.0449 402.2462 365.4046 311.458 266.7118 228.5887 179.0802 [2,] 442.8769 448.9551 417.7538 394.5954 362.542 302.2882 256.4113 211.9198 [,9] [,10] [,11] [1,] 114.3788 43.60688 NA [2,] 158.6212 87.39312 NA $out [1] 458 457 $group [1] 1 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(432, 433.5, 439, 442, 442, 402, 410.5, 435.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7e05j1280748511.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,] 263.10 270.00 279.0000 [2,] 275.05 287.25 293.0625 [3,] 292.75 310.75 305.0625 [4,] 310.10 335.00 319.1875 [5,] 322.70 345.00 335.1250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 276.7635 288.9709 293.1467 [2,] 308.7365 332.5291 316.9783 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > try(system("convert tmp/1t0ov1280748511.ps tmp/1t0ov1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/2t0ov1280748511.ps tmp/2t0ov1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/3t0ov1280748511.ps tmp/3t0ov1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/4395g1280748511.ps tmp/4395g1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/5395g1280748511.ps tmp/5395g1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/6395g1280748511.ps tmp/6395g1280748511.png",intern=TRUE)) character(0) > try(system("convert tmp/7e05j1280748511.ps tmp/7e05j1280748511.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.186 0.935 10.571