R version 2.6.0 (2007-10-03) Copyright (C) 2007 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(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) > par1 = '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2007), Mean Plot (v1.0.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanplot.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description > par1 <- as.numeric(par1) > (n <- length(x)) [1] 73 > (np <- floor(n / par1)) [1] 6 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 7 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 513 523 578 613 628 620 542 [2,] 503 519 565 611 612 588 NA [3,] 471 509 547 594 595 566 NA [4,] 471 512 555 595 597 557 NA [5,] 476 519 562 591 593 561 NA [6,] 475 517 561 589 590 549 NA [7,] 470 510 555 584 580 532 NA [8,] 461 509 544 573 574 526 NA [9,] 455 501 537 567 573 511 NA [10,] 456 507 543 569 573 499 NA [11,] 517 569 594 621 620 555 NA [12,] 525 580 611 629 626 565 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/1rx291194894791.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/2zjd11194894791.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/35gxh1194894791.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/4eirw1194894791.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,] 513.0 503.0 471.0 471 476.0 475 470.0 461 455 456 517.0 525.0 [2,] 532.5 519.0 509.0 512 519.0 517 510.0 509 501 499 555.0 565.0 [3,] 578.0 576.5 556.5 556 561.5 555 543.5 535 524 525 581.5 595.5 [4,] 616.5 611.0 594.0 595 591.0 589 580.0 573 567 569 620.0 626.0 [5,] 628.0 612.0 595.0 597 593.0 590 584.0 574 573 573 621.0 629.0 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 527.8366 517.157 501.6723 502.4623 515.0577 508.5577 498.3477 493.7179 [2,] 628.1634 635.843 611.3277 609.5377 607.9423 601.4423 588.6523 576.2821 [,9] [,10] [,11] [,12] [1,] 481.4279 479.8477 539.5729 556.153 [2,] 566.5721 570.1523 623.4271 634.847 $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(513, 532.5, 578, 616.5, 628, 503, 519, 576.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5pgev1194894791.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] [1,] 455.0 501.0 537.0 567.0 573 499.0 542 [2,] 465.5 509.0 545.5 578.5 577 529.0 542 [3,] 473.0 514.5 558.0 592.5 594 556.0 542 [4,] 508.0 521.0 571.5 612.0 616 565.5 542 [5,] 525.0 523.0 594.0 629.0 628 620.0 542 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 453.6155 509.0267 546.1412 577.2204 576.2118 539.3521 542 [2,] 492.3845 519.9733 569.8588 607.7796 611.7882 572.6479 542 $out [1] 569 580 611 $group [1] 2 2 3 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(455, 465.5, 473, 508, 525, 501, 509, 514.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6449t1194894791.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,] 524.0000 524.00 531.5000 [2,] 534.8333 539.25 542.0625 [3,] 547.4167 556.25 553.8125 [4,] 570.0952 577.25 571.1875 [5,] 589.3333 595.50 595.5000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 531.3335 538.9179 540.5284 [2,] 563.4999 573.5821 567.0966 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(524, 534.833333333333, 547.416666666667, 570.095238095238, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1rx291194894791.ps tmp/1rx291194894791.png") > system("convert tmp/2zjd11194894791.ps tmp/2zjd11194894791.png") > system("convert tmp/35gxh1194894791.ps tmp/35gxh1194894791.png") > system("convert tmp/4eirw1194894791.ps tmp/4eirw1194894791.png") > system("convert tmp/5pgev1194894791.ps tmp/5pgev1194894791.png") > system("convert tmp/6449t1194894791.ps tmp/6449t1194894791.png") > > > proc.time() user system elapsed 2.414 1.357 2.639