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(4.031636,3.702076,3.056167,3.280707,2.984728,3.693712,3.226317,2.190349,2.599515,3.080288,2.929672,2.922548,3.234943,2.983081,3.284389,3.806511,3.784579,2.645654,3.092081,3.204859,3.107225,3.466909,2.984404,3.218072,2.82731,3.182049,2.236319,2.033218,1.644804,1.627971,1.677559,2.330828,2.493615,2.257172,2.655517,2.298655,2.600402,3.04523,2.790583,3.227052,2.967479,2.938817,3.277961,3.423985,3.072646,2.754253,2.910431,3.174369,3.068387,3.089543,2.906654,2.931161,3.02566,2.939551,2.691019,3.19812,3.07639,2.863873,3.013802,3.053364,2.864753,3.057062,2.959365,3.252258,3.602988,3.497704,3.296867,3.602417,3.3001,3.40193,3.502591,3.402348,3.498551,3.199823,2.700064,2.801034,2.898628,2.800854,2.399942,2.402724,2.202331,2.102594,1.798293,1.202484,1.400201,1.200832,1.298083,1.099742,1.001377,0.8361743) > 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] 90 > (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] 8 8 8 8 8 8 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 4.031636 3.234943 2.827310 2.600402 3.068387 2.864753 3.498551 1.4002010 [2,] 3.702076 2.983081 3.182049 3.045230 3.089543 3.057062 3.199823 1.2008320 [3,] 3.056167 3.284389 2.236319 2.790583 2.906654 2.959365 2.700064 1.2980830 [4,] 3.280707 3.806511 2.033218 3.227052 2.931161 3.252258 2.801034 1.0997420 [5,] 2.984728 3.784579 1.644804 2.967479 3.025660 3.602988 2.898628 1.0013770 [6,] 3.693712 2.645654 1.627971 2.938817 2.939551 3.497704 2.800854 0.8361743 [7,] 3.226317 3.092081 1.677559 3.277961 2.691019 3.296867 2.399942 NA [8,] 2.190349 3.204859 2.330828 3.423985 3.198120 3.602417 2.402724 NA [9,] 2.599515 3.107225 2.493615 3.072646 3.076390 3.300100 2.202331 NA [10,] 3.080288 3.466909 2.257172 2.754253 2.863873 3.401930 2.102594 NA [11,] 2.929672 2.984404 2.655517 2.910431 3.013802 3.502591 1.798293 NA [12,] 2.922548 3.218072 2.298655 3.174369 3.053364 3.402348 1.202484 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -0.329560 -0.251862 0.354739 0.444828 0.021156 0.192309 -0.298728 [2,] -0.645909 0.301308 -0.945730 -0.254647 -0.182889 -0.097697 -0.499759 [3,] 0.224540 0.522122 -0.203101 0.436469 0.024507 0.292893 0.100970 [4,] -0.295979 -0.021932 -0.388414 -0.259573 0.094499 0.350730 0.097594 [5,] 0.708984 -1.138925 -0.016833 -0.028662 -0.086109 -0.105284 -0.097774 [6,] -0.467395 0.446427 0.049588 0.339144 -0.248532 -0.200837 -0.400912 [7,] -1.035968 0.112778 0.653269 0.146024 0.507101 0.305550 0.002782 [8,] 0.409166 -0.097634 0.162787 -0.351339 -0.121730 -0.302317 -0.200393 [9,] 0.480773 0.359684 -0.236443 -0.318393 -0.212517 0.101830 -0.099737 [10,] -0.150616 -0.482505 0.398345 0.156178 0.149929 0.100661 -0.304301 [11,] -0.007124 0.233668 -0.356862 0.263938 0.039562 -0.100243 -0.595809 [12,] 0.312395 -0.390762 0.301747 -0.105982 -0.188611 0.096203 0.197717 [,8] [1,] -0.1993690 [2,] 0.0972510 [3,] -0.1983410 [4,] -0.0983650 [5,] -0.1652027 [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] 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/1m3jb1293620372.ps",horizontal=F,onefile=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/2m3jb1293620372.ps",horizontal=F,onefile=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/3fu1v1293620372.ps",horizontal=F,onefile=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/4fu1v1293620372.ps",horizontal=F,onefile=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,] 2.600402 2.983081 2.236319 2.033218 1.001377 0.8361743 1.677559 2.190349 [2,] 2.713856 3.014156 2.468191 2.417126 2.271716 2.1368125 2.545481 2.366776 [3,] 2.966570 3.073303 2.848618 3.079106 2.976103 2.8698355 3.092081 3.198120 [4,] 3.366747 3.190936 3.007766 3.266482 3.314324 3.2186275 3.252139 3.314422 [5,] 4.031636 3.199823 3.284389 3.806511 3.784579 3.6937120 3.296867 3.602417 [,9] [,10] [,11] [,12] [1,] 2.202331 2.102594 2.655517 2.298655 [2,] 2.546565 2.505713 2.782974 2.610602 [3,] 3.072646 2.863873 2.929672 3.053364 [4,] 3.091807 3.241109 2.999103 3.196220 [5,] 3.300100 3.466909 3.013802 3.402348 $n [1] 8 8 8 8 8 8 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2.601856 2.974550 2.547204 2.604644 2.393688 2.265518 2.670076 2.632201 [2,] 3.331284 3.172055 3.150033 3.553569 3.558519 3.474153 3.514086 3.764039 [,9] [,10] [,11] [,12] [1,] 2.747036 2.424706 2.800603 2.703642 [2,] 3.398256 3.303040 3.058741 3.403086 $out [1] 1.400201 3.702076 1.200832 1.298083 1.099742 3.502591 1.798293 1.202484 $group [1] 1 2 2 3 4 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.600402, 2.713856, 2.96657, 3.366747, 4.031636, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5fu1v1293620372.ps",horizontal=F,onefile=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] [1,] -0.3295600 -0.945730 -0.203101 -0.3884140 -0.1652027 -0.467395 0.0027820 [2,] -0.2752950 -0.572834 -0.086917 -0.2777760 -0.1352433 -0.324722 0.0577800 [3,] -0.0891065 -0.218768 0.162755 -0.0601485 -0.0919415 -0.200837 0.1460240 [4,] 0.2735240 -0.000223 0.364681 0.0960465 -0.0227475 0.194366 0.4063255 [5,] 0.4448280 0.301308 0.522122 0.3507300 -0.0168330 0.446427 0.6532690 [,8] [,9] [,10] [,11] [,12] [1,] -0.3513390 -0.318393 -0.4825050 -0.5958090 -0.3907620 [2,] -0.2513550 -0.224480 -0.2274585 -0.2285525 -0.1472965 [3,] -0.1217300 -0.099737 0.1006610 -0.0071240 0.0962030 [4,] 0.0325765 0.230757 0.1530535 0.1366150 0.2497320 [5,] 0.4091660 0.480773 0.3983450 0.2639380 0.3123950 $n [1] 8 8 8 8 8 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.3956847 -0.5386367 -0.08951413 -0.2689711 -0.15478330 -0.5108280 [2,] 0.2174717 0.1011007 0.41502413 0.1486741 -0.02909970 0.1091540 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.06212175 -0.29128931 -0.3715972 -0.1265746 -0.2251961 -0.1408960 [2,] 0.35416975 0.04782931 0.1721232 0.3278966 0.2109481 0.3333020 $out [1] 0.708984 -1.138925 -1.035968 $group [1] 5 5 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.32956, -0.275295, -0.0891065000000001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6q3iz1293620372.ps",horizontal=F,onefile=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,] 2.190349 2.645654 1.627971 2.600402 2.691019 2.864753 1.202484 0.8361743 [2,] 2.926110 3.038242 1.855389 2.850507 2.918907 3.154660 2.152462 1.0013770 [3,] 3.068227 3.211466 2.277914 3.006355 3.019731 3.351015 2.551394 1.1502870 [4,] 3.487210 3.375649 2.574566 3.200711 3.072388 3.500147 2.849831 1.2980830 [5,] 4.031636 3.806511 3.182049 3.423985 3.198120 3.602988 3.498551 1.4002010 $n [1] 12 12 12 12 12 12 12 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2.812306 3.057572 1.949892 2.846624 2.949727 3.193436 2.233320 0.958902 [2,] 3.324149 3.365359 2.605935 3.166085 3.089735 3.508594 2.869468 1.341672 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(2.190349, 2.92611, 3.0682275, 3.4872095, 4.031636, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7q3iz1293620372.ps",horizontal=F,onefile=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,] 2.622555 2.848618 2.735161 [2,] 2.745950 2.899754 2.829893 [3,] 2.818318 3.014734 2.884310 [4,] 2.877164 3.076205 2.918818 [5,] 2.940773 3.198120 3.035714 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 2.758471 2.934253 2.843751 [2,] 2.878166 3.095214 2.924870 $out [1] 3.108093 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(2.6225546625, 2.7459501875, 2.81831828571429, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > try(system("convert tmp/1m3jb1293620372.ps tmp/1m3jb1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/2m3jb1293620372.ps tmp/2m3jb1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/3fu1v1293620372.ps tmp/3fu1v1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/4fu1v1293620372.ps tmp/4fu1v1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/5fu1v1293620372.ps tmp/5fu1v1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/6q3iz1293620372.ps tmp/6q3iz1293620372.png",intern=TRUE)) character(0) > try(system("convert tmp/7q3iz1293620372.ps tmp/7q3iz1293620372.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.233 0.968 3.782