Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
68.733 41.381 51.310 39.036 42.157 52.574 39.168 41.325 40.645 42.303 38.167 8.980 66.356 50.335 47.212 42.247 45.793 48.248 40.074 39.602 41.283 41.900 29.812 7.236 56.184 33.958 34.532 30.277 30.850 35.492 33.544 27.829 33.663 35.690 27.356 8.033 62.798 37.581 44.753 37.546 36.830 50.683 38.487 36.522 45.544 43.575 36.921 11.393 74.787 49.019 56.601 47.637 49.806 50.499 42.092 39.062 44.382 43.635 41.082 17.244 70.170 43.949 52.333 41.032 47.758 76.116 30.917 32.996 31.951 26.775 30.268 18.214 47.957 31.901 35.559 30.408 30.083 35.043 30.475 28.309 31.394 36.313 40.357 38.918 44.368 33.298 29.366 28.282 30.943 32.699 29.764 25.524 29.807 35.112 32.192 36.214
Sample Range:
(leave blank to include all observations)
From:
To:
# blockwidth
Chart options
R Code
par1 <- as.numeric(par1) x <- na.omit(x) (n <- length(x)) (np <- floor(n / par1)) 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 arr darr 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 bitmap(file='plot1.png') plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.mean,0) dev.off() bitmap(file='plot2.png') plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.median,0) dev.off() bitmap(file='plot3.png') plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') mtext(paste('#blocks = ',np)) abline(overall.midrange,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='plot4b.png') 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')) 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.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')) 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