Supposons que vous ayez une variable explicative où représente une coordonnée donnée. Vous avez également une variable de réponse . Maintenant, nous pouvons combiner les deux variables comme:
Dans ce cas, nous choisissons simplement et est une matrice de covariance qui décrit la relation entre et . Cela ne décrit que la valeur de et à . Puisque nous avons plus de points provenant d'autres emplacements pour et , nous pouvons décrire plus de valeurs de de la manière suivante:
Vous remarquerez que nous avons réorganisé les composants de et pour obtenir tous les dans une colonne et après cela, concaténer tous les ensemble. Chaque composant est une fonction de corrélation et est comme ci-dessus. La raison pour laquelle nous avons la covariance est parce que l' on suppose qu'il est possible de séparer la matrice de covariance comme .
Question 1: Lorsque je calcule le , ce que je fais en fait, c'est de générer un ensemble de valeurs de basé sur , correct? J'ai déjà donc je serais plus intéressé à prédire un nouveau point . Dans ce cas, je devrais avoir une matrice définie comme
dans lequel est un vecteur . Par conséquent, nous pouvons construire un vecteur (sans réarrangement):
Et maintenant je viens de réorganiser pour obtenir une distribution conjointe et obtenir le conditionnel .p(y(s0)∣x0,X,Y)
Est-ce correct?
Question 2: Pour prédire, l'article que je lis indique que je dois utiliser cette distribution conditionnelle et obtenir une postérieure distribution , mais je ne sais pas comment obtenir la distribution postérieure des paramètres. Je pourrais peut-être utiliser la distribution que je pense est exactement le même que , puis utilisez simplement le théorème de Bayes pour obtenirp ( μ , T , ϕ ∣ x ( s 0 ) , Y , X ) ( X x ( s 0 ) Y ) p ( X , x ( s 0 ) , Y ∣ μ , T , ϕ ) p
Question 3: À la fin du sous-chapitre, l'auteur dit ceci:
Pour la prédiction, nous n'avons pas . Cela ne crée pas de nouveaux problèmes car il peut être traité comme une variable latente et incorporé dans Cela ne fait que générer un tirage supplémentaire dans chaque itération de Gibbs et est un ajout trivial à la tâche de calcul. x ′
Que signifie ce paragraphe?
Soit dit en passant, cette procédure se trouve dans cet article (page 8), mais comme vous pouvez le voir, j'ai besoin d'un peu plus de détails.
Merci!
la source
Réponses:
Question 1: Étant donné votre modèle de probabilité conjoint la distribution conditionnelle de étant donné est également normal, avec une moyenne et matrice de variance-covariance Y X μ 2 + Σ 21 Σ - 1 11 ( X - μ
Question 2: Le prédictif est défini comme c'est-à-dire en intégrant les paramètres en utilisant la distribution postérieure de ces postérieurs, étant donné les données actuelles . Il y a donc un peu plus dans la réponse complète. Évidemment, si vous avez seulement besoin de simuler à partir du prédictif, votre notion de simulation conjointe à partir de puis à partir de est valide.p(y(s0)∣x(s0),X,Y)
Question 3: Dans le cas où n'est pas observé, la paire peut être prédite à partir d'un autre prédicteurx(s0) (x(s0),y(s0))
Lors de la simulation à partir de ce prédictif, car il n'est pas disponible sous une forme gérable, un échantillonneur Gibbs peut être exécuté qui simule de manière itérative
ou bien fusionner les étapes 4 et 5 en une seule étape
la source