Distribution de Burr

Bonsoir,
On sait que la fonction de distribution de Burr est:

F(x)=1 − (1 + x^−ρ/γ)^1/ρ, x ≥ 0, γ > 0, ρ < 0

Je voudrais générer un échantillon qui suit la loi de Burr(γ, ρ), j'ai utilisé le package actuar. Ici un exemple

Code:

---

1 2 3 4 5 6 7

library(actuar) n=10 gamma=0.6 rho=1/3 rate=1 x=rburr(n, shape1=1/gamma, shape2=rho ,scale=1/rate)

---

1. Est-ce que mon travail est correct ?
2. C'est quoi "rate" ?
3. Est-ce qu'on prend toujours rate=1.

Merci d'avance :(