ss.phiC: Murray...

In mpoll/scale: The Scaleable Langevin Exact Algorithm -- Reference Implementation

Usage

ss.phiC(p.mat.curr, l.bound = NULL, u.bound = NULL)

Arguments

p.mat.curr
l.bound
u.bound

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 (p.mat.curr, l.bound = NULL, u.bound = NULL) 
{
    l.bds <- pmin(0, l.bound)
    u.bds <- pmax(0, u.bound)
    distance <- pmax(abs(l.bds), abs(u.bds))
    alpha.bds <- alpha.prime.bds <- numeric(dimen)
    for (j in 1:dimen) {
        alpha.bds[j] <- dsz * sum(distance * dim.grad.max[j, 
            , 1])
        alpha.prime.bds[j] <- dsz * sum(distance * dim.grad.max[j, 
            , 2])
    }
    A.bds <- alpha.bds^2
    phi.bds <- (sum(2 * abs(alpha.cent) * alpha.bds) + sum(A.bds) + 
        sum(alpha.prime.bds))/2
    phiL <- phi.cent - phi.bds
    phiU <- phi.cent + phi.bds
    intensity <- phiU - phiL
    list(distance = distance, alpha.bds = alpha.bds, alpha.prime.bds = alpha.prime.bds, 
        A.bds = A.bds, phiL = phiL, phiU = phiU, intensity = intensity)
  }

mpoll/scale documentation built on Dec. 9, 2019, 7:15 a.m.