Dear R forum

>From the given data, I have estimated the parameters of Wakeby distribution using lmomco package as

library(lmomco)

(amounts <- read.csv("input_S.csv")$amount)

# ___________________________________________________________

# Wakeby distribution - Parameter estimation

N                      =
 length(amounts)
lmr                    = lmom.ub(amounts)
parameters_of_Wakeby   = parwak(lmr)


> parameters_of_Wakeby

$type
[1]
 "wak"

$para
                  xi                alpha 
1.18813927666405e+04 0.00000000000000e+00 
                beta                gamma 
0.00000000000000e+00 8.11391042554567e+04 
               delta 
9.57554297149062e-01 

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.

Thanking in advance

Katherine

 
	[[alternative HTML version deleted]]

