R version 2.11.1 (2010-05-31) Copyright (C) 2010 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(12.008,9.169,8.788,8.417,8.247,8.197,8.236,8.253,7.733,8.366,8.626,8.863,10.102,8.463,9.114,8.563,8.872,8.301,8.301,8.278,7.736,7.973,8.268,9.476,11.100,8.962,9.173,8.738,8.459,8.078,8.411,8.291,7.810,8.616,8.312,9.692,9.911,8.915,9.452,9.112,8.472,8.230,8.384,8.625,8.221,8.649,8.625,10.443,10.357,8.586,8.892,8.329,8.101,7.922,8.120,7.838,7.735,8.406,8.209,9.451,10.041,9.411,10.405,8.467,8.464,8.102,7.627,7.513,7.510,8.291,8.064,9.383,9.706,8.579,9.474,8.318,8.213,8.059,9.111,7.708,7.680,8.014,8.007,8.718,9.486,9.113,9.025,8.476,7.952,7.759,7.835,7.600,7.651,8.319,8.812,8.630) > par1 = '8' > #'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] 96 > (np <- floor(n / par1)) [1] 12 > 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] 12 12 12 12 12 12 12 12 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 12.008 7.733 8.872 11.100 7.810 8.472 10.357 7.735 8.464 9.706 7.680 [2,] 9.169 8.366 8.301 8.962 8.616 8.230 8.586 8.406 8.102 8.579 8.014 [3,] 8.788 8.626 8.301 9.173 8.312 8.384 8.892 8.209 7.627 9.474 8.007 [4,] 8.417 8.863 8.278 8.738 9.692 8.625 8.329 9.451 7.513 8.318 8.718 [5,] 8.247 10.102 7.736 8.459 9.911 8.221 8.101 10.041 7.510 8.213 9.486 [6,] 8.197 8.463 7.973 8.078 8.915 8.649 7.922 9.411 8.291 8.059 9.113 [7,] 8.236 9.114 8.268 8.411 9.452 8.625 8.120 10.405 8.064 9.111 9.025 [8,] 8.253 8.563 9.476 8.291 9.112 10.443 7.838 8.467 9.383 7.708 8.476 [,12] [,13] [1,] 7.952 NA [2,] 7.759 NA [3,] 7.835 NA [4,] 7.600 NA [5,] 7.651 NA [6,] 8.319 NA [7,] 8.812 NA [8,] 8.630 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -2.839 0.633 -0.571 -2.138 0.806 -0.242 -1.771 0.671 -0.362 -1.127 [2,] -0.381 0.260 0.000 0.211 -0.304 0.154 0.306 -0.197 -0.475 0.895 [3,] -0.371 0.237 -0.023 -0.435 1.380 0.241 -0.563 1.242 -0.114 -1.156 [4,] -0.170 1.239 -0.542 -0.279 0.219 -0.404 -0.228 0.590 -0.003 -0.105 [5,] -0.050 -1.639 0.237 -0.381 -0.996 0.428 -0.179 -0.630 0.781 -0.154 [6,] 0.039 0.651 0.295 0.333 0.537 -0.024 0.198 0.994 -0.227 1.052 [7,] 0.017 -0.551 1.208 -0.120 -0.340 1.818 -0.282 -1.938 1.319 -1.403 [8,] -0.520 0.309 1.624 -0.481 -0.640 -0.086 -0.103 -0.003 0.323 -0.028 [,11] [,12] [,13] [1,] 0.334 -0.193 NA [2,] -0.007 0.076 NA [3,] 0.711 -0.235 NA [4,] 0.768 0.051 NA [5,] -0.373 0.668 NA [6,] -0.088 0.493 NA [7,] -0.549 -0.182 NA [8,] -0.524 NA 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/rcomp/tmp/12ocj1290624860.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/rcomp/tmp/22ocj1290624860.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/rcomp/tmp/3vxu41290624860.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/rcomp/tmp/4vxu41290624860.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] [1,] 7.6800 7.759 7.627 7.6000 7.5100 7.9220 8.0640 7.7080 [2,] 7.7725 8.166 8.108 8.2980 7.9185 8.0685 8.2520 8.2720 [3,] 8.4680 8.386 8.348 8.5210 8.2340 8.3050 8.7185 8.5195 [4,] 10.0315 8.601 8.840 8.8005 9.6985 8.7820 9.1125 9.2475 [5,] 12.0080 9.169 9.474 9.4510 10.1020 9.4110 9.4520 10.4430 $n [1] 12 12 12 12 12 12 12 12 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 7.437655 8.187594 8.01413 8.291806 7.42213 7.979568 8.32602 8.074568 [2,] 9.498345 8.584406 8.68187 8.750194 9.04587 8.630432 9.11098 8.964432 $out [1] 9.692 7.513 10.405 $group [1] 4 4 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" Warning message: In bxp(list(stats = c(7.68, 7.7725, 8.468, 10.0315, 12.008, 7.759, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5vxu41290624860.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] [1,] -2.8390 -0.4750 -1.1560 -0.5420 -1.6390 -0.2270 -1.9380 -0.6400 [2,] -1.4490 -0.2505 -0.4030 -0.2535 -0.5055 0.0075 -0.5500 -0.5005 [3,] -0.3020 0.0380 -0.0685 -0.0540 -0.1665 0.3140 -0.2320 -0.0860 [4,] 0.4835 0.2355 0.4760 0.4045 0.3325 0.5940 0.6125 0.1530 [5,] 0.8060 0.8950 1.3800 1.2390 0.7810 1.0520 1.8180 0.3230 $n [1] 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1.1834262 -0.1836679 -0.4694178 -0.3541182 -0.5487174 0.04649341 [2,] 0.5794262 0.2596679 0.3324178 0.2461182 0.2157174 0.58150659 [,7] [,8] [1,] -0.7622241 -0.3973195 [2,] 0.2982241 0.2253195 $out [1] 1.624 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" Warning message: In bxp(list(stats = c(-2.839, -1.449, -0.302000000000000, 0.4835, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/66pb71290624860.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] [1,] 8.1970 7.7330 7.7360 8.0780 7.8100 8.2210 7.8380 7.7350 7.5100 7.7080 [2,] 8.2415 8.4145 8.1205 8.3510 8.4640 8.3070 8.0115 8.3075 7.5700 8.1360 [3,] 8.3350 8.5945 8.2895 8.5985 9.0135 8.5485 8.2245 8.9390 8.0830 8.4485 [4,] 8.9785 8.9885 8.5865 9.0675 9.5720 8.6370 8.7390 9.7460 8.3775 9.2925 [5,] 9.1690 9.1140 8.8720 9.1730 9.9110 8.6490 8.8920 10.4050 9.3830 9.7060 [,11] [,12] [,13] [1,] 7.6800 7.6000 NA [2,] 8.0105 7.7050 NA [3,] 8.5970 7.8935 NA [4,] 9.0690 8.4745 NA [5,] 9.4860 8.8120 NA $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 7.923301 8.273855 8.029186 8.198253 8.394555 8.364157 7.818108 8.135433 [2,] 8.746699 8.915145 8.549814 8.998747 9.632445 8.732843 8.630892 9.742567 [,9] [,10] [,11] [,12] [,13] [1,] 7.631919 7.802462 8.005707 7.463646 NA [2,] 8.534081 9.094538 9.188293 8.323354 NA $out [1] 12.008 10.102 9.476 11.100 10.443 10.357 $group [1] 1 2 3 4 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" NA Warning message: In bxp(list(stats = c(8.197, 8.2415, 8.335, 8.9785, 9.169, 7.733, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/76pb71290624860.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,] 8.424167 8.23400 8.394375 [2,] 8.459083 8.32650 8.441000 [3,] 8.592500 8.42700 8.612250 [4,] 8.761792 8.52025 8.765812 [5,] 8.990750 8.71850 8.830000 $n [1] 8 8 8 $conf [,1] [,2] [,3] [1,] 8.423403 8.318768 8.430805 [2,] 8.761597 8.535232 8.793695 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(8.42416666666667, 8.45908333333333, 8.5925, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/12ocj1290624860.ps tmp/12ocj1290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/22ocj1290624860.ps tmp/22ocj1290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/3vxu41290624860.ps tmp/3vxu41290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/4vxu41290624860.ps tmp/4vxu41290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/5vxu41290624860.ps tmp/5vxu41290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/66pb71290624860.ps tmp/66pb71290624860.png",intern=TRUE)) character(0) > try(system("convert tmp/76pb71290624860.ps tmp/76pb71290624860.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.010 1.300 3.342