R/gee_func.r

Defines functions weightf

weightf <- function(ehats, ahats, med) {
    w = ahats/(ehats - med)
    w[is.infinite(w)] <- 0
    w[is.nan(w)] <- 0
    w[w == 0] <- max(w)
    list(w = w)
}
herbps10/rlme documentation built on Jan. 9, 2018, 10:45 p.m.