R version 2.7.2 (2008-08-25) 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(2173,2363,2126,1905,2121,1983,1734,2074,2049,2406,2558,2251,2059,2397,1747,1707,2319,1631,1627,1791,2034,1997,2169,2028,2253,2218,1855,2187,1852,1570,1851,1954,1828,2251,2277,2085,2282,2266,1878,2267,2069,1746,2299,2360,2214,2825,2355,2333,3016,2155,2172,2150,2533,2058,2160,2259,2498,2695,2799,2945,2930,2318,2540,2570,2669,2450,2842,3439,2677,2979,2257,2842,2546,2455,2293,2379,2478,2054,2272,2351,2271,2542,2304,2194,2722,2395,2146,1894,2548,2087,2063,2481,2476,2212,2834,2148,2598) > par1 = '60' > #'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] 97 > (np <- floor(n / par1)) [1] 1 > 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] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 > arr [,1] [,2] [1,] 2173 2930 [2,] 2363 2318 [3,] 2126 2540 [4,] 1905 2570 [5,] 2121 2669 [6,] 1983 2450 [7,] 1734 2842 [8,] 2074 3439 [9,] 2049 2677 [10,] 2406 2979 [11,] 2558 2257 [12,] 2251 2842 [13,] 2059 2546 [14,] 2397 2455 [15,] 1747 2293 [16,] 1707 2379 [17,] 2319 2478 [18,] 1631 2054 [19,] 1627 2272 [20,] 1791 2351 [21,] 2034 2271 [22,] 1997 2542 [23,] 2169 2304 [24,] 2028 2194 [25,] 2253 2722 [26,] 2218 2395 [27,] 1855 2146 [28,] 2187 1894 [29,] 1852 2548 [30,] 1570 2087 [31,] 1851 2063 [32,] 1954 2481 [33,] 1828 2476 [34,] 2251 2212 [35,] 2277 2834 [36,] 2085 2148 [37,] 2282 2598 [38,] 2266 NA [39,] 1878 NA [40,] 2267 NA [41,] 2069 NA [42,] 1746 NA [43,] 2299 NA [44,] 2360 NA [45,] 2214 NA [46,] 2825 NA [47,] 2355 NA [48,] 2333 NA [49,] 3016 NA [50,] 2155 NA [51,] 2172 NA [52,] 2150 NA [53,] 2533 NA [54,] 2058 NA [55,] 2160 NA [56,] 2259 NA [57,] 2498 NA [58,] 2695 NA [59,] 2799 NA [60,] 2945 NA > darr [,1] [,2] [1,] 190 -612 [2,] -237 222 [3,] -221 30 [4,] 216 99 [5,] -138 -219 [6,] -249 392 [7,] 340 597 [8,] -25 -762 [9,] 357 302 [10,] 152 -722 [11,] -307 585 [12,] -192 -296 [13,] 338 -91 [14,] -650 -162 [15,] -40 86 [16,] 612 99 [17,] -688 -424 [18,] -4 218 [19,] 164 79 [20,] 243 -80 [21,] -37 271 [22,] 172 -238 [23,] -141 -110 [24,] 225 528 [25,] -35 -327 [26,] -363 -249 [27,] 332 -252 [28,] -335 654 [29,] -282 -461 [30,] 281 -24 [31,] 103 418 [32,] -126 -5 [33,] 423 -264 [34,] 26 622 [35,] -192 -686 [36,] 197 450 [37,] -16 NA [38,] -388 NA [39,] 389 NA [40,] -198 NA [41,] -323 NA [42,] 553 NA [43,] 61 NA [44,] -146 NA [45,] 611 NA [46,] -470 NA [47,] -22 NA [48,] 683 NA [49,] -861 NA [50,] 17 NA [51,] -22 NA [52,] 383 NA [53,] -475 NA [54,] 102 NA [55,] 99 NA [56,] 239 NA [57,] 197 NA [58,] 104 NA [59,] 146 NA [60,] -15 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/1xsx11225664577.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/29g3c1225664577.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/32w811225664577.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/43lpy1225664577.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,] 2173.0 2318.0 2126 1905.0 2121 1983.0 1734 2074.0 2049 2406.0 2257.0 [2,] 2173.0 2318.0 2126 1905.0 2121 1983.0 1734 2074.0 2049 2406.0 2257.0 [3,] 2551.5 2340.5 2333 2237.5 2395 2216.5 2288 2756.5 2363 2692.5 2407.5 [4,] 2930.0 2363.0 2540 2570.0 2669 2450.0 2842 3439.0 2677 2979.0 2558.0 [5,] 2930.0 2363.0 2540 2570.0 2669 2450.0 2842 3439.0 2677 2979.0 2558.0 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 2251.0 2059.0 2397 1747 1707 2319.0 1631.0 1627.0 1791 2034.0 1997.0 [2,] 2251.0 2059.0 2397 1747 1707 2319.0 1631.0 1627.0 1791 2034.0 1997.0 [3,] 2546.5 2302.5 2426 2020 2043 2398.5 1842.5 1949.5 2071 2152.5 2269.5 [4,] 2842.0 2546.0 2455 2293 2379 2478.0 2054.0 2272.0 2351 2271.0 2542.0 [5,] 2842.0 2546.0 2455 2293 2379 2478.0 2054.0 2272.0 2351 2271.0 2542.0 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 2169.0 2028 2253.0 2218.0 1855.0 1894.0 1852 1570.0 1851 1954.0 1828 [2,] 2169.0 2028 2253.0 2218.0 1855.0 1894.0 1852 1570.0 1851 1954.0 1828 [3,] 2236.5 2111 2487.5 2306.5 2000.5 2040.5 2200 1828.5 1957 2217.5 2152 [4,] 2304.0 2194 2722.0 2395.0 2146.0 2187.0 2548 2087.0 2063 2481.0 2476 [5,] 2304.0 2194 2722.0 2395.0 2146.0 2187.0 2548 2087.0 2063 2481.0 2476 [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [1,] 2212.0 2277.0 2085.0 2282 2266 1878 2267 2069 1746 2299 2360 2214 [2,] 2212.0 2277.0 2085.0 2282 2266 1878 2267 2069 1746 2299 2360 2214 [3,] 2231.5 2555.5 2116.5 2440 2266 1878 2267 2069 1746 2299 2360 2214 [4,] 2251.0 2834.0 2148.0 2598 2266 1878 2267 2069 1746 2299 2360 2214 [5,] 2251.0 2834.0 2148.0 2598 2266 1878 2267 2069 1746 2299 2360 2214 [,46] [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [1,] 2825 2355 2333 3016 2155 2172 2150 2533 2058 2160 2259 2498 [2,] 2825 2355 2333 3016 2155 2172 2150 2533 2058 2160 2259 2498 [3,] 2825 2355 2333 3016 2155 2172 2150 2533 2058 2160 2259 2498 [4,] 2825 2355 2333 3016 2155 2172 2150 2533 2058 2160 2259 2498 [5,] 2825 2355 2333 3016 2155 2172 2150 2533 2058 2160 2259 2498 [,58] [,59] [,60] [1,] 2695 2799 2945 [2,] 2695 2799 2945 [3,] 2695 2799 2945 [4,] 2695 2799 2945 [5,] 2695 2799 2945 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1705.758 2290.225 1870.467 1494.543 1782.759 1694.754 1050.111 1231.483 [2,] 3397.242 2390.775 2795.533 2980.457 3007.241 2738.246 3525.889 4281.517 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 1661.380 2052.328 2071.214 1886.218 1758.410 2361.201 1409.993 1292.222 [2,] 3064.620 3332.672 2743.786 3206.782 2846.590 2490.799 2630.007 2793.778 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 2220.861 1369.912 1228.887 1445.352 1887.717 1660.610 2085.674 1925.54 [2,] 2576.139 2315.088 2670.113 2696.648 2417.283 2878.390 2387.326 2296.46 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 1963.520 2108.751 1675.386 1713.152 1422.409 1250.893 1720.148 1628.720 [2,] 3011.480 2504.249 2325.614 2367.848 2977.591 2406.107 2193.852 2806.280 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [1,] 1428.036 2187.928 1933.204 2046.115 2086.956 2266 1878 2267 2069 1746 [2,] 2875.964 2275.072 3177.796 2186.885 2793.044 2266 1878 2267 2069 1746 [,43] [,44] [,45] [,46] [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [1,] 2299 2360 2214 2825 2355 2333 3016 2155 2172 2150 2533 2058 [2,] 2299 2360 2214 2825 2355 2333 3016 2155 2172 2150 2533 2058 [,55] [,56] [,57] [,58] [,59] [,60] [1,] 2160 2259 2498 2695 2799 2945 [2,] 2160 2259 2498 2695 2799 2945 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" Warning message: In bxp(list(stats = c(2173, 2173, 2551.5, 2930, 2930, 2318, 2318, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5voph1225664577.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] [,11] [1,] -612 -237.0 -221.0 99.0 -219.0 -249.0 340.0 -762.0 302.0 -722 -307 [2,] -612 -237.0 -221.0 99.0 -219.0 -249.0 340.0 -762.0 302.0 -722 -307 [3,] -211 -7.5 -95.5 157.5 -178.5 71.5 468.5 -393.5 329.5 -285 139 [4,] 190 222.0 30.0 216.0 -138.0 392.0 597.0 -25.0 357.0 152 585 [5,] 190 222.0 30.0 216.0 -138.0 392.0 597.0 -25.0 357.0 152 585 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [1,] -296 -91.0 -650 -40 99.0 -688 -4 79.0 -80.0 -37 -238 -141.0 [2,] -296 -91.0 -650 -40 99.0 -688 -4 79.0 -80.0 -37 -238 -141.0 [3,] -244 123.5 -406 23 355.5 -556 107 121.5 81.5 117 -33 -125.5 [4,] -192 338.0 -162 86 612.0 -424 218 164.0 243.0 271 172 -110.0 [5,] -192 338.0 -162 86 612.0 -424 218 164.0 243.0 271 172 -110.0 [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [1,] 225.0 -327 -363 -252 -335.0 -461.0 -24.0 103.0 -126.0 -264.0 26 [2,] 225.0 -327 -363 -252 -335.0 -461.0 -24.0 103.0 -126.0 -264.0 26 [3,] 376.5 -181 -306 40 159.5 -371.5 128.5 260.5 -65.5 79.5 324 [4,] 528.0 -35 -249 332 654.0 -282.0 281.0 418.0 -5.0 423.0 622 [5,] 528.0 -35 -249 332 654.0 -282.0 281.0 418.0 -5.0 423.0 622 [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] -686 197.0 -16 -388 389 -198 -323 553 61 -146 611 -470 [2,] -686 197.0 -16 -388 389 -198 -323 553 61 -146 611 -470 [3,] -439 323.5 -16 -388 389 -198 -323 553 61 -146 611 -470 [4,] -192 450.0 -16 -388 389 -198 -323 553 61 -146 611 -470 [5,] -192 450.0 -16 -388 389 -198 -323 553 61 -146 611 -470 [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [1,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [2,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [3,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [4,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [5,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [,59] [,60] [1,] 146 -15 [2,] 146 -15 [3,] 146 -15 [4,] 146 -15 [5,] 146 -15 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 1 [39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1107.0174 -520.308 -375.9244 26.78424 -268.99553 -644.6436 181.3722 [2,] 685.0174 505.308 184.9244 288.21576 -88.00447 787.6436 755.6278 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] -1216.8976 268.0524 -1261.4579 -857.568 -360.1918 -355.7911 -951.2076 [2,] 429.8976 390.9476 691.4579 1135.568 -127.8082 602.7911 139.2076 [,15] [,16] [,17] [,18] [,19] [,20] [,21] [1,] -117.7708 -217.6383 -850.9484 -141.0248 26.53556 -279.3649 -227.1064 [2,] 163.7708 928.6383 -261.0516 355.0248 216.46444 442.3649 461.1064 [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] -491.0638 -160.13409 37.9797 -507.2308 -433.3641 -612.4616 -945.4392 [2,] 425.0638 -90.86591 715.0203 145.2308 -178.6359 692.4616 1264.4392 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [1,] -571.4839 -212.2548 -91.42704 -200.68467 -688.0361 -341.8683 -990.9110 [2,] -171.5161 469.2548 612.42704 69.68467 847.0361 989.8683 112.9110 [,36] [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [1,] 40.84114 -16 -388 389 -198 -323 553 61 -146 611 -470 [2,] 606.15886 -16 -388 389 -198 -323 553 61 -146 611 -470 [,47] [,48] [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [1,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [2,] -22 683 -861 17 -22 383 -475 102 99 239 197 104 [,59] [,60] [1,] 146 -15 [2,] 146 -15 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" "43" "44" "45" [46] "46" "47" "48" "49" "50" "51" "52" "53" "54" "55" "56" "57" "58" "59" "60" Warning message: In bxp(list(stats = c(-612, -612, -211, 190, 190, -237, -237, -7.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6rqiq1225664577.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] [1,] 1570.0 1894 [2,] 1968.5 2271 [3,] 2164.5 2455 [4,] 2309.0 2598 [5,] 2799.0 2979 $n [1] 60 37 $conf [,1] [,2] [1,] 2095.046 2370.062 [2,] 2233.954 2539.938 $out [1] 2825 3016 2945 3439 $group [1] 1 1 1 2 $names [1] "1" NA > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7e12d1225664577.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,] 1746.00 1746.00 1746.00 [2,] 2133.25 2133.25 2133.25 [3,] 2266.50 2266.50 2266.50 [4,] 2416.75 2416.75 2416.75 [5,] 2825.00 2825.00 2825.00 $n [1] 60 60 60 $conf [,1] [,2] [,3] [1,] 2208.672 2208.672 2208.672 [2,] 2324.328 2324.328 2324.328 $out [1] 3016 2945 3016 2945 3016 2945 $group [1] 1 1 2 2 3 3 $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/1xsx11225664577.ps tmp/1xsx11225664577.png") > system("convert tmp/29g3c1225664577.ps tmp/29g3c1225664577.png") > system("convert tmp/32w811225664577.ps tmp/32w811225664577.png") > system("convert tmp/43lpy1225664577.ps tmp/43lpy1225664577.png") > system("convert tmp/5voph1225664577.ps tmp/5voph1225664577.png") > system("convert tmp/6rqiq1225664577.ps tmp/6rqiq1225664577.png") > system("convert tmp/7e12d1225664577.ps tmp/7e12d1225664577.png") > > > proc.time() user system elapsed 1.566 0.993 1.959