Rendering of the source text

<?xml version="1.0"?>
<model xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#" cmeta:id="wodarz_2003" name="wodarz_2003">

	<documentation xmlns="http://cellml.org/tmp-documentation">
		<article>
			<articleinfo>
				<title>Modelling the Evolution of Immunological Memory</title>
				<author>
					<firstname>Catherine</firstname>
					<surname>Lloyd</surname>
					<affiliation>
						<shortaffil>Auckland Bioengineering Institute, The University of Auckland</shortaffil>
					</affiliation>
				</author>
			</articleinfo>
			<section id="sec_status">
				<title>Model Status</title>
				<para>
            This model runs in Open Cell but due to the timescale being in days it is not able to be easily graphed in COR. The model can produce output that matches figure 2A in the published paper.  The units have been checked and they balance.
          </para>
			</section>
			<sect1 id="sec_structure">
				<title>Model Structure</title>
				<para>
ABSTRACT: Memory is a central characteristic of immune responses. It is defined as an elevated number of specific immune cells that remain after resolution of infection and can protect the host against reinfection. The evolution of immunological memory is subject to debate. The advantages of memory discussed so far include protection from reinfection, control of chronic infection, and the transfer of immune function to the next generation. Mathematical models are used to identify a new force that can drive the evolution of immunological memory: the duration of memory can regulate the degree of competition between different pathogens. While a long duration of memory provides lasting protection against reinfection, it may also allow an inferior pathogen species to persist. This can be detrimental for the host if the inferior pathogen is more virulent. On the other hand, a shorter duration of memory ensures that an inferior pathogen species is excluded. This can be beneficial for the host if the inferior pathogen is more virulent. Thus, while in the absence of pathogen diversity memory is always expected to evolve to a long duration, under specific circumstances, memory can evolve toward shorter durations in the presence of pathogen diversity.    
</para>
				
				<para>
The original paper reference is cited below:
</para>
				<para>
