Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
48.948 48.338 48.885 49.217 48.056 49.018 49.490 48.113 49.365 48.093 47.829 47.732 50.759 48.358 49.717 50.546 48.841 49.932 48.932 47.950 47.409 48.803 48.695 47.111 48.749 49.028 47.833 51.445 50.113 49.565 50.343 50.491 50.147 51.764 51.276 51.436 50.403 49.517 49.077 47.781 48.446 47.979 50.436 49.250 49.498 49.500 49.500 51.100 48.900 49.600 50.500 51.000 49.800 50.500 50.704 49.377 50.999 50.166 49.777 50.437 50.673 50.524 50.715 49.319 49.440 49.123 50.036 49.688 48.633 50.103 49.662 49.916
minimum value of shape parameter
maximum value of shape parameter
Chart options
R Code
library(MASS) PPCCGamma <- function(shape, rate, x) { x <- sort(x) pp <- ppoints(x) cor(qgamma(pp, shape=shape, rate=rate), x) } par1 <- as.numeric(par1) par2 <- as.numeric(par2) if (par1 < 0.1) par1 <- 0.1 if (par1 > 50) par1 <- 50 if (par2 < 0.1) par2 <- 0.1 if (par2 > 50) par2 <- 50 par1h <- par1*10 par2h <- par2*10 sortx <- sort(x) c <- array(NA,dim=c(par2h)) for (i in par1h:par2h) { c[i] <- cor(qgamma(ppoints(x), shape=i/10,rate=2),sortx) } bitmap(file='test1.png') plot((par1h:par2h)/10,c[par1h:par2h],xlab='shape',ylab='correlation',main='PPCC Plot - Gamma') dev.off() f<-fitdistr(x, 'gamma') f$estimate f$sd xlab <- paste('Gamma(shape=',round(f$estimate[[1]],2)) xlab <- paste(xlab,', rate=') xlab <- paste(xlab,round(f$estimate[[2]],2)) xlab <- paste(xlab,')') bitmap(file='test2.png') qqplot(qgamma(ppoints(x), shape=f$estimate[[1]], rate=f$estimate[[2]]), x, main='QQ plot (Gamma)', xlab=xlab ) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Parameter',1,TRUE) a<-table.element(a,'Estimated Value',1,TRUE) a<-table.element(a,'Standard Deviation',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'shape',header=TRUE) a<-table.element(a,f$estimate[1]) a<-table.element(a,f$sd[1]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'rate',header=TRUE) a<-table.element(a,f$estimate[2]) a<-table.element(a,f$sd[2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.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