Difference between revisions of "Team:Tokyo Tech/Model"
| Line 255: | Line 255: | ||
</div> | </div> | ||
<div class="mbr-figure pt-5"> | <div class="mbr-figure pt-5"> | ||
| − | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/2/2a/T--Tokyo_Tech--modeling12.png" alt="Mobirise" title="" style="width: 60%;"> |
</div> | </div> | ||
</div> | </div> | ||
| Line 440: | Line 440: | ||
</div> | </div> | ||
<div class="mbr-figure pt-5"> | <div class="mbr-figure pt-5"> | ||
| − | <img src="https://static.igem.org/mediawiki/2018/ | + | <img src="https://static.igem.org/mediawiki/2018/d/d2/T--Tokyo_Tech--modeling3.png" alt="Mobirise" title="" style="width: 60%;"> |
</div> | </div> | ||
</div> | </div> | ||
| Line 488: | Line 488: | ||
<figure class="mbr-figure container"> | <figure class="mbr-figure container"> | ||
<div class="image-block" style="width: 36%;"> | <div class="image-block" style="width: 36%;"> | ||
| − | <img src=" | + | <img src="https://static.igem.org/mediawiki/2018/6/63/T--Tokyo_Tech--modeling7.png" width="1400" alt="Mobirise" title=""> |
</div> | </div> | ||
| Line 500: | Line 500: | ||
<figure class="mbr-figure container"> | <figure class="mbr-figure container"> | ||
<div class="image-block" style="width: 36%;"> | <div class="image-block" style="width: 36%;"> | ||
| − | <img src=" | + | <img src="https://static.igem.org/mediawiki/2018/5/54/T--Tokyo_Tech--modeling8.png" width="1400" alt="Mobirise" title=""> |
</div> | </div> | ||
| Line 530: | Line 530: | ||
<figure class="mbr-figure container"> | <figure class="mbr-figure container"> | ||
<div class="image-block" style="width: 66%;"> | <div class="image-block" style="width: 66%;"> | ||
| − | <img src=" | + | <img src="https://static.igem.org/mediawiki/2018/9/9e/T--Tokyo_Tech--modeling1.png" width="1400" alt="Mobirise" title=""> |
</div> | </div> | ||
| Line 596: | Line 596: | ||
<figure class="mbr-figure container"> | <figure class="mbr-figure container"> | ||
<div class="image-block" style="width: 84%;"> | <div class="image-block" style="width: 84%;"> | ||
| − | <img src=" | + | <img src="https://static.igem.org/mediawiki/2018/2/28/T--Tokyo_Tech--modeling5.png" width="1400" alt="Mobirise" title=""> |
</div> | </div> | ||
| Line 644: | Line 644: | ||
<figure class="mbr-figure container"> | <figure class="mbr-figure container"> | ||
<div class="image-block" style="width: 39%;"> | <div class="image-block" style="width: 39%;"> | ||
| − | <a href="https://github.com/iGEMTokyoTech/Dengue-Prediction" target="_blank"><img src=" | + | <a href="https://github.com/iGEMTokyoTech/Dengue-Prediction" target="_blank"><img src="https://static.igem.org/mediawiki/2018/f/ff/T--Tokyo_Tech--modeling9.png" width="1400" alt="Mobirise" title=""></a> |
</div> | </div> | ||
Revision as of 18:08, 17 October 2018
<!DOCTYPE html>
Modeling
Establishment of dengue virus serotype prediction system
Abstract
In the past iGEM, the measurement theory model of artificial gene circuits has been centrally treated. We tried to apply a new model to predict the fluctuation in the number of infected people with dengue.
We used the data of yearly occurrence of dengue serotypes by primary or secondary dengue infection in children along with climate data in Bangkok, Thailand from 1973 to 1999 to create a model to predict change in the number and ratio of patients in each serotype.
This model can predict an epidemic of each serotype of dengue fever, thus can be used by health care institutions to prepare vaccine production distribution and production, furthermore allowing them to take special care of secondary dengue infection patients.
From the data we found that some region had distinct epidemic cycle of each serotype while some region had more chaotic pattern, therefore we decided to use simplex projection method for our prediction model, which is suitable for distinguishing chaotic time series from random noise.
If the dengue fever inspection method we developed in the Wet experiment spreads and the data about patients infected with dengue virus is fulfilled, it is possible to estimate the number of infected people more precisely. This can be expected to be useful for measures against abandoned tropical infections.
Introduction
Dengue is a mosquito-borne viral disease that has rapidly spread in all regions of World Health Organization in recent years. Dengue is widespread throughout the tropics, with local variations in risk influenced by rainfall, temperature and unplanned rapid urbanization.
There are 4 distinct, but closely related, serotypes of the virus that cause dengue (DEN-1, DEN-2, DEN-3 and DEN-4). Recovery from infection by one provides lifelong immunity against that particular serotype. However, cross-immunity to the other serotypes after recovery is only partial and temporary. Subsequent infections by other serotypes increase the risk of developing severe dengue. Severe dengue is a leading cause of serious illness and death among children in some Asian and Latin American countries.
Reference: http://www.who.int/en/news-room/fact-sheets/detail/dengue-and-severe-dengue
Overview of the model
We predicted how the number of people infected with dengue virus and the proportion of each serotype changes by using a mathematical method, Simplex Projection to the data about the yearly occurrence of dengue serotypes by primary or secondary dengue infection in children and climate data in Bangkok, Thailand from 1973 to 1999. We constructed a model that automatically predicts the number of infected individuals and the variation of each serotype from dengue fever dengue-infected subjects data and climate data.
Simplex projection
Simplex projection, originally proposed by George Sugihara and Bob May in 1990 is a powerful tool for detecting patterns (chaos) in what otherwise appears to be randomness. Moreover, it will also tell you how complicated this pattern is, and then predict the future, often with far greater accuracy than any other technique.
Ignoring ties in distances, minimal distances, minimal weights and other potential hazards, the following algorithm performs Simplex projection to predict Tp time-steps ahead.
Notation
| Description | |
|---|---|
| E | E denotes the embedding dimension. |
| K | k denotes the number of nearest neighbors we use. For the simplex method, the default is k = E + 1 but for the S-map method it can be much larger. |
| Tp | Tp denotes how many time-steps into the future we are trying to predict. |
| X | X ∈ R denotes a (potentially long) time series. |
| y | y ∈ R is a vector of lagged observations for which we want to make a prediction — in the simplest case where all components of the vector are single time step lags, y1 represents the current value, y2 is the value one time step prior and yE is the value E − 1 time steps prior. |
Figure1 Algorithm of Simplex projection
Figure2 Simple illustration of algorithm of Simplex projection
Estimation of distribution of dengue serotypes
We predicted how the number of people infected with dengue virus and the proportion of each serotype changes by using a mathematical method, Simplex Projection to the data about the yearly occurrence of dengue serotypes by primary or secondary dengue infection in children in Bangkok, Thailand from 1973 to 1999 and climate data in Bangkok, Thailand from 1973 to 1999.
In this analysis, we evaluate the prediction accuracy by rho (Pearson’s correlation coefficient between predictions and observations) values. The following formula is the definition of rho. It can take a range of values from +1 to -1. A value of 0 indicates that there is no association between the two variables. A value greater than 0 indicates a positive association; that is, as the value of one variable increases, so does the value of the other variable. A value less than 0 indicates a negative association; that is, as the value of one variable increases, the value of the other variable decreases.
Figure3 Definition of rho
The values of rho in each prediction are shown in the table below. If the data are increased more, there is no doubt that the accuracy of the prediction model will improve.
| Prediction type | Pearson’s correlation coefficient between predictions and observations |
|---|---|
| Change in scaled subjects of DENV1 | 0.691924 |
| Change in scaled subjects of DENV2 | 0.4864992 |
| Change in scaled subjects of DENV3 | 0.3443239 |
| Change in scaled subjects of DENV4 | 0.4620425 |
| Change in existence rate of DENV1 | 0.6788994 |
| Change in existence rate of DENV2 | 0.7765715 |
| Change in existence rate of DENV3 | 0.6569529 |
| Change in existence rate of DENV4 | 0.8263565 |
The results of the predictions are shown in the following figures.
Figure4 Prediction of scaled subjects number infected with each dengue virus in Bangkok, Thailand from 1973 to 2000
Figure5 Prediction of rate of each dengue virus in Bangkok, Thailand from 1973 to 2000
Co-prediction
The goal of co-prediction is to quantify dynamic similarity between two time series. Given two time series, x and y, we assume that the dynamics can be represented as:
Then co-prediction is a way to quantify how closely F and G resemble each other.
We thought that there might be something related to the temperature affecting the dengue-mediated mosquito activity and the total number of patients in each serotype of dengue.
In order to identify the relationship, co-prediction was performed using the total number of patients in each serotype and the temperature data from 1973 to 1993 in Bangkok, Thailand.
The time series x is total value of infected individuals of each serotype and the time series y is the temperature value.
The figure below shows how much accuracy of prediction was improved by co-prediction.
Figure6 Co-prediction by data of total infected individuals of each serotype (x) and temperature data (y)
Here, we see that the predictions of total infected individuals of each serotype (x) are better when we try to use the inferred dynamics from temperature data (y) to make forecasts. Likewise, the predictions of temperature data (y) have better forecast skill by using the inferred dynamics from total infected individuals of each serotype (x).
Since the time series come from a model simulation where x and y have coefficients, we can infer their dynamical maps are actually related.
Predicton model
In the above fluctuation prediction, prediction accuracy could be improved by using climate data for fluctuation prediction in addition to dengue virus infected patients’ data. It is expected that the prediction of the number of people infected with dengue virus of each serotype will be carried out accurately when the data of infected people with dengue fever increase. Here we constructed a prediction model for the fluctuation in the number of dengue infected people in the coming year from dengue fever infected person data and climate data. With this model we can see which serotype dengue virus will become popular next year from past data.
Figure7 Conceptual diagram of the prediction model
The GitHub repository about this prediction model is here.
Discussion
When using Simplex projection, it is said that time series data of at least 35 to 40 points is required. Despite very few data on dengue fever infected people, we were able to make predictions with some degree of accuracy. Furthermore, by adding climate data to infected data, it is possible to improve prediction accuracy.If we can gather data on dengue-infected individuals of each serotype from Wet's dengue infection detection kit created in this project, we can improve the accuracy of prediction and further help prevent infection with dengue.
Reference
http://www.who.int/en/news-room/fact-sheets/detail/dengue-and-severe-dengue
http://deepeco.ucsd.edu/simplex/
https://cran.r-project.org/web/packages/rEDM/vignettes/rEDM-tutorial.html
https://cran.r-project.org/web/packages/rEDM/vignettes/rEDM-coprediction.html
https://cran.r-project.org/web/packages/rEDM/vignettes/rEDM-algorithms.pdf
https://statistics.laerd.com/statistical-guides/pearson-correlation-coefficient-statistical-guide.php
http://www.ajtmh.org/content/journals/10.4269/ajtmh.2003.68.191?ijkey=a7c6f78b661b24d3734e97302ea1bc8e754b0b01&keytype2=tf_ipsecsha#dataandmedia
