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. 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(493,481,462,457,442,439,488,521,501,485,464,460,467,460,448,443,436,431,484,510,513,503,471,471,476,475,470,461,455,456,517,525,523,519,509,512,519,517,510,509,501,507,569,580,578,565,547,555,562,561,555,544,537,543,594,611,613,611,594,595,591,589,584,573,567,569,621,629,628,612,595,597,593,590,580,574,573,573,620,626,620,588,566,557,561,549,532,526,511,499,555,565,542,527,510,514,517,508,493,490,469,478) > 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] 102 > (np <- floor(n / par1)) [1] 8 > 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] 9 9 9 9 9 9 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 493 467 476 519 562 591 593 561 517 [2,] 481 460 475 517 561 589 590 549 508 [3,] 462 448 470 510 555 584 580 532 493 [4,] 457 443 461 509 544 573 574 526 490 [5,] 442 436 455 501 537 567 573 511 469 [6,] 439 431 456 507 543 569 573 499 478 [7,] 488 484 517 569 594 621 620 555 NA [8,] 521 510 525 580 611 629 626 565 NA [9,] 501 513 523 578 613 628 620 542 NA [10,] 485 503 519 565 611 612 588 527 NA [11,] 464 471 509 547 594 595 566 510 NA [12,] 460 471 512 555 595 597 557 514 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -12 -7 -1 -2 -1 -2 -3 -12 -9 [2,] -19 -12 -5 -7 -6 -5 -10 -17 -15 [3,] -5 -5 -9 -1 -11 -11 -6 -6 -3 [4,] -15 -7 -6 -8 -7 -6 -1 -15 -21 [5,] -3 -5 1 6 6 2 0 -12 9 [6,] 49 53 61 62 51 52 47 56 NA [7,] 33 26 8 11 17 8 6 10 NA [8,] -20 3 -2 -2 2 -1 -6 -23 NA [9,] -16 -10 -4 -13 -2 -16 -32 -15 NA [10,] -21 -32 -10 -18 -17 -17 -22 -17 NA [11,] -4 0 3 8 1 2 -9 4 NA [12,] 7 5 7 7 -4 -4 4 3 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/1q0sk1229259410.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/2dq951229259410.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/3p85k1229259410.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/4rjli1229259410.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,] 467 460 448 443 436 431 484.0 510.0 501.0 485.0 464.0 460.0 [2,] 493 481 470 461 455 456 502.5 523.0 518.0 511.0 490.0 491.5 [3,] 519 517 510 509 501 499 562.0 572.5 560.0 546.0 528.5 534.5 [4,] 562 561 555 544 537 543 607.0 618.5 616.5 599.5 580.0 576.0 [5,] 593 590 584 574 573 573 621.0 629.0 628.0 612.0 595.0 597.0 $n [1] 9 9 9 9 9 9 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 482.66 474.8667 465.2333 465.2867 457.8133 453.18 503.6248 519.1523 [2,] 555.34 559.1333 554.7667 552.7133 544.1867 544.82 620.3752 625.8477 [,9] [,10] [,11] [,12] [1,] 504.9765 496.5626 478.2247 487.2971 [2,] 615.0235 595.4374 578.7753 581.7029 $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(467, 493, 519, 562, 593, 460, 481, 517, 561, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5taw91229259410.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,] -12 -19 -11 -21 -12 47.0 6.0 -23.0 -16 -22.0 -9.0 -4.0 [2,] -9 -15 -9 -15 -3 50.0 8.0 -13.0 -16 -21.5 -2.0 -0.5 [3,] -3 -10 -6 -7 1 52.5 10.5 -2.0 -14 -17.5 1.5 4.5 [4,] -2 -6 -5 -6 6 58.5 21.5 0.5 -7 -17.0 3.5 7.0 [5,] -1 -5 -1 -1 9 62.0 33.0 3.0 -2 -17.0 8.0 7.0 $n [1] 9 9 9 9 9 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] -6.6866667 -14.74 -8.106667 -11.74 -3.74 47.75178 2.958706 -9.541294 [2,] 0.6866667 -5.26 -3.893333 -2.26 5.74 57.24822 18.041294 5.541294 [,9] [,10] [,11] [,12] [1,] -19.02753 -20.01376 -1.572379 0.3103923 [2,] -8.97247 -14.98624 4.572379 8.6896077 $out [1] -32 -32 -10 $group [1] 9 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-12, -9, -3, -2, -1, -19, -15, -10, -6, -5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/60ow91229259410.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] [1,] 439.0 431.0 455.0 501.0 537.0 567.0 557.0 499.0 469.0 [2,] 458.5 445.5 465.5 509.5 549.5 578.5 573.0 512.5 478.0 [3,] 472.5 469.0 492.5 533.0 578.0 593.0 584.0 529.5 491.5 [4,] 490.5 493.5 518.0 567.0 603.0 616.5 606.5 552.0 508.0 [5,] 521.0 513.0 525.0 580.0 613.0 629.0 626.0 565.0 517.0 $n [1] 12 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 457.9046 447.1069 468.5544 506.7739 553.5983 575.668 568.7204 511.4838 [2,] 487.0954 490.8931 516.4456 559.2261 602.4017 610.332 599.2796 547.5162 [,9] [1,] 472.149 [2,] 510.851 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(439, 458.5, 472.5, 490.5, 521, 431, 445.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7whwr1229259410.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,] 499.0000 499.00 496.0000 [2,] 511.7222 509.50 507.5000 [3,] 531.5000 523.75 530.8125 [4,] 553.6250 553.00 554.7500 [5,] 570.8750 572.50 569.3750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 512.3879 503.9094 509.2615 [2,] 550.6121 543.5906 552.3635 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(499, 511.722222222222, 531.5, 553.625, 570.875, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1q0sk1229259410.ps tmp/1q0sk1229259410.png") > system("convert tmp/2dq951229259410.ps tmp/2dq951229259410.png") > system("convert tmp/3p85k1229259410.ps tmp/3p85k1229259410.png") > system("convert tmp/4rjli1229259410.ps tmp/4rjli1229259410.png") > system("convert tmp/5taw91229259410.ps tmp/5taw91229259410.png") > system("convert tmp/60ow91229259410.ps tmp/60ow91229259410.png") > system("convert tmp/7whwr1229259410.ps tmp/7whwr1229259410.png") > > > proc.time() user system elapsed 1.201 0.933 1.450