“…Some programs developed in R for the L-L distribution fitting and estimation of its parameters are given as follows: library(MASS) # LHD data x=c (110,13,72,4,45,56,19,27,36,90,19,7,2,118,44,8,277,4,8,10,79,103,6,18,147,96,22,3,24,3,9,99,82,121,54,79,99,18,5,21,1,3,5,1,59,22,17,35,35,29) n=50 logL=function(par){ al=par [1] be=par [2] te=par [3] -sum(log((1+x)*(1+be*x)*teˆ2+al*be*(1+te+te*x))-log(te+1)-(al+1) *log(1+be*x)-te*x)} est=nlm(logL,c(1,.01,.01)) alpha=est $estimate…”