R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(100.09,100.09,100.06,100.11,100.08,100.08,100.08,100.08,99.94,99.79,99.98,99.98,99.98,99.98,99.93,99.93,99.93,99.93,99.93,99.93,99.76,99.48,99.55,99.56,99.57,99.56,99.65,99.5,99.5,99.49,99.49,99.46,99.52,99.49,99.55,99.57,99.57,99.57,99.57,99.57,99.57,99.57,99.57,99.53,100.38,100.32,100.46,100.47,100.47,100.47,100.51,100.5,100.51,100.51,100.51,100.51,101.65,102.13,102.2,102.13,102.13,102.12,102.13,102.05,102,102.01,102.01,102.02,102.78,103.39,103.41,103.5,103.5,103.49,103.38,103.24,103.25,103.25,103.25,103.25,103.83,104.33,104.36,104.48,104.5,104.48,104.35,104.48,104.48,104.47,104.47,104.86,105.22,105.96,106.03,106.03) > 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] 96 > (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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 100.09 99.98 99.57 99.57 100.47 102.13 103.50 104.50 NA [2,] 100.09 99.98 99.56 99.57 100.47 102.12 103.49 104.48 NA [3,] 100.06 99.93 99.65 99.57 100.51 102.13 103.38 104.35 NA [4,] 100.11 99.93 99.50 99.57 100.50 102.05 103.24 104.48 NA [5,] 100.08 99.93 99.50 99.57 100.51 102.00 103.25 104.48 NA [6,] 100.08 99.93 99.49 99.57 100.51 102.01 103.25 104.47 NA [7,] 100.08 99.93 99.49 99.57 100.51 102.01 103.25 104.47 NA [8,] 100.08 99.93 99.46 99.53 100.51 102.02 103.25 104.86 NA [9,] 99.94 99.76 99.52 100.38 101.65 102.78 103.83 105.22 NA [10,] 99.79 99.48 99.49 100.32 102.13 103.39 104.33 105.96 NA [11,] 99.98 99.55 99.55 100.46 102.20 103.41 104.36 106.03 NA [12,] 99.98 99.56 99.57 100.47 102.13 103.50 104.48 106.03 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 0.00 0.00 -0.01 0.00 0.00 -0.01 -0.01 -0.02 NA [2,] -0.03 -0.05 0.09 0.00 0.04 0.01 -0.11 -0.13 NA [3,] 0.05 0.00 -0.15 0.00 -0.01 -0.08 -0.14 0.13 NA [4,] -0.03 0.00 0.00 0.00 0.01 -0.05 0.01 0.00 NA [5,] 0.00 0.00 -0.01 0.00 0.00 0.01 0.00 -0.01 NA [6,] 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 NA [7,] 0.00 0.00 -0.03 -0.04 0.00 0.01 0.00 0.39 NA [8,] -0.14 -0.17 0.06 0.85 1.14 0.76 0.58 0.36 NA [9,] -0.15 -0.28 -0.03 -0.06 0.48 0.61 0.50 0.74 NA [10,] 0.19 0.07 0.06 0.14 0.07 0.02 0.03 0.07 NA [11,] 0.00 0.01 0.02 0.01 -0.07 0.09 0.12 0.00 NA [12,] 0.00 0.01 0.00 0.00 0.00 0.00 0.02 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/html/rcomp/tmp/13s191263845122.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/213gk1263845122.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/3r0x91263845122.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/4x3co1263845122.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] [1,] 99.570 99.560 99.570 99.500 99.500 99.490 99.490 99.460 99.520 [2,] 99.775 99.775 99.790 99.750 99.750 99.750 99.750 99.730 99.850 [3,] 100.280 100.280 100.285 100.305 100.295 100.295 100.295 100.295 101.015 [4,] 102.815 102.805 102.755 102.645 102.625 102.630 102.630 102.635 103.305 [5,] 104.500 104.480 104.350 104.480 104.480 104.470 104.470 104.860 105.220 [,10] [,11] [,12] [1,] 99.480 99.550 99.560 [2,] 99.640 99.765 99.775 [3,] 101.225 101.330 101.300 [4,] 103.860 103.885 103.990 [5,] 105.960 106.030 106.030 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 98.58181 98.5874 98.6287 98.68781 98.68898 98.6862 98.6862 [2,] 101.97819 101.9726 101.9413 101.92219 101.90102 101.9038 101.9038 [,8] [,9] [,10] [,11] [,12] [1,] 98.67223 99.08499 98.86765 99.0285 98.94544 [2,] 101.91777 102.94501 103.58235 103.6315 103.65456 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(99.57, 99.775, 100.28, 102.815, 104.5, 99.56, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5scgs1263845122.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,] -0.020 -0.130 -0.150 -0.030 -0.010 0 -0.040 -0.170 -0.280 0.020 -0.070 [2,] -0.010 -0.080 -0.110 -0.015 -0.005 0 -0.015 -0.040 -0.105 0.045 0.000 [3,] -0.005 -0.015 -0.005 0.000 0.000 0 0.000 0.470 0.225 0.070 0.010 [4,] 0.000 0.025 0.025 0.005 0.000 0 0.005 0.805 0.555 0.105 0.055 [5,] 0.000 0.090 0.130 0.010 0.000 0 0.010 1.140 0.740 0.190 0.120 [,12] [1,] 0.000 [2,] 0.000 [3,] 0.000 [4,] 0.005 [5,] 0.010 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.0105861436 -0.07365451 -0.08041294 -0.01117229 -0.002793072 0 [2,] 0.0005861436 0.04365451 0.07041294 0.01117229 0.002793072 0 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.01117229 -0.002029132 -0.1436855 0.03648314 -0.02072379 -0.002985919 [2,] 0.01117229 0.942029132 0.5936855 0.10351686 0.04072379 0.002985919 $out [1] -0.05 0.01 0.39 0.02 $group [1] 4 5 7 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.019999999999996, -0.009999999999998, -0.00499999999999545, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6suc71263845122.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,] 99.940 99.48 99.460 99.53 100.470 102.000 103.240 104.350 NA [2,] 99.980 99.66 99.490 99.57 100.505 102.015 103.250 104.475 NA [3,] 100.080 99.93 99.510 99.57 100.510 102.125 103.435 104.490 NA [4,] 100.085 99.93 99.565 100.35 101.890 103.085 104.080 105.590 NA [5,] 100.110 99.98 99.650 100.47 102.200 103.500 104.480 106.030 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0321 99.80685 99.4758 99.21424 99.8783 101.6370 103.0564 103.9814 [2,] 100.1279 100.05315 99.5442 99.92576 101.1417 102.6130 103.8136 104.9986 [,9] [1,] NA [2,] NA $out [1] 99.79 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(99.94, 99.98, 100.08, 100.085, 100.11, 99.48, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7zn5f1263845122.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,] 101.1637 100.280 101.0763 [2,] 101.1688 100.290 101.0800 [3,] 101.2125 100.295 101.1606 [4,] 101.7481 101.120 101.5694 [5,] 101.9650 101.330 101.8113 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 100.9482 99.91643 100.9374 [2,] 101.4768 100.67357 101.3838 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(101.16375, 101.16875, 101.2125, 101.748125, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/13s191263845122.ps tmp/13s191263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/213gk1263845122.ps tmp/213gk1263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/3r0x91263845122.ps tmp/3r0x91263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/4x3co1263845122.ps tmp/4x3co1263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/5scgs1263845122.ps tmp/5scgs1263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/6suc71263845122.ps tmp/6suc71263845122.png",intern=TRUE)) character(0) > try(system("convert tmp/7zn5f1263845122.ps tmp/7zn5f1263845122.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.242 1.007 14.081