[R] lmomco package - Random number generation using Wakeby distribution
dwinsemius at comcast.net
Mon Jan 21 20:46:01 CET 2013
On Jan 21, 2013, at 10:30 AM, Katherine Gobin wrote:
> Dear R forum
>> From the given data, I have estimated the parameters of Wakeby distribution using lmomco package as
> (amounts <- read.csv("input_S.csv")$amount)
> # ___________________________________________________________
> # Wakeby distribution - Parameter estimation
> N =
> lmr = lmom.ub(amounts)
> parameters_of_Wakeby = parwak(lmr)
It appears you have a) not included the code that produced that output and b) failed to read the Index page for that package
?rlmomco # Random Deviates of a Distribution
So on the assumption that you have an object in your workspace named "parameters_of_Wakeby" and it is an lmomco produced object like that returned by lmom2par() I would try:
> xi alpha
> 1.18813927666405e+04 0.00000000000000e+00
> beta gamma
> 0.00000000000000e+00 8.11391042554567e+04
> This means the scale parameters are 0.
> However, assuming, all the five parameters of Wakeby distribution (viz. location parameter m (xi), the scale parameters a, b, and shape parameters g and d are available.
> Then, how do I generate say 100 random no.s using Wakeby distribution w.r.t. these
> 5 available parameters.
> I couldn't find any information about this in lmomco. Kindly guide if random no.s can be generated or not and if yes, how it can be done in r.
You should have been able to find this with:
Alameda, CA, USA
More information about the R-help