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(67.643 + ,69.371 + ,66.294 + ,70.768 + ,71.774 + ,73.388 + ,74.040 + ,73.238 + ,78.121 + ,69.825 + ,71.099 + ,70.676 + ,69.515 + ,68.246 + ,68.594 + ,70.405 + ,61.223 + ,60.542 + ,61.952 + ,68.173 + ,67.240 + ,68.739 + ,69.234 + ,65.570 + ,67.408 + ,64.630 + ,68.848 + ,73.370 + ,74.292 + ,76.525 + ,74.368 + ,75.674 + ,74.868 + ,79.824 + ,80.022 + ,79.942 + ,80.622 + ,80.079 + ,79.212 + ,80.626 + ,83.551 + ,80.407 + ,85.053 + ,86.399 + ,88.536 + ,89.008 + ,89.652 + ,88.904 + ,87.472 + ,88.631 + ,87.221 + ,88.759 + ,90.127 + ,88.709 + ,90.030 + ,88.697 + ,88.762 + ,89.475 + ,88.936 + ,90.411 + ,90.004 + ,92.725 + ,90.252 + ,93.226 + ,92.575 + ,93.125 + ,95.987 + ,97.175 + ,97.321 + ,98.577 + ,99.026 + ,101.851 + ,99.958 + ,97.875 + ,97.927 + ,95.149 + ,94.551 + ,93.999 + ,93.297 + ,89.901 + ,89.742 + ,87.096 + ,86.863 + ,86.718 + ,80.020 + ,63.483 + ,51.289 + ,44.071 + ,43.654 + ,66.115 + ,84.518 + ,83.395 + ,78.307 + ,80.049 + ,78.346 + ,78.317 + ,75.918 + ,73.739 + ,74.530 + ,74.179 + ,76.974 + ,75.408 + ,73.336 + ,69.210 + ,67.286 + ,64.606 + ,64.159 + ,64.423 + ,66.411 + ,64.270 + ,56.521 + ,50.599 + ,54.751 + ,62.227 + ,63.932 + ,65.391 + ,75.744 + ,74.590 + ,76.035 + ,74.427 + ,73.354 + ,73.081 + ,75.309 + ,75.463 + ,75.910 + ,76.151 + ,76.882 + ,78.632 + ,80.137 + ,82.490 + ,79.896 + ,81.303 + ,79.344 + ,81.355 + ,82.328 + ,79.669 + ,77.249 + ,75.101 + ,72.520 + ,72.438 + ,72.653 + ,71.429 + ,69.189 + ,66.451 + ,63.354 + ,61.379 + ,61.880 + ,62.274 + ,62.429 + ,63.905 + ,63.917 + ,64.295 + ,61.930 + ,60.440 + ,59.353 + ,58.695 + ,60.569 + ,60.386 + ,60.938 + ,61.795 + ,63.304 + ,64.270 + ,63.492 + ,61.333 + ,59.341 + ,58.412 + ,58.725 + ,59.277 + ,58.562 + ,57.858 + ,58.790 + ,58.243 + ,57.044 + ,57.339 + ,59.429 + ,60.575 + ,61.950 + ,61.712) > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = '' > par2 = '' > par1 = '12' > ylab = '' > xlab = '' > main = '' > #'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] 178 > (np <- floor(n / par1)) [1] 14 > 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] 15 15 15 15 15 15 15 15 15 15 14 14 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 67.643 69.515 67.408 80.622 87.472 90.004 99.958 80.020 75.918 66.411 [2,] 69.371 68.246 64.630 80.079 88.631 92.725 97.875 63.483 73.739 64.270 [3,] 66.294 68.594 68.848 79.212 87.221 90.252 97.927 51.289 74.530 56.521 [4,] 70.768 70.405 73.370 80.626 88.759 93.226 95.149 44.071 74.179 50.599 [5,] 71.774 61.223 74.292 83.551 90.127 92.575 94.551 43.654 76.974 54.751 [6,] 73.388 60.542 76.525 80.407 88.709 93.125 93.999 66.115 75.408 62.227 [7,] 74.040 61.952 74.368 85.053 90.030 95.987 93.297 84.518 73.336 63.932 [8,] 73.238 68.173 75.674 86.399 88.697 97.175 89.901 83.395 69.210 65.391 [9,] 78.121 67.240 74.868 88.536 88.762 97.321 89.742 78.307 67.286 75.744 [10,] 69.825 68.739 79.824 89.008 89.475 98.577 87.096 80.049 64.606 74.590 [11,] 71.099 69.234 80.022 89.652 88.936 99.026 86.863 78.346 64.159 76.035 [12,] 70.676 65.570 79.942 88.904 90.411 101.851 86.718 78.317 64.423 74.427 [,11] [,12] [,13] [,14] [,15] [1,] 73.354 79.344 63.354 60.569 58.562 [2,] 73.081 81.355 61.379 60.386 57.858 [3,] 75.309 82.328 61.880 60.938 58.790 [4,] 75.463 79.669 62.274 61.795 58.243 [5,] 75.910 77.249 62.429 63.304 57.044 [6,] 76.151 75.101 63.905 64.270 57.339 [7,] 76.882 72.520 63.917 63.492 59.429 [8,] 78.632 72.438 64.295 61.333 60.575 [9,] 80.137 72.653 61.930 59.341 61.950 [10,] 82.490 71.429 60.440 58.412 61.712 [11,] 79.896 69.189 59.353 58.725 NA [12,] 81.303 66.451 58.695 59.277 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1.728 -1.269 -2.778 -0.543 1.159 2.721 -2.083 -16.537 -2.179 -2.141 [2,] -3.077 0.348 4.218 -0.867 -1.410 -2.473 0.052 -12.194 0.791 -7.749 [3,] 4.474 1.811 4.522 1.414 1.538 2.974 -2.778 -7.218 -0.351 -5.922 [4,] 1.006 -9.182 0.922 2.925 1.368 -0.651 -0.598 -0.417 2.795 4.152 [5,] 1.614 -0.681 2.233 -3.144 -1.418 0.550 -0.552 22.461 -1.566 7.476 [6,] 0.652 1.410 -2.157 4.646 1.321 2.862 -0.702 18.403 -2.072 1.705 [7,] -0.802 6.221 1.306 1.346 -1.333 1.188 -3.396 -1.123 -4.126 1.459 [8,] 4.883 -0.933 -0.806 2.137 0.065 0.146 -0.159 -5.088 -1.924 10.353 [9,] -8.296 1.499 4.956 0.472 0.713 1.256 -2.646 1.742 -2.680 -1.154 [10,] 1.274 0.495 0.198 0.644 -0.539 0.449 -0.233 -1.703 -0.447 1.445 [11,] -0.423 -3.664 -0.080 -0.748 1.475 2.825 -0.145 -0.029 0.264 -1.608 [12,] -1.161 1.838 0.680 -1.432 -0.407 -1.893 -6.698 -2.399 1.988 -1.073 [,11] [,12] [,13] [,14] [,15] [1,] -0.273 2.011 -1.975 -0.183 -0.704 [2,] 2.228 0.973 0.501 0.552 0.932 [3,] 0.154 -2.659 0.394 0.857 -0.547 [4,] 0.447 -2.420 0.155 1.509 -1.199 [5,] 0.241 -2.148 1.476 0.966 0.295 [6,] 0.731 -2.581 0.012 -0.778 2.090 [7,] 1.750 -0.082 0.378 -2.159 1.146 [8,] 1.505 0.215 -2.365 -1.992 1.375 [9,] 2.353 -1.224 -1.490 -0.929 -0.238 [10,] -2.594 -2.240 -1.087 0.313 NA [11,] 1.407 -2.738 -0.658 0.552 NA [12,] -1.959 -3.097 1.874 -0.715 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/19ic21290677660.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/22st51290677660.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/32st51290677660.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/42st51290677660.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,] 58.5620 57.8580 51.289 44.0710 43.654 57.3390 59.4290 60.575 59.3410 [2,] 66.9095 63.8765 61.409 62.0345 61.826 64.0875 63.9245 66.782 67.2630 [3,] 73.3540 69.3710 68.848 73.3700 74.292 75.1010 74.0400 73.238 75.7440 [4,] 80.3210 80.7170 80.770 80.1475 80.400 78.4660 84.7855 84.897 84.3365 [5,] 99.9580 97.8750 97.927 95.1490 94.551 93.9990 95.9870 97.175 97.3210 [,10] [,11] [,12] [1,] 58.4120 58.7250 58.695 [2,] 66.6725 69.1890 65.570 [3,] 74.5900 77.1905 76.372 [4,] 84.7930 86.8630 86.718 [5,] 98.5770 99.0260 101.851 $n [1] 15 15 15 15 15 15 15 15 15 15 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 67.88272 62.50085 60.9496 65.98072 66.71466 69.23523 65.52967 65.84791 [2,] 78.82528 76.24115 76.7464 80.75928 81.86934 80.96677 82.55033 80.62809 [,9] [,10] [,11] [,12] [1,] 68.7788 67.19767 69.72725 67.44178 [2,] 82.7092 81.98233 84.65375 85.30222 $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(58.562, 66.9095, 73.354, 80.321, 99.958, 57.858, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5u1a81290677660.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] [1,] -2.778 -3.0770 -5.9220 -2.4200 -3.1440 -2.5810 -4.126 -5.0880 -2.6800 [2,] -2.112 -1.9415 -1.6030 -0.6245 -1.0495 -0.7400 -1.228 -1.4285 -1.3570 [3,] -0.704 0.3480 0.3940 0.4470 0.2950 0.7310 0.378 0.0650 -0.2380 [4,] 0.488 0.8615 1.6745 1.4385 1.5450 1.8975 1.326 1.4400 1.3775 [5,] 2.721 4.2180 4.5220 4.1520 2.2330 4.6460 1.750 4.8830 4.9560 [,10] [,11] [,12] [1,] -2.5940 -1.6080 -3.097 [2,] -1.0870 -0.7480 -1.959 [3,] -0.0175 -0.1125 -1.117 [4,] 0.4950 0.5520 0.680 [5,] 1.4450 1.4750 1.988 $n [1] 15 15 15 15 15 15 15 15 15 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.7646810 -0.7954958 -0.94307 -0.3946096 -0.7634373 -0.3449793 -0.6639151 [2,] 0.3566810 1.4914958 1.73107 1.2886096 1.3534373 1.8069793 1.4199151 [,8] [,9] [,10] [,11] [,12] [1,] -1.105217 -1.3535509 -0.6855355 -0.6614546 -2.231377820 [2,] 1.235217 0.8775509 0.6505355 0.4364546 -0.002622180 $out [1] -16.537 -12.194 -7.749 -7.218 -9.182 22.461 7.476 18.403 6.221 [10] 10.353 -8.296 -3.664 2.825 -2.738 -6.698 $group [1] 1 2 2 3 4 5 5 6 7 8 9 11 11 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2.77800000000001, -2.112, -0.704, 0.488000000000003, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6u1a81290677660.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] [1,] 66.2940 60.5420 64.6300 79.2120 87.2210 90.0040 86.718 43.6540 64.1590 [2,] 69.5980 63.7610 71.1090 80.5145 88.6640 92.6500 88.419 57.3860 65.9460 [3,] 70.9335 68.2095 74.6180 84.3020 88.7605 94.6065 93.648 78.3120 73.5375 [4,] 73.3130 68.9865 78.1745 88.7200 89.7525 97.9490 96.512 80.0345 74.9690 [5,] 78.1210 70.4050 80.0220 89.6520 90.4110 101.8510 99.958 84.5180 76.9740 [,10] [,11] [,12] [,13] [,14] [,15] [1,] 50.5990 73.0810 66.4510 58.6950 58.4120 57.044 [2,] 59.3740 75.3860 71.9335 60.9095 59.3090 57.858 [3,] 64.8305 76.5165 73.8770 62.1020 60.7535 58.676 [4,] 74.5085 80.0165 79.5065 63.6295 62.5495 60.575 [5,] 76.0350 82.4900 82.3280 64.2950 64.2700 61.950 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 69.23906 65.82611 71.39538 80.55942 88.26403 92.18959 89.95673 67.98187 [2,] 72.62794 70.59289 77.84062 88.04458 89.25697 97.02341 97.33927 88.64213 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 69.42205 57.92755 74.4045 70.4229 60.86139 59.27549 57.31848 [2,] 77.65295 71.73345 78.6285 77.3311 63.34261 62.23151 60.03352 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" NA Warning message: In bxp(list(stats = c(66.294, 69.598, 70.9335, 73.313, 78.121, 60.542, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/75aat1290677660.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,] 71.90640 73.2380 71.08950 [2,] 72.56803 73.2960 71.19487 [3,] 75.09690 74.1660 73.98513 [4,] 75.94033 75.4225 75.76625 [5,] 76.46679 77.1905 77.17650 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 73.55877 73.19609 71.90009 [2,] 76.63503 75.13591 76.07016 $out [1] 69.371 68.848 $group [1] 2 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(71.9064, 72.5680333333333, 75.0969, 75.9403333333333, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/19ic21290677660.ps tmp/19ic21290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/22st51290677660.ps tmp/22st51290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/32st51290677660.ps tmp/32st51290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/42st51290677660.ps tmp/42st51290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/5u1a81290677660.ps tmp/5u1a81290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/6u1a81290677660.ps tmp/6u1a81290677660.png",intern=TRUE)) character(0) > try(system("convert tmp/75aat1290677660.ps tmp/75aat1290677660.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.940 1.446 2.229