Displaying 1 result from an estimated 1 matches for "hct0".
Did you mean:
hc0
2007 Nov 06
1
Integrating a function
...he results?
I have tried to solve this problem by adding a “while loop”. It works but
the expression becomes to complex if I want to perform further optimization.
Have any one a solution for this problem?
Thanks a lot in advance.
Example
# General equation
beta1 <- 0.001
beta2 <- 0
Hct0 <- 0.27
tau <- 70
k <- 0.001
alpha <- 0.3
T2 <-40
step <- function(x){
if(x>=0) step <- 1
else step <- 0
}
Nmax <- 250
n <- 1
while(n<Nmax){
surviv <- function(n){
surviv <- exp(-k*n)/(1+exp(alpha*(n...