Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
3744.5 4098.5 4094 4083 4081.5 4227.5 4246.5 4468 4418.5 4508.5 4836.5 4797 4798.5 4933 4793.5 5020 4987 4738 4784.5 4721.5 4756.5 4666 4535.5 4347 4272 4152 4096.5 4193.5 4244.5 4415.5 4365 4280 4288.5 4289 4166 4179.5 4389.5 4577.5 4322.5 4075.5 3785 3740.5 3755 3795.5 3698 3678.5 3760.5 3892 3844 3831 3957.5 3995 4150 4388 4416 4533.5 4789.5 5160.5 5041 5151.5 5334 5631.5 5620.5 5690 5713 5715.5 5881 6028.5 6018.5 6343 6347 6450.5 6337 6354 6246.5 6520.5 6618.5 6593.5 6481.5 6357 6274 5939 5825.5 5755.5 5764 5668 5760 5801 5814 5726 5668 5489.5 5498 5619.5 5630 5496.5 5318.5 5142.5 5097.5 5021 4989 5008.5 5011.5 5026.5 5080 4979 5133 5143.5 5062.5 5190.5 5116.5 5129 5057.5 5085 5140 5251.5
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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation