R version 2.7.0 (2008-04-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(492.865,480.961,461.935,456.608,441.977,439.148,488.180,520.564,501.492,485.025,464.196,460.170,467.037,460.070,447.988,442.867,436.087,431.328,484.015,509.673,512.927,502.831,470.984,471.067,476.049,474.605,470.439,461.251,454.724,455.626,516.847,525.192,522.975,518.585,509.239,512.238,519.164,517.009,509.933,509.127,500.857,506.971,569.323,579.714,577.992,565.464,547.344,554.788,562.325,560.854,555.332,543.599,536.662,542.722,593.530,610.763,612.613,611.324,594.167,595.454,590.865,589.379,584.428,573.100,567.456,569.028,620.735,628.884,628.232,612.117,595.404,597.141,593.408,590.072,579.799,574.205,572.775,572.942,619.567,625.809,619.916,587.625,565.742,557.274,560.576,548.854,531.673,525.919,511.038,498.662,555.362,564.591,541.657,527.070,509.846,514.258,516.922) > 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] 97 > (np <- floor(n / par1)) [1] 8 > 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] 9 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 492.865 467.037 476.049 519.164 562.325 590.865 593.408 560.576 516.922 [2,] 480.961 460.070 474.605 517.009 560.854 589.379 590.072 548.854 NA [3,] 461.935 447.988 470.439 509.933 555.332 584.428 579.799 531.673 NA [4,] 456.608 442.867 461.251 509.127 543.599 573.100 574.205 525.919 NA [5,] 441.977 436.087 454.724 500.857 536.662 567.456 572.775 511.038 NA [6,] 439.148 431.328 455.626 506.971 542.722 569.028 572.942 498.662 NA [7,] 488.180 484.015 516.847 569.323 593.530 620.735 619.567 555.362 NA [8,] 520.564 509.673 525.192 579.714 610.763 628.884 625.809 564.591 NA [9,] 501.492 512.927 522.975 577.992 612.613 628.232 619.916 541.657 NA [10,] 485.025 502.831 518.585 565.464 611.324 612.117 587.625 527.070 NA [11,] 464.196 470.984 509.239 547.344 594.167 595.404 565.742 509.846 NA [12,] 460.170 471.067 512.238 554.788 595.454 597.141 557.274 514.258 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] -11.904 -6.967 -1.444 -2.155 -1.471 -1.486 -3.336 -11.722 NA [2,] -19.026 -12.082 -4.166 -7.076 -5.522 -4.951 -10.273 -17.181 NA [3,] -5.327 -5.121 -9.188 -0.806 -11.733 -11.328 -5.594 -5.754 NA [4,] -14.631 -6.780 -6.527 -8.270 -6.937 -5.644 -1.430 -14.881 NA [5,] -2.829 -4.759 0.902 6.114 6.060 1.572 0.167 -12.376 NA [6,] 49.032 52.687 61.221 62.352 50.808 51.707 46.625 56.700 NA [7,] 32.384 25.658 8.345 10.391 17.233 8.149 6.242 9.229 NA [8,] -19.072 3.254 -2.217 -1.722 1.850 -0.652 -5.893 -22.934 NA [9,] -16.467 -10.096 -4.390 -12.528 -1.289 -16.115 -32.291 -14.587 NA [10,] -20.829 -31.847 -9.346 -18.120 -17.157 -16.713 -21.883 -17.224 NA [11,] -4.026 0.083 2.999 7.444 1.287 1.737 -8.468 4.412 NA [12,] 6.867 4.982 6.926 7.537 -4.589 -3.733 3.302 2.664 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/1beg41212254055.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/2mlhj1212254055.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/331pq1212254055.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/4tttz1212254055.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] [1,] 467.037 460.0700 447.9880 442.8670 436.0870 431.3280 484.0150 509.6730 [2,] 492.865 477.7830 466.1870 458.9295 448.3505 447.3870 502.5135 522.8780 [3,] 519.164 532.9315 520.8030 517.5230 505.9475 502.8165 562.3425 572.1525 [4,] 562.325 575.1165 567.5655 558.3495 552.0590 555.8750 606.5485 618.2860 [5,] 593.408 590.0720 584.4280 574.2050 572.7750 572.9420 620.7350 628.8840 [,9] [,10] [,11] [,12] [1,] 501.4920 485.0250 464.1960 460.1700 [2,] 517.9510 510.7080 490.1115 491.6525 [3,] 559.8245 546.2670 528.5950 534.5230 [4,] 616.2645 599.4745 579.9545 576.3640 [5,] 628.2320 612.1170 595.4040 597.1410 $n [1] 9 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 482.5817 478.5596 464.1715 461.9856 448.0144 442.2135 504.2271 518.8562 [2,] 555.7463 587.3034 577.4345 573.0604 563.8806 563.4195 620.4579 625.4488 [,9] [,10] [,11] [,12] [1,] 504.9052 496.6808 478.4074 487.2019 [2,] 614.7438 595.8532 578.7826 581.8441 $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(467.037, 492.865, 519.164, 562.325, 593.408, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5k1of1212254055.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] [1,] -11.9040 -19.0260 -11.733 -14.8810 -12.3760 46.6250 6.2420 -22.9340 [2,] -9.3445 -14.6315 -10.258 -11.4505 -3.7940 49.9200 8.2470 -12.4825 [3,] -2.7455 -8.6745 -5.674 -6.8585 0.5345 52.1970 9.8100 -1.9695 [4,] -1.4785 -5.2365 -5.224 -6.0855 3.8160 58.9605 21.4455 0.5990 [5,] -1.4440 -4.1660 -0.806 -1.4300 6.1140 62.3520 32.3840 3.2540 [,9] [,10] [,11] [,12] [1,] -16.4670 -21.883 -8.4680 -4.5890 [2,] -16.2910 -21.356 -1.9715 -0.5345 [3,] -13.5575 -17.672 1.5120 4.1420 [4,] -7.2430 -16.935 3.7055 6.8965 [5,] -1.2890 -16.713 7.4440 7.5370 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -7.139561 -13.922682 -8.486065 -9.855466 -3.716555 47.14685 2.437128 [2,] 1.648561 -3.426318 -2.861935 -3.861534 4.785555 57.24715 17.182872 [,8] [,9] [,10] [,11] [,12] [1,] -9.277014 -18.611843 -20.14163 -1.659254 -0.009063288 [2,] 5.338014 -8.503157 -15.20237 4.683254 8.293063288 $out [1] -32.291 -31.847 -9.346 $group [1] 9 10 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-11.904, -9.34449999999998, -2.74549999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68iwa1212254055.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] [1,] 439.1480 431.3280 454.724 500.8570 536.6620 567.4560 557.2740 498.6620 [2,] 458.3890 445.4275 465.845 509.5300 549.4655 578.7640 572.8585 512.6480 [3,] 472.5785 469.0105 492.644 533.2540 577.9275 593.1345 583.7120 529.3715 [4,] 490.5225 493.4230 517.716 567.3935 603.1085 616.4260 606.4875 552.1080 [5,] 520.5640 512.9270 525.192 579.7140 612.6130 628.8840 625.8090 564.5910 [,9] [1,] 516.922 [2,] 516.922 [3,] 516.922 [4,] 516.922 [5,] 516.922 $n [1] 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 457.9222 447.1194 468.9853 506.8621 553.4606 575.9566 568.3736 511.3735 [2,] 487.2348 490.9016 516.3027 559.6459 602.3944 610.3124 599.0504 547.3695 [,9] [1,] 516.922 [2,] 516.922 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(439.148, 458.389, 472.5785, 490.5225, 520.564, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7vrrw1212254055.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,] 502.0534 502.8165 497.9489 [2,] 514.2627 518.3435 510.2066 [3,] 531.5693 530.7632 530.9886 [4,] 553.6000 553.0457 554.4789 [5,] 570.6488 572.1525 569.2797 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 513.6273 514.9353 510.7956 [2,] 549.5114 546.5912 551.1815 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(502.053375, 514.2626875, 531.569347222222, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1beg41212254055.ps tmp/1beg41212254055.png") > system("convert tmp/2mlhj1212254055.ps tmp/2mlhj1212254055.png") > system("convert tmp/331pq1212254055.ps tmp/331pq1212254055.png") > system("convert tmp/4tttz1212254055.ps tmp/4tttz1212254055.png") > system("convert tmp/5k1of1212254055.ps tmp/5k1of1212254055.png") > system("convert tmp/68iwa1212254055.ps tmp/68iwa1212254055.png") > system("convert tmp/7vrrw1212254055.ps tmp/7vrrw1212254055.png") > > > proc.time() user system elapsed 4.779 2.787 5.131