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(831581,808744,899237,929532,883165,908232,955613,937590,849396,978630,868513,1156102,1505713,1415151,1545021,1681193,1457973,1638575,1688972,1563924,1596359,1722061,1549332,2264959,1420268,1415099,1597279,1605693,1575400,1654752,1553966,1570959,1642414,1664774,1551560,2304365,1644081,1425600,1569344,1456489,1610786,1601519,1496600,1486452,1637939,1605759,1504221,1993384,1507620,1477037,1679184,1504731,1570141,1734191,1657498,1652164,1610941,1813765,1711573,2165466,1492778,1385488,1470589,1514657,1641395,1606185,1581162,1517847,1630080,1604623,1548973,2125558) > 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,] 831581 1505713 1420268 1644081 1507620 1492778 NA [2,] 808744 1415151 1415099 1425600 1477037 1385488 NA [3,] 899237 1545021 1597279 1569344 1679184 1470589 NA [4,] 929532 1681193 1605693 1456489 1504731 1514657 NA [5,] 883165 1457973 1575400 1610786 1570141 1641395 NA [6,] 908232 1638575 1654752 1601519 1734191 1606185 NA [7,] 955613 1688972 1553966 1496600 1657498 1581162 NA [8,] 937590 1563924 1570959 1486452 1652164 1517847 NA [9,] 849396 1596359 1642414 1637939 1610941 1630080 NA [10,] 978630 1722061 1664774 1605759 1813765 1604623 NA [11,] 868513 1549332 1551560 1504221 1711573 1548973 NA [12,] 1156102 2264959 2304365 1993384 2165466 2125558 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22837 -90562 -5169 -218481 -30583 -107290 NA [2,] 90493 129870 182180 143744 202147 85101 NA [3,] 30295 136172 8414 -112855 -174453 44068 NA [4,] -46367 -223220 -30293 154297 65410 126738 NA [5,] 25067 180602 79352 -9267 164050 -35210 NA [6,] 47381 50397 -100786 -104919 -76693 -25023 NA [7,] -18023 -125048 16993 -10148 -5334 -63315 NA [8,] -88194 32435 71455 151487 -41223 112233 NA [9,] 129234 125702 22360 -32180 202824 -25457 NA [10,] -110117 -172729 -113214 -101538 -102192 -55650 NA [11,] 287589 715627 752805 489163 453893 576585 NA [12,] 349611 -844691 -660284 -485764 -672688 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/1pufn1242130420.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/2c4tg1242130420.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/3c09c1242130420.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/4ldol1242130420.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,] 1420268 1385488 1470589 1456489 1457973 1601519 1496600 1486452 1596359 [2,] 1420268 1385488 1470589 1456489 1457973 1601519 1496600 1486452 1596359 [3,] 1499246 1415125 1557182 1509694 1572770 1622380 1567564 1540886 1620510 [4,] 1507620 1425600 1597279 1605693 1610786 1654752 1657498 1570959 1637939 [5,] 1507620 1477037 1679184 1681193 1641395 1734191 1688972 1652164 1642414 [,10] [,11] [,12] [1,] 1604623 1504221 1993384 [2,] 1604623 1504221 1993384 [3,] 1635266 1549152 2145512 [4,] 1722061 1551560 2264959 [5,] 1813765 1551560 2304365 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 1442901 1389251 1475463 1413453 1474201 1588043 1463780 1486376 1593690 [2,] 1555590 1440999 1638902 1605935 1671340 1656717 1671348 1595395 1647331 [,10] [,11] [,12] [1,] 1559515 1518617 1970337 [2,] 1711018 1579688 2320687 $out [1] 831581 1644081 808744 899237 929532 883165 908232 955613 937590 [10] 849396 978630 868513 1711573 1156102 $group [1] 1 1 2 3 4 5 6 7 8 9 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(1420268, 1420268, 1499245.5, 1507620, 1507620, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5k4fj1242130420.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,] -218481.0 85101 -174453.0 -223220.0 -35210.0 -104919 -125048.0 -88194 [2,] -107290.0 90493 -112855.0 -46367.0 -9267.0 -100786 -63315.0 -41223 [3,] -60572.5 136807 19354.5 17558.5 52209.5 -50858 -14085.5 51945 [4,] -22837.0 182180 44068.0 126738.0 164050.0 47381 -5334.0 112233 [5,] -5169.0 202147 136172.0 154297.0 180602.0 50397 16993.0 151487 [,9] [,10] [,11] [,12] [1,] -32180 -113214.0 287589 -844691 [2,] -25457 -113214.0 453893 -672688 [3,] 74031 -106154.5 532874 -660284 [4,] 129234 -101538.0 715627 -485764 [5,] 202824 -101538.0 752805 -485764 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -115047.413 77665.92 -81865.9 -94099.82 -59585.56 -146430.5 -51485.12 [2,] -6097.587 195948.08 120574.9 129216.82 164004.56 44714.5 23314.12 [,8] [,9] [,10] [,11] [,12] [1,] -47039.08 -25749.69 -113685.9 364047.1 -792364 [2,] 150929.08 173811.69 -98623.1 701700.9 -528204 $out [1] -172729 -55650 349611 $group [1] 10 10 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-218481, -107290, -60572.5, -22837, -5169, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6xxxt1242130420.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,] 808744.0 1415151 1415099 1425600 1477037 1385488 NA [2,] 858954.5 1525367 1552763 1491526 1538880 1503718 NA [3,] 903734.5 1580142 1586340 1585432 1654831 1565068 NA [4,] 946601.5 1685082 1648583 1624362 1722882 1618132 NA [5,] 978630.0 1722061 1664774 1644081 1813765 1641395 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 863758.1 1507294 1542635 1524844 1570907 1512882 NA [2,] 943710.9 1652989 1630044 1646019 1738755 1617253 NA $out [1] 1156102 2264959 2304365 1993384 2165466 2125558 $group [1] 1 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(808744, 858954.5, 903734.5, 946601.5, 978630, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7hn3n1242130420.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,] 1400340 1415125 1407939 [2,] 1451769 1525290 1528746 [3,] 1458293 1562373 1541862 [4,] 1509215 1621445 1622343 [5,] 1564935 1635266 1656323 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1432091 1518516 1499172 [2,] 1484494 1606230 1584552 $out [1] 1321186 2001639 2145512 2133257 $group [1] 1 1 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1400340.16666667, 1451769.25, 1458292.83333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1pufn1242130420.ps tmp/1pufn1242130420.png") > system("convert tmp/2c4tg1242130420.ps tmp/2c4tg1242130420.png") > system("convert tmp/3c09c1242130420.ps tmp/3c09c1242130420.png") > system("convert tmp/4ldol1242130420.ps tmp/4ldol1242130420.png") > system("convert tmp/5k4fj1242130420.ps tmp/5k4fj1242130420.png") > system("convert tmp/6xxxt1242130420.ps tmp/6xxxt1242130420.png") > system("convert tmp/7hn3n1242130420.ps tmp/7hn3n1242130420.png") > > > proc.time() user system elapsed 1.233 0.963 1.566