R/gmDA.R

In ctmcd: Estimating the Parameters of a Continuous-Time Markov Chain from Discrete-Time Data

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