Evolution of Immunological Memory and the Regulation of Competition between Pathogens, Dominik Wodarz, 2003, <emphasis>Current Biology</emphasis>, 13, 1648-1652.  <ulink url="http://www.ncbi.nlm.nih.gov/pubmed/13678598">PubMed ID: 13678598</ulink>
				</para>
				<informalfigure float="0" id="fig_cell_diagram1">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell diagram1</title>
							</objectinfo>
							<imagedata fileref="wodarz_2003a.png"/>
						</imageobject>
					</mediaobject>
					<caption>Interactions between a single pathogen population and a host population.</caption>
				</informalfigure>
				<informalfigure float="0" id="fig_cell_diagram2">
					<mediaobject>
						<imageobject>
							<objectinfo>
								<title>cell diagram2</title>
							</objectinfo>
							<imagedata fileref="wodarz_2003b.png"/>
						</imageobject>
					</mediaobject>
					<caption>Interactions between two pathogen populations and a host population.</caption>
				</informalfigure>
			</sect1>
		</article>
	</documentation>
	
	
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#S">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">S</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">uninfected and susceptible hosts</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#I1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">I1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">infected with pathogen 1 host cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#I2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">I2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">infected with pathogen 2 host cells</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#I12">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">I12</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">hosts recovered from infection with pathogen 2 then infected with pathogen 1</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#I21">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">I21</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">hosts recovered from infection with pathogen 1 then infected with pathogen 2</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#R1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">recovered and immune hosts after infection with pathogen 1</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#R2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">recovered and immune hosts after infection with pathogen 2</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#R12">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">R12</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">recovered and immune host after infection with pathogen 1 and pathogen 2</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#P1">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P1</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">population of pathogens 1</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
		<rdf:Description rdf:about="#P2">
			<dc:title xmlns:dc="http://purl.org/dc/elements/1.1/">P2</dc:title>
			<dcterms:alternative xmlns:dcterms="http://purl.org/dc/terms/">population of pathogens 2</dcterms:alternative>
		</rdf:Description>
	</rdf:RDF>
	
	
	   <units name="day">
      <unit multiplier="86400" units="second"/>
   </units>
   <units name="first_order_rate_constant">
      <unit exponent="-1" units="day"/>
   </units>
   <component name="environment">
      <variable cmeta:id="environment_time" name="time" public_interface="out" units="day"/>
   </component>
   <component name="S">
      <variable cmeta:id="S_S" initial_value="100" name="S" public_interface="out" units="dimensionless"/>
      <variable initial_value="0.5" name="r" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="epsilon" units="dimensionless"/>
      <variable name="H" units="dimensionless"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="g" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="R_1" public_interface="in" units="dimensionless"/>
      <variable name="R_2" public_interface="in" units="dimensionless"/>
      <variable name="R_12" public_interface="in" units="dimensionless"/>
      <variable name="P_1" public_interface="in" units="dimensionless"/>
      <variable name="P_2" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <variable name="I_1" public_interface="in" units="dimensionless"/>
      <variable name="I_2" public_interface="in" units="dimensionless"/>
      <variable name="I_12" public_interface="in" units="dimensionless"/>
      <variable name="I_21" public_interface="in" units="dimensionless"/>
      <variable name="P" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>P</ci>
            <apply>
               <plus/>
               <ci>P_1</ci>
               <ci>P_2</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>H</ci>
            <apply>
               <plus/>
               <ci>S</ci>
               <ci>I_1</ci>
               <ci>R_1</ci>
               <ci>I_2</ci>
               <ci>R_2</ci>
               <ci>I_12</ci>
               <ci>I_21</ci>
               <ci>R_12</ci>
            </apply>
         </apply>
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>S</ci>
            </apply>
            <apply>
               <plus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <minus/>
                        <apply>
                           <divide/>
                           <apply>
                              <times/>
                              <ci>r</ci>
                              <ci>H</ci>
                           </apply>
                           <apply>
                              <plus/>
                              <apply>
                                 <times/>
                                 <ci>epsilon</ci>
                                 <ci>H</ci>
                              </apply>
                              <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="dimensionless">1</cn>
                           </apply>
                        </apply>
                        <apply>
                           <times/>
                           <ci>d</ci>
                           <ci>S</ci>
                        </apply>
                     </apply>
                     <apply>
                        <times/>
                        <ci>beta_1</ci>
                        <ci>S</ci>
                        <ci>P_1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>beta_2</ci>
                     <ci>S</ci>
                     <ci>P_2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>g</ci>
                  <apply>
                     <plus/>
                     <ci>R_1</ci>
                     <ci>R_2</ci>
                     <ci>R_12</ci>
                  </apply>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I1">
      <variable cmeta:id="I1_I_1" initial_value="0" name="I_1" public_interface="out" units="dimensionless"/>
      <variable name="a_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="alpha_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="S" public_interface="in" units="dimensionless"/>
      <variable name="P_1" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta_1</ci>
                     <ci>S</ci>
                     <ci>P_1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a_1</ci>
                     <ci>I_1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_1</ci>
                  <ci>I_1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I2">
      <variable cmeta:id="I2_I_2" initial_value="0" name="I_2" public_interface="out" units="dimensionless"/>
      <variable name="a_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="alpha_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="S" public_interface="in" units="dimensionless"/>
      <variable name="P_2" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta_2</ci>
                     <ci>S</ci>
                     <ci>P_2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a_2</ci>
                     <ci>I_2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_2</ci>
                  <ci>I_2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I12">
      <variable cmeta:id="I12_I_12" initial_value="0" name="I_12" public_interface="out" units="dimensionless"/>
      <variable name="a_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="alpha_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="R_1" public_interface="in" units="dimensionless"/>
      <variable name="P_2" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_12</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta_2</ci>
                     <ci>R_1</ci>
                     <ci>P_2</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a_2</ci>
                     <ci>I_12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_2</ci>
                  <ci>I_12</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="I21">
      <variable cmeta:id="I21_I_21" initial_value="0" name="I_21" public_interface="out" units="dimensionless"/>
      <variable name="a_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="alpha_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="R_2" public_interface="in" units="dimensionless"/>
      <variable name="P_1" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>I_21</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <times/>
                     <ci>beta_1</ci>
                     <ci>R_2</ci>
                     <ci>P_1</ci>
                  </apply>
                  <apply>
                     <times/>
                     <ci>a_1</ci>
                     <ci>I_21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>alpha_1</ci>
                  <ci>I_21</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="R1">
      <variable cmeta:id="R1_R_1" initial_value="0" name="R_1" public_interface="out" units="dimensionless"/>
      <variable name="alpha_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="g" public_interface="in" units="first_order_rate_constant"/>
      <variable name="I_1" public_interface="in" units="dimensionless"/>
      <variable name="P_2" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>R_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>alpha_1</ci>
                        <ci>I_1</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>d</ci>
                        <ci>R_1</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>g</ci>
                     <ci>R_1</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_2</ci>
                  <ci>R_1</ci>
                  <ci>P_2</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="R2">
      <variable cmeta:id="R2_R_2" initial_value="0" name="R_2" public_interface="out" units="dimensionless"/>
      <variable name="alpha_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="beta_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="g" public_interface="in" units="first_order_rate_constant"/>
      <variable name="I_2" public_interface="in" units="dimensionless"/>
      <variable name="P_1" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>R_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <minus/>
                     <apply>
                        <times/>
                        <ci>alpha_2</ci>
                        <ci>I_2</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>d</ci>
                        <ci>R_2</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>g</ci>
                     <ci>R_2</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>beta_1</ci>
                  <ci>R_2</ci>
                  <ci>P_1</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="R12">
      <variable initial_value="0" name="R_12" public_interface="out" units="dimensionless"/>
      <variable name="alpha_2" public_interface="in" units="first_order_rate_constant"/>
      <variable name="alpha_1" public_interface="in" units="first_order_rate_constant"/>
      <variable name="d" public_interface="in" units="first_order_rate_constant"/>
      <variable name="g" public_interface="in" units="first_order_rate_constant"/>
      <variable name="I_21" public_interface="in" units="dimensionless"/>
      <variable name="I_12" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>R_12</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <minus/>
                  <apply>
                     <plus/>
                     <apply>
                        <times/>
                        <ci>alpha_2</ci>
                        <ci>I_12</ci>
                     </apply>
                     <apply>
                        <times/>
                        <ci>alpha_1</ci>
                        <ci>I_21</ci>
                     </apply>
                  </apply>
                  <apply>
                     <times/>
                     <ci>d</ci>
                     <ci>R_12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>g</ci>
                  <ci>R_12</ci>
               </apply>
            </apply>
         </apply>
      </math>
   </component>
   <component name="P1">
      <variable cmeta:id="P1_P_1" initial_value="1" name="P_1" public_interface="out" units="dimensionless"/>
      <variable initial_value="1" name="k_1" units="first_order_rate_constant"/>
      <variable name="u" public_interface="in" units="first_order_rate_constant"/>
      <variable name="I_1" public_interface="in" units="dimensionless"/>
      <variable name="I_21" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <variable cmeta:id="P1_log_P1" name="log_P1" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_1</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_1</ci>
                  <apply>
                     <plus/>
                     <ci>I_1</ci>
                     <ci>I_21</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>u</ci>
                  <ci>P_1</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>log_P1</ci>
            <apply>
               <log/>
               <ci>P_1</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="P2">
      <variable cmeta:id="P2_P_2" initial_value="1" name="P_2" public_interface="out" units="dimensionless"/>
      <variable initial_value="1" name="k_2" units="first_order_rate_constant"/>
      <variable name="u" public_interface="in" units="first_order_rate_constant"/>
      <variable name="I_2" public_interface="in" units="dimensionless"/>
      <variable name="I_12" public_interface="in" units="dimensionless"/>
      <variable name="time" public_interface="in" units="day"/>
      <variable cmeta:id="P2_log_P2" name="log_P2" units="dimensionless"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <apply>
               <diff/>
               <bvar>
                  <ci>time</ci>
               </bvar>
               <ci>P_2</ci>
            </apply>
            <apply>
               <minus/>
               <apply>
                  <times/>
                  <ci>k_2</ci>
                  <apply>
                     <plus/>
                     <ci>I_2</ci>
                     <ci>I_12</ci>
                  </apply>
               </apply>
               <apply>
                  <times/>
                  <ci>u</ci>
                  <ci>P_2</ci>
               </apply>
            </apply>
         </apply>
         <apply>
            <eq/>
            <ci>log_P2</ci>
            <apply>
               <log/>
               <ci>P_2</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="memory_duration">
      <variable name="G" units="dimensionless"/>
      <variable initial_value="0.01" name="g" public_interface="out" units="first_order_rate_constant"/>
      <math xmlns="http://www.w3.org/1998/Math/MathML">
         <apply>
            <eq/>
            <ci>G</ci>
            <apply>
               <divide/>
               <cn xmlns:cellml="http://www.cellml.org/cellml/1.0#" cellml:units="first_order_rate_constant">1</cn>
               <ci>g</ci>
            </apply>
         </apply>
      </math>
   </component>
   <component name="kinetic_parameters">
      <variable initial_value="1" name="beta_1" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="alpha_1" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="beta_2" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.1" name="alpha_2" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.03" name="a_1" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="1" name="a_2" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.01" name="d" public_interface="out" units="first_order_rate_constant"/>
      <variable initial_value="0.5" name="u" public_interface="out" units="first_order_rate_constant"/>
   </component>
   <connection>
      <map_components component_1="S" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I12" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="I21" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="R12" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="P1" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="P2" component_2="environment"/>
      <map_variables variable_1="time" variable_2="time"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="P1"/>
      <map_variables variable_1="P_1" variable_2="P_1"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="P2"/>
      <map_variables variable_1="P_2" variable_2="P_2"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="R1"/>
      <map_variables variable_1="R_1" variable_2="R_1"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="R2"/>
      <map_variables variable_1="R_2" variable_2="R_2"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="R12"/>
      <map_variables variable_1="R_12" variable_2="R_12"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="S"/>
      <map_variables variable_1="S" variable_2="S"/>
      <map_variables variable_1="I_1" variable_2="I_1"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="S"/>
      <map_variables variable_1="S" variable_2="S"/>
      <map_variables variable_1="I_2" variable_2="I_2"/>
   </connection>
   <connection>
      <map_components component_1="I12" component_2="S"/>
      <map_variables variable_1="I_12" variable_2="I_12"/>
   </connection>
   <connection>
      <map_components component_1="I21" component_2="S"/>
      <map_variables variable_1="I_21" variable_2="I_21"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="P1"/>
      <map_variables variable_1="I_1" variable_2="I_1"/>
      <map_variables variable_1="P_1" variable_2="P_1"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="P2"/>
      <map_variables variable_1="I_2" variable_2="I_2"/>
      <map_variables variable_1="P_2" variable_2="P_2"/>
   </connection>
   <connection>
      <map_components component_1="I21" component_2="P1"/>
      <map_variables variable_1="I_21" variable_2="I_21"/>
      <map_variables variable_1="P_1" variable_2="P_1"/>
   </connection>
   <connection>
      <map_components component_1="I12" component_2="P2"/>
      <map_variables variable_1="I_12" variable_2="I_12"/>
      <map_variables variable_1="P_2" variable_2="P_2"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="I1"/>
      <map_variables variable_1="I_1" variable_2="I_1"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="I12"/>
      <map_variables variable_1="R_1" variable_2="R_1"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="I21"/>
      <map_variables variable_1="R_2" variable_2="R_2"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="I2"/>
      <map_variables variable_1="I_2" variable_2="I_2"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="P2"/>
      <map_variables variable_1="P_2" variable_2="P_2"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="P1"/>
      <map_variables variable_1="P_1" variable_2="P_1"/>
   </connection>
   <connection>
      <map_components component_1="R12" component_2="I21"/>
      <map_variables variable_1="I_21" variable_2="I_21"/>
   </connection>
   <connection>
      <map_components component_1="R12" component_2="I12"/>
      <map_variables variable_1="I_12" variable_2="I_12"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="kinetic_parameters"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
   </connection>
   <connection>
      <map_components component_1="I1" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="a_1" variable_2="a_1"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
   </connection>
   <connection>
      <map_components component_1="I2" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="a_2" variable_2="a_2"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
   </connection>
   <connection>
      <map_components component_1="I21" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="a_1" variable_2="a_1"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
   </connection>
   <connection>
      <map_components component_1="I12" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="a_2" variable_2="a_2"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="d" variable_2="d"/>
      <map_variables variable_1="beta_2" variable_2="beta_2"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="beta_1" variable_2="beta_1"/>
      <map_variables variable_1="d" variable_2="d"/>
   </connection>
   <connection>
      <map_components component_1="R12" component_2="kinetic_parameters"/>
      <map_variables variable_1="alpha_1" variable_2="alpha_1"/>
      <map_variables variable_1="alpha_2" variable_2="alpha_2"/>
      <map_variables variable_1="d" variable_2="d"/>
   </connection>
   <connection>
      <map_components component_1="P1" component_2="kinetic_parameters"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="P2" component_2="kinetic_parameters"/>
      <map_variables variable_1="u" variable_2="u"/>
   </connection>
   <connection>
      <map_components component_1="S" component_2="memory_duration"/>
      <map_variables variable_1="g" variable_2="g"/>
   </connection>
   <connection>
      <map_components component_1="R1" component_2="memory_duration"/>
      <map_variables variable_1="g" variable_2="g"/>
   </connection>
   <connection>
      <map_components component_1="R2" component_2="memory_duration"/>
      <map_variables variable_1="g" variable_2="g"/>
   </connection>
   <connection>
      <map_components component_1="R12" component_2="memory_duration"/>
      <map_variables variable_1="g" variable_2="g"/>
   </connection>
<rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#"><rdf:Description rdf:about="rdf:#c719dd67-0bf3-441e-8dcc-9bd7be04bf67"><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#0cd404ed-f8a6-4537-9ff9-f4d1cf378fe9"/></N><ORG xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#4c5a7ebe-9bfa-4ffd-aa26-3939043c4175"/></ORG><EMAIL xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#7994b0f0-d145-482c-860c-7e20a57ae564"/></EMAIL></rdf:Description><rdf:Description rdf:about="#I12"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">I12</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">infected hosts</alternative></rdf:Description><rdf:Description rdf:about="rdf:#c2ffd07b-2fd1-41b0-bc90-b086171bc24c"><subject xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#8194e217-9709-4ea3-8220-045c3ec7b106"/></subject></rdf:Description><rdf:Description rdf:nodeID="n1"><rdf:_1 xml:lang="en">immunology</rdf:_1><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/></rdf:type><rdf:_2 xml:lang="en">interspecific pathogenic competition</rdf:_2></rdf:Description><rdf:Description rdf:about="rdf:#7994b0f0-d145-482c-860c-7e20a57ae564"><rdf:type><rdf:Description rdf:about="http://imc.org/vCard/3.0#internet"/></rdf:type><rdf:value xml:lang="en">c.lloyd@auckland.ac.nz</rdf:value></rdf:Description><rdf:Description rdf:about="#I21"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">I21</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">infected hosts</alternative></rdf:Description><rdf:Description rdf:nodeID="n2"><boundIntervals xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n3"/></boundIntervals></rdf:Description><rdf:Description rdf:about="rdf:#f39648d7-22c1-44ff-8cd1-615f9026e25c"><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#citationAuthorsSeq"/></creator><issued xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#34be0937-46e5-4c7f-bf21-cf0d501445be"/></issued><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Evolution of Immunological Memory and the Regulation of Competition between Pathogens</title><Journal xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#4d81b96d-ac9a-4c9e-b606-53bf6a3e2bec"/></Journal><last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">1652</last_page><volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">13</volume><first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">1648</first_page></rdf:Description><rdf:Description rdf:about="rdf:#34be0937-46e5-4c7f-bf21-cf0d501445be"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2003-09-16</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#4c5a7ebe-9bfa-4ffd-aa26-3939043c4175"><Orgunit xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The Bioengineering Institute</Orgunit><Orgname xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">The University of Auckland</Orgname></rdf:Description><rdf:Description rdf:about="rdf:#citationAuthorsSeq"><rdf:_1><rdf:Description rdf:about="rdf:#author1Vcard"/></rdf:_1><rdf:type><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/></rdf:type></rdf:Description><rdf:Description rdf:about="rdf:#ac74765f-535e-48f3-9de6-4ae538fbec5c"><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#9025aa44-9e18-4f72-929a-c86b0e27f9fa"/></creator><rdf:value xml:lang="en">This is the CellML description of Wodarz's 2nd 2003 mathematical model of the evolution of immunological memory and the regulation of competition between pathogens.</rdf:value></rdf:Description><rdf:Description rdf:about="rdf:#e5516720-c1e2-4ae2-bbe6-925204067934"><Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">13678598</Pubmed_id><JournalArticle xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#f39648d7-22c1-44ff-8cd1-615f9026e25c"/></JournalArticle></rdf:Description><rdf:Description rdf:about="rdf:#ac4ccc76-4569-48be-8324-7c378243f2cf"><W3CDTF xmlns="http://purl.org/dc/terms/" xml:lang="en">2003-12-17</W3CDTF></rdf:Description><rdf:Description rdf:about="rdf:#0cd404ed-f8a6-4537-9ff9-f4d1cf378fe9"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Lloyd</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine</Given><Other xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">May</Other></rdf:Description><rdf:Description rdf:about="#S"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">S</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">uninfected and susceptible hosts</alternative></rdf:Description><rdf:Description rdf:nodeID="n4"><endingValue xmlns="http://www.cellml.org/metadata/simulation/1.0#" xml:lang="en">100.0</endingValue></rdf:Description><rdf:Description rdf:about="rdf:#9025aa44-9e18-4f72-929a-c86b0e27f9fa"><FN xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Catherine Lloyd</FN></rdf:Description><rdf:Description rdf:about="#I2"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">I2</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">infected hosts</alternative></rdf:Description><rdf:Description rdf:about="rdf:#author1VcardN"><Family xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Wodarz</Family><Given xmlns="http://www.w3.org/2001/vcard-rdf/3.0#" xml:lang="en">Dominik</Given></rdf:Description><rdf:Description rdf:about="rdf:#4d81b96d-ac9a-4c9e-b606-53bf6a3e2bec"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">Current Biology</title></rdf:Description><rdf:Description rdf:about="rdf:#author1Vcard"><rdf:type><rdf:Description rdf:about="http://www.cellml.org/bqs/1.0#Person"/></rdf:type><N xmlns="http://www.w3.org/2001/vcard-rdf/3.0#"><rdf:Description rdf:about="rdf:#author1VcardN"/></N></rdf:Description><rdf:Description rdf:nodeID="n3"><rdf:rest><rdf:Description rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#nil"/></rdf:rest><rdf:first><rdf:Description rdf:nodeID="n4"/></rdf:first></rdf:Description><rdf:Description rdf:about=""><creator xmlns="http://purl.org/dc/elements/1.1/"><rdf:Description rdf:about="rdf:#c719dd67-0bf3-441e-8dcc-9bd7be04bf67"/></creator><created xmlns="http://purl.org/dc/terms/"><rdf:Description rdf:about="rdf:#ac4ccc76-4569-48be-8324-7c378243f2cf"/></created><publisher xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
        The University of Auckland, Bioengineering Institute
      </publisher></rdf:Description><rdf:Description rdf:about="rdf:#8194e217-9709-4ea3-8220-045c3ec7b106"><rdf:value><rdf:Description rdf:nodeID="n1"/></rdf:value><subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en">keyword</subject_type></rdf:Description><rdf:Description rdf:about="#I1"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">I1</title><alternative xmlns="http://purl.org/dc/terms/" xml:lang="en">infected hosts</alternative></rdf:Description><rdf:Description rdf:about="#wodarz_2003"><title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en">
        Wodarz's 2nd 2003 mathematical model of the evolution of immunological 
        memory and the regulation of competition between pathogens.
      </title><cmeta:comment><rdf:Description rdf:about="rdf:#ac74765f-535e-48f3-9de6-4ae538fbec5c"/></cmeta:comment><simulation xmlns="http://www.cellml.org/metadata/simulation/1.0#"><rdf:Description rdf:nodeID="n2"/></simulation><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#e5516720-c1e2-4ae2-bbe6-925204067934"/></reference><reference xmlns="http://www.cellml.org/bqs/1.0#"><rdf:Description rdf:about="rdf:#c2ffd07b-2fd1-41b0-bc90-b086171bc24c"/></reference></rdf:Description></rdf:RDF></model>