R version 2.6.1 (2007-11-26) Copyright (C) 2007 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(88900,87280,85519,83647,81616,80100,94027,102327,104296,101593,94816,93535,93618,92330,90751,88576,86102,85494,103432,108870,109713,106960,103195,102348,102158,100431,97649,95611,93035,93579,111777,116065,116609,112934,107660,107965,107772,106201,102288,99217,96511,96456,113021,117836,118492,113922,109317,107496,105524,103824,101833,99436,96915,96072,111941,116008,117557,113445,108762,106661,102824,101912,99005,97894,96256,95606,108948,111223,113142,106078,100992,97413) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 88900 93618 102158 107772 105524 102824 NA [2,] 87280 92330 100431 106201 103824 101912 NA [3,] 85519 90751 97649 102288 101833 99005 NA [4,] 83647 88576 95611 99217 99436 97894 NA [5,] 81616 86102 93035 96511 96915 96256 NA [6,] 80100 85494 93579 96456 96072 95606 NA [7,] 94027 103432 111777 113021 111941 108948 NA [8,] 102327 108870 116065 117836 116008 111223 NA [9,] 104296 109713 116609 118492 117557 113142 NA [10,] 101593 106960 112934 113922 113445 106078 NA [11,] 94816 103195 107660 109317 108762 100992 NA [12,] 93535 102348 107965 107496 106661 97413 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1620 -1288 -1727 -1571 -1700 -912 NA [2,] -1761 -1579 -2782 -3913 -1991 -2907 NA [3,] -1872 -2175 -2038 -3071 -2397 -1111 NA [4,] -2031 -2474 -2576 -2706 -2521 -1638 NA [5,] -1516 -608 544 -55 -843 -650 NA [6,] 13927 17938 18198 16565 15869 13342 NA [7,] 8300 5438 4288 4815 4067 2275 NA [8,] 1969 843 544 656 1549 1919 NA [9,] -2703 -2753 -3675 -4570 -4112 -7064 NA [10,] -6777 -3765 -5274 -4605 -4683 -5086 NA [11,] -1281 -847 305 -1821 -2101 -3579 NA [12,] 83 -190 -193 -1972 -3837 NA 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/1ww7t1197740032.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/245j91197740032.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/36si81197740032.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/4kj6a1197740032.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] [1,] 88900 87280.0 85519 83647.0 81616.0 80100.0 94027.0 102327.0 104296.0 [2,] 93618 92330.0 90751 88576.0 86102.0 85494.0 103432.0 108870.0 109713.0 [3,] 102491 101171.5 98327 96752.5 94645.5 94592.5 110362.5 113615.5 114875.5 [4,] 105524 103824.0 101833 99217.0 96511.0 96072.0 111941.0 116065.0 117557.0 [5,] 107772 106201.0 102288 99436.0 96915.0 96456.0 113021.0 117836.0 118492.0 [,10] [,11] [,12] [1,] 101593 94816.0 93535.0 [2,] 106078 100992.0 97413.0 [3,] 109947 105427.5 104504.5 [4,] 113445 108762.0 107496.0 [5,] 113922 109317.0 107965.0 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 94811.25 93757.5 91178.75 89888.71 87931.36 87769.35 104873.9 [2,] 110170.75 108585.5 105475.25 103616.29 101359.64 101415.65 115851.1 [,8] [,9] [,10] [,11] [,12] [1,] 108974.5 109815.9 105195.0 100415.6 98000.64 [2,] 118256.5 119935.1 114699.0 110439.4 111008.36 $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(88900, 93618, 102491, 105524, 107772, 87280, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/593l91197740032.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] [1,] -1727.0 -3913.0 -3071.0 -2706.0 -1516 13342 2275.0 544 -7064.0 -5274.0 [2,] -1700.0 -2907.0 -2397.0 -2576.0 -843 13927 4067.0 656 -4570.0 -5274.0 [3,] -1595.5 -2386.5 -2106.5 -2497.5 -629 16217 4551.5 1196 -3893.5 -4884.5 [4,] -1288.0 -1761.0 -1872.0 -2031.0 -55 17938 5438.0 1919 -2753.0 -4605.0 [5,] -912.0 -1579.0 -1111.0 -1638.0 544 18198 5438.0 1969 -2703.0 -3765.0 [,11] [,12] [1,] -3579 -3837 [2,] -2101 -1972 [3,] -1551 -193 [4,] -847 -190 [5,] 305 83 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1861.253 -3125.707 -2445.142 -2849.043 -1137.2855 13629.78 3667.161 [2,] -1329.747 -1647.293 -1767.858 -2145.957 -120.7145 18804.22 5435.839 [,8] [,9] [,10] [,11] [,12] [1,] 381.3242 -5065.524 -5316.027 -2359.8705 -1452.157 [2,] 2010.6758 -2721.476 -4452.973 -742.1295 1066.157 $out [1] 8300 -6777 $group [1] 7 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1727, -1700, -1595.5, -1288, -912, -3913, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6mj4c1197740032.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] [1,] 80100.0 85494.0 93035.0 96456.0 96072.0 95606.0 NA [2,] 84583.0 89663.5 96630.0 100752.5 100634.5 97653.5 NA [3,] 91217.5 97983.0 104909.0 107634.0 106092.5 101452.0 NA [4,] 98204.5 105196.0 112355.5 113471.5 112693.0 107513.0 NA [5,] 104296.0 109713.0 116609.0 118492.0 117557.0 113142.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 85004.64 90898.52 97736.5 101832.8 100592.5 96955.02 NA [2,] 97430.36 105067.48 112081.5 113435.2 111592.5 105948.98 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(80100, 84583, 91217.5, 98204.5, 104296, 85494, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7slhm1197740032.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,] 91217.83 94592.50 91735.38 [2,] 95118.83 97539.75 95705.62 [3,] 101351.17 103497.75 101634.00 [4,] 108173.17 110154.75 109081.69 [5,] 113301.50 114875.50 113945.12 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 95397.0 97743.96 95533.09 [2,] 107305.3 109251.54 107734.91 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(91217.8333333333, 95118.8333333333, 101351.166666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1ww7t1197740032.ps tmp/1ww7t1197740032.png") > system("convert tmp/245j91197740032.ps tmp/245j91197740032.png") > system("convert tmp/36si81197740032.ps tmp/36si81197740032.png") > system("convert tmp/4kj6a1197740032.ps tmp/4kj6a1197740032.png") > system("convert tmp/593l91197740032.ps tmp/593l91197740032.png") > system("convert tmp/6mj4c1197740032.ps tmp/6mj4c1197740032.png") > system("convert tmp/7slhm1197740032.ps tmp/7slhm1197740032.png") > > > proc.time() user system elapsed 2.568 1.592 2.795