maxbgrRatio: ~~function to do ... ~~

In jservadio/TrialistNof1: Bayesian N-of-1 calculations

Usage

maxbgrRatio(x)

Arguments

x

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x) 
{
    t1 = mcmc(x[, 1, drop = F])
    t2 = mcmc(x[, 2, drop = F])
    t3 = mcmc(x[, 3, drop = F])
    cols.to.drop = seq(dim(t1)[2])[apply(t1, 2, var) == 0]
    if (length(cols.to.drop > 0)) {
        t1 = t1[, -cols.to.drop, drop = F]
        t2 = t2[, -cols.to.drop, drop = F]
        t3 = t3[, -cols.to.drop, drop = F]
    }
    max.bgrRatio = max(gelman.diag(list(t1, t2, t3), autoburnin = F)[[1]][, 
        1])
    return(max.bgrRatio)
  }

jservadio/TrialistNof1 documentation built on May 20, 2019, 2:08 a.m.