AW: Conversion of micro constants to macro constants
From: Lutz.Harnisch@aventis.com
Subject: AW: Conversion of micro constants to macro constants
Date: Wed, 15 Aug 2001 10:23:52 +0200
The lambdas (alpha, beta, gamma) can numerically be calculated by evaluating
the Eigenvalues from the Eigenmatrix. The Eigenmatrix is constructed from
the linear DES. See Splus guide on eigen() for references.
The special case for an example of a mammilary model is given below as Splus
code, but any other linear DES could be handled in equivalently.
____________________________________
k10 <- 1.486
k12 <- 0.351
k21 <- 0.080
k13 <- 1.703
k31 <- 0.633
print(c(k10, k12, k21, k13, k31))
dg <- matrix(NA,3,3)
dg[1, ] <- c( - k10 - k12 - k13, k21, k31)
dg[2, ] <- c( k12 , - k21, 0)
dg[3, ] <- c( k13, 0, - k31)
print(dg)
print(-eigen(dg)$values)
_____________________________________
Lutz Harnisch
DMPK PopKin
Aventis Pharma
Frankfurt
Germany