Discussion :

[the_reward]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[chautho]

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.

[juliatheric]

Bonsoir,

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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