R version 3.6.3 (2020-02-29) -- "Holding the Windsock" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > source('/home/pw/wessanet/cretab') > > > > myrfcuid = '' > > x <- c(5.44,18.78,35.97,53.71,65.86,78.08,89.77,94.08,96.07,97.42) > par2 = '100' > par1 = '1' > par2 <- '100' > par1 <- '1' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Tue, 16 Aug 2016 15:25:14 +0200) > #Author: root > #To cite this work: Wessa P., (2013), Maximum-likelihood Weibull Distribution Fitting (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_fitdistrweibull.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > library(MASS) > PPCCWeibull <- function(shape, scale, x) + { + x <- sort(x) + pp <- ppoints(x) + cor(qweibull(pp, shape=shape, scale=scale), 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(qweibull(ppoints(x), shape=i/10,scale=2),sortx) + } > postscript(file="/home/pw/wessanet/rcomp/tmp/1tscc1587151873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot((par1h:par2h)/10,c[par1h:par2h],xlab='shape',ylab='correlation',main='PPCC Plot - Weibull') > dev.off() null device 1 > f<-fitdistr(x, 'weibull') > f$estimate shape scale 1.903455 70.660091 > f$sd shape scale 0.5363491 12.1765345 > xlab <- paste('Weibull(shape=',round(f$estimate[[1]],2)) > xlab <- paste(xlab,', scale=') > xlab <- paste(xlab,round(f$estimate[[2]],2)) > xlab <- paste(xlab,')') > postscript(file="/home/pw/wessanet/rcomp/tmp/2tb541587151873.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqplot(qweibull(ppoints(x), shape=f$estimate[[1]], scale=f$estimate[[2]]), x, main='QQ plot (Weibull)', xlab=xlab ) > grid() > dev.off() null device 1 > > 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,'scale',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="/home/pw/wessanet/rcomp/tmp/3yt0a1587151873.tab") > > try(system("convert /home/pw/wessanet/rcomp/tmp/1tscc1587151873.ps /home/pw/wessanet/rcomp/tmp/1tscc1587151873.png",intern=TRUE)) convert-im6.q16: profile 'icc': 'RGB ': RGB color space not permitted on grayscale PNG `/home/pw/wessanet/rcomp/tmp/1tscc1587151873.png' @ warning/png.c/MagickPNGWarningHandler/1654. character(0) > try(system("convert /home/pw/wessanet/rcomp/tmp/2tb541587151873.ps /home/pw/wessanet/rcomp/tmp/2tb541587151873.png",intern=TRUE)) convert-im6.q16: profile 'icc': 'RGB ': RGB color space not permitted on grayscale PNG `/home/pw/wessanet/rcomp/tmp/2tb541587151873.png' @ warning/png.c/MagickPNGWarningHandler/1654. character(0) > > proc.time() user system elapsed 0.480 0.088 0.564