Loi de Laplace

Bonjour tous le monde,
je voudrais faire une fonction qui génère un tirage suivant la loi de Laplace de paramètre (a,b), j'ai fait

Code:

---

1 2 3 4 5

laplace=function(a,b){ x=seq(-30,30,by=1) g=exp(-abs(x-a)/b)/2*b return(g) }

---

ma question est la suivante: est ce que c'est juste de définir x comme une séquence sachant que je dois crée une autre fonction qui génère un vecteur v=(c,d,e,f) tels que les quatre composantes sont tirés indépendamment suivant un loi de Laplace ?

Merci d'avance.

Bonjour,

Tu peux aussi faire en sorte que x soit tiré aléatoirement avec la fonction runif() (uniforme) afin d'avoir que de l'aléatoire.

Bonsoir,

Citation:

---

Envoyé par the_reward

ma question est la suivante: est ce que c'est juste de définir x comme une séquence sachant que je dois crée une autre fonction qui génère un vecteur v=(c,d,e,f) tels que les quatre composantes sont tirés indépendamment suivant un loi de Laplace ?
Merci d'avance.

---

Non, générer une séquence n'est pas la bonne façon de faire. Voir plutôt la bibliothèque (package) Runuran. En utilisant un exemple tiré de sa documentation, ta fonction s'écrirait :

Code:

---

1 2 3 4 5 6

> library("Runuran") > laplace <- function(position, echelle, nbre_composantes) { distr <- udlaplace(position, echelle) ur(pinvd.new(distr), nbre_composantes) } > laplace(0, 1, 4)

---

[1]  1.2193826  0.5357662 -0.4097137  1.1565536