File: gmDA.R

package info (click to toggle)
r-cran-ctmcd 1.4.1-2
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 1,384 kB
  • sloc: cpp: 183; ansic: 19; makefile: 2
file content (8 lines) | stat: -rw-r--r-- 194 bytes parent folder | download
1
2
3
4
5
6
7
8
gmDA <-
function (tmrel, te, logmethod = "Eigen") 
{
    gmest = logm(tmrel, method = logmethod)/te
    gmest[which(gmest < 0)] = 0
    diag(gmest) = -rowSums(gmest)
    return(gmest)
}