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(284.4,212.8,226.9,308.4,262,227.9,236.1,320.4,271.9,232.8,237,313.4,261.4,226.8,249.9,314.3,286.1,226.5,260.4,311.4,294.7,232.6,257.2,339.2,279.1,249.8,269.8,345.7,293.8,254.7,277.5,363.4,313.4,272.8,300.1,369.5,330.8,287.8,305.9,386.1,335.2,288,308.3,402.3,352.8,316.1,324.9,404.8,393,318.9,327,442.3,383.1,331.6,361.4,445.9,386.6,357.2,373.6,466.2,409.6,369.8,378.6,487,419.2,376.7,392.8,506.1,458.4,387.4,426.9,565,464.8,444.5,449.5,556.1,499.6,451.9,434.9,553.8,510,432.9,453.2,547.6,485.8,452.6,456.6,565.7,514.8,464.3,430.9,588.3,503.1,442.6,448,554.5,504.5,427.3,473.1,526.2,547.5,440.2,468.7,574.5,492.6,432.6,479.8,575.7,474.6,405.3,434.6,535.1,452.6,429.5,417.2,551.8,464,416.6,422.9,553.6,458.6,427.6,429.2,534.2,481.7,416,440.2,538.7,473.8,439.9,446.8,597.5,467.2,439.4,447.4,568.5,485.9,442.1,430.5,600,464.5,423.6,437,574,443,410,420,532,432,420,411,512) > par1 = '4' > #'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] 152 > (np <- floor(n / par1)) [1] 38 > 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 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 284.4 262.0 271.9 261.4 286.1 294.7 279.1 293.8 313.4 330.8 335.2 352.8 [2,] 212.8 227.9 232.8 226.8 226.5 232.6 249.8 254.7 272.8 287.8 288.0 316.1 [3,] 226.9 236.1 237.0 249.9 260.4 257.2 269.8 277.5 300.1 305.9 308.3 324.9 [4,] 308.4 320.4 313.4 314.3 311.4 339.2 345.7 363.4 369.5 386.1 402.3 404.8 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 393.0 383.1 386.6 409.6 419.2 458.4 464.8 499.6 510.0 485.8 514.8 503.1 [2,] 318.9 331.6 357.2 369.8 376.7 387.4 444.5 451.9 432.9 452.6 464.3 442.6 [3,] 327.0 361.4 373.6 378.6 392.8 426.9 449.5 434.9 453.2 456.6 430.9 448.0 [4,] 442.3 445.9 466.2 487.0 506.1 565.0 556.1 553.8 547.6 565.7 588.3 554.5 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 504.5 547.5 492.6 474.6 452.6 464.0 458.6 481.7 473.8 467.2 485.9 464.5 [2,] 427.3 440.2 432.6 405.3 429.5 416.6 427.6 416.0 439.9 439.4 442.1 423.6 [3,] 473.1 468.7 479.8 434.6 417.2 422.9 429.2 440.2 446.8 447.4 430.5 437.0 [4,] 526.2 574.5 575.7 535.1 551.8 553.6 534.2 538.7 597.5 568.5 600.0 574.0 [,37] [,38] [,39] [1,] 443 432 NA [2,] 410 420 NA [3,] 420 411 NA [4,] 532 512 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/1nxqw1231413547.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/2irqo1231413547.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/3eta51231413547.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/4nn9j1231413547.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] [1,] 261.4 212.80 226.9 308.4 [2,] 330.8 287.80 305.9 386.1 [3,] 447.8 407.65 418.6 529.1 [4,] 481.7 432.90 440.2 556.1 [5,] 547.5 464.30 479.8 600.0 $n [1] 38 38 38 38 $conf [,1] [,2] [,3] [,4] [1,] 409.1228 370.4594 384.1776 485.5273 [2,] 486.4772 444.8406 453.0224 572.6727 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" Warning message: In bxp(list(stats = c(261.4, 330.8, 447.8, 481.7, 547.5, 212.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5lqd81231413547.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,] 212.80 227.90 232.80 226.80 226.50 232.60 249.80 254.70 272.80 287.80 [2,] 219.85 232.00 234.90 238.35 243.45 244.90 259.80 266.10 286.45 296.85 [3,] 255.65 249.05 254.45 255.65 273.25 275.95 274.45 285.65 306.75 318.35 [4,] 296.40 291.20 292.65 287.85 298.75 316.95 312.40 328.60 341.45 358.45 [5,] 308.40 320.40 313.40 314.30 311.40 339.20 345.70 363.40 369.50 386.10 [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] 288.00 316.10 318.90 331.60 357.2 369.8 376.70 387.40 444.50 434.90 432.90 [2,] 298.15 320.50 322.95 346.50 365.4 374.2 384.75 407.15 447.00 443.40 443.05 [3,] 321.75 338.85 360.00 372.25 380.1 394.1 406.00 442.65 457.15 475.75 481.60 [4,] 368.75 378.80 417.65 414.50 426.4 448.3 462.65 511.70 510.45 526.70 528.80 [5,] 402.30 404.80 442.30 445.90 466.2 487.0 506.10 565.00 556.10 553.80 547.60 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [1,] 452.60 430.90 442.60 427.30 440.20 432.60 405.30 417.20 416.60 427.6 [2,] 454.60 447.60 445.30 450.20 454.45 456.20 419.95 423.35 419.75 428.4 [3,] 471.20 489.55 475.55 488.80 508.10 486.20 454.60 441.05 443.45 443.9 [4,] 525.75 551.55 528.80 515.35 561.00 534.15 504.85 502.20 508.80 496.4 [5,] 565.70 588.30 554.50 526.20 574.50 575.70 535.10 551.80 553.60 534.2 [,32] [,33] [,34] [,35] [,36] [,37] [,38] [,39] [1,] 416.00 439.90 439.40 430.50 423.60 410.0 411.0 NA [2,] 428.10 443.35 443.40 436.30 430.30 415.0 415.5 NA [3,] 460.95 460.30 457.30 464.00 450.75 431.5 426.0 NA [4,] 510.20 535.65 517.85 542.95 519.25 487.5 472.0 NA [5,] 538.70 597.50 568.50 600.00 574.00 532.0 512.0 NA $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 [39] 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 195.1755 202.282 208.8275 216.545 229.563 219.0305 232.896 236.275 263.3 [2,] 316.1245 295.818 300.0725 294.755 316.937 332.8695 316.004 335.025 350.2 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [1,] 269.686 265.976 292.793 285.187 318.53 331.91 335.561 344.459 360.0555 [2,] 367.014 377.524 384.907 434.813 425.97 428.29 452.639 467.541 525.2445 [,19] [,20] [,21] [,22] [,23] [,24] [,25] [,26] [1,] 407.0245 409.943 413.8575 414.9915 407.4295 409.585 437.3315 423.9255 [2,] 507.2755 541.557 549.3425 527.4085 571.6705 541.515 540.2685 592.2745 [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [1,] 424.6195 387.529 378.7585 373.1005 390.18 396.091 387.383 398.4845 [2,] 547.7805 521.671 503.3415 513.7995 497.62 525.809 533.217 516.1155 [,35] [,36] [,37] [,38] [,39] [1,] 379.7465 380.4795 374.225 381.365 NA [2,] 548.2535 521.0205 488.775 470.635 NA $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" NA Warning message: In bxp(list(stats = c(212.8, 219.85, 255.65, 296.4, 308.4, 227.9, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6uzes1231413547.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,] 80.76178 251.50 132.9000 [2,] 82.62147 252.20 138.9375 [3,] 85.86297 269.50 146.5000 [4,] 93.12732 288.85 156.7875 [5,] 99.00985 291.60 165.5500 $n [1] 4 4 4 $conf [,1] [,2] [,3] [1,] 77.56334 240.5465 132.3985 [2,] 94.16259 298.4535 160.6015 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(80.7617790183783, 82.6214656758718, 85.8629675689209, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1nxqw1231413547.ps tmp/1nxqw1231413547.png") > system("convert tmp/2irqo1231413547.ps tmp/2irqo1231413547.png") > system("convert tmp/3eta51231413547.ps tmp/3eta51231413547.png") > system("convert tmp/4nn9j1231413547.ps tmp/4nn9j1231413547.png") > system("convert tmp/5lqd81231413547.ps tmp/5lqd81231413547.png") > system("convert tmp/6uzes1231413547.ps tmp/6uzes1231413547.png") > > > proc.time() user system elapsed 1.201 0.967 4.692