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 = '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] 97 > (np <- floor(n / par1)) [1] 8 > 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] 9 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 2173 2059 2253 2282 3016 2930 2546 2722 2598 [2,] 2363 2397 2218 2266 2155 2318 2455 2395 NA [3,] 2126 1747 1855 1878 2172 2540 2293 2146 NA [4,] 1905 1707 2187 2267 2150 2570 2379 1894 NA [5,] 2121 2319 1852 2069 2533 2669 2478 2548 NA [6,] 1983 1631 1570 1746 2058 2450 2054 2087 NA [7,] 1734 1627 1851 2299 2160 2842 2272 2063 NA [8,] 2074 1791 1954 2360 2259 3439 2351 2481 NA [9,] 2049 2034 1828 2214 2498 2677 2271 2476 NA [10,] 2406 1997 2251 2825 2695 2979 2542 2212 NA [11,] 2558 2169 2277 2355 2799 2257 2304 2834 NA [12,] 2251 2028 2085 2333 2945 2842 2194 2148 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 190 338 -35 -16 -861 -612 -91 -327 NA [2,] -237 -650 -363 -388 17 222 -162 -249 NA [3,] -221 -40 332 389 -22 30 86 -252 NA [4,] 216 612 -335 -198 383 99 99 654 NA [5,] -138 -688 -282 -323 -475 -219 -424 -461 NA [6,] -249 -4 281 553 102 392 218 -24 NA [7,] 340 164 103 61 99 597 79 418 NA [8,] -25 243 -126 -146 239 -762 -80 -5 NA [9,] 357 -37 423 611 197 302 271 -264 NA [10,] 152 172 26 -470 104 -722 -238 622 NA [11,] -307 -141 -192 -22 146 585 -110 -686 NA [12,] -192 225 197 683 -15 -296 528 450 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/1egnq1225755836.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/2w5uu1225755836.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/3vtvn1225755836.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/4hz371225755836.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,] 2059 2155.0 1747.0 1707.0 1852.0 1570.0 1627.0 1791.0 1828.0 1997.0 2169.0 [2,] 2253 2242.0 1866.5 1899.5 2095.0 1688.5 1792.5 2014.0 2041.5 2231.5 2267.0 [3,] 2546 2340.5 2136.0 2168.5 2398.5 2018.5 2111.5 2305.0 2242.5 2474.0 2329.5 [4,] 2722 2396.0 2232.5 2323.0 2540.5 2072.5 2285.5 2420.5 2487.0 2760.0 2678.5 [5,] 3016 2455.0 2540.0 2570.0 2669.0 2450.0 2842.0 2481.0 2677.0 2979.0 2834.0 [,12] [1,] 2028.0 [2,] 2116.5 [3,] 2222.5 [4,] 2587.5 [5,] 2945.0 $n [1] 9 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2298.993 2254.473 1931.547 1931.927 2149.637 1803.992 1836.103 2077.923 [2,] 2793.007 2426.527 2340.453 2405.073 2647.363 2233.008 2386.897 2532.077 [,9] [,10] [,11] [,12] [1,] 1993.637 2178.772 2099.63 1959.393 [2,] 2491.363 2769.228 2559.37 2485.607 $out [1] 3439 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2059, 2253, 2546, 2722, 3016, 2155, 2242, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5xazc1225755836.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,] -861.0 -650.0 -252.0 -335.0 -688.0 -249.0 61.0 -146.0 -264.0 -722 -686.0 [2,] -469.5 -375.5 -130.5 -49.5 -468.0 -14.0 89.0 -136.0 80.0 -354 -249.5 [3,] -63.0 -243.0 4.0 157.5 -373.5 160.0 133.5 -52.5 286.5 65 -125.5 [4,] 87.0 -72.5 209.0 497.5 -250.5 336.5 379.0 117.0 390.0 162 62.0 [5,] 338.0 222.0 389.0 654.0 -138.0 553.0 597.0 243.0 611.0 622 146.0 [,12] [1,] -296.0 [2,] -103.5 [3,] 211.0 [4,] 489.0 [5,] 683.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -373.8689 -412.26015 -185.6496 -148.0621 -494.9986 -35.79433 -28.49816 [2,] 247.8689 -73.73985 193.6496 463.0621 -252.0014 355.79433 295.49816 [,8] [,9] [,10] [,11] [,12] [1,] -193.82943 113.3295 -223.245 -299.50837 -119.979 [2,] 88.82943 459.6705 353.245 48.50837 541.979 $out [1] -762 585 $group [1] 8 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-861, -469.5, -63, 87, 338, -650, -375.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/60zqv1225755836.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,] 1734.0 1627.0 1570.0 1878.0 2058.0 2257 2054.0 1894.0 2598 [2,] 2016.0 1727.0 1851.5 2141.5 2157.5 2495 2271.5 2116.5 2598 [3,] 2123.5 2012.5 2019.5 2274.5 2378.5 2673 2327.5 2303.5 2598 [4,] 2307.0 2114.0 2234.5 2344.0 2747.0 2886 2466.5 2514.5 2598 [5,] 2558.0 2397.0 2277.0 2360.0 3016.0 3439 2546.0 2834.0 2598 $n [1] 12 12 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1990.773 1835.987 1844.811 2182.138 2109.625 2494.662 2238.559 2121.970 [2,] 2256.227 2189.013 2194.189 2366.862 2647.375 2851.338 2416.441 2485.030 [,9] [1,] 2598 [2,] 2598 $out [1] 1746 2825 $group [1] 4 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(1734, 2016, 2123.5, 2307, 2558, 1627, 1727, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7aj2m1225755836.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,] 1947.375 2018.50 1891.250 [2,] 2119.188 2152.25 2074.438 [3,] 2322.250 2273.75 2279.812 [4,] 2398.688 2369.50 2384.938 [5,] 2508.778 2546.00 2487.500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2194.768 2174.661 2138.191 [2,] 2449.732 2372.839 2421.434 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(1947.375, 2119.1875, 2322.25, 2398.6875, 2508.77777777778, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1egnq1225755836.ps tmp/1egnq1225755836.png") > system("convert tmp/2w5uu1225755836.ps tmp/2w5uu1225755836.png") > system("convert tmp/3vtvn1225755836.ps tmp/3vtvn1225755836.png") > system("convert tmp/4hz371225755836.ps tmp/4hz371225755836.png") > system("convert tmp/5xazc1225755836.ps tmp/5xazc1225755836.png") > system("convert tmp/60zqv1225755836.ps tmp/60zqv1225755836.png") > system("convert tmp/7aj2m1225755836.ps tmp/7aj2m1225755836.png") > > > proc.time() user system elapsed 1.413 0.964 1.659