R version 2.8.0 (2008-10-20) 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(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/freestat/rcomp/tmp/1f2o31275505439.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/freestat/rcomp/tmp/2quno1275505439.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/freestat/rcomp/tmp/3quno1275505439.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/freestat/rcomp/tmp/4quno1275505439.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/freestat/rcomp/tmp/5jlmr1275505439.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/freestat/rcomp/tmp/6jlmr1275505439.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 > > try(system("convert tmp/1f2o31275505439.ps tmp/1f2o31275505439.png",intern=TRUE)) character(0) > try(system("convert tmp/2quno1275505439.ps tmp/2quno1275505439.png",intern=TRUE)) character(0) > try(system("convert tmp/3quno1275505439.ps tmp/3quno1275505439.png",intern=TRUE)) character(0) > try(system("convert tmp/4quno1275505439.ps tmp/4quno1275505439.png",intern=TRUE)) character(0) > try(system("convert tmp/5jlmr1275505439.ps tmp/5jlmr1275505439.png",intern=TRUE)) character(0) > try(system("convert tmp/6jlmr1275505439.ps tmp/6jlmr1275505439.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.641 1.252 1.833