Wessa.net

Send output to:


Data:
1066 1315 4447 2513 1154 2167 2497 4182 3053 1690 1102 3386 1112 2687 773 1701 2592 3255 852 3308 2229 1410 4241 1613 1098 922 4766 1088 1626 854 1327 871 1930 2562 2405 873 3878 2195 1628 1577 3997 893 1597 1049 2499 703 2448 1122 2049 668 1354 665 1430 825 2359 925 1769 505 943 2965 2141 1183 1880 935 1286 2398 1313 2661 1714 722 1116 1868 2193 645 3674 3541 2613 1153 2890 3889 1758 541 888 3223 2890 690 1478 1527 4306 646 1318 1886 1712 988 3587 2692 2464 500 4270 3313 1997 675 2709 719 2137 693 2218 2398 1393 580 1156 976 1268 908 1897 3189 580 1196 763 4117 828 1087 1400 1443

minimum value of sdlog parameter


maximum value of sdlog parameter


Chart options


Title:
Label y-axis:
Label x-axis:

R Code
library(MASS) PPCC <- function(meanlog, sdlog, x) { x <- sort(x) pp <- ppoints(x) cor(qlnorm(pp, meanlog=meanlog, sdlog=sdlog), 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 <- par110 par2h <- par210 sortx <- sort(x) c <- array(NA,dim=c(par2h)) for (i in par1h:par2h) { c[i] <- cor(qlnorm(ppoints(x), meanlog=0,sdlog=i/10),sortx) } bitmap(file='test1.png') plot((par1h:par2h)/10,c[par1h:par2h],xlab='sdlog',ylab='correlation',main='PPCC Plot - Lognormal') dev.off() (f<-fitdistr(x, 'lognormal')) xlab <- paste('Lognormal(meanlog=',round(f$estimate[[1]],2)) xlab <- paste(xlab,', sdlog=') xlab <- paste(xlab,round(f$estimate[[2]],2)) xlab <- paste(xlab,')') bitmap(file='test2.png') qqplot(qlnorm(ppoints(x), meanlog=f$estimate[[1]], sdlog=f$estimate[[2]]), x, main='QQ plot (Lognormal)', 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,'meanlog',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,'sdlog',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')

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