// Created by libAntimony v2.13.0
model *restif06()

  // Compartments and Species:
  compartment env;
  substanceOnly species N in env, S in env, I_1 in env, I_2 in env, R_1 in env;
  substanceOnly species R_2 in env, I_1p in env, I_2p in env, R_p in env;

  // Assignment Rules:
  mu := 1/l_e;
  beta_1 := R0_1*gamma_1;
  gamma_1 := 365/tInf_1;
  beta_2 := R0_2*gamma_2;
  gamma_2 := 365/tInf_2;
  sigma := 1/tImm;
  Lambda_1 := beta_1*(I_1 + I_1p)/N;
  Lambda_2 := beta_2*(I_2 + I_2p)/N;
  I_1_frac := (I_1 + I_1p)/N;
  I_2_frac := (I_2 + I_2p)/N;
  S_frac := S/N;
  R1_frac := (R_1 + R_p)/N;
  R2_frac := (R_2 + R_p)/N;
  Rp_frac := R_p/N;

  // Reactions:
  r1:  => S; mu*N;
  r2: S => ; mu*S;
  r3: I_1 => ; mu*I_1;
  r4: I_2 => ; mu*I_2;
  r5: R_1 => ; mu*R_1;
  r6: R_2 => ; mu*R_2;
  r7: I_1p => ; mu*I_1p;
  r8: I_2p => ; mu*I_2p;
  r9: R_p => ; mu*R_p;
  r10: S => I_1; beta_1*((I_1 + I_1p)/N)*S;
  r11: S => I_2; beta_2*((I_2 + I_2p)/N)*S;
  r12: R_2 => I_1p; ((1 - psi)*beta_1*(I_1 + I_1p)/N)*R_2;
  r13: R_1 => I_2p; ((1 - psi)*beta_2*(I_2 + I_2p)/N)*R_1;
  r14: I_1 => R_1; gamma_1*I_1;
  r15: I_2 => R_2; gamma_2*I_2;
  r16: I_1p => R_p; gamma_1*I_1p;
  r17: I_2p => R_p; gamma_2*I_2p;
  r18: R_1 => S; sigma*R_1;
  r19: R_2 => S; sigma*R_2;
  r20: R_p => S; sigma*R_p;

  // Species initializations:
  N = 1;
  S = 0.0588912;
  I_1 = 0.003775;
  I_2 = 1e-06;
  R_1 = 0.93733;
  R_2 = 0;
  I_1p = 0;
  I_2p = 0;
  R_p = 0;

  // Compartment initializations:
  env = 1;
  mu has per_year;

  // Variable initializations:
  l_e = 50;
  l_e has time_unit;
  beta_1 has per_year;
  R0_1 = 17;
  R0_1 has dimensionless;
  gamma_1 has per_year;
  tInf_1 = 21;
  tInf_1 has day;
  beta_2 has per_year;
  R0_2 = 17;
  R0_2 has dimensionless;
  gamma_2 has per_year;
  tInf_2 = 21;
  tInf_2 has day;
  sigma has per_year;
  tImm = 20;
  tImm has time_unit;
  Lambda_1 has per_year;
  Lambda_2 has per_year;
  I_1_frac has dimensionless;
  I_2_frac has dimensionless;
  S_frac has dimensionless;
  R1_frac has dimensionless;
  R2_frac has dimensionless;
  Rp_frac has dimensionless;
  psi = 0.2;
  psi has dimensionless;

  // Other declarations:
  var mu, beta_1, gamma_1, beta_2, gamma_2, sigma, Lambda_1, Lambda_2, I_1_frac;
  var I_2_frac, S_frac, R1_frac, R2_frac, Rp_frac;
  const env, l_e, R0_1, tInf_1, R0_2, tInf_2, tImm, psi;

  // Unit definitions:
  unit substance = item;
  unit time_unit = 3.1536e+07 second;
  unit day = 86400 second;
  unit per_year = 1 / 3.1536e+07 second;

  // Display Names:
  substance is "individuals";
  time_unit is "years";
  day is "days";
  env is "environment";
  l_e is "life expectancy";
  tInf_1 is "infectious period 1";
  tInf_2 is "infectious period 2";
  tImm is "immune period";
  r1 is "Birth";
  r2 is "Death in S";
  r3 is "Death in I_1";
  r4 is "Death in I_2";
  r5 is "Death in R_1";
  r6 is "Death in R_2";
  r7 is "Death in I_1p";
  r8 is "Death in I_2p";
  r9 is "Death in R_p";
  r10 is "Primary Infection with strain 1";
  r11 is "Primary Infection with strain 2";
  r12 is "Secondary Infection with strain 1";
  r13 is "Secondary Infection with strain 2";
  r14 is "Recovery (I_1)";
  r15 is "Recovery (I_2)";
  r16 is "Recovery (I_1p)";
  r17 is "Recovery (I_2p)";
  r18 is "Loss of Immunity (R_1)";
  r19 is "Loss of Immunity (R_2)";
  r20 is "Loss of Immunity (R_p)";

  // SBO terms:
  env.sboTerm = 290
  N.sboTerm = 240
  S.sboTerm = 240
  I_1.sboTerm = 240
  I_2.sboTerm = 240
  R_1.sboTerm = 240
  R_2.sboTerm = 240
  I_1p.sboTerm = 240
  I_2p.sboTerm = 240
  R_p.sboTerm = 240
  r1.sboTerm = 375
  r2.sboTerm = 375
  r3.sboTerm = 375
  r4.sboTerm = 375
  r5.sboTerm = 375
  r6.sboTerm = 375
  r7.sboTerm = 375
  r8.sboTerm = 375
  r9.sboTerm = 375
  r10.sboTerm = 375
  r11.sboTerm = 375
  r12.sboTerm = 375
  r13.sboTerm = 375
  r14.sboTerm = 375
  r15.sboTerm = 375
  r16.sboTerm = 375
  r17.sboTerm = 375
  r18.sboTerm = 375
  r19.sboTerm = 375
  r20.sboTerm = 375

  // CV terms:
  env identity "http://identifiers.org/go/GO:0043245"
  N hypernym "http://identifiers.org/taxonomy/9606"
  S hypernym "http://identifiers.org/taxonomy/9606"
  I_1 version "http://identifiers.org/taxonomy/9606",
              "http://identifiers.org/taxonomy/520"
  I_2 version "http://identifiers.org/taxonomy/9606",
              "http://identifiers.org/taxonomy/520"
  R_1 hypernym "http://identifiers.org/taxonomy/9606"
  R_2 hypernym "http://identifiers.org/taxonomy/9606"
  I_1p version "http://identifiers.org/taxonomy/9606",
               "http://identifiers.org/taxonomy/520"
  I_2p version "http://identifiers.org/taxonomy/9606",
               "http://identifiers.org/taxonomy/520"
  R_p hypernym "http://identifiers.org/taxonomy/9606"
  r1 hypernym "http://identifiers.org/go/GO:0007567"
  r2 hypernym "http://identifiers.org/go/GO:0016265"
  r3 hypernym "http://identifiers.org/go/GO:0016265"
  r4 hypernym "http://identifiers.org/go/GO:0016265"
  r5 hypernym "http://identifiers.org/go/GO:0016265"
  r6 hypernym "http://identifiers.org/go/GO:0016265"
  r7 hypernym "http://identifiers.org/go/GO:0016265"
  r8 hypernym "http://identifiers.org/go/GO:0016265"
  r9 hypernym "http://identifiers.org/go/GO:0016265"
end

restif06 is "Restif2006 - Whooping cough"

restif06 model_entity_is "http://identifiers.org/biomodels.db/MODEL1003290000"
restif06 model_entity_is "http://identifiers.org/biomodels.db/BIOMD0000000249"
restif06 description "http://identifiers.org/pubmed/16615206"
restif06 origin "http://identifiers.org/pubmed/460412",
                "http://identifiers.org/pubmed/460424"
restif06 taxon "http://identifiers.org/taxonomy/9606",
               "http://identifiers.org/taxonomy/520"
restif06 property "http://identifiers.org/doid/DOID:1116"
restif06 hypernym "http://identifiers.org/go/GO:0044403",
                  "http://identifiers.org/go/GO:0009814"
restif06 property "http://identifiers.org/mamo/MAMO_0000046"
