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(196.9 + ,192.1 + ,201.8 + ,186.9 + ,218 + ,214.4 + ,227.5 + ,204.1 + ,225.8 + ,223.7 + ,244.7 + ,243.9 + ,257.3 + ,234.5 + ,251.4 + ,243.8 + ,247.4 + ,245.3 + ,262.5 + ,270 + ,259.9 + ,262.2 + ,244.9 + ,249.3 + ,268.2 + ,231.2 + ,264.3 + ,252.7 + ,275.5 + ,261.5 + ,275.5 + ,272.3 + ,268.6 + ,270.4 + ,267.7 + ,275 + ,272.6 + ,248.6 + ,279.4 + ,270.5 + ,292.8 + ,297.8 + ,296.8 + ,290.9 + ,282.8 + ,312.8 + ,303.2 + ,301.4 + ,289.8 + ,279.6 + ,302.2 + ,299.1 + ,319.7 + ,310.9 + ,315.2 + ,338.5 + ,315.6 + ,321.2 + ,318.5 + ,342.7 + ,261.4 + ,287 + ,331.5 + ,326.9 + ,338.6 + ,337 + ,358.4 + ,344.5 + ,345.7 + ,344.1 + ,317.4 + ,354.5 + ,345.2 + ,314.1 + ,352.5 + ,361.2 + ,365.9 + ,332.5 + ,364 + ,359.1 + ,345.6 + ,366.9 + ,370.2 + ,359.9 + ,366.6 + ,336.3 + ,368.5 + ,374.2 + ,384.3 + ,358.9 + ,407.7 + ,433.3 + ,404.7 + ,392.7 + ,409.7 + ,416.5 + ,414.3 + ,404.3 + ,421.4 + ,372.6 + ,404.7 + ,420.2 + ,438.4 + ,449.1 + ,445.8 + ,413.8 + ,420.5 + ,442.3 + ,438.9 + ,394.5 + ,416.8 + ,402.9 + ,424.5 + ,432.3 + ,484.1 + ,492.7 + ,496.3 + ,471.9 + ,491.2 + ,512.9 + ,482.4 + ,407.9 + ,448.5 + ,431.1 + ,498.8 + ,497.1 + ,517.1 + ,487.7 + ,512.5 + ,550.1 + ,532.5 + ,524.1 + ,515.7 + ,461 + ,529.3 + ,467.4 + ,559.8 + ,536.5 + ,531.9 + ,546.5 + ,547.4 + ,536.1 + ,482.8 + ,551 + ,532.9 + ,484.1 + ,554.8 + ,537 + ,558 + ,511.4 + ,502.9 + ,558.6 + ,545.1 + ,574.3 + ,542.2 + ,600 + ,588.6 + ,524.4 + ,618.5 + ,580.9 + ,557.2 + ,571.2 + ,597.5 + ,601.7 + ,558.9 + ,600.9 + ,601 + ,615.7 + ,578.1 + ,495.9 + ,526.8 + ,522.1 + ,605.1 + ,574.4 + ,609.7 + ,580.7 + ,565.1 + ,590.7 + ,571.5 + ,601.3 + ,567.3 + ,467.9 + ,588.9 + ,579.4 + ,502.6 + ,568.7 + ,616 + ,586.2 + ,575.5 + ,599.9 + ,568.2 + ,516 + ,493.4 + ,496.8 + ,529.9 + ,491.7 + ,543.2 + ,490.8 + ,554.7 + ,625.7 + ,605 + ,645.2 + ,645.2 + ,611.8 + ,600.3 + ,549.8 + ,635.5 + ,617.7 + ,643.5 + ,485.7 + ,689.5 + ,692 + ,677.3 + ,704.7 + ,668.6 + ,717.8 + ,689.8 + ,640.4 + ,675.2 + ,528.1 + ,538 + ,527.2 + ,655.6 + ,650.6 + ,623.7 + ,748.4 + ,727.4 + ,750.5 + ,678.9 + ,659.5 + ,691.9 + ,639.8 + ,663.8 + ,572.9 + ,592.5 + ,734.8 + ,696.1 + ,589.2 + ,662.9 + ,661.2 + ,672.1 + ,583.7 + ,705.5 + ,631 + ,733.3 + ,674.9 + ,695.5 + ,634.1 + ,630.6 + ,635.2 + ,554.1 + ,623.9 + ,679.3 + ,565.6 + ,564.1 + ,637.2 + ,650.8 + ,602.7 + ,587.5 + ,619.2 + ,616.5 + ,637.9 + ,557.9 + ,594 + ,668.7 + ,603.3 + ,674.5 + ,573.4 + ,706 + ,693.7 + ,627.5 + ,550.7 + ,592.3 + ,660.2 + ,597.3 + ,641 + ,663.6 + ,595.9 + ,638.4 + ,665.4 + ,671.4 + ,637 + ,685.7 + ,705.8 + ,704.8 + ,734.4 + ,674.2 + ,748.6 + ,763.4 + ,658 + ,627.5 + ,528.9 + ,488.3 + ,575.5 + ,735.6 + ,685.3 + ,613.6 + ,629.5 + ,634.7 + ,652.6 + ,728.3 + ,634.3 + ,690.7 + ,676.3 + ,675.4 + ,595.6 + ,712.4 + ,735.8 + ,544.4 + ,567 + ,510 + ,564 + ,630.7 + ,496.7 + ,660.9 + ,601.2 + ,655.2 + ,591.6 + ,606.1 + ,560.7 + ,368.3 + ,371.6 + ,413.9 + ,413.9 + ,389 + ,399.2 + ,429.8 + ,395.6 + ,472 + ,486 + ,525 + ,396 + ,511 + ,525 + ,492 + ,517 + ,525 + ,474 + ,539 + ,468 + ,543 + ,532 + ,565 + ,535 + ,534 + ,546 + ,494 + ,552 + ,511 + ,451 + ,537 + ,494 + ,549 + ,544 + ,598 + ,583 + ,582 + ,589 + ,578 + ,561 + ,592 + ,504 + ,545 + ,547 + ,585 + ,562 + ,520 + ,581 + ,590 + ,562 + ,548 + ,567 + ,542 + ,473 + ,531 + ,462 + ,479 + ,533 + ,552 + ,547 + ,562 + ,524 + ,479 + ,445 + ,406 + ,475 + ,589 + ,495 + ,484 + ,536 + ,555 + ,565 + ,564 + ,573 + ,569 + ,588 + ,546 + ,508 + ,560 + ,558 + ,516 + ,549 + ,595 + ,586 + ,597 + ,592 + ,538 + ,590 + ,576 + ,451 + ,538 + ,555 + ,532 + ,530 + ,553 + ,626 + ,601 + ,573 + ,569 + ,562 + ,468 + ,483 + ,460 + ,411 + ,458 + ,455 + ,600 + ,605 + ,545 + ,549 + ,415 + ,568 + ,577 + ,517 + ,558 + ,518 + ,489 + ,502 + ,569 + ,540 + ,550 + ,557 + ,542 + ,542 + ,582 + ,525 + ,584 + ,562 + ,639 + ,613 + ,604 + ,613 + ,625 + ,654 + ,638) > 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] 455 > (np <- floor(n / par1)) [1] 37 > 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] 38 38 38 38 38 38 38 38 38 38 38 37 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 196.9 257.3 268.2 272.6 289.8 261.4 345.2 366.6 414.3 438.9 482.4 515.7 [2,] 192.1 234.5 231.2 248.6 279.6 287.0 314.1 336.3 404.3 394.5 407.9 461.0 [3,] 201.8 251.4 264.3 279.4 302.2 331.5 352.5 368.5 421.4 416.8 448.5 529.3 [4,] 186.9 243.8 252.7 270.5 299.1 326.9 361.2 374.2 372.6 402.9 431.1 467.4 [5,] 218.0 247.4 275.5 292.8 319.7 338.6 365.9 384.3 404.7 424.5 498.8 559.8 [6,] 214.4 245.3 261.5 297.8 310.9 337.0 332.5 358.9 420.2 432.3 497.1 536.5 [7,] 227.5 262.5 275.5 296.8 315.2 358.4 364.0 407.7 438.4 484.1 517.1 531.9 [8,] 204.1 270.0 272.3 290.9 338.5 344.5 359.1 433.3 449.1 492.7 487.7 546.5 [9,] 225.8 259.9 268.6 282.8 315.6 345.7 345.6 404.7 445.8 496.3 512.5 547.4 [10,] 223.7 262.2 270.4 312.8 321.2 344.1 366.9 392.7 413.8 471.9 550.1 536.1 [11,] 244.7 244.9 267.7 303.2 318.5 317.4 370.2 409.7 420.5 491.2 532.5 482.8 [12,] 243.9 249.3 275.0 301.4 342.7 354.5 359.9 416.5 442.3 512.9 524.1 551.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 532.9 588.6 578.1 567.3 493.4 600.3 689.8 678.9 672.1 679.3 668.7 663.6 [2,] 484.1 524.4 495.9 467.9 496.8 549.8 640.4 659.5 583.7 565.6 603.3 595.9 [3,] 554.8 618.5 526.8 588.9 529.9 635.5 675.2 691.9 705.5 564.1 674.5 638.4 [4,] 537.0 580.9 522.1 579.4 491.7 617.7 528.1 639.8 631.0 637.2 573.4 665.4 [5,] 558.0 557.2 605.1 502.6 543.2 643.5 538.0 663.8 733.3 650.8 706.0 671.4 [6,] 511.4 571.2 574.4 568.7 490.8 485.7 527.2 572.9 674.9 602.7 693.7 637.0 [7,] 502.9 597.5 609.7 616.0 554.7 689.5 655.6 592.5 695.5 587.5 627.5 685.7 [8,] 558.6 601.7 580.7 586.2 625.7 692.0 650.6 734.8 634.1 619.2 550.7 705.8 [9,] 545.1 558.9 565.1 575.5 605.0 677.3 623.7 696.1 630.6 616.5 592.3 704.8 [10,] 574.3 600.9 590.7 599.9 645.2 704.7 748.4 589.2 635.2 637.9 660.2 734.4 [11,] 542.2 601.0 571.5 568.2 645.2 668.6 727.4 662.9 554.1 557.9 597.3 674.2 [12,] 600.0 615.7 601.3 516.0 611.8 717.8 750.5 661.2 623.9 594.0 641.0 748.6 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 763.4 728.3 630.7 389.0 525 511 592 542 406 546 576 468 [2,] 658.0 634.3 496.7 399.2 474 451 504 473 475 508 451 483 [3,] 627.5 690.7 660.9 429.8 539 537 545 531 589 560 538 460 [4,] 528.9 676.3 601.2 395.6 468 494 547 462 495 558 555 411 [5,] 488.3 675.4 655.2 472.0 543 549 585 479 484 516 532 458 [6,] 575.5 595.6 591.6 486.0 532 544 562 533 536 549 530 455 [7,] 735.6 712.4 606.1 525.0 565 598 520 552 555 595 553 600 [8,] 685.3 735.8 560.7 396.0 535 583 581 547 565 586 626 605 [9,] 613.6 544.4 368.3 511.0 534 582 590 562 564 597 601 545 [10,] 629.5 567.0 371.6 525.0 546 589 562 524 573 592 573 549 [11,] 634.7 510.0 413.9 492.0 494 578 548 479 569 538 569 415 [12,] 652.6 564.0 413.9 517.0 552 561 567 445 588 590 562 568 [,37] [,38] [1,] 577 582 [2,] 517 525 [3,] 558 584 [4,] 518 562 [5,] 489 639 [6,] 502 613 [7,] 569 604 [8,] 540 613 [9,] 550 625 [10,] 557 654 [11,] 542 638 [12,] 542 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/www/html/rcomp/tmp/1n1tu1274089485.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2n1tu1274089485.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3n1tu1274089485.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4gsaf1274089485.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,] 196.90 231.2 201.8 186.9 218.0 245.3 315.2 204.10 225.80 223.7 244.7 243.9 [2,] 406.00 399.2 421.4 395.6 424.5 432.3 484.1 449.10 445.80 413.8 415.0 442.3 [3,] 537.45 479.0 538.5 506.5 524.0 531.0 560.0 562.85 554.45 564.5 540.0 561.0 [4,] 600.30 525.0 618.5 573.4 605.1 572.9 606.1 619.20 601.00 600.9 578.0 601.3 [5,] 763.40 659.5 705.5 676.3 733.3 693.7 735.6 735.80 704.80 748.4 727.4 750.5 $n [1] 38 38 38 38 38 38 38 38 38 38 38 37 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 487.649 446.7562 487.9813 460.9281 477.7104 494.9628 528.7302 519.2517 [2,] 587.251 511.2438 589.0187 552.0719 570.2896 567.0372 591.2698 606.4483 [,9] [,10] [,11] [,12] [1,] 514.6707 516.5444 498.2215 519.6997 [2,] 594.2293 612.4556 581.7785 602.3003 $out [1] 192.1 214.4 227.5 262.5 275.5 296.8 $group [1] 2 6 7 7 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(196.9, 406, 537.45, 600.3, 763.4, 231.2, 399.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5gsaf1274089485.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] [1,] 186.90 234.50 252.70 248.60 279.60 317.40 332.50 336.30 372.60 394.50 [2,] 199.35 245.10 262.90 276.00 300.65 322.15 345.40 367.55 409.25 420.65 [3,] 216.20 250.35 268.40 291.85 315.40 337.80 359.50 388.50 420.35 455.40 [4,] 226.65 261.05 273.65 299.60 320.45 345.10 364.95 408.70 440.35 491.95 [5,] 244.70 270.00 275.50 312.80 342.70 358.40 370.20 433.30 449.10 512.90 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [1,] 407.90 461.00 484.10 524.40 495.90 502.60 490.80 549.80 527.20 572.90 [2,] 465.45 499.25 522.15 565.05 545.95 541.65 495.10 609.00 580.85 616.15 [3,] 497.95 534.00 543.65 593.05 576.25 572.10 548.95 656.05 653.10 662.05 [4,] 520.60 546.95 558.30 601.35 596.00 587.55 618.75 690.75 708.60 685.40 [5,] 550.10 559.80 600.00 618.50 609.70 616.00 645.20 717.80 750.50 734.80 [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] 554.10 557.90 550.70 595.9 488.30 510.00 368.30 389.0 468.0 451.0 504 [2,] 627.25 576.55 594.80 651.0 594.55 565.50 413.90 397.6 509.5 524.0 546 [3,] 634.65 609.60 634.25 672.8 632.10 654.85 576.15 479.0 534.5 555.0 562 [4,] 685.20 637.55 671.60 705.3 671.65 701.55 618.40 514.0 544.5 582.5 583 [5,] 733.30 679.30 706.00 748.6 763.40 735.80 660.90 525.0 565.0 598.0 592 [,32] [,33] [,34] [,35] [,36] [,37] [,38] [1,] 445.0 406.0 508 530.0 411.0 489.0 525.0 [2,] 476.0 489.5 542 535.0 456.5 517.5 583.0 [3,] 527.5 559.5 559 558.5 475.5 542.0 613.0 [4,] 544.5 571.0 591 574.5 558.5 557.5 631.5 [5,] 562.0 589.0 597 626.0 605.0 577.0 654.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 [26] 12 12 12 12 12 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 203.7483 243.0751 263.4969 281.0859 306.3691 327.3324 350.5831 369.7312 [2,] 228.6517 257.6249 273.3031 302.6141 324.4309 348.2676 368.4169 407.2688 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 406.1651 422.8796 472.7957 512.2437 527.1617 576.4933 553.4219 551.1647 [2,] 434.5349 487.9204 523.1043 555.7563 560.1383 609.6067 599.0781 593.0353 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 492.5524 618.7633 594.8324 630.4646 608.2186 581.7775 599.221 648.0334 [2,] 605.3476 693.3367 711.3676 693.6354 661.0814 637.4225 669.279 697.5666 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 596.9342 592.7967 482.8762 425.9092 518.5363 528.3178 545.124 496.2567 [2,] 667.2658 716.9033 669.4238 532.0908 550.4637 581.6822 578.876 558.7433 [,33] [,34] [,35] [,36] [,37] [,38] [1,] 522.3273 536.6508 540.4838 428.9771 523.7557 589.8952 [2,] 596.6727 581.3492 576.5162 522.0229 560.2443 636.1048 $out [1] 231.2 261.4 287.0 314.1 467.9 485.7 451.0 $group [1] 3 6 6 7 16 18 35 $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" NA Warning message: In bxp(list(stats = c(186.9, 199.35, 216.2, 226.65, 244.7, 234.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6gsaf1274089485.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 120.2711 467.40 116.7750 [2,] 126.5070 481.00 138.0375 [3,] 132.3440 505.15 159.5000 [4,] 135.8473 520.00 172.7250 [5,] 148.3517 566.50 190.1500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 128.0838 487.3618 143.6788 [2,] 136.6042 522.9382 175.3212 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(120.271113636022, 126.506977092194, 132.344027209236, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1n1tu1274089485.ps tmp/1n1tu1274089485.png",intern=TRUE)) character(0) > try(system("convert tmp/2n1tu1274089485.ps tmp/2n1tu1274089485.png",intern=TRUE)) character(0) > try(system("convert tmp/3n1tu1274089485.ps tmp/3n1tu1274089485.png",intern=TRUE)) character(0) > try(system("convert tmp/4gsaf1274089485.ps tmp/4gsaf1274089485.png",intern=TRUE)) character(0) > try(system("convert tmp/5gsaf1274089485.ps tmp/5gsaf1274089485.png",intern=TRUE)) character(0) > try(system("convert tmp/6gsaf1274089485.ps tmp/6gsaf1274089485.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.140 0.834 1.458