R/posteriors.R

In bfp: Bayesian Fractional Polynomials

#####################################################################################
## Author: Daniel Sabanes Bove [daniel *.* sabanesbove *a*t* ifspm *.* uzh *.* ch]
## Project: Bayesian FPs
## 
## Time-stamp: <[posteriors.R] by DSB Mit 16/06/2010 11:32 (CEST)>
##
## Description:
## Extract posterior model probability estimates from BayesMfp object.
##
## History:
## 04/07/2008   copy from thesis function collection.
## 02/10/2009   checks, in particular if there are two estimates for the models when ind=2
## 04/12/2009   update to the new data mode that there is always a frequency estimate,
##              but this may be NA when exhaustive search was used.
#####################################################################################

`posteriors` <-
function (BayesMfpObject,
          ind = 1)         # ind = 1 means normalized posteriors, ind = 2 sampling freqs          
{
    stopifnot(ind %in% c(1, 2))
    
    if(ind == 2){
        ## does the first model have two probability estimates?
        stopifnot(! is.na(BayesMfpObject[[1]]$posterior[ind]))
    }

    ## now extract the required one from all models
    sapply (BayesMfpObject, function (one) one[["posterior"]] [ind])
}