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(-0.350799688370536 + ,-0.140769858461987 + ,-0.176994002962568 + ,0.562290817956577 + ,0.144201690134944 + ,-0.111030533067628 + ,-0.243045038133344 + ,-0.151508019282699 + ,-0.396513842133135 + ,0.0329753101728757 + ,-0.0833088349918484 + ,0.156956311022783 + ,0.523621600794454 + ,-0.216145977751341 + ,-0.208014642793551 + ,-0.341520041804527 + ,-0.114593780635957 + ,-0.558823753124076 + ,0.155427432046773 + ,0.138052779163656 + ,-0.4090756185683 + ,0.110316654295769 + ,-0.250998024914370 + ,0.329481615878391 + ,0.346878506478169 + ,0.323299830417798 + ,-0.0310312875854553 + ,-0.84926656008146 + ,-0.77680329705359 + ,0.250909338924582 + ,0.0831309127438593 + ,0.435693250316565 + ,-0.130263895758751 + ,-0.232129238451711 + ,0.165342569504054 + ,0.142248180166828 + ,-0.280455237088636 + ,0.155402490351535 + ,-0.142784384867873 + ,0.28524670566859 + ,0.157418555269204 + ,-0.172706911536110 + ,-0.0149928813013760 + ,0.306296873369522 + ,-0.345207720981392 + ,0.722409122191829 + ,-0.158636568880998 + ,-0.431031704131918 + ,-0.308095207623409 + ,0.447320732827413 + ,0.447799821943335 + ,0.0214235496574372 + ,0.152471875346584 + ,0.163942423990300 + ,-0.00498382593447193 + ,0.308356992481182 + ,-0.230794586105762 + ,-0.0542129596699259 + ,-0.137751473869932 + ,0.137775171895457 + ,-0.208611181469940 + ,0.530877388565287 + ,-0.307691326359942 + ,0.287643632777539 + ,0.239698181689028 + ,-0.0324235808931527 + ,-0.110760321621385 + ,0.187857496084540 + ,-0.629717422499043 + ,-0.464897710438974 + ,0.381446921054368 + ,0.388355710424064 + ,-0.59380461268784 + ,0.456215309938507 + ,-0.35419239371042 + ,0.278101258603134 + ,-0.32386448817015 + ,-0.115433965509661 + ,-0.0483653154188056 + ,-0.0424966476262035 + ,-0.179962353470623 + ,0.634520247547924 + ,1.07894182435992 + ,0.395987866990447 + ,-0.219758520860607 + ,0.361998505460662 + ,0.660014172078838 + ,-0.331005583090054 + ,0.695781527796754 + ,0.646124350414131 + ,-0.0835481139905882 + ,-0.763493350940545 + ,0.0627068862169823 + ,-0.0788278540379063 + ,-1.01076620067800 + ,-0.208872220642831 + ,-0.515267639977992 + ,0.319253920126987 + ,-0.739240050461461 + ,-0.0632470686715278 + ,-0.182648154761123 + ,-0.242076757818195 + ,-0.704413262583499 + ,0.851791451893186 + ,-0.183365626294880) > 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] 105 > (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] 9 9 9 9 9 9 9 9 9 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.35079969 0.5236216 0.34687851 -0.28045524 -0.308095208 -0.20861118 [2,] -0.14076986 -0.2161460 0.32329983 0.15540249 0.447320733 0.53087739 [3,] -0.17699400 -0.2080146 -0.03103129 -0.14278438 0.447799822 -0.30769133 [4,] 0.56229082 -0.3415200 -0.84926656 0.28524671 0.021423550 0.28764363 [5,] 0.14420169 -0.1145938 -0.77680330 0.15741856 0.152471875 0.23969818 [6,] -0.11103053 -0.5588238 0.25090934 -0.17270691 0.163942424 -0.03242358 [7,] -0.24304504 0.1554274 0.08313091 -0.01499288 -0.004983826 -0.11076032 [8,] -0.15150802 0.1380528 0.43569325 0.30629687 0.308356992 0.18785750 [9,] -0.39651384 -0.4090756 -0.13026390 -0.34520772 -0.230794586 -0.62971742 [10,] 0.03297531 0.1103167 -0.23212924 0.72240912 -0.054212960 -0.46489771 [11,] -0.08330883 -0.2509980 0.16534257 -0.15863657 -0.137751474 0.38144692 [12,] 0.15695631 0.3294816 0.14224818 -0.43103170 0.137775172 0.38835571 [,7] [,8] [,9] [1,] -0.59380461 -0.21975852 -0.51526764 [2,] 0.45621531 0.36199851 0.31925392 [3,] -0.35419239 0.66001417 -0.73924005 [4,] 0.27810126 -0.33100558 -0.06324707 [5,] -0.32386449 0.69578153 -0.18264815 [6,] -0.11543397 0.64612435 -0.24207676 [7,] -0.04836532 -0.08354811 -0.70441326 [8,] -0.04249665 -0.76349335 0.85179145 [9,] -0.17996235 0.06270689 -0.18336563 [10,] 0.63452025 -0.07882785 NA [11,] 1.07894182 -1.01076620 NA [12,] 0.39598787 -0.20887222 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.21002983 -0.739767579 -0.02357868 0.4358577 0.7554159405 0.73948857 [2,] -0.03622414 0.008131335 -0.35433112 -0.2981869 0.0004790891 -0.83856871 [3,] 0.73928482 -0.133505399 -0.81823527 0.4280311 -0.4263762723 0.59533496 [4,] -0.41808913 0.226926261 0.07246326 -0.1278282 0.1310483257 -0.04794545 [5,] -0.25523222 -0.444229972 1.02771264 -0.3301255 0.0114705486 -0.27212176 [6,] -0.13201451 0.714251185 -0.16777843 0.1577140 -0.1689262499 -0.07833674 [7,] 0.09153702 -0.017374653 0.35256234 0.3212898 0.3133408184 0.29861782 [8,] -0.24500582 -0.547128398 -0.56595715 -0.6515046 -0.5391515786 -0.81757492 [9,] 0.42948915 0.519392273 -0.10186534 1.0676168 0.1765816264 0.16481971 [10,] -0.11628415 -0.361314679 0.39747181 -0.8810457 -0.0835385142 0.84634463 [11,] 0.24026515 0.580479641 -0.02309439 -0.2723951 0.2755266458 0.00690879 [12,] 0.36666529 0.017396891 -0.42270342 0.1229365 -0.3463863534 -0.98216032 [,7] [,8] [,9] [1,] 1.050019923 0.58175703 0.8345216 [2,] -0.810407704 0.29801567 -1.0584940 [3,] 0.632293652 -0.99101976 0.6759930 [4,] -0.601965747 1.02678711 -0.1194011 [5,] 0.208430523 -0.04965718 -0.0594286 [6,] 0.067068650 -0.72967246 -0.4623365 [7,] 0.005868668 -0.67994524 1.5562047 [8,] -0.137465706 0.82620024 -1.0351571 [9,] 0.814482601 -0.14153474 NA [10,] 0.444421577 -0.93193835 NA [11,] -0.682953957 0.80189398 NA [12,] -0.615746388 -0.30639542 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/1b0dt1260910975.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/2gwab1260910975.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/3elpx1260910975.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/47cg91260910975.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] [1,] -0.5152676 -0.2161460 -0.35419239 -0.84926656 -0.3238645 -0.5588238 [2,] -0.3507997 0.1554025 -0.30769133 -0.33100558 -0.1826482 -0.1727069 [3,] -0.2804552 0.3232998 -0.17699400 0.02142355 0.1442017 -0.1110305 [4,] -0.2086112 0.4473207 -0.03103129 0.28524671 0.1574186 0.1639424 [5,] -0.2086112 0.5308774 -0.03103129 0.56229082 0.2396982 0.6461244 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.243045038 -0.15150802 -0.62971742 -0.46489771 -0.2509980 -0.43103170 [2,] -0.110760322 -0.04249665 -0.39651384 -0.15547855 -0.2048173 -0.03554852 [3,] -0.048365315 0.18785750 -0.23079459 -0.01061882 -0.1105302 0.14960225 [4,] -0.004983826 0.30835699 -0.17996235 0.37241845 0.2733947 0.35891866 [5,] 0.083130913 0.43569325 0.06270689 0.72240912 0.3814469 0.39598787 $n [1] 9 9 9 9 9 9 9 9 9 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.3553412 0.1695562 -0.32270162 -0.3031360 -0.03490011 -0.28833252 [2,] -0.2055693 0.4770434 -0.03128638 0.3459831 0.32330349 0.06627145 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.104074270 0.003074579 -0.3448450 -0.3055097 -0.3776663 -0.07075279 [2,] 0.007343639 0.372640413 -0.1167441 0.2842720 0.1566060 0.36995728 $out [1] 0.5236216 0.3468785 -0.5938046 0.4477998 0.6600142 -0.7392401 [7] -0.7768033 0.6957815 0.1554274 -0.7044133 -0.7634934 0.8517915 [13] 1.0789418 -1.0107662 $group [1] 1 1 1 3 3 3 5 5 7 7 8 8 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.515267639977992, -0.350799688370536, -0.280455237088636, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5prj91260910975.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] [1,] -0.02357868 -1.0584939706 -0.9910198 -0.41808913 -0.44422997 -0.46233650 [2,] 0.21002983 -0.8104077036 -0.4263763 -0.12782815 -0.27212176 -0.16892625 [3,] 0.58175703 -0.2981868752 0.4280311 -0.04794545 -0.05942860 -0.13201451 [4,] 0.75541594 0.0004790891 0.6322937 0.13104833 0.01147055 0.06706865 [5,] 1.05001992 0.2980156666 0.7392848 0.22692626 0.20843052 0.15771403 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.017374653 -1.0351571 -0.14153474 -0.93193835 -0.6829540 -0.9821603 [2,] 0.005868668 -0.6515046 0.03147718 -0.62118019 -0.1477448 -0.5192249 [3,] 0.298617818 -0.5471284 0.30303539 -0.09991133 0.1235870 -0.3263909 [4,] 0.321289755 -0.2450058 0.66693744 0.42094669 0.4280031 0.0701667 [5,] 0.352562338 -0.1374657 1.06761684 0.84634463 0.8018940 0.3666653 $n [1] 9 9 9 9 9 9 9 9 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.2945203 -0.7252539 -0.1295351 -0.18428706 -0.20878722 -0.256305152 [2,] 0.8689937 0.1288802 0.9855973 0.08839616 0.08993001 -0.007723858 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.1324960 -0.7612178 -0.05194183 -0.6820584 -0.1980341 -0.655633494 [2,] 0.4647396 -0.3330390 0.65801261 0.4822357 0.4452080 0.002851721 $out [1] -0.7397676 -0.6019657 1.0267871 1.0277126 0.7142512 -0.7296725 -0.6799452 [8] 1.5562047 0.8262002 $group [1] 1 4 4 5 6 6 7 7 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.023578676060371, 0.210029829908549, 0.581757026321269, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/607qk1260910975.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] [1,] -0.3965138 -0.5588238 -0.8492666 -0.43103170 -0.30809521 -0.62971742 [2,] -0.2100195 -0.2962590 -0.1811966 -0.22658107 -0.09598222 -0.25815125 [3,] -0.1259002 -0.1613042 0.1126895 -0.07888863 0.07959936 0.07771696 [4,] 0.0885885 0.1467401 0.2871046 0.22133263 0.23614971 0.33454528 [5,] 0.1569563 0.5236216 0.4356933 0.72240912 0.44779982 0.53087739 [,7] [,8] [,9] [1,] -0.59380461 -1.01076620 -0.73924005 [2,] -0.25191342 -0.27538205 -0.51526764 [3,] -0.04543098 -0.08118798 -0.18336563 [4,] 0.42610159 0.50406143 -0.06324707 [5,] 1.07894182 0.69578153 0.31925392 $n [1] 12 12 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.26209732 -0.36335909 -0.1009058 -0.2831851 -0.07188824 -0.1926159 [2,] 0.01029693 0.04075067 0.3262848 0.1254078 0.23108696 0.3480498 [,7] [,8] [,9] [1,] -0.3546782 -0.4366974 -0.42142979 [2,] 0.2638162 0.2743214 0.05469854 $out [1] 0.5622908 0.8517915 $group [1] 1 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(-0.396513842133135, -0.210019520547956, -0.125900195764807, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7xdcu1260910975.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,] -0.271354909 -0.28045524 -0.288238098 [2,] -0.101315806 -0.14401227 -0.113616690 [3,] -0.009334961 -0.02949207 -0.008498522 [4,] 0.098815906 0.14690197 0.097518674 [5,] 0.248605816 0.32329983 0.301361612 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.10061638 -0.1621800 -0.10479878 [2,] 0.08194646 0.1031959 0.08780174 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.271354908843878, -0.101315806050663, -0.0093349611640666, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > try(system("convert tmp/1b0dt1260910975.ps tmp/1b0dt1260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/2gwab1260910975.ps tmp/2gwab1260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/3elpx1260910975.ps tmp/3elpx1260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/47cg91260910975.ps tmp/47cg91260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/5prj91260910975.ps tmp/5prj91260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/607qk1260910975.ps tmp/607qk1260910975.png",intern=TRUE)) character(0) > try(system("convert tmp/7xdcu1260910975.ps tmp/7xdcu1260910975.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.223 0.972 1.666