| <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> |
| |
| <html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns="http://www.w3.org/TR/REC-html40"> |
| <head> |
| <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> |
| |
| <title>help</title> |
| </head> |
| |
| <body> |
| <h1>Lyapunov View</h1> |
| <p>The <em>Lyapunov View</em> allows a user to apply tools developed to analyze |
| dynamical systems. The RMS comparison view is a useful measure of the average |
| difference between two scenarios. However, if a model and reference scenario |
| each describe an epidemic that begins and end in the same state (zero |
| infectious), the RMS error will eventually fall to zero as a function of |
| time, even in case of a “bad” model. In addition to measuring the average error, |
| it is useful to look for other measures that might provide a “fingerprint” for |
| the spatiotemporal dynamics of an infectious disease. Like many dynamical |
| systems, infectious disease is a process of many variables. However, it is often |
| possible to capture the essential dynamics by looking at just a few system |
| variables in an appropriate phase space. In its most general formalism, any |
| dynamical system is defined by a fixed set of equations that govern the time |
| dependence of the system’s state variables. The state variables at some instant |
| in time define a point in phase space. A SEIR model defines a four dimensional |
| phase space. An SI model defines a two dimensional space. Examining a reduced |
| set of dimension may be thought of as taking slice through phase space (for |
| example in the SI plane). |
| <p class="MsoNormal" style="text-indent:11.5pt">At the state of the system |
| changes with time, the point (S(t), I(t)) in phase space defines a <i>trajectory</i> |
| in the SI plane. Consider an epidemic that begins with one infectious person and |
| virtually the entire population susceptible at t=0, S(0) ~ 1. The trajectory |
| will begin at time zero along the S axis near 1. As the disease spreads, the |
| susceptible population (S) will decrease and the infectious population (I) will |
| increase. The detailed shape of the this trajectory will depend on the time it |
| takes for the disease to spread to different population centers, as well as the |
| (susceptible) population density function. The peaks and valleys along the |
| trajectory in SI phase space proved a signature or fingerprint for an epidemic |
| the shape of which depends on the disease, the disease vectors, the population |
| distribution, etc. The mathematics of dynamical systems provide us with a |
| formalism to compare trajectories in a phase space. Given a single set of rules |
| (e.g., a disease model), two simulations that begin infinitesimally close |
| together in phase space may evolve different in time and space. This separation |
| in phase space can be measure quantitatively. </p> |
| <p class="MsoNormal" style="text-indent:11.5pt">The Vector</p> |
| <p class="MsoNormal" style="text-indent:11.5pt" align="center"> <span style="font-size: |
| 9.0pt"><img border="0" src="img/lyapEq2.jpg" width="147" height="39"></span> </p> |
| <p class="MsoNormal" style="text-indent:11.5pt">defines a trajectory in SI space. The initial separation at |
| time zero be defined as </p> |
| <p class="MsoNormal" style="text-indent:11.5pt" align="center"> |
| <img border="0" src="img/lyapEq3.jpg" width="66" height="38">. |
| </p> |
| <p class="MsoNormal" style="text-indent:11.5pt">The rate of separation of two trajectories in phase space will often obey the |
| equation</p> |
| <p class="MsoNormal" align="center" style="text-align:center;text-indent:11.5pt"> |
| <img border="0" src="img/lyapEq1.jpg" width="234" height="48"></p> |
| <p class="MsoNormal" style="text-indent:0in"> where |
| <span style="font-family:Symbol">l </span>is the Lyapunov Exponent. This |
| exponent is a characteristic of the dynamical system that defines the rate of |
| separation of infinitesimally close trajectories in phase space. </p> |
| <p>To use the Lyapunov view:<ol> |
| <li>Enter the <a href="../perspectives/analysis.html">Analysis Perspective</a></li> |
| <li>Click on the Lyapunov Tab</li> |
| <li>Use the Select Folder buttons to chose the folders containing the data |
| you wish to compare. The files should have the following |
| <a href="csvloggerview.html"> format.</a></li> |
| <li>Click "Compute Lyapunov Exponent"</li> |
| <li>Two charts will appear, the left hand chart will show the trajectories |
| in phase space (I vs. S) for the two data sets. The right hand chart will |
| show the Log of the integrated<br> |
| difference between the two trajectories as a function of time. The exponent |
| is the initial slope of this time varying difference plotted on a semi-log |
| plot.</li> |
| </ol> |
| </p> |
| <p><img border="0" src="img/Lyapunov.jpg" width="1000" height="608"> </p> |
| </body> |
| </html> |