Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
5813 7359 8242 8321 8463 8613 10088 10088 10230 10529 10601 10821 10821 10893 11113 11113 11263 11705 11997 12367 12517 12809 12959 13030 13251 13322 13543 13543 13621 13984 14284 14355 14355 14355 15089 15089 15238 15238 15388 15830 16193 16193 16342 16635 16635 16784 16855 16855 16934 17375 17375 17447 17518 17596 17739 18330 18330 18401 18401 18480 18551 18622 18622 18701 18701 18772 18843 18843 18843 18921 18921 18921 18921 18921 18993 19064 19064 19064 19064 19142 19142 19213 19285 19434 19434 19434 19434 19506 19506 19655 19726 19805 19876 19876 20026 20097 20168 20246 20246 20246 20246 20246 20318 20318 20318 20389 20389 20389 20467 20467 20538 20538 20538 20610 20688 20759 20759 21422 21493 21792
Sample Range:
(leave blank to include all observations)
From:
To:
# periods
Chart options
R Code
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()
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation