R version 2.8.1 (2008-12-22) 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(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,528,534,518,506,502,516,528) > par1 = '93' > #'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] 97 > (np <- floor(n / par1)) [1] 1 > 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] 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [77] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 > arr [,1] [,2] [1,] 467 506 [2,] 460 502 [3,] 448 516 [4,] 443 528 [5,] 436 NA [6,] 431 NA [7,] 484 NA [8,] 510 NA [9,] 513 NA [10,] 503 NA [11,] 471 NA [12,] 471 NA [13,] 476 NA [14,] 475 NA [15,] 470 NA [16,] 461 NA [17,] 455 NA [18,] 456 NA [19,] 517 NA [20,] 525 NA [21,] 523 NA [22,] 519 NA [23,] 509 NA [24,] 512 NA [25,] 519 NA [26,] 517 NA [27,] 510 NA [28,] 509 NA [29,] 501 NA [30,] 507 NA [31,] 569 NA [32,] 580 NA [33,] 578 NA [34,] 565 NA [35,] 547 NA [36,] 555 NA [37,] 562 NA [38,] 561 NA [39,] 555 NA [40,] 544 NA [41,] 537 NA [42,] 543 NA [43,] 594 NA [44,] 611 NA [45,] 613 NA [46,] 611 NA [47,] 594 NA [48,] 595 NA [49,] 591 NA [50,] 589 NA [51,] 584 NA [52,] 573 NA [53,] 567 NA [54,] 569 NA [55,] 621 NA [56,] 629 NA [57,] 628 NA [58,] 612 NA [59,] 595 NA [60,] 597 NA [61,] 593 NA [62,] 590 NA [63,] 580 NA [64,] 574 NA [65,] 573 NA [66,] 573 NA [67,] 620 NA [68,] 626 NA [69,] 620 NA [70,] 588 NA [71,] 566 NA [72,] 557 NA [73,] 561 NA [74,] 549 NA [75,] 532 NA [76,] 526 NA [77,] 511 NA [78,] 499 NA [79,] 555 NA [80,] 565 NA [81,] 542 NA [82,] 527 NA [83,] 510 NA [84,] 514 NA [85,] 517 NA [86,] 508 NA [87,] 493 NA [88,] 490 NA [89,] 469 NA [90,] 478 NA [91,] 528 NA [92,] 534 NA [93,] 518 NA > darr [,1] [,2] [1,] -7 -4 [2,] -12 14 [3,] -5 12 [4,] -7 NA [5,] -5 NA [6,] 53 NA [7,] 26 NA [8,] 3 NA [9,] -10 NA [10,] -32 NA [11,] 0 NA [12,] 5 NA [13,] -1 NA [14,] -5 NA [15,] -9 NA [16,] -6 NA [17,] 1 NA [18,] 61 NA [19,] 8 NA [20,] -2 NA [21,] -4 NA [22,] -10 NA [23,] 3 NA [24,] 7 NA [25,] -2 NA [26,] -7 NA [27,] -1 NA [28,] -8 NA [29,] 6 NA [30,] 62 NA [31,] 11 NA [32,] -2 NA [33,] -13 NA [34,] -18 NA [35,] 8 NA [36,] 7 NA [37,] -1 NA [38,] -6 NA [39,] -11 NA [40,] -7 NA [41,] 6 NA [42,] 51 NA [43,] 17 NA [44,] 2 NA [45,] -2 NA [46,] -17 NA [47,] 1 NA [48,] -4 NA [49,] -2 NA [50,] -5 NA [51,] -11 NA [52,] -6 NA [53,] 2 NA [54,] 52 NA [55,] 8 NA [56,] -1 NA [57,] -16 NA [58,] -17 NA [59,] 2 NA [60,] -4 NA [61,] -3 NA [62,] -10 NA [63,] -6 NA [64,] -1 NA [65,] 0 NA [66,] 47 NA [67,] 6 NA [68,] -6 NA [69,] -32 NA [70,] -22 NA [71,] -9 NA [72,] 4 NA [73,] -12 NA [74,] -17 NA [75,] -6 NA [76,] -15 NA [77,] -12 NA [78,] 56 NA [79,] 10 NA [80,] -23 NA [81,] -15 NA [82,] -17 NA [83,] 4 NA [84,] 3 NA [85,] -9 NA [86,] -15 NA [87,] -3 NA [88,] -21 NA [89,] 9 NA [90,] 50 NA [91,] 6 NA [92,] -16 NA [93,] -12 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/1r6op1241642078.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/2yypd1241642078.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/33k6x1241642078.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/4ezdc1241642078.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] [,13] [1,] 467.0 460 448 443.0 436 431 484 510 513 503 471 471 476 [2,] 467.0 460 448 443.0 436 431 484 510 513 503 471 471 476 [3,] 486.5 481 482 485.5 436 431 484 510 513 503 471 471 476 [4,] 506.0 502 516 528.0 436 431 484 510 513 503 471 471 476 [5,] 506.0 502 516 528.0 436 431 484 510 513 503 471 471 476 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] 475 470 461 455 456 517 525 523 519 509 512 519 [2,] 475 470 461 455 456 517 525 523 519 509 512 519 [3,] 475 470 461 455 456 517 525 523 519 509 512 519 [4,] 475 470 461 455 456 517 525 523 519 509 512 519 [5,] 475 470 461 455 456 517 525 523 519 509 512 519 [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 517 510 509 501 507 569 580 578 565 547 555 562 [2,] 517 510 509 501 507 569 580 578 565 547 555 562 [3,] 517 510 509 501 507 569 580 578 565 547 555 562 [4,] 517 510 509 501 507 569 580 578 565 547 555 562 [5,] 517 510 509 501 507 569 580 578 565 547 555 562 [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [1,] 561 555 544 537 543 594 611 613 611 594 595 591 [2,] 561 555 544 537 543 594 611 613 611 594 595 591 [3,] 561 555 544 537 543 594 611 613 611 594 595 591 [4,] 561 555 544 537 543 594 611 613 611 594 595 591 [5,] 561 555 544 537 543 594 611 613 611 594 595 591 [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [1,] 589 584 573 567 569 621 629 628 612 595 597 593 [2,] 589 584 573 567 569 621 629 628 612 595 597 593 [3,] 589 584 573 567 569 621 629 628 612 595 597 593 [4,] 589 584 573 567 569 621 629 628 612 595 597 593 [5,] 589 584 573 567 569 621 629 628 612 595 597 593 [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [1,] 590 580 574 573 573 620 626 620 588 566 557 561 [2,] 590 580 574 573 573 620 626 620 588 566 557 561 [3,] 590 580 574 573 573 620 626 620 588 566 557 561 [4,] 590 580 574 573 573 620 626 620 588 566 557 561 [5,] 590 580 574 573 573 620 626 620 588 566 557 561 [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [1,] 549 532 526 511 499 555 565 542 527 510 514 517 [2,] 549 532 526 511 499 555 565 542 527 510 514 517 [3,] 549 532 526 511 499 555 565 542 527 510 514 517 [4,] 549 532 526 511 499 555 565 542 527 510 514 517 [5,] 549 532 526 511 499 555 565 542 527 510 514 517 [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [1,] 508 493 490 469 478 528 534 518 [2,] 508 493 490 469 478 528 534 518 [3,] 508 493 490 469 478 528 534 518 [4,] 508 493 490 469 478 528 534 518 [5,] 508 493 490 469 478 528 534 518 $n [1] 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [77] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 442.9281 434.0764 406.0284 390.5356 436 431 484 510 513 503 471 [2,] 530.0719 527.9236 557.9716 580.4644 436 431 484 510 513 503 471 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 471 476 475 470 461 455 456 517 525 523 519 509 [2,] 471 476 475 470 461 455 456 517 525 523 519 509 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] 512 519 517 510 509 501 507 569 580 578 565 547 [2,] 512 519 517 510 509 501 507 569 580 578 565 547 [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [1,] 555 562 561 555 544 537 543 594 611 613 611 594 [2,] 555 562 561 555 544 537 543 594 611 613 611 594 [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [1,] 595 591 589 584 573 567 569 621 629 628 612 595 [2,] 595 591 589 584 573 567 569 621 629 628 612 595 [,60] [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [1,] 597 593 590 580 574 573 573 620 626 620 588 566 [2,] 597 593 590 580 574 573 573 620 626 620 588 566 [,72] [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [1,] 557 561 549 532 526 511 499 555 565 542 527 510 [2,] 557 561 549 532 526 511 499 555 565 542 527 510 [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [1,] 514 517 508 493 490 469 478 528 534 518 [2,] 514 517 508 493 490 469 478 528 534 518 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" "73" "74" "75" [76] "76" "77" "78" "79" "80" "81" "82" "83" "84" "85" "86" "87" "88" "89" "90" [91] "91" "92" "93" Warning message: In bxp(list(stats = c(467, 467, 486.5, 506, 506, 460, 460, 481, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5z3h51241642078.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] [,13] [,14] [1,] -7.0 -12 -5.0 -7 -5 53 26 3 -10 -32 0 5 -1 -5 [2,] -7.0 -12 -5.0 -7 -5 53 26 3 -10 -32 0 5 -1 -5 [3,] -5.5 1 3.5 -7 -5 53 26 3 -10 -32 0 5 -1 -5 [4,] -4.0 14 12.0 -7 -5 53 26 3 -10 -32 0 5 -1 -5 [5,] -4.0 14 12.0 -7 -5 53 26 3 -10 -32 0 5 -1 -5 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] -9 -6 1 61 8 -2 -4 -10 3 7 -2 -7 [2,] -9 -6 1 61 8 -2 -4 -10 3 7 -2 -7 [3,] -9 -6 1 61 8 -2 -4 -10 3 7 -2 -7 [4,] -9 -6 1 61 8 -2 -4 -10 3 7 -2 -7 [5,] -9 -6 1 61 8 -2 -4 -10 3 7 -2 -7 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] -1 -8 6 62 11 -2 -13 -18 8 7 -1 -6 [2,] -1 -8 6 62 11 -2 -13 -18 8 7 -1 -6 [3,] -1 -8 6 62 11 -2 -13 -18 8 7 -1 -6 [4,] -1 -8 6 62 11 -2 -13 -18 8 7 -1 -6 [5,] -1 -8 6 62 11 -2 -13 -18 8 7 -1 -6 [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [1,] -11 -7 6 51 17 2 -2 -17 1 -4 -2 -5 [2,] -11 -7 6 51 17 2 -2 -17 1 -4 -2 -5 [3,] -11 -7 6 51 17 2 -2 -17 1 -4 -2 -5 [4,] -11 -7 6 51 17 2 -2 -17 1 -4 -2 -5 [5,] -11 -7 6 51 17 2 -2 -17 1 -4 -2 -5 [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [,61] [,62] [1,] -11 -6 2 52 8 -1 -16 -17 2 -4 -3 -10 [2,] -11 -6 2 52 8 -1 -16 -17 2 -4 -3 -10 [3,] -11 -6 2 52 8 -1 -16 -17 2 -4 -3 -10 [4,] -11 -6 2 52 8 -1 -16 -17 2 -4 -3 -10 [5,] -11 -6 2 52 8 -1 -16 -17 2 -4 -3 -10 [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [,72] [,73] [,74] [1,] -6 -1 0 47 6 -6 -32 -22 -9 4 -12 -17 [2,] -6 -1 0 47 6 -6 -32 -22 -9 4 -12 -17 [3,] -6 -1 0 47 6 -6 -32 -22 -9 4 -12 -17 [4,] -6 -1 0 47 6 -6 -32 -22 -9 4 -12 -17 [5,] -6 -1 0 47 6 -6 -32 -22 -9 4 -12 -17 [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [,84] [,85] [,86] [1,] -6 -15 -12 56 10 -23 -15 -17 4 3 -9 -15 [2,] -6 -15 -12 56 10 -23 -15 -17 4 3 -9 -15 [3,] -6 -15 -12 56 10 -23 -15 -17 4 3 -9 -15 [4,] -6 -15 -12 56 10 -23 -15 -17 4 3 -9 -15 [5,] -6 -15 -12 56 10 -23 -15 -17 4 3 -9 -15 [,87] [,88] [,89] [,90] [,91] [,92] [,93] [1,] -3 -21 9 50 6 -16 -12 [2,] -3 -21 9 50 6 -16 -12 [3,] -3 -21 9 50 6 -16 -12 [4,] -3 -21 9 50 6 -16 -12 [5,] -3 -21 9 50 6 -16 -12 $n [1] 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [77] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -8.851686 -28.04795 -15.49289 -7 -5 53 26 3 -10 -32 0 [2,] -2.148314 30.04795 22.49289 -7 -5 53 26 3 -10 -32 0 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 5 -1 -5 -9 -6 1 61 8 -2 -4 -10 3 [2,] 5 -1 -5 -9 -6 1 61 8 -2 -4 -10 3 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] 7 -2 -7 -1 -8 6 62 11 -2 -13 -18 8 [2,] 7 -2 -7 -1 -8 6 62 11 -2 -13 -18 8 [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [1,] 7 -1 -6 -11 -7 6 51 17 2 -2 -17 1 [2,] 7 -1 -6 -11 -7 6 51 17 2 -2 -17 1 [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [1,] -4 -2 -5 -11 -6 2 52 8 -1 -16 -17 2 [2,] -4 -2 -5 -11 -6 2 52 8 -1 -16 -17 2 [,60] [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [1,] -4 -3 -10 -6 -1 0 47 6 -6 -32 -22 -9 [2,] -4 -3 -10 -6 -1 0 47 6 -6 -32 -22 -9 [,72] [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [,82] [,83] [1,] 4 -12 -17 -6 -15 -12 56 10 -23 -15 -17 4 [2,] 4 -12 -17 -6 -15 -12 56 10 -23 -15 -17 4 [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [,93] [1,] 3 -9 -15 -3 -21 9 50 6 -16 -12 [2,] 3 -9 -15 -3 -21 9 50 6 -16 -12 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" [61] "61" "62" "63" "64" "65" "66" "67" "68" "69" "70" "71" "72" "73" "74" "75" [76] "76" "77" "78" "79" "80" "81" "82" "83" "84" "85" "86" "87" "88" "89" "90" [91] "91" "92" "93" Warning message: In bxp(list(stats = c(-7, -7, -5.5, -4, -4, -12, -12, 1, 14, 14, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6m5ak1241642078.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] [1,] 431 502 [2,] 508 504 [3,] 537 511 [4,] 578 522 [5,] 629 528 $n [1] 93 4 $conf [,1] [,2] [1,] 525.5313 496.78 [2,] 548.4687 525.22 $out numeric(0) $group numeric(0) $names [1] "1" NA Warning message: In bxp(list(stats = c(431, 508, 537, 578, 629, 502, 504, 511, 522, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/74uxe1241642078.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,] 431 431 431 [2,] 508 508 508 [3,] 537 537 537 [4,] 578 578 578 [5,] 629 629 629 $n [1] 93 93 93 $conf [,1] [,2] [,3] [1,] 525.5313 525.5313 525.5313 [2,] 548.4687 548.4687 548.4687 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/1r6op1241642078.ps tmp/1r6op1241642078.png") > system("convert tmp/2yypd1241642078.ps tmp/2yypd1241642078.png") > system("convert tmp/33k6x1241642078.ps tmp/33k6x1241642078.png") > system("convert tmp/4ezdc1241642078.ps tmp/4ezdc1241642078.png") > system("convert tmp/5z3h51241642078.ps tmp/5z3h51241642078.png") > system("convert tmp/6m5ak1241642078.ps tmp/6m5ak1241642078.png") > system("convert tmp/74uxe1241642078.ps tmp/74uxe1241642078.png") > > > proc.time() user system elapsed 1.497 1.022 2.307