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(374 + ,572 + ,402 + ,589 + ,507 + ,628 + ,698 + ,451 + ,694 + ,0 + ,488 + ,526 + ,343 + ,494 + ,447 + ,0 + ,470 + ,366 + ,517 + ,483 + ,485 + ,530 + ,308 + ,481 + ,437 + ,468 + ,502 + ,408 + ,479 + ,436 + ,410 + ,451 + ,344 + ,411 + ,0 + ,427 + ,454 + ,365 + ,499 + ,416 + ,430 + ,470 + ,325 + ,452 + ,442 + ,488 + ,446 + ,523 + ,594 + ,439 + ,588 + ,503 + ,444 + ,525 + ,375 + ,472 + ,436 + ,458 + ,514 + ,0 + ,472 + ,360 + ,450 + ,549 + ,361 + ,466 + ,387 + ,457 + ,470 + ,396 + ,471 + ,422 + ,404 + ,414 + ,342 + ,459 + ,379 + ,0 + ,410 + ,319 + ,411 + ,371 + ,365 + ,429 + ,333 + ,392 + ,469 + ,432 + ,534 + ,379 + ,436 + ,448 + ,358 + ,492 + ,387 + ,529 + ,475 + ,439 + ,459 + ,361 + ,0 + ,0 + ,394 + ,425 + ,341 + ,455 + ,403 + ,471 + ,523 + ,389 + ,531 + ,468 + ,398 + ,446 + ,355 + ,435 + ,353 + ,0 + ,400 + ,332 + ,389 + ,355 + ,384 + ,406 + ,356 + ,336 + ,351 + ,278 + ,265 + ,229 + ,387 + ,435 + ,317 + ,490 + ,472 + ,440 + ,429 + ,350 + ,489 + ,494 + ,436 + ,436 + ,375 + ,429 + ,0 + ,434 + ,472 + ,362 + ,440 + ,433 + ,400 + ,442 + ,316 + ,432 + ,401 + ,434 + ,488 + ,377 + ,484 + ,377 + ,0 + ,0 + ,300 + ,389 + ,337 + ,376 + ,377 + ,331 + ,339 + ,356 + ,280 + ,249 + ,196 + ,268 + ,379 + ,401 + ,404 + ,397 + ,419 + ,421 + ,407 + ,296 + ,468 + ,475 + ,422 + ,456 + ,339 + ,446 + ,419 + ,346 + ,327 + ,326 + ,403 + ,359 + ,358 + ,421 + ,322 + ,367 + ,394 + ,356 + ,418 + ,344 + ,372 + ,358 + ,373 + ,379 + ,0 + ,348 + ,369 + ,341 + ,390 + ,279 + ,325 + ,354 + ,346 + ,358 + ,296 + ,356 + ,337 + ,360 + ,474 + ,362 + ,440 + ,443 + ,435 + ,429 + ,341 + ,434 + ,329 + ,416 + ,430 + ,307 + ,408 + ,322 + ,0 + ,324 + ,303 + ,369 + ,328 + ,258 + ,372 + ,298 + ,376 + ,306 + ,359 + ,418 + ,311 + ,355 + ,335 + ,345 + ,318 + ,291 + ,340 + ,0 + ,356 + ,419 + ,296 + ,361 + ,371 + ,392 + ,383 + ,286 + ,362 + ,358) > 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] 264 > (np <- floor(n / par1)) [1] 22 > 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] 22 22 22 22 22 22 22 22 22 22 22 22 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 374 343 437 454 594 472 404 333 475 523 389 317 0 [2,] 572 494 468 365 439 360 414 392 439 389 355 490 434 [3,] 402 447 502 499 588 450 342 469 459 531 384 472 472 [4,] 589 0 408 416 503 549 459 432 361 468 406 440 362 [5,] 507 470 479 430 444 361 379 534 0 398 356 429 440 [6,] 628 366 436 470 525 466 0 379 0 446 336 350 433 [7,] 698 517 410 325 375 387 410 436 394 355 351 489 400 [8,] 451 483 451 452 472 457 319 448 425 435 278 494 442 [9,] 694 485 344 442 436 470 411 358 341 353 265 436 316 [10,] 0 530 411 488 458 396 371 492 455 0 229 436 432 [11,] 488 308 0 446 514 471 365 387 403 400 387 375 401 [12,] 526 481 427 523 0 422 429 529 471 332 435 429 434 [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 488 339 407 403 373 296 329 372 340 NA [2,] 377 356 296 359 379 356 416 298 0 NA [3,] 484 280 468 358 0 337 430 376 356 NA [4,] 377 249 475 421 348 360 307 306 419 NA [5,] 0 196 422 322 369 474 408 359 296 NA [6,] 0 268 456 367 341 362 322 418 361 NA [7,] 300 379 339 394 390 440 0 311 371 NA [8,] 389 401 446 356 279 443 324 355 392 NA [9,] 337 404 419 418 325 435 303 335 383 NA [10,] 376 397 346 344 354 429 369 345 286 NA [11,] 377 419 327 372 346 341 328 318 362 NA [12,] 331 421 326 358 358 434 258 291 358 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/1br511243890473.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/29nau1243890473.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/3o41b1243890473.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/4cj7t1243890473.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,] 296.0 296 280.0 249 296 268.0 300.0 278.0 265.0 229 308 258.0 [2,] 339.0 356 358.0 360 356 336.0 351.0 356.0 337.0 345 341 332.0 [3,] 381.5 384 448.5 412 403 366.5 388.5 438.5 393.5 386 376 424.5 [4,] 454.0 439 472.0 459 444 446.0 410.0 451.0 436.0 436 403 435.0 [5,] 594.0 494 588.0 589 534 525.0 489.0 494.0 485.0 530 488 529.0 $n [1] 22 22 22 22 22 22 22 22 22 22 22 22 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 342.7614 356.0409 410.0983 378.6511 373.3566 329.4457 368.6254 406.4986 [2,] 420.2386 411.9591 486.9017 445.3489 432.6434 403.5543 408.3746 470.5014 [,9] [,10] [,11] [,12] [1,] 360.1511 355.346 355.1149 389.8037 [2,] 426.8489 416.654 396.8851 459.1963 $out [1] 0 572 0 0 0 0 0 196 628 0 0 0 698 517 0 694 0 0 0 [20] 514 0 $group [1] 1 2 2 3 4 5 5 5 6 6 6 6 7 7 7 9 10 10 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(296, 339, 381.5, 454, 594, 296, 356, 384, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/50skv1243890473.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] [,12] [1,] 374.0 308.0 344.0 365 375.0 360.0 319.0 333.0 341 332 229.0 317.0 [2,] 426.5 354.5 409.0 423 437.5 391.5 353.5 383.0 351 354 307.0 402.0 [3,] 516.5 475.5 431.5 449 465.0 453.5 391.5 434.0 414 399 355.5 436.0 [4,] 608.5 489.5 459.5 479 519.5 470.5 412.5 480.5 457 457 388.0 480.5 [5,] 698.0 530.0 502.0 523 594.0 549.0 459.0 534.0 475 531 435.0 494.0 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 316.0 300.0 196.0 296 322.0 279 296.0 258.0 291.0 286.0 NA [2,] 381.0 315.5 274.0 333 357.0 333 348.5 305.0 308.5 318.0 NA [3,] 432.5 376.5 367.5 413 363.0 351 395.5 326.0 340.0 359.5 NA [4,] 437.0 383.0 402.5 451 398.5 371 437.5 388.5 365.5 377.0 NA [5,] 472.0 484.0 421.0 475 421.0 390 474.0 430.0 418.0 419.0 NA $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 433.4886 413.9256 408.4666 423.458 427.5992 417.4676 364.5897 389.5296 [2,] 599.5114 537.0744 454.5334 474.542 502.4008 489.5324 418.4103 478.4704 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 365.6527 352.021 318.5554 400.1956 406.958 345.7128 308.8903 359.1794 [2,] 462.3473 445.979 392.4446 471.8044 458.042 407.2872 426.1097 466.8206 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] 344.0716 333.6679 354.9065 287.9151 314.0019 332.5897 NA [2,] 381.9284 368.3321 436.0935 364.0849 365.9981 386.4103 NA $out [1] 0 0 0 325 0 0 0 0 0 0 488 0 0 0 0 0 $group [1] 1 2 3 4 5 7 9 9 10 13 14 14 14 18 20 22 $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" NA Warning message: In bxp(list(stats = c(374, 426.5, 516.5, 608.5, 698, 308, 354.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6qs0z1243890473.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,] 88.36715 429.0 58.000 [2,] 103.20510 521.5 83.625 [3,] 116.13866 553.0 91.875 [4,] 127.74019 591.5 102.000 [5,] 162.13179 628.0 110.500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 104.9480 521.0725 83.49404 [2,] 127.3293 584.9275 100.25596 $out [1] 64.09827 698.00000 216.00000 $group [1] 1 2 2 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(88.3671502695048, 103.205102996635, 116.138663158088, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1br511243890473.ps tmp/1br511243890473.png") > system("convert tmp/29nau1243890473.ps tmp/29nau1243890473.png") > system("convert tmp/3o41b1243890473.ps tmp/3o41b1243890473.png") > system("convert tmp/4cj7t1243890473.ps tmp/4cj7t1243890473.png") > system("convert tmp/50skv1243890473.ps tmp/50skv1243890473.png") > system("convert tmp/6qs0z1243890473.ps tmp/6qs0z1243890473.png") > > > proc.time() user system elapsed 1.104 0.836 1.394