R version 2.10.1 (2009-12-14) 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(19097.1,19304.6,19601.7,16006.9,17681.2,19790.4,17014.2,17424.5,18908.9,15692.1,15160,15794.3,16032.1,16065,16236.8,12521,14762.1,15446.9,13635,14212.6,15021.7,14134.3,13721.4,14384.5,15638.6,19711.6,20359.8,16141.4,20056.9,20605.5,19325.8,20547.7,19211.2,19009.5,18746.8,16471.5,18957.2,20515.2,18374.4,16192.9,18147.5,19301.4,18344.7,17183.6,19630,17167.2,17428.5,16016.5,18466.5,18406.6,18174.1,14851.9,16260.7,18329.6,18003.8,15903.8,19554.2,16554.2,16198.9,16571.8,17535.2,16198.1,17487.5,13768,14915.8,17160.9,15607.4,16181.5,17413.2,15116.3,14544.5,15050.6,15535.4,15919.3,15853.1,12336.4,14355.5,16040.8,13867.7,14656.6) > 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] 80 > (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] 7 7 7 7 7 7 7 7 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 19097.1 16032.1 15638.6 18957.2 18466.5 17535.2 15535.4 [2,] 19304.6 16065.0 19711.6 20515.2 18406.6 16198.1 15919.3 [3,] 19601.7 16236.8 20359.8 18374.4 18174.1 17487.5 15853.1 [4,] 16006.9 12521.0 16141.4 16192.9 14851.9 13768.0 12336.4 [5,] 17681.2 14762.1 20056.9 18147.5 16260.7 14915.8 14355.5 [6,] 19790.4 15446.9 20605.5 19301.4 18329.6 17160.9 16040.8 [7,] 17014.2 13635.0 19325.8 18344.7 18003.8 15607.4 13867.7 [8,] 17424.5 14212.6 20547.7 17183.6 15903.8 16181.5 14656.6 [9,] 18908.9 15021.7 19211.2 19630.0 19554.2 17413.2 NA [10,] 15692.1 14134.3 19009.5 17167.2 16554.2 15116.3 NA [11,] 15160.0 13721.4 18746.8 17428.5 16198.9 14544.5 NA [12,] 15794.3 14384.5 16471.5 16016.5 16571.8 15050.6 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 207.5 32.9 4073.0 1558.0 -59.9 -1337.1 383.9 [2,] 297.1 171.8 648.2 -2140.8 -232.5 1289.4 -66.2 [3,] -3594.8 -3715.8 -4218.4 -2181.5 -3322.2 -3719.5 -3516.7 [4,] 1674.3 2241.1 3915.5 1954.6 1408.8 1147.8 2019.1 [5,] 2109.2 684.8 548.6 1153.9 2068.9 2245.1 1685.3 [6,] -2776.2 -1811.9 -1279.7 -956.7 -325.8 -1553.5 -2173.1 [7,] 410.3 577.6 1221.9 -1161.1 -2100.0 574.1 788.9 [8,] 1484.4 809.1 -1336.5 2446.4 3650.4 1231.7 NA [9,] -3216.8 -887.4 -201.7 -2462.8 -3000.0 -2296.9 NA [10,] -532.1 -412.9 -262.7 261.3 -355.3 -571.8 NA [11,] 634.3 663.1 -2275.3 -1412.0 372.9 506.1 NA [12,] 237.8 1254.1 2485.7 2450.0 963.4 484.8 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/yougetitorg/rcomp/tmp/1w0kp1305823673.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/yougetitorg/rcomp/tmp/2tafx1305823673.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/yougetitorg/rcomp/tmp/37ru91305823673.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/yougetitorg/rcomp/tmp/4tqos1305823673.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] [1,] 15535.40 15919.30 15853.10 12336.40 14355.50 15446.90 13635.00 14212.60 [2,] 15835.35 16131.55 16862.15 13144.50 14838.95 16600.85 14737.55 15280.20 [3,] 17535.20 18406.60 18174.10 14851.90 16260.70 18329.60 17014.20 16181.50 [4,] 18711.85 19508.10 18988.05 16074.15 17914.35 19545.90 18174.25 17304.05 [5,] 19097.10 20515.20 20359.80 16192.90 20056.90 20605.50 19325.80 17424.50 [,9] [,10] [,11] [,12] [1,] 15021.70 14134.30 13721.40 14384.5 [2,] 17413.20 15116.30 14544.50 15050.6 [3,] 19060.05 16123.15 15679.45 15905.4 [4,] 19554.20 17167.20 17428.50 16471.5 [5,] 19630.00 19009.50 18746.80 16571.8 $n [1] 7 7 7 7 7 7 7 7 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15817.4 16390.18 16904.55 13102.36 14424.12 16570.86 14961.86 14972.89 [2,] 19253.0 20423.02 19443.65 16601.44 18097.28 20088.34 19066.54 17390.11 [,9] [,10] [,11] [,12] [1,] 17679.04 14800.25 13819.18 14988.87 [2,] 20441.06 17446.05 17539.72 16821.93 $out [1] 20547.7 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(15535.4, 15835.35, 17535.2, 18711.85, 19097.1, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/5q91f1305823673.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] [1,] -1337.10 -232.50 -3719.50 1147.80 548.60 -2776.2 -1161.10 -1336.50 [2,] -13.50 -149.35 -3717.65 1541.55 919.35 -1992.5 -375.40 809.10 [3,] 207.50 171.80 -3594.80 1954.60 1685.30 -1553.5 574.10 1358.05 [4,] 970.95 472.65 -3419.45 2130.10 2089.05 -1118.2 683.25 2446.40 [5,] 1558.00 1289.40 -3322.20 2241.10 2245.10 -325.8 1221.90 3650.40 [,9] [,10] [,11] [,12] [1,] -3216.80 -571.8 -2275.3 237.80 [2,] -3000.00 -532.1 -1412.0 484.80 [3,] -2379.85 -384.1 439.5 1108.75 [4,] -887.40 -262.7 634.3 2450.00 [5,] -201.70 -262.7 663.1 2485.70 $n [1] 7 7 7 7 7 7 7 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -380.3977 -199.6484 -3772.88 1603.127 986.774 -2075.618 -58.1087 [2,] 795.3977 543.2484 -3416.72 2306.073 2383.826 -1031.382 1206.3087 [,8] [,9] [,10] [,11] [,12] [1,] 301.9386 -3742.545 -557.8717 -880.4296 -158.8675 [2,] 2414.1614 -1017.155 -210.3283 1759.4296 2376.3675 $out [1] 4073.0 -2140.8 -4218.4 -2181.5 3915.5 -2100.0 261.3 $group [1] 1 2 3 3 4 7 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1337.1, -13.5000000000009, 207.5, 970.95, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/6t6bq1305823673.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,] 15160.00 12521.00 15638.60 16016.5 14851.9 13768.00 12336.4 [2,] 15900.60 13927.85 17609.15 17175.4 16229.8 14983.20 14111.6 [3,] 17552.85 14573.30 19268.50 18246.1 17287.8 15894.45 15096.0 [4,] 19200.85 15739.50 20208.35 19129.3 18368.1 17287.05 15886.2 [5,] 19790.40 16236.80 20605.50 20515.2 19554.2 17535.20 16040.8 $n [1] 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 16047.58 13746.99 18082.99 17354.91 16312.51 14843.65 14104.68 [2,] 19058.12 15399.61 20454.01 19137.29 18263.09 16945.25 16087.32 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(15160, 15900.6, 17552.85, 19200.85, 19790.4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/787r41305823673.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,] 14545.50 14851.90 14609.33 [2,] 16122.81 16014.28 16023.42 [3,] 16592.14 16637.45 16416.28 [4,] 18014.84 18251.85 17872.46 [5,] 18289.87 19060.05 18627.79 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 15729.17 15616.88 15572.91 [2,] 17455.11 17658.02 17259.64 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(14545.5, 16122.8083333333, 16592.1428571429, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1w0kp1305823673.ps tmp/1w0kp1305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/2tafx1305823673.ps tmp/2tafx1305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/37ru91305823673.ps tmp/37ru91305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/4tqos1305823673.ps tmp/4tqos1305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/5q91f1305823673.ps tmp/5q91f1305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/6t6bq1305823673.ps tmp/6t6bq1305823673.png",intern=TRUE)) character(0) > try(system("convert tmp/787r41305823673.ps tmp/787r41305823673.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.460 1.140 2.115