logistic.centering: Compute centering value for logistic regression

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

Usage

logistic.centering(datum.idx.start = 1, length.idxs = dsz)

Arguments

datum.idx.start
length.idxs

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 (datum.idx.start = 1, length.idxs = dsz) 
{
    call.data <- data.subset(datum.idx.start, length.idxs)
    subset.design.mat <- as.matrix(call.data$subset.design[, 
        -1, drop = FALSE])
    subset.data <- call.data$subset.data
    subset.glm <- glm(formula = subset.data ~ subset.design.mat, 
        family = binomial(link = "logit"))
    center <- subset.glm$coef
    precon <- vcov(subset.glm)
    list(center = center, precon = precon, precon.inv = solve(precon), 
        logistic.glm = subset.glm, det.metric = det(precon))
  }

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