R/gee_func.r

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)
}

Try the rlme package in your browser

Any scripts or data that you put into this service are public.

rlme documentation built on May 2, 2019, 3:47 p.m.