Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
209.704 208.923 208.131 206.492 222.706 221.848 209.704 201.630 202.411 202.411 203.280 204.842 207.273 207.273 205.711 201.630 222.706 225.918 221.067 209.704 214.566 207.273 210.562 212.135 213.774 209.704 210.562 204.842 222.706 228.349 223.498 214.566 224.279 213.774 223.498 222.706 225.137 216.205 225.918 225.137 239.712 236.423 223.498 216.986 225.918 213.774 222.706 224.279 227.568 220.286 224.279 226.710 235.642 228.349 218.636 208.131 217.855 191.125 204.061 211.343 218.636 208.131 208.131 208.131 213.774 205.711 195.129 186.274 192.698 167.618 182.985 191.917 193.556 184.624 185.405 182.985 191.125 185.405 174.130 165.979 179.762 149.831 169.268 178.123 178.123 167.618 157.905 157.124 165.979 157.905 142.549 131.967 143.330 116.611 140.899 153.824 157.905 148.973 137.687 145.761 148.973 146.542 122.243 110.968 119.031 94.743 119.823 128.755 136.037 123.893 112.530 119.031 122.243 115.819 91.531 80.949 90.662 63.943 93.093 110.968
Sample Range:
(leave blank to include all observations)
From:
To:
# blockwidth
Chart options
R Code
par1 <- as.numeric(par1) (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