Home » date » 2009 » Dec » 23 »

*Unverified author* R Software Module: /rwasp_sdplot.wasp (opens new window with default values) Title produced by software: Standard Deviation Plot Date of computation: Wed, 23 Dec 2009 07:01:40 -0700   Cite this page as follows: Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw.htm/, Retrieved Mon, 20 May 2013 07:30:55 +0000   Original text written by user:   IsPrivate? No (this computation is public)   User-defined keywords: KDGP2W83   System-generated keywords (parent): (pk = 0) Estimated Impact 48   Dataseries X: » Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability « 0.62 0.62 0.62 0.63 0.62 0.63 0.63 0.62 0.63 0.63 0.63 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.63 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.64 0.65 0.65 0.65 0.65 0.65 0.65 0.66 0.66 0.66 0.66 0.68 0.69 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.71 0.7 0.71 0.7 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.7 0.71 0.71 0.7 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.7 0.7 0.68 0.68 0.69 0.69 0.7 0.7 0.7 0.7 0.7 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.71 0.76 0.77 0.78 0.85 0.89 0.9 0.91 0.91 0.91 0.9 0.89 0.88 0.87 0.86 0.87 0.87 0.87 0.85   Output produced by software: Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 3 seconds R Server 'Gwilym Jenkins' @ 72.249.127.135   Charts produced by software: http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/16sge1261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/16sge1261576897.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/20n001261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/20n001261576897.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/3hmj01261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/3hmj01261576897.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/4xnuy1261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/4xnuy1261576897.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/5fr5n1261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/5fr5n1261576897.ps (opens in new window) Click here to open pdf file. http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/6g4tc1261576897.png (opens in new window) http://www.freestatistics.org/blog/date/2009/Dec/23/t1261577067wils9nlb2vy77nw/6g4tc1261576897.ps (opens in new window) Click here to open pdf file.   Parameters (Session): par1 = 12 ;   Parameters (R input): par1 = 12 ;   R code (references can be found in the software module): par1 <- as.numeric(par1) (n <- length(x)) (np <- floor(n / par1)) arr <- array(NA,dim=c(par1,np+1)) ari <- array(0,dim=par1) j <- 0 for (i in 1:n) { j = j + 1 ari[j] = ari[j] + 1 arr[j,ari[j]] <- x[i] if (j == par1) j = 0 } ari arr arr.sd <- array(NA,dim=par1) arr.range <- array(NA,dim=par1) arr.iqr <- array(NA,dim=par1) for (j in 1:par1) { arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) } overall.sd <- sqrt(var(x)) overall.range <- max(x) - min(x) overall.iqr <- quantile(x,0.75) - quantile(x,0.25) bitmap(file='plot1.png') plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.sd,0) dev.off() bitmap(file='plot2.png') plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.range,0) dev.off() bitmap(file='plot3.png') plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') mtext(paste('# blocks = ',np)) abline(overall.iqr,0) dev.off() bitmap(file='plot4.png') 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')) dev.off() bitmap(file='plot5.png') 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')) dev.off() bitmap(file='plot6.png') z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) names(z) <- list('S.D.','Range','IQR') (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) dev.off()