Displaying 1 result from an estimated 1 matches for "gradient2".
Did you mean:
gradient
2006 Apr 20
2
nlminb( ) : one compartment open PK model
...)*S)*exp((-1/(2*(S^2)))*
(log.conc-(log(dose*Ke*Ka/(Cl*(Ka-Ke)))
+log(exp(-Ke*time)-exp(-Ka*time))))^2)),
namevec = c("Ke","Ka","Cl","S"),
function.arg = function(Ke, Ka, Cl, S, dose, time, log.conc) NULL )
gradient2.1compart <- function(parms, dose, time, log.conc)
{
Ke <- parms[1]; Ka <- parms[2]; Cl <- parms[3]; S <- parms[4]
colSums(attr(deriv2.1compart(Ke, Ka, Cl, S, dose, time, log.conc), "gradient"))
}
attach(foo.frame)
inits <- c(Ke = .5,
Ka = .5,
Cl = 1,
S = 1)...