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.0349102355942728,0.00807721497901564,-0.121638609033658,0.271140172031791,0.115362521983019,0.208637719955046,0.00736333218421053,0.364107986611877,-0.401950824320542,0.208115628227064,-0.231184614403349,0.0886065933683892,-0.115016757667903,0.0105278025279736,0.238703049625312,0.0833760938194638,-0.244263809620781,0.0300469072378363,0.160888921911825,0.136128315686798,-0.178471109158176,-0.156191416101216,-0.123271475512903,0.287855689322172,0.114229814421261,-0.373596176543768,0.0534370488169519,-0.0094243656563928,-0.196264153183377,0.0544183083541583,0.458801613285757,-0.481978233278928,0.231921304283570,-0.120496663818209,-0.00860173725973467,-0.101427264965322,-0.096296702560096,0.0893328589991427,-0.338915896651276,-0.579111335739363,-0.240498095063384,-0.379196305811148,-0.198420345132734,0.206304690993777,-0.00161455539305067,-0.152996595141273,0.0272720847132301,0.322114095042851) > 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] 48 > (np <- floor(n / par1)) [1] 4 > 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] 4 4 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -0.034910236 -0.11501676 0.114229814 -0.096296703 NA [2,] 0.008077215 0.01052780 -0.373596177 0.089332859 NA [3,] -0.121638609 0.23870305 0.053437049 -0.338915897 NA [4,] 0.271140172 0.08337609 -0.009424366 -0.579111336 NA [5,] 0.115362522 -0.24426381 -0.196264153 -0.240498095 NA [6,] 0.208637720 0.03004691 0.054418308 -0.379196306 NA [7,] 0.007363332 0.16088892 0.458801613 -0.198420345 NA [8,] 0.364107987 0.13612832 -0.481978233 0.206304691 NA [9,] -0.401950824 -0.17847111 0.231921304 -0.001614555 NA [10,] 0.208115628 -0.15619142 -0.120496664 -0.152996595 NA [11,] -0.231184614 -0.12327148 -0.008601737 0.027272085 NA [12,] 0.088606593 0.28785569 -0.101427265 0.322114095 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 0.04298745 0.12554456 -0.487825991 0.1856296 NA [2,] -0.12971582 0.22817525 0.427033225 -0.4282488 NA [3,] 0.39277878 -0.15532696 -0.062861414 -0.2401954 NA [4,] -0.15577765 -0.32763990 -0.186839788 0.3386132 NA [5,] 0.09327520 0.27431072 0.250682462 -0.1386982 NA [6,] -0.20127439 0.13084201 0.404383305 0.1807760 NA [7,] 0.35674465 -0.02476061 -0.940779847 0.4047250 NA [8,] -0.76605881 -0.31459942 0.713899538 -0.2079192 NA [9,] 0.61006645 0.02227969 -0.352417968 -0.1513820 NA [10,] -0.43930024 0.03291994 0.111894927 0.1802687 NA [11,] 0.31979121 0.41112716 -0.092825528 0.2948420 NA [12,] -0.20362335 -0.17362587 0.005130562 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/1o01g1259700193.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/2zibu1259700193.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/30idu1259700193.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/4ivtq1259700193.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.11501676 -0.373596177 -0.33891590 -0.57911134 -0.24426381 -0.37919631 [2,] -0.10565673 -0.182759481 -0.23027725 -0.29426785 -0.24238095 -0.17457470 [3,] -0.06560347 0.009302509 -0.03410078 0.03697586 -0.21838112 0.04223261 [4,] 0.03965979 0.049930331 0.14607005 0.17725813 -0.04045082 0.13152801 [5,] 0.11422981 0.089332859 0.23870305 0.27114017 0.11536252 0.20863772 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.19842035 -0.4819782 -0.40195082 -0.15619142 -0.231184614 -0.101427265 [2,] -0.09552851 -0.1729250 -0.29021097 -0.15459401 -0.177228045 -0.006410336 [3,] 0.08412613 0.1712165 -0.09004283 -0.13674663 -0.065936606 0.188231141 [4,] 0.30984527 0.2852063 0.11515337 0.04380948 0.009335174 0.304984892 [5,] 0.45880161 0.3641080 0.23192130 0.20811563 0.027272085 0.322114095 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.18040352 -0.1745224 -0.3314151 -0.3355297 -0.37790593 -0.1995885 [2,] 0.04919658 0.1931275 0.2632136 0.4094814 -0.05885632 0.2840538 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.2361192 -0.1907072 -0.4102807 -0.29348538 -0.21332155 -0.05777109 [2,] 0.4043714 0.5331402 0.2301950 0.01999213 0.08144834 0.43423337 $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(-0.115016757667903, -0.105656730113999, -0.0656034690771844, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5w67z1259700193.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.4878260 -0.42824876 -0.2401954 -0.3276399 -0.13869821 -0.20127439 [2,] -0.2224193 -0.27898229 -0.1977612 -0.2572398 -0.02271151 -0.03521619 [3,] 0.0842660 0.04922971 -0.1090942 -0.1713087 0.17197883 0.15580899 [4,] 0.1555871 0.32760424 0.1649587 0.0914178 0.26249659 0.29257963 [5,] 0.1856296 0.42703323 0.3927788 0.3386132 0.27431072 0.40438330 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.9407798 -0.7660588 -0.35241797 -0.43930024 -0.09282553 -0.203623351 [2,] -0.4827702 -0.5403291 -0.25190000 -0.20319015 0.10100824 -0.188624613 [3,] 0.1659920 -0.2612593 -0.06455117 0.07240743 0.30731661 -0.173625875 [4,] 0.3807348 0.2529901 0.31617307 0.14608180 0.36545919 -0.084247656 [5,] 0.4047250 0.7138995 0.61006645 0.18026868 0.41112716 0.005130562 $n [1] 4 4 4 4 4 4 4 4 4 4 4 3 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.214359 -0.4299736 -0.3956429 -0.4467483 -0.05333557 -0.1031497 [2,] 0.382891 0.5284331 0.1774545 0.1041308 0.39729323 0.4147677 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.516177 -0.8879816 -0.5133289 -0.2035174 0.09840036 -0.26883994 [2,] 0.848161 0.3654629 0.3842266 0.3483323 0.51623286 -0.07841181 $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(-0.487825990965029, -0.222419270195870, 0.0842660053845825, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/69vpj1259700193.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] [1,] -0.40195082 -0.24426381 -0.481978233 -0.57911134 NA [2,] -0.07827442 -0.13973145 -0.158380409 -0.28970700 NA [3,] 0.04834190 0.02028735 -0.009013051 -0.12464665 NA [4,] 0.20837667 0.14850862 0.084324061 0.05830247 NA [5,] 0.36410799 0.28785569 0.231921304 0.32211410 NA $n [1] 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [1,] -0.08240159 -0.1111809 -0.1197122 -0.28337610 NA [2,] 0.17908539 0.1517556 0.1016861 0.03408281 NA $out [1] 0.4588016 $group [1] 3 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-0.401950824320542, -0.0782744223139654, 0.0483419041737024, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7rxgj1259700193.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.14141588 -0.21838112 -0.17989850 [2,] -0.07518051 -0.07798972 -0.08186367 [3,] -0.04874793 -0.01239914 -0.03332911 [4,] 0.01730867 0.06317937 0.05299844 [5,] 0.14928728 0.18823114 0.16875921 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.090932867 -0.07678730 -0.09484063 [2,] -0.006562996 0.05198903 0.02818240 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.141415883971131, -0.0751805053125491, -0.048747931759538, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1o01g1259700193.ps tmp/1o01g1259700193.png") > system("convert tmp/2zibu1259700193.ps tmp/2zibu1259700193.png") > system("convert tmp/30idu1259700193.ps tmp/30idu1259700193.png") > system("convert tmp/4ivtq1259700193.ps tmp/4ivtq1259700193.png") > system("convert tmp/5w67z1259700193.ps tmp/5w67z1259700193.png") > system("convert tmp/69vpj1259700193.ps tmp/69vpj1259700193.png") > system("convert tmp/7rxgj1259700193.ps tmp/7rxgj1259700193.png") > > > proc.time() user system elapsed 1.162 0.920 1.466