Discussion :

[the_reward]

Bonsoir tout le monde,

Je dois simuler un tirage suivant une loi de Laplace, puis faire une fonction qui calcule le log de la densité d'un loi de Laplace. Voilà ce que j'ai fait :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
Nsim=1000
rlaplace = function(mu,b){
  u=runif(Nsim,min=-0.5,max=0.5)
  x = mu-b*sign(u)*log(1-2*abs(u))
  return(x) 
}
res = rlaplace(-5,4)
dloglaplace = function(mu,b){
  g = log(rlaplace(mu,b))
  return(g)
}
res1 = dloglaplace(-5,4)

Mon problème, c'est que quand je passe au log, je trouve des NAN. J'ai pris plusieurs paramètres mais il y a toujours des NAN !
Quelqu'un peut-il me dire où est le problème svp ?

Merci d'avance

[juliatheric]

Bonjour,
Qu'en est-il de la discussion suivante que vous aviez vous-même initiée ?
http://www.developpez.net/forums/d14...r/loi-laplace/

[dev_ggy]

Bonjour,

Ta fonction rlaplace est négative, donc ton passage au logarithme est impossible.

Pourquoi passes-tu au logarithme d'une loi de probabilité qui peut être négative parfois ? Ne fais-tu pas une confusion avec le maximum de vraisemblance ?

Cordialement.

[the_reward]

Bonjour juliatheric,

Je vous remercie pour votre réponse. En fait, j'ai cherché sur Internet et j'ai trouvé que je pouvais simuler une loi de Laplace à partir d'une variable aléatoire qui suit une loi uniforme et en utilisant l'inverse de la fonction de répartition d'une densité de Laplace (c'est ce que j'ai fait dans la fonction rlaplace).

Ensuite, je dois calculer le log d'une densité de Laplace, mais le résultat de rlaplace donne des valeurs négatives, donc le log n'est pas défini !

Je ne vois pas comment calculer le log...

[dev_ggy]

Ta fonction de Laplace définie à partir d'une loi uniforme semble correcte.

Le log de la densité, c'est le logarithme de la probabilité et non de la variable aléatoire. ( Ce n'est pas log(x) mais log(1/2bexp(abs(x-mu)/b))

Elle est proche de -log 2 - log b - abs(x-mu)/b.

Cordialement.

[the_reward]

bonjour dev_ggy,

en fait j'ai pensé à ça, faire le log de la densité -log2-log b-abd(x-mu)/b, mais comment je vais gérer le x?

Cordialement.