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. Natural language support but running in an English locale 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(97,100.7,101.4,101.5,101.8,101.5,102.2,101.8,98.5,98.4,97.5,97.7,98.3,99.6,99.4,96.7,96.9,96.1,97.9,99.2,97.8,94.9,93.3,91.5,89.1,92.3,91.8,92.1,94.4,92.8,92.6,92.3,92.1,89.8,87.4,87.7,86.3,89.1,90.4,87.1,86.7,84.4,88.4,88.9,88.5,87.2,86.2,83.4,87.5,85.7,87.4,86.8,87.9,85.9,87.7,87,86.8,86.2,86.1,87.5,85.7,88.9,89.8,91.4,95.2,94.1,96.8,96.1,96.6,94.2,93.9,96.5,93.4,95,95.2,94,97,96.9,96.3,96.3,97.3,95.7,96.4,95.1,94.6,95.9,96.2,94.3,98.3,95.9,92.1,94.6,94.7,96.7,97.5,96.2,97.1,95.9,94.5,99.4,101.3,101.4,100.9,101.4,103.1,102.4,101.1,102,103.9,101.7,101.2,101.9,101.1,103.1,103.3,101.4,102.8,103,102.6,102.2) > 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] 120 > (np <- floor(n / par1)) [1] 10 > 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] 10 10 10 10 10 10 10 10 10 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 97.0 98.3 89.1 86.3 87.5 85.7 93.4 94.6 97.1 103.9 NA [2,] 100.7 99.6 92.3 89.1 85.7 88.9 95.0 95.9 95.9 101.7 NA [3,] 101.4 99.4 91.8 90.4 87.4 89.8 95.2 96.2 94.5 101.2 NA [4,] 101.5 96.7 92.1 87.1 86.8 91.4 94.0 94.3 99.4 101.9 NA [5,] 101.8 96.9 94.4 86.7 87.9 95.2 97.0 98.3 101.3 101.1 NA [6,] 101.5 96.1 92.8 84.4 85.9 94.1 96.9 95.9 101.4 103.1 NA [7,] 102.2 97.9 92.6 88.4 87.7 96.8 96.3 92.1 100.9 103.3 NA [8,] 101.8 99.2 92.3 88.9 87.0 96.1 96.3 94.6 101.4 101.4 NA [9,] 98.5 97.8 92.1 88.5 86.8 96.6 97.3 94.7 103.1 102.8 NA [10,] 98.4 94.9 89.8 87.2 86.2 94.2 95.7 96.7 102.4 103.0 NA [11,] 97.5 93.3 87.4 86.2 86.1 93.9 96.4 97.5 101.1 102.6 NA [12,] 97.7 91.5 87.7 83.4 87.5 96.5 95.1 96.2 102.0 102.2 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 3.7 1.3 3.2 2.8 -1.8 3.2 1.6 1.3 -1.2 -2.2 NA [2,] 0.7 -0.2 -0.5 1.3 1.7 0.9 0.2 0.3 -1.4 -0.5 NA [3,] 0.1 -2.7 0.3 -3.3 -0.6 1.6 -1.2 -1.9 4.9 0.7 NA [4,] 0.3 0.2 2.3 -0.4 1.1 3.8 3.0 4.0 1.9 -0.8 NA [5,] -0.3 -0.8 -1.6 -2.3 -2.0 -1.1 -0.1 -2.4 0.1 2.0 NA [6,] 0.7 1.8 -0.2 4.0 1.8 2.7 -0.6 -3.8 -0.5 0.2 NA [7,] -0.4 1.3 -0.3 0.5 -0.7 -0.7 0.0 2.5 0.5 -1.9 NA [8,] -3.3 -1.4 -0.2 -0.4 -0.2 0.5 1.0 0.1 1.7 1.4 NA [9,] -0.1 -2.9 -2.3 -1.3 -0.6 -2.4 -1.6 2.0 -0.7 0.2 NA [10,] -0.9 -1.6 -2.4 -1.0 -0.1 -0.3 0.7 0.8 -1.3 -0.4 NA [11,] 0.2 -1.8 0.3 -2.8 1.4 2.6 -1.3 -1.3 0.9 -0.4 NA [12,] 0.6 -2.4 -1.4 4.1 -1.8 -3.1 -0.5 0.9 1.9 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/freestat/rcomp/tmp/1cdbl1289390135.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/freestat/rcomp/tmp/2cdbl1289390135.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/freestat/rcomp/tmp/3cdbl1289390135.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/freestat/rcomp/tmp/44mto1289390135.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,] 85.7 85.70 87.40 86.80 86.70 84.4 87.70 87.0 86.80 86.2 86.10 [2,] 87.5 89.10 90.40 91.40 94.40 92.8 92.10 92.3 92.10 89.8 87.40 [3,] 94.0 95.45 94.85 94.15 96.95 96.0 96.55 96.2 96.95 95.3 95.15 [4,] 97.1 99.60 99.40 99.40 101.10 101.4 100.90 101.4 98.50 98.4 97.50 [5,] 103.9 101.70 101.40 101.90 101.80 103.1 103.30 101.8 103.10 103.0 102.60 [,12] [1,] 83.40 [2,] 87.70 [3,] 95.65 [4,] 97.70 [5,] 102.20 $n [1] 10 10 10 10 10 10 10 10 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 89.20346 90.20378 90.35324 90.15288 93.60241 91.7031 92.15317 [2,] 98.79654 100.69622 99.34676 98.14712 100.29759 100.2969 100.94683 [,8] [,9] [,10] [,11] [,12] [1,] 91.65328 93.7523 91.0031 90.10364 90.6536 [2,] 100.74672 100.1477 99.5969 100.19636 100.6464 $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(85.7, 87.5, 94, 97.1, 103.9, 85.7, 89.1, 95.45, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/54mto1289390135.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] [,12] [1,] -2.20 -1.40 -3.30 -0.8 -2.40 -3.80 -1.90 -1.40 -2.9 -2.40 -2.8 -3.1 [2,] -1.20 -0.50 -1.90 0.2 -2.00 -0.50 -0.70 -0.40 -2.3 -1.30 -1.3 -1.8 [3,] 1.45 0.25 -0.25 1.5 -0.95 0.45 -0.15 -0.05 -1.0 -0.65 -0.1 -0.5 [4,] 3.20 0.90 0.70 3.0 -0.10 1.80 0.50 1.00 -0.1 -0.10 0.9 0.9 [5,] 3.70 1.70 1.60 4.0 2.00 4.00 1.30 1.70 2.0 0.80 2.6 4.1 $n [1] 10 10 10 10 10 10 10 10 10 10 10 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.7484154 -0.4494958 -1.549064 0.1010084 -1.8993157536 -0.6991717 [2,] 3.6484154 0.9494958 1.049064 2.8989916 -0.0006842464 1.5991717 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.7495678 -0.7494958 -2.09920771 -1.24956784 -1.1992077 -1.922 [2,] 0.4495678 0.6494958 0.09920771 -0.05043216 0.9992077 0.922 $out [1] 4.9 2.5 -3.3 $group [1] 3 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2.20000000000000, -1.19999999999999, 1.45000000000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6xvar1289390135.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,] 97.00 91.50 87.40 83.40 85.70 85.70 93.40 92.10 94.50 101.10 NA [2,] 98.05 95.50 89.45 86.25 86.15 90.60 95.05 94.60 98.25 101.55 NA [3,] 101.05 97.35 92.10 87.15 86.90 94.15 96.00 95.90 101.20 102.40 NA [4,] 101.65 98.75 92.45 88.70 87.50 96.30 96.65 96.45 101.70 103.05 NA [5,] 102.20 99.60 94.40 90.40 87.90 96.80 97.30 98.30 103.10 103.90 NA $n [1] 12 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 99.40802 95.86765 90.73168 86.03254 86.28426 91.5502 95.27023 95.0562 [2,] 102.69198 98.83235 93.46832 88.26746 87.51574 96.7498 96.72977 96.7438 [,9] [,10] [,11] [1,] 99.62643 101.7158 NA [2,] 102.77357 103.0842 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(97, 98.05, 101.05, 101.65, 102.2, 91.5, 95.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7xvar1289390135.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,] 93.29 94.000 92.48750 [2,] 94.34 95.000 93.73750 [3,] 94.79 95.550 94.91250 [4,] 95.82 96.375 96.44375 [5,] 96.06 96.950 97.50000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 94.11496 94.92285 93.67816 [2,] 95.46504 96.17715 96.14684 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(93.29, 94.34, 94.79, 95.82, 96.06, 94, 95, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1cdbl1289390135.ps tmp/1cdbl1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/2cdbl1289390135.ps tmp/2cdbl1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/3cdbl1289390135.ps tmp/3cdbl1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/44mto1289390135.ps tmp/44mto1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/54mto1289390135.ps tmp/54mto1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/6xvar1289390135.ps tmp/6xvar1289390135.png",intern=TRUE)) character(0) > try(system("convert tmp/7xvar1289390135.ps tmp/7xvar1289390135.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.947 1.465 2.246