Niedawno przystąpiłem do dopasowywania modeli mieszanych z regresją w ramach Bayesa, używając algorytmu MCMC (właściwie funkcja MCMCglmm w R).
Wydaje mi się, że zrozumiałem, jak zdiagnozować zbieżność procesu szacowania (ślad, wykres Geweke'a, autokorelacja, rozkład tylny ...).
Jedną z rzeczy, która uderza mnie w ramy Bayesa, jest to, że wiele wysiłku wydaje się poświęcać na wykonanie tych diagnostyki, podczas gdy wydaje się, że niewiele robi się w zakresie sprawdzania resztek dopasowanego modelu. Na przykład w MCMCglmm funkcja residual.mcmc () istnieje, ale tak naprawdę nie jest jeszcze zaimplementowana (tzn. Zwraca: „reszty jeszcze nie zaimplementowane dla obiektów MCMCglmm”; ta sama historia dla przewidywania.mcmc ()). Wydaje się, że brakuje go również w innych pakietach, a ogólnie rzecz biorąc, mało jest omawiana w literaturze, którą znalazłem (poza DIC, który jest również dość intensywnie omawiany).
Czy ktoś mógłby wskazać mi kilka przydatnych odniesień, a najlepiej R kodu, z którym mógłbym grać lub modyfikować?
Wielkie dzięki.