Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
9.026 9.787 9.536 9.49 9.736 9.694 9.647 9.753 10.07 10.137 9.984 9.732 9.103 9.155 9.308 9.394 9.948 10.177 10.002 9.728 10.002 10.063 10.018 9.96 10.236 10.893 10.756 10.94 10.997 10.827 10.166 10.186 10.457 10.368 10.244 10.511 10.812 10.738 10.171 9.721 9.897 9.828 9.924 10.371 10.846 10.413 10.709 10.662 10.57 10.297 10.635 10.872 10.296 10.383 10.431 10.574 10.653 10.805 10.872 10.625 10.407 10.463 10.556 10.646 10.702 11.353 11.346 11.451 11.964 12.574 13.031 13.812 14.544 14.931 14.886 16.005 17.064 15.168 16.05 15.839 15.137 14.954 15.648 15.305
Sample Range:
(leave blank to include all observations)
From:
To:
Seasonal Period
12
12
4
6
12
Chart options
R Code
par1 <- as.numeric(par1) (n <- length(x)) (np <- floor(n / par1)) arr <- array(NA,dim=c(par1,np)) j <- 0 k <- 1 for (i in 1:(np*par1)) { j = j + 1 arr[j,k] <- x[i] if (j == par1) { j = 0 k=k+1 } } arr arr.mean <- array(NA,dim=np) arr.sd <- array(NA,dim=np) arr.range <- array(NA,dim=np) for (j in 1:np) { arr.mean[j] <- mean(arr[,j],na.rm=TRUE) arr.sd[j] <- sd(arr[,j],na.rm=TRUE) arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE) } arr.mean arr.sd arr.range (lm1 <- lm(arr.sd~arr.mean)) (lnlm1 <- lm(log(arr.sd)~log(arr.mean))) (lm2 <- lm(arr.range~arr.mean)) bitmap(file='test1.png') plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation') dev.off() bitmap(file='test2.png') plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range') dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Section',header=TRUE) a<-table.element(a,'Mean',header=TRUE) a<-table.element(a,'Standard Deviation',header=TRUE) a<-table.element(a,'Range',header=TRUE) a<-table.row.end(a) for (j in 1:np) { a<-table.row.start(a) a<-table.element(a,j,header=TRUE) a<-table.element(a,arr.mean[j]) a<-table.element(a,arr.sd[j] ) a<-table.element(a,arr.range[j] ) a<-table.row.end(a) } a<-table.end(a) table.save(a,file='mytable.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'alpha',header=TRUE) a<-table.element(a,lm1$coefficients[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'beta',header=TRUE) a<-table.element(a,lm1$coefficients[[2]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'S.D.',header=TRUE) a<-table.element(a,summary(lm1)$coefficients[2,2]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'T-STAT',header=TRUE) a<-table.element(a,summary(lm1)$coefficients[2,3]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,summary(lm1)$coefficients[2,4]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable1.tab') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'alpha',header=TRUE) a<-table.element(a,lnlm1$coefficients[[1]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'beta',header=TRUE) a<-table.element(a,lnlm1$coefficients[[2]]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'S.D.',header=TRUE) a<-table.element(a,summary(lnlm1)$coefficients[2,2]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'T-STAT',header=TRUE) a<-table.element(a,summary(lnlm1)$coefficients[2,3]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'p-value',header=TRUE) a<-table.element(a,summary(lnlm1)$coefficients[2,4]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Lambda',header=TRUE) a<-table.element(a,1-lnlm1$coefficients[[2]]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable2.tab')
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