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(27.00,26.88,27.38,26.82,27.00,26.15,25.85,26.17,25.66,25.72,25.42,25.10,26.20,26.39,26.27,26.63,21.10,20.30,20.61,21.05,20.45,20.91,21.22,20.85,21.90,22.71,22.40,22.81,23.96,23.37,23.55,23.01,22.63,22.63,22.00,22.15,22.00,22.00,21.84,22.10,22.37,21.83,21.77,21.89,20.76,20.21,20.19,20.01,19.16,18.50,17.41,18.14,18.60,18.32,18.40,18.16,17.29,16.65,16.36,16.32,17.37,17.30,18.10,19.00,18.38,18.41,18.10,17.87,18.70,18.81,18.88,19.44,18.60,18.80,18.62,18.24,17.84,17.85,17.67,17.99,18.15,18.39,18.07,18.39) > par1 = '42' > #'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] 2 > 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] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [39] 2 2 2 2 > arr [,1] [,2] [,3] [1,] 27.00 21.77 NA [2,] 26.88 21.89 NA [3,] 27.38 20.76 NA [4,] 26.82 20.21 NA [5,] 27.00 20.19 NA [6,] 26.15 20.01 NA [7,] 25.85 19.16 NA [8,] 26.17 18.50 NA [9,] 25.66 17.41 NA [10,] 25.72 18.14 NA [11,] 25.42 18.60 NA [12,] 25.10 18.32 NA [13,] 26.20 18.40 NA [14,] 26.39 18.16 NA [15,] 26.27 17.29 NA [16,] 26.63 16.65 NA [17,] 21.10 16.36 NA [18,] 20.30 16.32 NA [19,] 20.61 17.37 NA [20,] 21.05 17.30 NA [21,] 20.45 18.10 NA [22,] 20.91 19.00 NA [23,] 21.22 18.38 NA [24,] 20.85 18.41 NA [25,] 21.90 18.10 NA [26,] 22.71 17.87 NA [27,] 22.40 18.70 NA [28,] 22.81 18.81 NA [29,] 23.96 18.88 NA [30,] 23.37 19.44 NA [31,] 23.55 18.60 NA [32,] 23.01 18.80 NA [33,] 22.63 18.62 NA [34,] 22.63 18.24 NA [35,] 22.00 17.84 NA [36,] 22.15 17.85 NA [37,] 22.00 17.67 NA [38,] 22.00 17.99 NA [39,] 21.84 18.15 NA [40,] 22.10 18.39 NA [41,] 22.37 18.07 NA [42,] 21.83 18.39 NA > darr [,1] [,2] [,3] [1,] -0.12 0.12 NA [2,] 0.50 -1.13 NA [3,] -0.56 -0.55 NA [4,] 0.18 -0.02 NA [5,] -0.85 -0.18 NA [6,] -0.30 -0.85 NA [7,] 0.32 -0.66 NA [8,] -0.51 -1.09 NA [9,] 0.06 0.73 NA [10,] -0.30 0.46 NA [11,] -0.32 -0.28 NA [12,] 1.10 0.08 NA [13,] 0.19 -0.24 NA [14,] -0.12 -0.87 NA [15,] 0.36 -0.64 NA [16,] -5.53 -0.29 NA [17,] -0.80 -0.04 NA [18,] 0.31 1.05 NA [19,] 0.44 -0.07 NA [20,] -0.60 0.80 NA [21,] 0.46 0.90 NA [22,] 0.31 -0.62 NA [23,] -0.37 0.03 NA [24,] 1.05 -0.31 NA [25,] 0.81 -0.23 NA [26,] -0.31 0.83 NA [27,] 0.41 0.11 NA [28,] 1.15 0.07 NA [29,] -0.59 0.56 NA [30,] 0.18 -0.84 NA [31,] -0.54 0.20 NA [32,] -0.38 -0.18 NA [33,] 0.00 -0.38 NA [34,] -0.63 -0.40 NA [35,] 0.15 0.01 NA [36,] -0.15 -0.18 NA [37,] 0.00 0.32 NA [38,] -0.16 0.16 NA [39,] 0.26 0.24 NA [40,] 0.27 -0.32 NA [41,] -0.54 0.32 NA [42,] -0.06 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/17b8r1243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/2zi681243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/3kiq31243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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/4frb61243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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,] 21.770 21.890 20.76 20.210 20.190 20.01 19.160 18.500 17.410 18.14 18.60 [2,] 21.770 21.890 20.76 20.210 20.190 20.01 19.160 18.500 17.410 18.14 18.60 [3,] 24.385 24.385 24.07 23.515 23.595 23.08 22.505 22.335 21.535 21.93 22.01 [4,] 27.000 26.880 27.38 26.820 27.000 26.15 25.850 26.170 25.660 25.72 25.42 [5,] 27.000 26.880 27.38 26.820 27.000 26.15 25.850 26.170 25.660 25.72 25.42 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 18.32 18.4 18.160 17.29 16.65 16.36 16.32 17.37 17.300 18.100 19.000 [2,] 18.32 18.4 18.160 17.29 16.65 16.36 16.32 17.37 17.300 18.100 19.000 [3,] 21.71 22.3 22.275 21.78 21.64 18.73 18.31 18.99 19.175 19.275 19.955 [4,] 25.10 26.2 26.390 26.27 26.63 21.10 20.30 20.61 21.050 20.450 20.910 [5,] 25.10 26.2 26.390 26.27 26.63 21.10 20.30 20.61 21.050 20.450 20.910 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [1,] 18.38 18.41 18.1 17.87 18.70 18.81 18.88 19.440 18.600 18.800 18.620 [2,] 18.38 18.41 18.1 17.87 18.70 18.81 18.88 19.440 18.600 18.800 18.620 [3,] 19.80 19.63 20.0 20.29 20.55 20.81 21.42 21.405 21.075 20.905 20.625 [4,] 21.22 20.85 21.9 22.71 22.40 22.81 23.96 23.370 23.550 23.010 22.630 [5,] 21.22 20.85 21.9 22.71 22.40 22.81 23.96 23.370 23.550 23.010 22.630 [,34] [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [1,] 18.240 17.84 17.85 17.670 17.990 18.150 18.390 18.07 18.39 [2,] 18.240 17.84 17.85 17.670 17.990 18.150 18.390 18.07 18.39 [3,] 20.435 19.92 20.00 19.835 19.995 19.995 20.245 20.22 20.11 [4,] 22.630 22.00 22.15 22.000 22.000 21.840 22.100 22.37 21.83 [5,] 22.630 22.00 22.15 22.000 22.000 21.840 22.100 22.37 21.83 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [39] 2 2 2 2 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 18.54189 18.81003 16.67395 16.13012 15.98667 16.22022 15.03074 13.76586 [2,] 30.22811 29.95997 31.46605 30.89988 31.20333 29.93978 29.97926 30.90414 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 12.31786 13.46141 14.3905 14.13519 13.58562 13.08021 11.74729 10.49006 [2,] 30.75214 30.39859 29.6295 29.28481 31.01438 31.46979 31.81271 32.78994 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 13.43434 13.86343 15.37018 14.98539 16.64951 17.82109 16.62707 16.90396 [2,] 24.02566 22.75657 22.60982 23.36461 21.90049 22.08891 22.97293 22.35604 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 15.75453 14.88261 16.41625 16.34109 15.74448 17.01429 15.54472 16.20147 [2,] 24.24547 25.69739 24.68375 25.27891 27.09552 25.79571 26.60528 25.60853 [,33] [,34] [,35] [,36] [,37] [,38] [,39] [,40] [1,] 16.14491 15.53037 15.27233 15.19592 14.9974 15.51491 15.87243 16.10008 [2,] 25.10509 25.33963 24.56767 24.80408 24.6726 24.47509 24.11757 24.38992 [,41] [,42] [1,] 15.41592 16.26673 [2,] 25.02408 23.95327 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" Warning message: In bxp(list(stats = c(21.77, 21.77, 24.385, 27, 27, 21.89, 21.89, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5b54x1243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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.12 -1.130 -0.560 -0.02 -0.850 -0.850 -0.66 -1.09 0.060 -0.30 -0.32 [2,] -0.12 -1.130 -0.560 -0.02 -0.850 -0.850 -0.66 -1.09 0.060 -0.30 -0.32 [3,] 0.00 -0.315 -0.555 0.08 -0.515 -0.575 -0.17 -0.80 0.395 0.08 -0.30 [4,] 0.12 0.500 -0.550 0.18 -0.180 -0.300 0.32 -0.51 0.730 0.46 -0.28 [5,] 0.12 0.500 -0.550 0.18 -0.180 -0.300 0.32 -0.51 0.730 0.46 -0.28 [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [1,] 0.08 -0.240 -0.870 -0.64 -5.53 -0.80 0.31 -0.070 -0.6 0.46 -0.620 [2,] 0.08 -0.240 -0.870 -0.64 -5.53 -0.80 0.31 -0.070 -0.6 0.46 -0.620 [3,] 0.59 -0.025 -0.495 -0.14 -2.91 -0.42 0.68 0.185 0.1 0.68 -0.155 [4,] 1.10 0.190 -0.120 0.36 -0.29 -0.04 1.05 0.440 0.8 0.90 0.310 [5,] 1.10 0.190 -0.120 0.36 -0.29 -0.04 1.05 0.440 0.8 0.90 0.310 [,23] [,24] [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [1,] -0.37 -0.31 -0.23 -0.31 0.11 0.07 -0.590 -0.84 -0.54 -0.38 -0.38 -0.630 [2,] -0.37 -0.31 -0.23 -0.31 0.11 0.07 -0.590 -0.84 -0.54 -0.38 -0.38 -0.630 [3,] -0.17 0.37 0.29 0.26 0.26 0.61 -0.015 -0.33 -0.17 -0.28 -0.19 -0.515 [4,] 0.03 1.05 0.81 0.83 0.41 1.15 0.560 0.18 0.20 -0.18 0.00 -0.400 [5,] 0.03 1.05 0.81 0.83 0.41 1.15 0.560 0.18 0.20 -0.18 0.00 -0.400 [,35] [,36] [,37] [,38] [,39] [,40] [,41] [,42] [1,] 0.01 -0.180 0.00 -0.16 0.24 -0.320 -0.54 -0.06 [2,] 0.01 -0.180 0.00 -0.16 0.24 -0.320 -0.54 -0.06 [3,] 0.08 -0.165 0.16 0.00 0.25 -0.025 -0.11 -0.06 [4,] 0.15 -0.150 0.32 0.16 0.26 0.270 0.32 -0.06 [5,] 0.15 -0.150 0.32 0.16 0.26 0.270 0.32 -0.06 $n [1] 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 [39] 2 2 2 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.2681349 -2.136083 -0.5661723 -0.1434457 -1.2635432 -1.18947579 [2,] 0.2681349 1.506083 -0.5438277 0.3034457 0.2335432 0.03947579 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -1.2648841 -1.4479927 -0.3535432 -0.7690938 -0.3446891 -0.5495733 [2,] 0.9248841 -0.1520073 1.1435432 0.9290938 -0.2553109 1.7295733 [,13] [,14] [,15] [,16] [,17] [,18] [1,] -0.5054083 -1.3329215 -1.2572287 -8.764278 -1.2690938 -0.1467492 [2,] 0.4554083 0.3429215 0.9772287 2.944278 0.4290938 1.5067492 [,19] [,20] [,21] [,22] [,23] [,24] [,25] [1,] -0.3847866 -1.464120 0.1884194 -1.1940227 -0.6168915 -1.149431 -0.8719179 [2,] 0.7547866 1.664120 1.1715806 0.8840227 0.2768915 1.889431 1.4519179 [,26] [,27] [,28] [,29] [,30] [,31] [1,] -1.013641 -0.07516861 -0.596607 -1.299813 -1.4695733 -0.9967492 [2,] 1.533641 0.59516861 1.816607 1.269813 0.8095733 0.6567492 [,32] [,33] [,34] [,35] [,36] [,37] [1,] -0.50344574 -0.6145469 -0.7719626 -0.07641202 -0.1985169 -0.1975132 [2,] -0.05655426 0.2345469 -0.2580374 0.23641202 -0.1314831 0.5175132 [,38] [,39] [,40] [,41] [,42] [1,] -0.3575132 0.2276554 -0.684165 -1.0708167 -0.06 [2,] 0.3575132 0.2723446 0.634165 0.8508167 -0.06 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" "38" "39" "40" "41" "42" Warning message: In bxp(list(stats = c(-0.120000000000001, -0.120000000000001, 0, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/65d751243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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] [1,] 20.30 16.650 NA [2,] 21.90 17.870 NA [3,] 22.76 18.385 NA [4,] 26.15 18.810 NA [5,] 27.38 20.210 NA $n [1] 42 42 0 $conf [,1] [,2] [,3] [1,] 21.72385 18.15583 NA [2,] 23.79615 18.61417 NA $out [1] 21.77 21.89 20.76 16.36 16.32 $group [1] 2 2 2 2 2 $names [1] "1" "2" NA Warning message: In bxp(list(stats = c(20.3, 21.9, 22.76, 26.15, 27.38, 16.65, 17.87, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7j07a1243855507.ps",horizontal=F,pagecentre=F,paper="special",width=11.111111111111,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,] 18.3100 18.3100 18.3100 [2,] 19.9950 19.9950 19.9950 [3,] 20.7175 20.7175 20.7175 [4,] 22.0100 22.0100 22.0100 [5,] 24.3850 24.3850 24.3850 $n [1] 42 42 42 $conf [,1] [,2] [,3] [1,] 20.22624 20.22624 20.22624 [2,] 21.20876 21.20876 21.20876 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" > dev.off() null device 1 > > system("convert tmp/17b8r1243855507.ps tmp/17b8r1243855507.png") > system("convert tmp/2zi681243855507.ps tmp/2zi681243855507.png") > system("convert tmp/3kiq31243855507.ps tmp/3kiq31243855507.png") > system("convert tmp/4frb61243855507.ps tmp/4frb61243855507.png") > system("convert tmp/5b54x1243855507.ps tmp/5b54x1243855507.png") > system("convert tmp/65d751243855507.ps tmp/65d751243855507.png") > system("convert tmp/7j07a1243855507.ps tmp/7j07a1243855507.png") > > > proc.time() user system elapsed 2.277 1.755 3.060