R version 2.6.2 (2008-02-08) 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. Natural language support but running in an English locale 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(464,675,703,887,1139,1077,1318,1260,1120,963,996,960,530,883,894,1045,1199,1287,1565,1577,1076,918,1008,1063,544,635,804,980,1018,1064,1404,1286,1104,999,996,1015,615,722,832,977,1270,1437,1520,1708,1151,934,1159,1209,699,830,996,1124,1458,1270,1753,2258,1208,1241,1265,1828,809,997,1164,1205,1538,1513,1378,2083,1357,1536,1526,1376,779,1005,1193,1522,1539,1546,2116,2326,1596,1356,1553,1613,814,1150,1225,1691,1759,1754,2100,2062,2012,1897,1964,2186,966,1549,1538,1612,2078,2137,2907,2249,1883,1739,1828,1868,1138,1430,1809,1763,2200,2067,2503,2141,2103,1972,2181,2344,970,1199,1718,1683,2025,2051,2439,2353,2230,1852,2147,2286,1007,1665,1642,1518,1831,2207,2822,2393,2306,1785,2047,2171,1212,1335,2011,1860,1954,2152,2835,2224,2182,1992,2389,2724,891,1247,2017,2257,2255,2255,3057,3330,1896,2096,2374,2535,1041,1728,2201,2455,2204,2660,3670,2665,2639,2226,2586,2684,1185,1749,2459,2618,2585,3310,3923) > 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] 187 > (np <- floor(n / par1)) [1] 15 > 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] 16 16 16 16 16 16 16 15 15 15 15 15 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 464 530 544 615 699 809 779 814 966 1138 970 1007 1212 [2,] 675 883 635 722 830 997 1005 1150 1549 1430 1199 1665 1335 [3,] 703 894 804 832 996 1164 1193 1225 1538 1809 1718 1642 2011 [4,] 887 1045 980 977 1124 1205 1522 1691 1612 1763 1683 1518 1860 [5,] 1139 1199 1018 1270 1458 1538 1539 1759 2078 2200 2025 1831 1954 [6,] 1077 1287 1064 1437 1270 1513 1546 1754 2137 2067 2051 2207 2152 [7,] 1318 1565 1404 1520 1753 1378 2116 2100 2907 2503 2439 2822 2835 [8,] 1260 1577 1286 1708 2258 2083 2326 2062 2249 2141 2353 2393 2224 [9,] 1120 1076 1104 1151 1208 1357 1596 2012 1883 2103 2230 2306 2182 [10,] 963 918 999 934 1241 1536 1356 1897 1739 1972 1852 1785 1992 [11,] 996 1008 996 1159 1265 1526 1553 1964 1828 2181 2147 2047 2389 [12,] 960 1063 1015 1209 1828 1376 1613 2186 1868 2344 2286 2171 2724 [,14] [,15] [,16] [1,] 891 1041 1185 [2,] 1247 1728 1749 [3,] 2017 2201 2459 [4,] 2257 2455 2618 [5,] 2255 2204 2585 [6,] 2255 2660 3310 [7,] 3057 3670 3923 [8,] 3330 2665 NA [9,] 1896 2639 NA [10,] 2096 2226 NA [11,] 2374 2586 NA [12,] 2535 2684 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/1xzrk1210602790.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/28gua1210602790.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/3al6v1210602790.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/4xc2c1210602791.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] [,10] [,11] [1,] 464.0 635.0 703.0 887.0 1018 1064.0 1318.0 1260.0 1076.0 918.0 996 [2,] 657.0 856.5 945.0 1084.5 1364 1362.0 1542.5 1885.0 1179.5 1120.0 1212 [3,] 852.5 1174.5 1381.5 1567.0 1795 1902.5 2277.5 2224.0 1883.0 1739.0 1828 [4,] 1024.0 1489.5 1910.0 1811.5 2139 2179.5 2871.0 2339.5 2142.5 1934.5 2164 [5,] 1212.0 1749.0 2459.0 2618.0 2585 3310.0 3923.0 2665.0 2639.0 2226.0 2586 [,12] [1,] 960.0 [2,] 1292.5 [3,] 1868.0 [4,] 2315.0 [5,] 2724.0 $n [1] 16 16 16 16 16 16 16 15 15 15 15 15 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 707.535 924.465 1000.325 1279.835 1488.875 1579.588 1752.742 2038.585 [2,] 997.465 1424.535 1762.675 1854.165 2101.125 2225.412 2802.258 2409.415 [,9] [,10] [,11] [,12] [1,] 1490.14 1406.721 1439.628 1450.867 [2,] 2275.86 2071.279 2216.372 2285.133 $out [1] 3330 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(464, 657, 852.5, 1024, 1212, 635, 856.5, 1174.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5ajxo1210602791.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] [,11] [1,] 464.0 530 635 615.0 699.0 809.0 1005.0 814 966.0 1430.0 970.0 [2,] 795.0 906 892 883.0 1060.0 1184.5 1274.5 1458 1580.5 1786.0 1700.5 [3,] 979.5 1054 1007 1155.0 1253.0 1377.0 1542.5 1828 1848.0 2085.0 2038.0 [4,] 1129.5 1243 1084 1353.5 1605.5 1531.0 1604.5 2037 2107.5 2190.5 2258.0 [5,] 1318.0 1577 1286 1708.0 2258.0 1538.0 1613.0 2186 2249.0 2503.0 2439.0 [,12] [,13] [,14] [,15] [,16] [1,] 1007.0 1335.0 1247.0 1728.0 1185 [2,] 1653.5 1907.0 1956.5 2202.5 2104 [3,] 1939.0 2081.5 2255.0 2520.5 2585 [4,] 2256.5 2306.5 2454.5 2662.5 2964 [5,] 2822.0 2835.0 3057.0 2684.0 3923 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 826.9323 900.292 919.4275 940.4018 1004.194 1218.959 1391.985 1563.914 [2,] 1132.0677 1207.708 1094.5725 1369.5982 1501.806 1535.041 1693.015 2092.086 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1607.632 1900.505 1783.721 1663.968 1899.285 2027.859 2310.691 2071.422 [2,] 2088.368 2269.495 2292.279 2214.032 2263.715 2482.141 2730.309 3098.578 $out [1] 544 1404 2083 779 2116 2326 2907 1138 1212 891 3330 1041 3670 $group [1] 3 3 6 7 7 7 9 10 13 14 14 15 15 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] NA Warning message: In bxp(list(stats = c(464, 795, 979.5, 1129.5, 1318, 530, 906, 1054, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68ajy1210602791.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,] 376.2936 748.0 337.500 [2,] 462.2587 1435.5 636.500 [3,] 530.4869 1660.5 790.375 [4,] 583.6115 1917.0 957.500 [5,] 619.8570 2605.0 1299.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 475.1371 1440.885 643.9647 [2,] 585.8368 1880.115 936.7853 $out [1] 237.7253 824.3957 $group [1] 1 1 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(376.293585515353, 462.258662562769, 530.48693759086, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1xzrk1210602790.ps tmp/1xzrk1210602790.png") > system("convert tmp/28gua1210602790.ps tmp/28gua1210602790.png") > system("convert tmp/3al6v1210602790.ps tmp/3al6v1210602790.png") > system("convert tmp/4xc2c1210602791.ps tmp/4xc2c1210602791.png") > system("convert tmp/5ajxo1210602791.ps tmp/5ajxo1210602791.png") > system("convert tmp/68ajy1210602791.ps tmp/68ajy1210602791.png") > > > proc.time() user system elapsed 2.40 1.35 2.64