R version 2.9.0 (2009-04-17) 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(100 + ,101 + ,102.01 + ,103.0301 + ,104.060401 + ,105.101005 + ,106.1520151 + ,107.2135352 + ,108.2856706 + ,109.3685273 + ,110.4622125 + ,111.5668347 + ,112.682503 + ,113.809328 + ,114.9474213 + ,116.0968955 + ,117.2578645 + ,118.4304431 + ,119.6147476 + ,120.810895 + ,122.019004 + ,123.239194 + ,124.471586 + ,125.7163018 + ,126.9734649 + ,128.2431995 + ,129.5256315 + ,130.8208878 + ,132.1290967 + ,133.4503877 + ,134.7848915 + ,136.1327404 + ,137.4940679 + ,138.8690085 + ,140.2576986 + ,141.6602756 + ,143.0768784 + ,144.5076471 + ,145.9527236 + ,147.4122509 + ,148.8863734 + ,150.3752371 + ,151.8789895 + ,153.3977794 + ,154.9317572 + ,156.4810747 + ,158.0458855 + ,159.6263443 + ,161.2226078 + ,162.8348338 + ,164.4631822 + ,166.107814 + ,167.7688921 + ,169.4465811 + ,171.1410469 + ,172.8524573 + ,174.5809819 + ,176.3267917 + ,178.0900597 + ,179.8709603 + ,181.6696699 + ,183.4863666 + ,185.3212302 + ,187.1744425 + ,189.0461869 + ,190.9366488 + ,192.8460153 + ,194.7744755 + ,196.7222202 + ,198.6894424 + ,200.6763368 + ,202.6831002 + ,204.7099312 + ,800 + ,206.7570305 + ,208.8246008 + ,210.9128468 + ,213.0219753 + ,215.1521951 + ,217.303717 + ,219.4767542 + ,221.6715217 + ,223.8882369 + ,226.1271193 + ,228.3883905 + ,230.6722744 + ,232.9789971 + ,235.3087871 + ,237.661875 + ,240.0384937 + ,242.4388787 + ,244.8632675 + ,247.3119001 + ,249.7850191 + ,252.2828693 + ,254.805698) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0000 112.6825 126.9735 143.0769 161.2226 181.6697 204.7099 228.3884 [2,] 101.0000 113.8093 128.2432 144.5076 162.8348 183.4864 800.0000 230.6723 [3,] 102.0100 114.9474 129.5256 145.9527 164.4632 185.3212 206.7570 232.9790 [4,] 103.0301 116.0969 130.8209 147.4123 166.1078 187.1744 208.8246 235.3088 [5,] 104.0604 117.2579 132.1291 148.8864 167.7689 189.0462 210.9128 237.6619 [6,] 105.1010 118.4304 133.4504 150.3752 169.4466 190.9366 213.0220 240.0385 [7,] 106.1520 119.6147 134.7849 151.8790 171.1410 192.8460 215.1522 242.4389 [8,] 107.2135 120.8109 136.1327 153.3978 172.8525 194.7745 217.3037 244.8633 [9,] 108.2857 122.0190 137.4941 154.9318 174.5810 196.7222 219.4768 247.3119 [10,] 109.3685 123.2392 138.8690 156.4811 176.3268 198.6894 221.6715 249.7850 [11,] 110.4622 124.4716 140.2577 158.0459 178.0901 200.6763 223.8882 252.2829 [12,] 111.5668 125.7163 141.6603 159.6263 179.8710 202.6831 226.1271 254.8057 [,9] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.000000 1.126825 1.269735 1.430769 1.612226 1.816697 595.290069 [2,] 1.010000 1.138093 1.282432 1.445076 1.628348 1.834864 -593.242969 [3,] 1.020100 1.149474 1.295256 1.459527 1.644632 1.853212 2.067570 [4,] 1.030301 1.160969 1.308209 1.474122 1.661078 1.871744 2.088246 [5,] 1.040604 1.172579 1.321291 1.488864 1.677689 1.890462 2.109128 [6,] 1.051010 1.184304 1.334504 1.503752 1.694466 1.909367 2.130220 [7,] 1.061520 1.196147 1.347849 1.518790 1.711410 1.928460 2.151522 [8,] 1.072135 1.208109 1.361328 1.533978 1.728525 1.947745 2.173037 [9,] 1.082857 1.220190 1.374941 1.549318 1.745810 1.967222 2.194768 [10,] 1.093685 1.232392 1.388690 1.564811 1.763268 1.986894 2.216715 [11,] 1.104622 1.244716 1.402577 1.580459 1.780901 2.006763 2.238882 [12,] 1.115668 1.257163 1.416603 1.596263 1.798710 2.026831 2.261271 [,8] [,9] [1,] 2.283884 NA [2,] 2.306723 NA [3,] 2.329790 NA [4,] 2.353088 NA [5,] 2.376619 NA [6,] 2.400385 NA [7,] 2.424389 NA [8,] 2.448633 NA [9,] 2.473119 NA [10,] 2.497850 NA [11,] 2.522829 NA [12,] 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/1ozrj1289796734.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/2ozrj1289796734.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/3z98m1289796734.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/4z98m1289796734.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] [1,] 100.0000 101.0000 102.0100 103.0301 104.0604 105.1010 106.1520 107.2135 [2,] 119.8280 121.0263 122.2365 123.4589 124.6935 125.9404 127.1998 128.4718 [3,] 152.1497 153.6712 155.2080 156.7600 158.3276 159.9109 161.5100 163.1251 [4,] 193.1898 207.0793 196.0391 197.9995 199.9795 201.9793 203.9991 206.0391 [5,] 228.3884 230.6723 232.9790 235.3088 237.6619 240.0385 242.4389 244.8633 [,9] [,10] [,11] [,12] [1,] 108.2857 109.3685 110.4622 111.5668 [2,] 129.7565 131.0541 132.3646 133.6883 [3,] 164.7564 166.4039 168.0680 169.7487 [4,] 208.0995 210.1805 212.2823 214.4051 [5,] 247.3119 249.7850 252.2829 254.8057 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 111.1688 105.6008 113.9808 115.1206 116.2718 117.4345 118.6088 119.7949 [2,] 193.1307 201.7417 196.4351 198.3995 200.3835 202.3873 204.4112 206.4553 [,9] [,10] [,11] [,12] [1,] 120.9929 122.2028 123.4248 124.6591 [2,] 208.5199 210.6051 212.7111 214.8382 $out [1] 800 $group [1] 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100, 119.82798395, 152.1497431, 193.18980055, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5z98m1289796734.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] [1,] 1.000000 1.010000 1.020100 1.030301 1.040604 1.051010 1.061520 1.072135 [2,] 1.198280 1.074047 1.222365 1.234589 1.246935 1.259404 1.271998 1.284718 [3,] 1.521497 1.363754 1.552080 1.567600 1.583276 1.599109 1.615100 1.631251 [4,] 2.050290 1.731606 1.960391 1.979995 1.999795 2.019793 2.039991 2.060391 [5,] 2.283884 2.306723 2.329790 2.353088 2.376619 2.400385 2.424389 2.448633 [,9] [,10] [,11] [,12] [1,] 1.082857 1.093685 1.104622 1.115668 [2,] 1.297565 1.310541 1.323646 1.336883 [3,] 1.647564 1.664039 1.680680 1.596263 [4,] 2.080995 2.101805 2.122823 1.912770 [5,] 2.473119 2.497850 2.522829 2.261271 $n [1] 8 8 8 8 8 8 8 8 8 8 8 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1.045552 0.9964322 1.139808 1.151206 1.162718 1.174345 1.186088 1.197949 [2,] 1.997443 1.7310763 1.964351 1.983995 2.003835 2.023873 2.044112 2.064553 [,9] [,10] [,11] [,12] [1,] 1.209929 1.222028 1.234248 1.252353 [2,] 2.085199 2.106051 2.127111 1.940174 $out [1] 595.2901 -593.2430 $group [1] 1 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1, 1.19827980000000, 1.52149735000002, 2.05029029999999, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ai771289796734.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,] 100.0000 112.6825 126.9735 143.0769 161.2226 181.6697 204.7099 228.3884 [2,] 102.5200 115.5222 130.1733 146.6825 165.2855 186.2478 209.8687 234.1439 [3,] 105.6265 119.0226 134.1176 151.1271 170.2938 191.8913 216.2280 241.2387 [4,] 108.8271 122.6291 138.1815 155.7064 175.4539 197.7058 222.7799 248.5485 [5,] 111.5668 125.7163 141.6603 159.6263 179.8710 202.6831 226.1271 254.8057 [,9] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 102.7498 115.7811 130.4650 147.0112 165.6559 186.6653 210.3391 234.6687 [2,] 108.5032 122.2641 137.7703 155.2430 174.9317 197.1174 222.1168 247.8087 [,9] [1,] NA [2,] NA $out [1] 800 $group [1] 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7ai771289796734.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,] 157.3404 152.1497 155.4152 [2,] 162.6562 155.9840 159.9111 [3,] 167.5849 160.7105 163.8921 [4,] 172.6628 165.5802 168.8582 [5,] 175.2571 169.7487 173.1092 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 163.0208 156.3336 159.8113 [2,] 172.1489 165.0873 167.9729 $out [1] 233.0692 $group [1] 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(157.3404307125, 162.6562071875, 167.5848529375, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ozrj1289796734.ps tmp/1ozrj1289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/2ozrj1289796734.ps tmp/2ozrj1289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/3z98m1289796734.ps tmp/3z98m1289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/4z98m1289796734.ps tmp/4z98m1289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/5z98m1289796734.ps tmp/5z98m1289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/6ai771289796734.ps tmp/6ai771289796734.png",intern=TRUE)) character(0) > try(system("convert tmp/7ai771289796734.ps tmp/7ai771289796734.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.205 0.950 1.691