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(286.1,307,358.1,341.8,378.8,375.2,295.6,362.7,409.6,336.8,389.1,389.3,355.9,542,648.4,452,582.4,506.5,555.5,530.4,609.4,543.9,616.2,634.6,541.7,549.8,627.6,797.4,689.8,1576.6,1572.1,1626.4,1972.4,1509.6,1584.9,1880,1324,1777.7,2172.4,1780.3,2134.9,1838.4,1557,1755.2,1702,1577.5,1485.9,2179.1,1740.9,1724.5,2328.1,1774.1,2224.2,1536.3,1521.2,2051.8,2483.1,1929.8,1808.6,2584.9,1997.9,1639.9,2379.1,1715,2750.9,1865.4,1647.4,2180.4,2593,2057.2,2635.8,2315.4,1863.6,2038,2235.8,2222.1,2636.9,2076.8,1935.5,2086.3,2470.9,1854.6,2041.3,2170.8,1905.5,2130.2,2791.2,2539.7,2661.3,1764.9,2176.9,2458.5,2179,2242.5,2089.6,2661.6,2112,2367.3,2543,2603.9,3146.7,1789.2,2114.8,2236.3,2288.1,2173.2,1877.7,2807.4,2357.4,2107.7,2856.8,2510.8,2875,2229.7,2055.1,2545.4,2775.1,2252.2,2091.7,2433) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 286.1 355.9 541.7 1324.0 1740.9 1997.9 1863.6 1905.5 2112.0 2357.4 NA [2,] 307.0 542.0 549.8 1777.7 1724.5 1639.9 2038.0 2130.2 2367.3 2107.7 NA [3,] 358.1 648.4 627.6 2172.4 2328.1 2379.1 2235.8 2791.2 2543.0 2856.8 NA [4,] 341.8 452.0 797.4 1780.3 1774.1 1715.0 2222.1 2539.7 2603.9 2510.8 NA [5,] 378.8 582.4 689.8 2134.9 2224.2 2750.9 2636.9 2661.3 3146.7 2875.0 NA [6,] 375.2 506.5 1576.6 1838.4 1536.3 1865.4 2076.8 1764.9 1789.2 2229.7 NA [7,] 295.6 555.5 1572.1 1557.0 1521.2 1647.4 1935.5 2176.9 2114.8 2055.1 NA [8,] 362.7 530.4 1626.4 1755.2 2051.8 2180.4 2086.3 2458.5 2236.3 2545.4 NA [9,] 409.6 609.4 1972.4 1702.0 2483.1 2593.0 2470.9 2179.0 2288.1 2775.1 NA [10,] 336.8 543.9 1509.6 1577.5 1929.8 2057.2 1854.6 2242.5 2173.2 2252.2 NA [11,] 389.1 616.2 1584.9 1485.9 1808.6 2635.8 2041.3 2089.6 1877.7 2091.7 NA [12,] 389.3 634.6 1880.0 2179.1 2584.9 2315.4 2170.8 2661.6 2807.4 2433.0 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 20.9 186.1 8.1 453.7 -16.4 -358.0 174.4 224.7 255.3 -249.7 [2,] 51.1 106.4 77.8 394.7 603.6 739.2 197.8 661.0 175.7 749.1 [3,] -16.3 -196.4 169.8 -392.1 -554.0 -664.1 -13.7 -251.5 60.9 -346.0 [4,] 37.0 130.4 -107.6 354.6 450.1 1035.9 414.8 121.6 542.8 364.2 [5,] -3.6 -75.9 886.8 -296.5 -687.9 -885.5 -560.1 -896.4 -1357.5 -645.3 [6,] -79.6 49.0 -4.5 -281.4 -15.1 -218.0 -141.3 412.0 325.6 -174.6 [7,] 67.1 -25.1 54.3 198.2 530.6 533.0 150.8 281.6 121.5 490.3 [8,] 46.9 79.0 346.0 -53.2 431.3 412.6 384.6 -279.5 51.8 229.7 [9,] -72.8 -65.5 -462.8 -124.5 -553.3 -535.8 -616.3 63.5 -114.9 -522.9 [10,] 52.3 72.3 75.3 -91.6 -121.2 578.6 186.7 -152.9 -295.5 -160.5 [11,] 0.2 18.4 295.1 693.2 776.3 -320.4 129.5 572.0 929.7 341.3 [12,] -33.4 -92.9 -556.0 -438.2 -587.0 -451.8 -265.3 -549.6 -450.0 NA [,11] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] 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/16p081257891021.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/2d79t1257891021.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/3yjgg1257891021.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/4hlcb1257891021.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,] 286.10 307.0 358.10 341.8 378.80 1536.30 1521.20 1626.40 609.40 [2,] 541.70 549.8 648.40 797.4 689.80 1536.30 1521.20 1626.40 1702.00 [3,] 1802.25 1751.1 2281.95 1777.2 2430.55 1777.05 1609.75 2069.05 2233.55 [4,] 1997.90 2107.7 2543.00 2510.8 2750.90 1865.40 2055.10 2236.30 2483.10 [5,] 2357.40 2367.3 2856.80 2603.9 3146.70 2229.70 2176.90 2545.40 2775.10 [,10] [,11] [,12] [1,] 543.9 616.20 1880.00 [2,] 1509.6 1485.90 1880.00 [3,] 1892.2 1843.15 2247.25 [4,] 2173.2 2089.60 2584.90 [5,] 2252.2 2635.80 2807.40 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1074.674 972.711 1335.332 921.117 1400.742 1612.619 1342.992 1764.320 [2,] 2529.826 2529.489 3228.568 2633.283 3460.358 1941.481 1876.508 2373.780 [,9] [,10] [,11] [,12] [1,] 1843.281 1560.639 1541.517 1895.054 [2,] 2623.819 2223.761 2144.783 2599.446 $out [1] 375.2 506.5 295.6 555.5 362.7 530.4 409.6 336.8 389.1 389.3 634.6 $group [1] 6 6 7 7 8 8 9 10 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(286.1, 541.7, 1802.25, 1997.9, 2357.4, 307, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/52bnq1257891021.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] [1,] -358.00 51.10 -664.10 -107.6 -1357.5 -281.40 -25.1 -279.50 -616.30 [2,] -16.40 106.40 -392.10 121.6 -885.5 -174.60 67.1 46.90 -535.80 [3,] 97.65 296.25 -223.95 359.4 -602.7 -47.35 174.5 154.35 -293.65 [4,] 224.70 661.00 -13.70 450.1 -75.9 49.00 490.3 384.60 -72.80 [5,] 453.70 749.10 169.80 542.8 886.8 325.60 533.0 431.30 63.50 [,10] [,11] [,12] [1,] -295.50 -320.4 -587.0 [2,] -152.90 18.4 -549.6 [3,] -19.65 318.2 -450.0 [4,] 75.30 693.2 -265.3 [5,] 186.70 929.7 -33.4 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -22.81317 19.14973 -413.01373 195.2683 -1007.2084 -159.06948 -36.94759 [2,] 218.11317 573.35027 -34.88627 523.5317 -198.1916 64.36948 385.94759 [,8] [,9] [,10] [,11] [,12] [1,] -14.37838 -524.98326 -133.66782 -18.95698 -599.7313 [2,] 323.07838 -62.31674 94.36782 655.35698 -300.2687 $out [1] 1035.9 412.0 578.6 $group [1] 4 6 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-358, -16.4000000000001, 97.65, 224.7, 453.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ryvu1257891021.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] [1,] 286.10 452.00 541.70 1324.00 1521.20 1639.90 1854.60 1764.90 1789.20 [2,] 321.90 518.45 658.70 1567.25 1732.70 1790.20 1986.75 2109.90 2113.40 [3,] 360.40 549.70 1540.85 1766.45 1869.20 2118.80 2081.55 2210.75 2262.20 [4,] 383.95 612.80 1605.65 1986.65 2276.15 2486.05 2228.95 2600.50 2573.45 [5,] 409.60 648.40 1972.40 2179.10 2584.90 2750.90 2470.90 2791.20 3146.70 [,10] [,11] [1,] 2055.10 NA [2,] 2168.70 NA [3,] 2395.20 NA [4,] 2660.25 NA [5,] 2875.00 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 332.0986 506.6663 1108.940 1575.159 1621.329 1801.418 1971.081 1986.984 [2,] 388.7014 592.7337 1972.760 1957.741 2117.071 2436.182 2192.019 2434.516 [,9] [,10] [,11] [1,] 2052.368 2171.001 NA [2,] 2472.032 2619.399 NA $out [1] 355.9 2636.9 $group [1] 2 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(286.1, 321.9, 360.4, 383.95, 409.6, 452, 518.45, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/79o721257891021.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,] 1448.500 1609.750 1456.300 [2,] 1549.505 1777.125 1717.613 [3,] 1667.895 1867.675 1785.881 [4,] 1921.155 2240.400 1915.125 [5,] 2008.090 2430.550 2124.825 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1498.383 1656.372 1695.794 [2,] 1837.407 2078.978 1875.968 $out [1] 1356.038 2249.812 $group [1] 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1448.5, 1549.505, 1667.895, 1921.155, 2008.09, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/16p081257891021.ps tmp/16p081257891021.png") > system("convert tmp/2d79t1257891021.ps tmp/2d79t1257891021.png") > system("convert tmp/3yjgg1257891021.ps tmp/3yjgg1257891021.png") > system("convert tmp/4hlcb1257891021.ps tmp/4hlcb1257891021.png") > system("convert tmp/52bnq1257891021.ps tmp/52bnq1257891021.png") > system("convert tmp/6ryvu1257891021.ps tmp/6ryvu1257891021.png") > system("convert tmp/79o721257891021.ps tmp/79o721257891021.png") > > > proc.time() user system elapsed 1.172 0.938 1.906