R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(121,117.7,115.4,114.3,109.5,108.1,108.2,99.1,101.2,98.1,95.5,97.9,98.2,98.7,95.6,95.8,94.4,96.5,103.3,104.3,104.5,102.3,103.8,103.1,102.2,106.3,102.1,94,102.6,102.6,106.7,107.9,109.3,105.9,109.1,108.5,111.7,109.8,109.1,108.5,108.5,106.2,117.1,109.8,115.2,115.9,119.2,121,118.6,117.6,114.6,110.6,102.5,101.6,107.4,105.8,102.8,104,100.4,100.6,107.9,106.9,106.5,103,90.5,90.6,94.4,89.4,92.5,94.4,91.7,93.3) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 121.0 98.2 102.2 111.7 118.6 107.9 NA [2,] 117.7 98.7 106.3 109.8 117.6 106.9 NA [3,] 115.4 95.6 102.1 109.1 114.6 106.5 NA [4,] 114.3 95.8 94.0 108.5 110.6 103.0 NA [5,] 109.5 94.4 102.6 108.5 102.5 90.5 NA [6,] 108.1 96.5 102.6 106.2 101.6 90.6 NA [7,] 108.2 103.3 106.7 117.1 107.4 94.4 NA [8,] 99.1 104.3 107.9 109.8 105.8 89.4 NA [9,] 101.2 104.5 109.3 115.2 102.8 92.5 NA [10,] 98.1 102.3 105.9 115.9 104.0 94.4 NA [11,] 95.5 103.8 109.1 119.2 100.4 91.7 NA [12,] 97.9 103.1 108.5 121.0 100.6 93.3 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3.3 0.5 4.1 -1.9 -1.0 -1.0 NA [2,] -2.3 -3.1 -4.2 -0.7 -3.0 -0.4 NA [3,] -1.1 0.2 -8.1 -0.6 -4.0 -3.5 NA [4,] -4.8 -1.4 8.6 0.0 -8.1 -12.5 NA [5,] -1.4 2.1 0.0 -2.3 -0.9 0.1 NA [6,] 0.1 6.8 4.1 10.9 5.8 3.8 NA [7,] -9.1 1.0 1.2 -7.3 -1.6 -5.0 NA [8,] 2.1 0.2 1.4 5.4 -3.0 3.1 NA [9,] -3.1 -2.2 -3.4 0.7 1.2 1.9 NA [10,] -2.6 1.5 3.2 3.3 -3.6 -2.7 NA [11,] 2.4 -0.7 -0.6 1.8 0.2 1.6 NA [12,] 0.3 -0.9 3.2 -2.4 7.3 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/15e4w1319657986.ps",horizontal=F,onefile=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/2k3me1319657986.ps",horizontal=F,onefile=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/3d1n51319657986.ps",horizontal=F,onefile=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/4u4y61319657986.ps",horizontal=F,onefile=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,] 98.2 98.70 95.6 94.00 90.50 90.6 103.30 89.40 92.50 94.40 91.7 [2,] 102.2 106.30 102.1 95.80 94.40 96.5 103.30 99.10 101.20 98.10 95.5 [3,] 109.8 108.35 107.8 105.75 102.55 102.1 107.05 105.05 103.65 103.15 102.1 [4,] 118.6 117.60 114.6 110.60 108.50 106.2 108.20 107.90 109.30 105.90 109.1 [5,] 121.0 117.70 115.4 114.30 109.50 108.1 108.20 109.80 115.20 115.90 119.2 [,12] [1,] 93.30 [2,] 97.90 [3,] 101.85 [4,] 108.50 [5,] 121.00 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 99.22147 101.0611 99.7371 96.20352 93.45504 95.84319 103.8893 [2,] 120.37853 115.6389 115.8629 115.29648 111.64496 108.35681 110.2107 [,8] [,9] [,10] [,11] [,12] [1,] 99.37372 98.42524 98.11875 93.32756 95.01266 [2,] 110.72628 108.87476 108.18125 110.87244 108.68734 $out [1] 117.1 94.4 $group [1] 7 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(98.2, 102.2, 109.8, 118.6, 121, 98.7, 106.3, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5y50r1319657986.ps",horizontal=F,onefile=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,] -3.3 -4.20 -8.1 -12.5 -2.30 0.10 -9.1 -3.00 -3.40 -3.60 -0.7 -2.4 [2,] -1.9 -3.10 -4.0 -8.1 -1.40 3.80 -7.3 0.20 -3.10 -2.70 -0.6 -0.9 [3,] -1.0 -2.65 -2.3 -3.1 -0.45 4.95 -3.3 1.75 -0.75 -0.55 0.9 0.3 [4,] 0.5 -0.70 -0.6 0.0 0.10 6.80 1.0 3.10 1.20 3.20 1.8 3.2 [5,] 4.1 -0.40 0.2 8.6 2.10 10.90 1.2 5.40 1.90 3.30 2.4 7.3 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2.5480775 -4.198078 -4.4931098 -8.324762 -1.4175484 3.014903 -8.653768 [2,] 0.5480775 -1.101922 -0.1068902 2.124762 0.5175484 6.885097 2.053768 [,8] [,9] [,10] [,11] [,12] [1,] -0.1205937 -3.523639 -4.355691 -0.6480775 -2.597050 [2,] 3.6205937 2.023639 3.255691 2.4480775 3.197050 $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(-3.3, -1.90000000000001, -1, 0.5, 4.09999999999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/617f91319657986.ps",horizontal=F,onefile=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] [1,] 95.50 94.40 94.0 106.20 100.40 89.40 NA [2,] 98.60 96.15 102.4 108.80 102.05 91.15 NA [3,] 108.15 100.50 106.1 110.75 104.90 93.85 NA [4,] 114.85 103.55 108.2 116.50 112.60 104.75 NA [5,] 121.00 104.50 109.3 121.00 118.60 107.90 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 100.7383 97.12481 103.4546 107.2380 100.0881 87.64695 NA [2,] 115.5617 103.87519 108.7454 114.2620 109.7119 100.05305 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(95.5, 98.6, 108.15, 114.85, 121, 94.4, 96.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7ifga1319657986.ps",horizontal=F,onefile=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,] 100.9333 101.850 101.5375 [2,] 103.0000 102.325 102.2687 [3,] 104.1583 104.350 103.8625 [4,] 106.7000 107.425 107.1437 [5,] 109.9333 109.800 111.0500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 102.4707 102.0239 101.6390 [2,] 105.8459 106.6761 106.0860 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(100.933333333333, 103, 104.158333333333, 106.7, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/15e4w1319657986.ps tmp/15e4w1319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/2k3me1319657986.ps tmp/2k3me1319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/3d1n51319657986.ps tmp/3d1n51319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/4u4y61319657986.ps tmp/4u4y61319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/5y50r1319657986.ps tmp/5y50r1319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/617f91319657986.ps tmp/617f91319657986.png",intern=TRUE)) character(0) > try(system("convert tmp/7ifga1319657986.ps tmp/7ifga1319657986.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.410 0.260 2.667