I am doing a simulation study based on a paper. In the simulation, x(data) values are generated from the quantile function of record-based transmuted Weibull distribution which is defined as

where p is a constant (p = 0.5), and q is derived from runif(2000), gamma=0.5 and theta = 0.05.

Does anyone know how to write this Lambert function in R? Is there any library that provides this function? If so, how to use it?

CodePudding user response：

library(sos); findFn("{lambert W}") finds R implementations of the Lambert W function in packages: emdbook, lamW, LambertW, condmixt, VGAM, pracma, spatstat.core ...