utilisation de la fonction InvGam
Bonjour,
Je souhaite tracer la courbe d'une fonction suivant une loi InvGam(46,1528).
J'ai du trouver cette loi (ce posterior) à partir des données :
- prior phi suit la loi InvGam(38,1110)
- likelihood y|phi suit la loi Normal(163,sigma) (variance inconnue)
- et 16 valeurs pour les yi
Dans l'aide de matlab, je vois que la fonction est définie par gaminv(P,A,B) où p est une proba et a et b des vecteurs. Je ne comprends pas ce que contiennent ces vecteurs.
Je pense ne pas avoir les données nécessaires pour utiliser cette fonction. Y-en-a-til d'autres?
Je vois 2 façons de tracer ma courbe :
1) en me servant du résultat, ma fonction suit une loi InvGam(46,1528) mais je sais pas sur quel intervalle. Je me demande si je peux faire :
Code:
1
2
3
| x=[161 164 181 174 163 162 164 174 160 168 172 156 172 161 174 164];
y = gaminv(1./x,46,1/1528);
plot(x,y,'-') |
Ceci n'est pas trop un problème de programmation, mais si quelqu'un sait ...
Ainsi, j'obtiens une courbe, mais je suis pas sûre qu'elle soit correcte
2) faire sans le résultat:
Code:
1
2
3
| phi = gaminv( [161 164 181 174 163 162 164 174 160 168 172 156 172 161 174 164] ,38,1110);
y = normpdf(phi,1164,sigma);
plot(x,y) |
Dans ce cas, je ne sais pas quoi faire de sigma.
Si quelqu'un des réponses à mes questions, merci de m'éclairer :D
