Last data update: 2014.03.03
|
R: Create a SISe3_sp model
Create a SISe3_sp model
Description
Create a SISe3_sp model to be used by the simulation framework.
Usage
SISe3_sp(u0, tspan, events = NULL, phi = NULL, upsilon_1 = NULL,
upsilon_2 = NULL, upsilon_3 = NULL, gamma_1 = NULL, gamma_2 = NULL,
gamma_3 = NULL, alpha = NULL, beta_t1 = NULL, beta_t2 = NULL,
beta_t3 = NULL, beta_t4 = NULL, end_t1 = NULL, end_t2 = NULL,
end_t3 = NULL, end_t4 = NULL, epsilon = NULL, distance = NULL,
coupling = NULL)
Arguments
u0 |
A data.frame with the initial state in each
node, see details.
|
tspan |
An increasing sequence of points in time where the
state of the system is to be returned.
|
events |
a data.frame with the scheduled events, see
siminf_model .
|
phi |
A numeric vector with the initial environmental
infectious pressure in each node. Default NULL which gives 0 in
each node.
|
upsilon_1 |
Indirect transmission rate of the environmental
infectious pressure in age category 1
|
upsilon_2 |
Indirect transmission rate of the environmental
infectious pressure in age category 2
|
upsilon_3 |
Indirect transmission rate of the environmental
infectious pressure in age category 3
|
gamma_1 |
The recovery rate from infected to susceptible for
age category 1
|
gamma_2 |
The recovery rate from infected to susceptible for
age category 2
|
gamma_3 |
The recovery rate from infected to susceptible for
age category 3
|
alpha |
Shed rate from infected individuals
|
beta_t1 |
The decay of the environmental infectious pressure
in the first interval of the year.
|
beta_t2 |
The decay of the environmental infectious pressure
in the second interval of the year.
|
beta_t3 |
The decay of the environmental infectious pressure
in the third interval of the year.
|
beta_t4 |
The decay of the environmental infectious pressure
in the fourth interval of the year.
|
end_t1 |
The non-inclusive day that ends interval 1.
|
end_t2 |
The non-inclusive day that ends interval 2.
|
end_t3 |
The non-inclusive day that ends interval 3.
|
end_t4 |
The non-inclusive day that ends interval 4.
|
epsilon |
The background infectious pressure
|
distance |
The distance matrix between neighboring nodes
|
coupling |
The coupling between neighboring nodes
|
Details
The argument u0 must be a data.frame with one row for
each node with the following columns:
- S_1
The number of sucsceptible in age category 1
- I_1
The number of infected in age category 1
- S_2
The number of sucsceptible in age category 2
- I_2
The number of infected in age category 2
- S_3
The number of sucsceptible in age category 3
- I_3
The number of infected in age category 3
Value
SISe3_sp
Beta
The time dependent beta is divided into four intervals of the year
where 0 <= day < 365 and
0 <= end_t1 < end_t2 < end_t3 < end_t4 <= 365
or
0 <= end_t4 < end_t1 < end_t2 < end_t3 < 364
Case 1: end_t1 < end_t4
INTERVAL_1 INTERVAL_2 INTERVAL_3 INTERVAL_4 INTERVAL_1
[0, end_t1) [end_t1, end_t2) [end_t2, end_t3) [end_t3, end_t4) [end_t4, 365]
Case 2: end_t4 < end_t1
INTERVAL_4 INTERVAL_1 INTERVAL_2 INTERVAL_3 INTERVAL_4
[0, end_t4) [end_t4, end_t1) [end_t1, end_t2) [end_t2, end_t3) [end_t3, 365)
Results
|