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 Nov. 25, 2022, 1:38 p.m.