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. 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(2648.9,2669.6,3042.3,2604.2,2732.1,2621.7,2483.7,2479.3,2684.6,2834.7,2566.1,2251.2,2350,2299.8,2542.8,2530.2,2508.1,2616.8,2534.1,2181.8,2578.9,2841.9,2529.9,2103.2,2326.2,2452.6,2782.1,2727.3,2648.2,2760.7,2613,2225.4,2713.9,2923.3,2707,2473.9,2521,2531.8,3068.8,2826.9,2674.2,2966.6,2798.8,2629.6,3124.6,3115.7,3083,2863.9,2728.7,2789.4,3225.7,3148.2,2836.5,3153.5,2656.9,2834.7,3172.5,2998.8,3103.1,2735.6,2818.1,2874.4,3438.5,2949.1,3306.8,3530,3003.8,3206.4,3514.6,3522.6,3525.5,2996.2,3231.1,3030,3541.7,3113.2,3390.8,3424.2,3079.8,3123.4,3317.1,3579.9,3317.9,2668.1) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2648.9 2350.0 2326.2 2521.0 2728.7 2818.1 3231.1 NA [2,] 2669.6 2299.8 2452.6 2531.8 2789.4 2874.4 3030.0 NA [3,] 3042.3 2542.8 2782.1 3068.8 3225.7 3438.5 3541.7 NA [4,] 2604.2 2530.2 2727.3 2826.9 3148.2 2949.1 3113.2 NA [5,] 2732.1 2508.1 2648.2 2674.2 2836.5 3306.8 3390.8 NA [6,] 2621.7 2616.8 2760.7 2966.6 3153.5 3530.0 3424.2 NA [7,] 2483.7 2534.1 2613.0 2798.8 2656.9 3003.8 3079.8 NA [8,] 2479.3 2181.8 2225.4 2629.6 2834.7 3206.4 3123.4 NA [9,] 2684.6 2578.9 2713.9 3124.6 3172.5 3514.6 3317.1 NA [10,] 2834.7 2841.9 2923.3 3115.7 2998.8 3522.6 3579.9 NA [11,] 2566.1 2529.9 2707.0 3083.0 3103.1 3525.5 3317.9 NA [12,] 2251.2 2103.2 2473.9 2863.9 2735.6 2996.2 2668.1 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 20.7 -50.2 126.4 10.8 60.7 56.3 -201.1 NA [2,] 372.7 243.0 329.5 537.0 436.3 564.1 511.7 NA [3,] -438.1 -12.6 -54.8 -241.9 -77.5 -489.4 -428.5 NA [4,] 127.9 -22.1 -79.1 -152.7 -311.7 357.7 277.6 NA [5,] -110.4 108.7 112.5 292.4 317.0 223.2 33.4 NA [6,] -138.0 -82.7 -147.7 -167.8 -496.6 -526.2 -344.4 NA [7,] -4.4 -352.3 -387.6 -169.2 177.8 202.6 43.6 NA [8,] 205.3 397.1 488.5 495.0 337.8 308.2 193.7 NA [9,] 150.1 263.0 209.4 -8.9 -173.7 8.0 262.8 NA [10,] -268.6 -312.0 -216.3 -32.7 104.3 2.9 -262.0 NA [11,] -314.9 -426.7 -233.1 -219.1 -367.5 -529.3 -649.8 NA [12,] 98.8 223.0 47.1 -135.2 82.5 234.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/rcomp/tmp/1477i1229624489.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/2frqo1229624489.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/3so5k1229624489.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/4wcio1229624489.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,] 2326.2 2299.8 2542.8 2530.20 2508.10 2616.80 2483.70 2181.80 2578.90 [2,] 2435.5 2492.2 2912.2 2665.75 2661.20 2691.20 2573.55 2352.35 2699.25 [3,] 2648.9 2669.6 3068.8 2826.90 2732.10 2966.60 2656.90 2629.60 3124.60 [4,] 2773.4 2831.9 3332.1 3031.15 3071.65 3288.85 2901.30 2979.05 3244.80 [5,] 3231.1 3030.0 3541.7 3148.20 3390.80 3530.00 3079.80 3206.40 3514.60 [,10] [,11] [,12] [1,] 2834.70 2529.90 2103.20 [2,] 2882.60 2636.55 2362.55 [3,] 2998.80 3083.00 2668.10 [4,] 3319.15 3210.50 2799.75 [5,] 3579.90 3525.50 2996.20 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2447.112 2466.737 2818.042 2608.689 2486.986 2609.693 2461.173 2255.345 [2,] 2850.688 2872.463 3319.558 3045.111 2977.214 3323.507 2852.627 3003.855 [,9] [,10] [,11] [,12] [1,] 2798.806 2738.099 2740.246 2407.011 [2,] 3450.394 3259.501 3425.754 2929.189 $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(2326.2, 2435.5, 2648.9, 2773.4, 3231.1, 2299.8, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5g2z31229624489.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] [1,] -50.2 243.00 -489.40 -311.70 -110.40 -526.20 -387.60 193.70 -173.70 -312.0 [2,] -19.7 351.10 -433.30 -115.90 71.05 -420.50 -260.75 256.75 -0.45 -265.3 [3,] 20.7 436.30 -241.90 -22.10 112.50 -167.80 -4.40 337.80 150.10 -216.3 [4,] 58.5 524.35 -66.15 202.75 257.80 -142.85 110.70 442.80 236.10 -14.9 [5,] 126.4 564.10 -12.60 357.70 317.00 -82.70 202.60 495.00 263.00 104.3 [,11] [,12] [1,] -649.8 -135.20 [2,] -478.0 47.10 [3,] -367.5 90.65 [4,] -274.0 223.00 [5,] -219.1 234.90 $n [1] 7 7 7 7 7 7 7 7 7 7 7 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -25.99978 332.8379 -461.15606 -212.3926 0.9759128 -333.608101 -226.2239 [2,] 67.39978 539.7621 -22.64394 168.1926 224.0240872 -1.991899 217.4239 [,8] [,9] [,10] [,11] [,12] [1,] 226.6939 8.836156 -365.83484 -489.3255 -22.81118 [2,] 448.9061 291.363844 -66.76516 -245.6745 204.11118 $out [1] -201.1 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-50.1999999999998, -19.6999999999998, 20.6999999999998, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6yz801229624489.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] [1,] 2251.20 2103.20 2225.40 2521.0 2656.90 2818.10 2668.1 NA [2,] 2524.90 2324.90 2463.25 2651.9 2762.50 2972.65 3096.5 NA [3,] 2635.30 2530.05 2677.60 2845.4 2917.65 3256.60 3274.1 NA [4,] 2708.35 2560.85 2744.00 3075.9 3150.85 3518.60 3407.5 NA [5,] 2834.70 2841.90 2923.30 3124.6 3225.70 3530.00 3579.9 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2551.627 2422.432 2549.548 2652.011 2740.521 3007.589 3132.251 NA [2,] 2718.973 2637.668 2805.652 3038.789 3094.779 3505.611 3415.949 NA $out [1] 3042.3 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(2251.2, 2524.9, 2635.3, 2708.35, 2834.7, 2103.2, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7uozo1229624489.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,] 2584.586 2629.6 2581.150 [2,] 2666.300 2662.5 2663.875 [3,] 2856.843 2779.5 2857.438 [4,] 3012.836 3033.8 2981.025 [5,] 3116.700 3124.6 3122.150 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2698.786 2610.148 2712.783 [2,] 3014.900 2948.852 3002.092 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(2584.58571428571, 2666.3, 2856.84285714286, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1477i1229624489.ps tmp/1477i1229624489.png") > system("convert tmp/2frqo1229624489.ps tmp/2frqo1229624489.png") > system("convert tmp/3so5k1229624489.ps tmp/3so5k1229624489.png") > system("convert tmp/4wcio1229624489.ps tmp/4wcio1229624489.png") > system("convert tmp/5g2z31229624489.ps tmp/5g2z31229624489.png") > system("convert tmp/6yz801229624489.ps tmp/6yz801229624489.png") > system("convert tmp/7uozo1229624489.ps tmp/7uozo1229624489.png") > > > proc.time() user system elapsed 1.234 0.993 2.022