R/gmDA.R

gmDA <-
function (tmrel, te, logmethod = "Eigen") 
{
    gmest = logm(tmrel, method = logmethod)/te
    gmest[which(gmest < 0)] = 0
    diag(gmest) = -rowSums(gmest)
    return(gmest)
}

Try the ctmcd package in your browser

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

ctmcd documentation built on May 29, 2024, 5:23 a.m.