Computation of standard errors for Bayes factors
Description
Compute the standard errors for the Bayes factors estimates.
Usage
1 2 
Arguments
pargrid 
A data frame with components "linkp", "phi", "omg", "kappa". Each row gives a combination of the parameters to compute the new standard errors. 
runs 
A list with outputs from the function

nbatches 
An integer scalar or vector of the same length as runs indicating the number of batches to create for computing the variance using the samples of the first stage. 
bfsize1 
A scalar or vector of the same length as

method 
Which method to use to calculate the Bayes factors: Reverse logistic or MengWong. 
reference 
Which model goes in the denominator. 
transf 
Whether to use the transformed sample mu for the computations. Otherwise it uses z. 
bmmcse.size 
Size for computing the batch means MonteCarlo variance using the samples of the second stage. 
Details
Uses the batch means method to compute the standard errors for Bayes factors.
Value
A list with components

pargrid
The inputted pargrid augmented with the computed standard errors. 
bfEstimate
The estimates of the Bayes factors 
bfSigma
The covariance matrix for the Bayes factors estimates.
References
Roy, V., Tan, A. and Flegal, J. (2015). Estimating standard errors for importance sampling estimators with multiple Markov chains. Technical report, Iowa State University. http://lib.dr.iastate.edu/stat_las_preprints/34
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker. Vote for new features on Trello.