R/marginalize.R

Defines functions marginalize

marginalize <-
function(members, beta, weights) {  
  w <- weights[members]/sum(weights[members])
  pvec <- colSums(beta[members,,drop=FALSE]*w)
  words <- order(pvec, decreasing=TRUE)
  return(list(beta=pvec, indices=words))
}
AntonioCoppola/stmCorrViz documentation built on May 5, 2019, 7:04 a.m.