llik_st: log-likelihood of a set of trees

In franciscorichter/dmea: Diversification and Macro-Evolution Analysis

Usage

llik_st(pars, setoftrees, impsam = F)

Arguments

pars	parameters the set of trees
setoftrees	list with a set of trees
impsam	if impsam=TRUE, then the likelihood will be adjusted by importance scalings weights, included on the list of trees

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 (pars, setoftrees, impsam = F)
{
    m = length(setoftrees)
    l = NULL
    for (i in 1:m) {
        s = setoftrees[[i]]
        if (impsam) {
            weight = s$weight
            if (weight == 0) {
                l[i] = 0
            }
            else {
                l[i] = llik(b = pars, n = s$n, E = s$E, t = s$wt) *
                  weight
            }
        }
        else {
            l[i] = llik(b = pars, n = s$n, E = s$E, t = s$wt)
        }
    }
    L = sum(l)
    return(L)
  }

franciscorichter/dmea documentation built on May 16, 2019, 1:54 p.m.