# A critical review of mathematical models and data used in diabetology

- A Boutayeb†
^{1}Email author and - A Chetouani†
^{1}

**5**:43

**DOI: **10.1186/1475-925X-5-43

© Boutayeb and Chetouani; licensee BioMed Central Ltd. 2006

**Received: **08 May 2006

**Accepted: **29 June 2006

**Published: **29 June 2006

## Abstract

The literature dealing with mathematical modelling for diabetes is abundant. During the last decades, a variety of models have been devoted to different aspects of diabetes, including glucose and insulin dynamics, management and complications prevention, cost and cost-effectiveness of strategies and epidemiology of diabetes in general. Several reviews are published regularly on mathematical models used for specific aspects of diabetes. In the present paper we propose a global overview of mathematical models dealing with many aspects of diabetes and using various tools. The review includes, side by side, models which are simple and/or comprehensive; deterministic and/or stochastic; continuous and/or discrete; using ordinary differential equations, partial differential equations, optimal control theory, integral equations, matrix analysis and computer algorithms.

## Introduction

It is now commonly admitted that diabetes is sweeping the globe as a silent epidemic largely contributing to the growing burden of non-communicable diseases and mainly encouraged by decreasing levels of activity and increasing prevalence of obesity [1–5]. Dramatic increase has occurred in both prevalence and incidence of diabetes, especially with the new threshold proposed by the Expert Committee on the diagnosis and classification of diabetes mellitus in 1997 [6] and adopted by the World Health Organization. During the last decades, a huge number of papers were published on different aspects of diabetes and its complications. In particular, an interesting literature has been devoted to studies collecting, analyzing and validating data concerning diabetes populations. A variety of mathematical models, statistical methods and computer algorithms have been proposed in order to understand different aspects of diabetes such as: glucose-insulin dynamics, epidemiology of diabetes and its complications, cost of diabetes and cost-effectiveness of strategies dealing with diabetes. Several reviews have been devoted to mathematical models and diabetes. In their majority, these reviews concentrated on specific aspects of diabetes such as glucose-insulin dynamics [7–13], computer algorithms and devices [14–16], sensors and control [17, 18], mathematical and software aspects [19], glycemic index [20], burden and cost of diabetes [21, 22]. On the one side, it is difficult to include in one review mathematical models published on different aspects of diabetes. On the other side, it would be very useful for researchers in this field to get a panorama of the models published so-far and their references. The benefit will be greater if such a review is published with free online access, especially for physicians and researchers in developing countries. The challenge is worth trying and the present paper is proposing a non-exhaustive overview including model structure, computer algorithms, data used as input and model validation.

## Glucose-insulin dynamics

### Mathematical models

The majority of mathematical models proposed in the literature were devoted to the dynamics of glucose-insulin, including Intra Venous Glucose Tolerance Test (IVGTT), Oral Glucose Test (OGTT), Frequently Sampled Intravenous Glucose Tolerance Test (FSIGT). In 1939, Himsworth and Ker [23] introduced the first approach to measure the insulin sensitivity *in vivo*. Mathematical models have been used to estimate the glucose disappearance and insulin-glucose dynamics in general. Bolie is among the pioneers in this field. In 1961, using ordinary differential equations, he proposed the following simple model [24]:

where *G* = *G*(*t*) represents the glucose concentration, *I = I*(*t*) represents the insulin and *p*, *a*
_{1}, *a*
_{2}, *a*
_{3}, *a*
_{4} are parameters.

Although various models (simple and comprehensive) were proposed by different authors [25–28] (see also the 244 references in the review by Bergman et al. (1985)) [7], especially those dealing with insulin sensitivity, the real start of modeling the glucose-insulin dynamics is thought to begin with the so-called minimal model proposed by the team of Bergman and Cobelli in the early eighties [29, 30]. The model was formulated as follows:

where, (*G*(*t*) *- p*
_{5})^{+} = *G*(*t*) *- p*
_{5} if *G*(*t*) > *p*
_{5} and 0 otherwise.

*X*(*t*) denotes an auxiliary function representing insulin-excitable tissue glucose uptake activity, *G*
_{
b
}, and *I*
_{
b
}are the subject's baseline glyceamia and insulinimia. *b*
_{0} – *b*
_{7} are parameters. It should be stressed that, although equations (3)–(5) were developed for describing FSIGT data, these equations were presented into two parts. Part 1 with equations (3)–(4) and part 2 using equation (5). A large number of papers have been published, using modified versions of the glucose minimal model (3)–(4) for describing OGTT and meal tests, while insulin minimal models derived from (5) are still limited to IVGTT. It should also be stated that the major contribution of the glucose minimal model (3)–(4) has been to provide means of estimating insulin sensitivity *S*
_{
I
}= *p* 3/*p* 2, avoiding the glucose clamp.

The same authors have further published papers, completing, testing or validating the results of the minimal model [7, 31–33]. An indication of the importance of the minimal model and subsequent research for diabetes understanding is given by the 2006 Banting medal awarded by the American Diabetes Association to Professor Bergman for his achievements.

Variant versions based on the minimal model were considered by different authors. An example of this category was proposed by Derouich and Boutayeb who used a modified version of the minimal model to introduce parameters related to physical exercise [34]:

where *q*
_{1}, *q*
_{2}, *q*
_{3} are parameters related to physical activity and defined as follows:

*q*
_{1} : the effect of physical exercise in accelerating the utilization of glucose by muscles and the liver insulin

*q*
_{2} : the effect of physical exercise in increasing the muscular and liver sensibility to the action of insulin.

*q*
_{3} : the effect of physical exercise in increasing the utilization of insulin. In other words, *q*
_{3} increases insulin effectiveness in enhancing glucose disposal and consequently improving insulin sensitivity to become: *S*
_{
I
}= (*p*
_{3} + *q*
_{3})(1 + *p*
_{2})*/P*
_{2}.

According to a paper published by Bergman in 2002 [35], more than 500 studies related to the minimal model can be found in the literature. More information on this historic model and related models can be found in literature [8, 9].

*p*

_{0}-

*p*

_{7}), and has been indisputably useful in physiological research, it has the following drawbacks:

- 1.
The model, as originally proposed, is to be regarded as composed of two separate parts. The first part uses equations (3) and (4) and the second part uses equation (5). For the last part, plasma glucose concentration is to be regarded as a known forcing function. In other words, the model parameter fitting has to be conducted in two steps: first, using the recorded insulin concentration as input data in order to derive the parameters in the two first equations, then using the recorded glucose as input data to derive the parameters in the third equation.

- 2.
Some of the mathematical results produced by this model are not realistic(problems of positive equilibrium and solutions not bounded).

- 3.
The artificial non-observable variable

*X*(*t*) is introduced to take account of the delay in the action of insulin.

Taking into account these remarks and stressing that the glucose-insulin system is an integrated physiologic dynamical system which should be dealt with as a whole, De Gaetano and Arino [36, 37] proposed an aggregated delay differential model called a dynamical model:

After renaming the parameters, the dynamical model takes the form [36]:

with *G*(*t*) = *G*
_{
b
}for -*b*
_{5} ≤ *t* < 0

Mukhopadhyay et al. (2004) recalled that this model has been shown to allow simultaneous estimation of both insulin secretion and glucose uptake parameters, to have positive, bounded solutions, and to be globally asymptotically stable around the pre-injection equilibrium blood glucose and insulin concentrations. They proposed an extension by introducing a generic weight function *ω* in the delay integral kernel for the pancreatic response to glucose. The new model obtained is as follows:

with *G*(*t*) = *G*
_{
b
}for *t* < 0

A more general model was proposed by Li et al. (2001). The authors noted that, while the dynamical model solves the problems of the minimal model, it implicitly or explicitly made a few assumptions that may not be necessary or realistic. Specifically, some of the interaction terms are too special and thus too restrictive. For example, the term *b*
_{4}
*I*(*t*)*G*(*t*) assumes that mass action law applies here while a more popular, general and realistic alternative is to replace this term by *b*
_{4}
*I*(*t*)*G*(*t*)*/*(*αG*(*t*) + 1). The way the delay is introduced is also restrictive. Consequently, the model proposed is the following:

with *G*(*t*) = *G*
_{
b
}for *-b*
_{5}
*≤ t <* 0 and *G*
_{
t
}(*θ*) = *G*(*t* + *θ*), *t* > 0, -*b*
_{5} ≤ *θ* < 0.

Other models of the glucose-insulin dynamics, using optimal control or partial differential equations were proposed by different authors [18, 40–44]. Cobelli and Tomaseth [41] discussed the optimal input design in a model of glucose kinetics. They proposed the following model:

*y*(*t*) = *c*(*p*)*x*(*t*) (15)

*z*(*t*) = *y*(*t*) + *e*(*t*) (16)

*h* [*x*(*t*), *u*(*t*), *p*] ≥ 0 (17)

Lam et al. [44] used a slightly modified version of the minimal model for the assessment of insulin sensitivity. An interesting survey of mathematical models using control for glucose-insulin and management of diabetes is given by Palerm in his Ph. D thesis, with some 350 references [18]. The author focussed on the Direct Model Reference Adaptive Control (DMRAC) and its reformulation.

The formulation of the general DMRAC algorithm is based on the following system

or more generally,

where *x*(*t*) is the (*nx* 1) state vector, *u*(*t*) is the (*mx* 1) control vector, *y*(*t*) is the (*qx* 1) output vector, *A*, *B*, *C* and *D* are matrices with appropriate dimensions. The objective is to find, without explicit knowledge of *A* and *B*, a control *u*(*t*) such that the output vector *y*(*t*) fellows a reference model.

### Computer algorithms

The regulation of blood glucose concentration is mainly achieved by acting on three control variables: insulin, meals and physical exercise. However, as stressed by Bellazzi et al. [15], the quantitative evaluation of meals and physical effort still represents a major problem in home monitoring. Consequently, the quasi totality of proposed control systems have focused on insulin therapy strategies. A number of devices-microsystems and computer approaches- have been reported in the literature with open, closed and partially closed algorithms, Selam and Charles [45], Lehman and Deutsch [46, 47]. Two reviews on the intravenous route to blood glucose control and subcutaneous route to insulin dependent diabetes therapy were recently published by Bellazzi et al. [15] and Parker et al. [16]. A partial list of software packages, commercially and freely available is given in a recent review by Makroglou et al. [19].

### The closed-loop strategy

An example of models with closed-loop strategy is a wearable artificial pancreas as proposed by Shimoda et al. [18, 48], based on the assumption that the relationship between plasma insulin and blood glucose concentration in a normal subject during an oral glucose bolus is as follows:

where *I* denotes the plasma insulin concentration; *G* is the blood glucose concentration; and *a*, *b* and *c* are parameters that can be estimated by nonlinear least squares method [49].

The insulin dynamics is described by the following ordinary differential system:

where *V* is the plasma volume, *IIR* represents the Insulin Infusion Rate; *X*, *Y* and *Z* are the insulin masses in the two subcutaneous compartments and in plasma respectively.

### The open-loop approach

According to the review by Bergman et al. [7], two distinct classes of methodologies were applied to open the glucose/insulin feedback relationship *in vivo*. The first method pioneered by Reaven and his colleagues [50], and labelled as the pancreatic or insulin suppression test, utilized pharmacologic means to render the pancreas blind to plasma glucose concentration. The second approach, labelled as the glucose clamp, was proposed by Andres [51]. The method uses a variable glucose infusion to establish a relatively constant plasma glucose concentration with or without exogenous insulin. Mathematical models are also proposed to deal with the deterioration of beta-cell [52–54].

### Epidemiological models applied to diabetes

Historically, since the first model of smallpox formulated by Bernoulli in 1760, an abundant literature was devoted to mathematical models dealing with communicable diseases such as measles, rubella, malaria, influenza, AIDS, dengue and others [56]. As indicated by a review published by Hethcote in 2000 [57], a tremendous variety of models have been formulated, mathematically analyzed, and applied to infectious diseases. Modelling has thus become an interesting tool providing conceptual results such as thresholds, basic reproduction numbers, contact numbers, and replacement numbers. Application of similar models for non communicable diseases is rather unusual. In this way, few authors have proposed epidemiological models for diabetes and obesity [58–66]. In [59], Boutayeb and Derouich considered two discrete models for the evolution from diabetes without complications to the stage of diabetes with complications.

In [60], using partial differential equations, the authors proposed an age structured continuous model for complications of diabetes. Supposing that *C* = *C*(*a*, *t*) and *D* = *D*(*a*, *t*) represent the numbers of diabetics with and without complications aged *a* at time *t*, respectively, and *n*(*a*, *t*) = *C*(*a*, *t*) + *D*(*a*, *t*) the size of the population of diabetics aged *a* at time *t*, different scenarios with different values were used for the following parameters: natural death rate (*d*(*a*, *t*)), death rate due to complications (*δ*(*a*, *t*)), incidence of diabetes with and without complications (*I*
_{1}(*a*, *t*)), (*I*
_{2}(*a*, *t*)) rate at which complications are developed (*p*(*a*, *t*)) and rate at which complications could eventually be cured (*q*(*a*, *t*)). The main objective of the authors was to show that, although diabetes is not curable at the moment, prevention of its complications (which is possible) would improve peoples quality of life and reduce costs of the national health and social services. Assuming that the number of males is equal to the number of females and that diabetes affects the people of the two sexes equally, the continuous age structured model is formalized by the following partial differential equations:

adding equations (25) and (26) and writing :

*C*(*a*, *t*) = *r*(*a*, *t*)*n*(*a*, *t*) (27)

leads to

In the same spirit, Boutayeb and colleagues [62, 63] proposed linear and non-linear population models of diabetes mellitus, using ordinary differential equations and numerical implementation.

## Data analysis, parameters estimation and validation

### Studies and trials

Worldwide, different studies were devoted to diabetes and its complications. These studies have been used directly or indirectly for data analysis, mathematical modeling and parameters validation. Among the most cited studies, Diabetes Control Complications Trial (DCCT) [67] and UK Diabetes Prevention Study (UKPDS) [68]. The first trial involved 1,441 volunteers with type 1 diabetes and 29 medical centers in the United States and Canada, and have shown that, diabetes complications can be reduced or at least delayed by a good regular glycemic control through intensive insulin therapy consisting of three or more insulin injections per day or in the use of insulin pumps. The main DCCT Study Findings were the following: Lowering blood glucose reduces the risks of eye disease, kidney failure and nerve disease by 76%, 50% and 60% respectively [69–72]. The second trial concerned over 5000 non insulin-dependent patients from 23 centres from all parts of England, Scotland and Northern Irland, showing that complications of diabetes can be prevented by a better control of blood glucose and blood pressure. Among the 70 papers published by the UKPDS group, we cite here some of those using mathematical models. The UKPDS risk engine was used as a model for the risk of coronary heart disease, myocardial infarctus and stroke in Type II diabetes [73–75]. Modeling glucose exposure as a risk factor for photocoagulation in Type II diabetes was considered in [76], whereas the UKPDS Outcomes Model was proposed to estimate the lifetime health outcomes of patients with Type 2 [77]. Data collection, parameters estimation and validation concerned a multitude of other studies such as: the Wisconsin epidemiologic study of diabetic retinopathy (WESDR) [78], Framingham Heart Study (FHS) [79], Diabetes Prevention Program (DPP) [80], Health Outcomes Prevention Evaluation (HOPE) [81, 82], The Health Plan Employer Data and Information Set (HEDIS) [83], Echantillon national temoin représentant des personnes diabétiques (ENTERED) [84], Multiple Risk Factor Intervention Trial (MRFIT) [85], Heart Protection Study (HPS) [86], Cholesterol and Recurrent Events (CARE) [87], ACE Inhibitors and Diabetic Nephropathy Trial (Lewis) [88], the IRMA-2 trial [89], Irbesartan Diabetic Nephropathy Trial (IDNT) [90], the Collaborative AtoRvastatin Diabetes Study (CARDS) [91].

### The archimedes model

- 1.
Biological modeling,

- 2.
Clinical Medicine,

- 3.
Operations research,

- 4.
Economic/system resources.

Thus, the authors present Archimedes as a new type of mathematical model which includes all four components. Archimedes is a very detailed, comprehensive, continuous simulation model. A person-by-person, object-by-object simulation, spanning from biological details to the care processes, logistics, resources, and costs of health care systems [92–94]. The model is written in differential equations for which different levels of detail may be considered [95]. The equations, assumptions, and sources are summarised in an online appendix (available at http://care.diabetesjournals.org). A validation of the Archimedes model of diabetes and its complications or a variety of populations, organ systems, treatments, and outcomes is given by Eddy & Schlessinger [96, 97]. The model was validated against 18 trials of which ten trials explicitly dealing with diabetes. Namely: the Diabetes Control and Complications Trial(DCCT) [67], the U.K. Prospective Diabetes Study (UKPDS) [68], the Diabetes Prevention Program (DPP) [80], the Health Outcomes Prevention Evaluation (HOPE) [81], the diabetes sub-population of the HOPE Trial (Micro-HOPE) [82], the Heart Protection Study (HPS) [86], Cholesterol and Recurrent Events (CARE) [87], the ACE Inhibitors and Diabetic Nephropathy Trial (Lewis) [88], the IRMA-2 trial [89], and the Irbesartan Diabetic Nephropathy Trial (IDNT) [90]. In general, between 10 and 30 equations are needed to represent the pathophysiology of the disease and calculate the effect of a specific treatment on a specific outcome in a specific population (not included equations for behaviors, care processes, logistics, and other nonbiological aspects of the model. As stressed by the Editorial of Diabetes Care [98], functional forms of the equation are given but values of the variables and parts of the model that describe micro- and macro-vascular complications are not provided. However, beyond these limitations, the model was used to predict 74 major outcomes, giving astounding results: In 71 out of the 74 clinical outcomes, the differences between the results calculated by the model and the observed ones were statistically not significant. More information on Archimedes model can be found in [92].

### Other models

Other models and computer algorithms were devoted to the burden, cost and cost-effectiveness of diabetes [99–101], telemedicine and home management of diabetes [102–104]. In a case study paper [104], Wu proposed the following model for self-management of Type 2 diabetes:

where *x* represents blood glucose level over the baseline at time *t* and *ω*
_{0} is the system natural frequency. Finally, in some papers and letters, mathematical models and guidelines for computer modeling of diabetes were subject to debate and criticism [105–107].

## Discussion

- 1.
to be physiologically based,

- 2.
having parameters that can be estimated with a reasonable precision,

- 3.
parameters with values that are reasonable and have physiological interpretation,

- 4.
best able to simulate the dynamics of the system with smallest number of identifiable parameters.

In section 2, variant versions of the minimal model were considered by different authors. For instance, in Derouich and Boutayeb [34], physical exercise was seen to be an interesting tool that improves insulin sensitivity (*S*
_{
I
}= (*p*
_{3} + *q*
_{3})(l + *p*
_{2})/*p*
_{2})). The authors stressed that new control strategies take a long time before they become affordable on a large scale, especially in developing countries where the majority of diabetics are struggling just to get insulin doses and where the price of a blood strip exceeds the individual income.

As indicated earlier, the minimal model has been indisputably useful in physiological research and served as starting point for many other models. The drawbacks raised by De Gaetano and Arino [36] were mainly based on the mathematical formulation (specifically problems of positive equilibrium and solutions not bounded). But they stressed that no criticism is implied regarding the practical usefulness of the minimal model. By the way, they acknowledge that their group uses the minimal model in the routine evaluation of insulin sensitivity in clinical patients. The so-called "dynamic model" proposed by the previous authors has not escaped from criticism since the assumptions made were judged to be not necessary or realistic. This judgement was made by Li et al. [38] who also criticised the restrictive way of using the delay and proposed a more general model. According to these authors, their general model was constructed for the study of IVGTT which focuses on the metabolism of glucose. However, except simulation and mathematical aspects (steady state, oscillatory glucose and insulin levels), no evidence is given on the real performance. For computer algorithms, as stressed earlier, the quasi totality of proposed control systems have focused on insulin therapy strategies. In many cases, the proposition of models are dictated by commercial purposes and the accessibility to their pragmatic use remains restrictive. Finally, among complex comprehensive models, Archimedes model can be seen as the most illustrative. The huge arsenal of computer and mathematical tools used by this model seems to be justified by the first validated results. However, the model has been published recently and more time is needed to get sufficient information for a critical discussion.

## Conclusion

During the last decades, an interesting number of papers have been published on mathematical models and computer algorithms. In the present review, the authors have tried to give a non-exhaustive panorama of the papers which have used mathematical modeling for different aspects of diabetes, including glucose-insulin dynamics, beta-cell function, epidemiology of diabetes, management and the burden of diabetes and its complications.

As indicated in section two, the award of 2006 banting medal by the American Diabetes Association to Professor Bergman for his achievements in diabetes research among which, the famous minimal model, is an indication of the importance of mathematical models for the understanding of diabetes and its management.

## Notes

## Declarations

### Acknowledgements

The paper is dedicated to the students : Wiam, Sara, Rajae, Oumnia, Mariam, Hanae and Fatima-Zohra for their support to diabetic patients.

The authors are grateful to Prof Wald Editor In Chief of the Biomedical Engineering Online for his comments and suggestions and also for granting a waiver for this publication. The authors wish also to thank anonymous reviewers for their constructive comments and suggestions.

## Authors’ Affiliations

## References

- The world health report 2003:
**Today's challenges.**Geneva, World Health Organization; [http://www.who.int/whr/2003/en] - International Diabetes Federation:
**IFD report 2003.**[http://www.idf.org/home/index.cfm] - Boutayeb A, Boutayeb S:
*The burden of non communicable diseases in developing countries.**Volume 4*. International Journal of Equity in Health; 2005:2.Google Scholar - Ratner RE:
**Type 2 diabetes mellitus: the grand overview.***Diabet Med*1998,**14:**54–57.MathSciNetGoogle Scholar - Jiwa F:
**Diabetes in the 1990s -an overview.***Stat Bull Metrop Co*1997,**78**(1):2–8.Google Scholar - The Expert Committee on the Diagnosis and Classification of Diabetes Mellitus:
**Report of the Expert Committee on the diagnosis and Classification of Diabetes Mellitus.***Diabetes Care*1997,**20:**1183–1197.Google Scholar - Bergman RN, Finegood DT, Ader M:
**Assessment of Insulin Sensitivity**in Vivo**.***Endicrine Reviews*1985,**6**(1):45–86.Google Scholar - Bergman RN:
**The minimal model: yesterday, today and tomorrow.**In*The minmal model Approach and Determination of Glucose Tolerance*.*Volume 7*. Edited by: Bergman RN, Lovejoy JC. Boston : Louisiana state university Press; 1997:3–50.Google Scholar - Bergman RN:
**The minimal model of glucose regulation: a biography.**In*Mathematical Modeling in Nutrition and Health*. Edited by: Novotny, Green, Boston. Kluwer Academic/Plenum; 2001.Google Scholar - Nucci G, Cobelli C:
**Models of subcatuneous insulin kinetics: a critical review.***Computer Methods and Programs in Biomedicine*2000,**62:**249–257. 10.1016/S0169-2607(00)00071-7Google Scholar - Sundell J, Knuuti J:
**Insulin and myocardial blood flow.***Cardiovascular Research*2003,**57**(2):312–319. 10.1016/S0008-6363(02)00718-6Google Scholar - Finegood DT:
**Application of the minimal model of glucose kinetics.**In*The minmal model Approach and Determination of Glucose Tolerance*.*Volume 7*. Edited by: Bergman RN, Lovejoy JC. Boston : Louisiana state university Press; 1997:51–122.Google Scholar - Mari A:
**Mathematical modelling in glucose metabolism and insulin secretion.***Current Opinion Clinical Nutrition Metabolism Care*2002,**5:**495–501. 10.1097/00075197-200209000-00007Google Scholar - Lehmann ED, Deutsch T:
**Application of computers in diabetes care: a review. II, Computers for decision support and education.***Med Inform*1995,**20:**303–329.Google Scholar - Bellazzi R, Nucci G, Cobelli C:
**The Subcutaneous Route to Insulin Dependent Diabetes Therapy: Closed-Loop and Partially Closed-Loop Control Strategies for insulin Delivery and Measuring Glucose Concentration.***IEEE Engrg Medicine Biol*2001,**20**(1):54–64. 10.1109/51.897828Google Scholar - Parker RS, DoyleIII FJ, Peppas NA:
**The Intraveneous Route to Blood Glucose Control: A Review of Control Algorithms for Noninvasive Monitoring and Regulation in Type 1 Diabetic Patients.***IEEE Engineering in Medicine and Biologie*2001,**20**(1):65–73. 10.1109/51.897829Google Scholar - Koschinsky T, Heinemann :
**Sensors for glucose monitoring: technical and clinical aspects.***Diabetes/Metabolism Research and Reviews*2001,**17:**113–123. 10.1002/dmrr.188Google Scholar - Palerm CCR:
**Drug Infusion Control: An Extended Direct Model Reference Adaptive Control Strategy.**In*Ph.D Thesis*.*Volume chap 3–4*. Rensselear Polytechnic Institute, Troy, New York; 2003:45–60.Google Scholar - Makroglou A, Li J, Kuang Y:
**Mathematical models and software tools for the glucose-insulin regulatory system and diabetes: an overview.***Applied Numerical Mathematics*2006,**56:**559–573. 10.1016/j.apnum.2005.04.023MathSciNetGoogle Scholar - Kalergis M, De Grandpre E, Andersons C:
**The Role of the Glycemic Index in the Prevention and Management of Diabetes: A Review and Discussion.***Canadian Journal of Diabetes*2005,**29**(1):27–38.Google Scholar - Atun R, Gurol-Urganci I:
*Health expenditure: an 'investment' rather than a cost ?*International Economics Programme; 2005.Google Scholar - Boutayeb A:
**The double burden of communicable and non communicable diseases in developing countries : A review.***Transactions of the Royal Society of Tropical Medicine and Hygiene*2006,**100:**191–199. 10.1016/j.trstmh.2005.07.021Google Scholar - Himsworth HP, Ker RB:
**Insulin-sensitive and insulin insensitive types of diabetes millitus.***Cli Sci*1939,**4:**119–122.Google Scholar - Bolie VW:
**Coefficients of normal blood glucose regulation.***J Appl Physiol*1961,**16:**783–788.Google Scholar - Ackerman E, Gatewood LC, Rosevear JW, Molnar GD:
**Model studies of blood glucose regulation.***Bull Math Biophys*1965,**27:**21–24.Google Scholar - Della C, Romano MR, Voehhelin MR, Seriam E:
**On a mathematical model for the analysis of the glucose tolerance curve.***Diabetes*1970,**19:**145–148.Google Scholar - Serge G, Turcogl M, Varcellone G:
**Modelling blood glucose and insulin kinetics in normal diabetic and obese subjects.***Diabetes*1973,**22:**94–97.Google Scholar - Srinivasan R, Kadish AH, Sridhar R:
**A mathematical model for the control mechanism of free-fatty acid and glucose metabolism in normal humans.***Comp Biomed Res*1970,**3:**146–149. 10.1016/0010-4809(70)90021-2Google Scholar - Bergman RN, Ider YZ, Bowden CR, Cobelli C:
**Quantitative Estimation of Insulin Sensitivity.***Am J Physiol*1979,**23**(6):E667-E677.Google Scholar - Toffolo G, Bergman RN, Finegood DT, Bowden CR, Cobelli C:
**Quantitaive Estimation of beta cell sensitivity to glucose in the intact organism: a Minimal Model of Insulin kinetics in the Dog.***Diabetes*1980,**29:**979–990.Google Scholar - Bergman RN, Bowden CR, Cobelli C:
**The Minimal Model approach to quantification of factors controlling glucose disposal in man.**In*Carbohydrate Metabolism*.*Volume chap 13*. Edited by: Cobelli, Bergman. John Wiley & Sons Ltd; 1981:269–293.Google Scholar - Cobelli C, Mari A:
**Validation of mathematical models complex endocrine-metabomism systems. A case study on a model of glucose regulation.***Med & Biot Eng & Comput*1983,**21:**390–399.Google Scholar - Cobelli C, Nucci G, del Prato S:
**A physiological simulation model of the glucose-insulin system in type I diabetes.***Diabetes, Nutrition & Metabolism*1998,**11**(1):78–80.Google Scholar - Derouich M, Boutayeb A:
**The effect of physical exercise on the dynamics of glucose and insulin.***Journal of Biomechanics*2002,**35:**911–917. 10.1016/S0021-9290(02)00055-6Google Scholar - Bergman RN:
**Pathogenesis and prediction of diabetes mellitus: lessons from integrative physiology.***Mount Sinai J Medicine*2002,**60:**280–290.Google Scholar - De Gaetano A, Arino O:
**Mathematical Modelling of the Intravenous Glucose Tolerance Test.***J Math Biol*2000,**40:**136–168. 10.1007/s002850050007MathSciNetGoogle Scholar - De Gaetano A, Arino O:
**A statistical approach to the determination of stability for dynamical systems modelling physiological processes.***Math Comput Modelling*2000,**31:**41–51. 10.1016/S0895-7177(00)00020-0Google Scholar - Li J, Kuang Y, Li B:
**Analysis of IVGTT Glucose-Insulin Interaction Models with time delay.***Discrete and Continous Dynamical Systems Series B*2000,**1**(1):103–124.MathSciNetGoogle Scholar - Mukhopadhyay A, De Gaetano A, Arino O:
**Modelling the intra-venous glucose tolerance test: A global study for single-distributed-delay model.***Discrete and Continous Dynamical Systems Series B*2004,**4**(2):407–417.MathSciNetGoogle Scholar - Cobelli C, Thomaseth K:
**Optimal input design for identification of compartimental models : theory and applications to a model of glucose kinetics.***Mathematical Biosciences*1985,**77:**267–270. 10.1016/0025-5564(85)90101-4MathSciNetGoogle Scholar - Cobelli C, Thomaseth K:
**The minimal model of glucose disappearence : optimal input studies.***Math Biosciences*1987,**83:**127–130. 10.1016/0025-5564(87)90107-6MathSciNetGoogle Scholar - Cobelli C, Thomaseth K:
**An optimalality of the impulse input for linear system identification.***Math Biosciences*1988,**89:**127–129. 10.1016/0025-5564(88)90094-6MathSciNetGoogle Scholar - Cobelli C, Thomaseth K:
**Optimal equidose inputs and role of measurement error for estimating the parameters of a compartmental model of glucose kinetics from continuous and discrete time optimal examples.***Math Biosciences***89:**135–137. 10.1016/0025-5564(88)90095-8 - Lam ZH, Hwang KS, Lee JY, Chase JG, Walker GC:
**Active insulin infusion using optimal and derivative weighted control.***Medical engineering physics*2002,**24:**663–672. 10.1016/S1350-4533(02)00147-9Google Scholar - Salam J, Charles MA:
**Devices for insulin administration.***Diabetes care*1999,**13:**955–979.Google Scholar - Lehman ED, Deutsch T:
**A physiological model of glucose insulin interaction in type 1 diabetes mellitus.***J Biomed Eng*1992,**14:**235–242.Google Scholar - Lehman ED, Deutsch T, Rondsari AV, Carson ER, Sonksen PH:
*Validation of a metabolic prototype to assist in the treatment of insulin dependent diabetes mellitus.**Volume 18*. Med Inform London; 1993:83–101.Google Scholar - Shimoda S, Nishida K, Sakakida M, Konno Y, ishinose K, Uehara M, Nowak T, Shichiri M:
**Closed-loop subcutaneous isulin infusion algorithm with a short acting insulin analog for long-term clinical application of a wearable artificial endocrine pancreas.***Frontiers Med Biol Eng*1997,**8:**197–211.Google Scholar - Hashiguchi Y, Sakakida M, Nishida K, Uemura T, Kajiwara K, Shichiri M:
**Development of a miniaturized glucose monitoring system by combining a needle-type glucose sensor with microdialysis sampling method. Long-term subcutaneous tissue glucose monitoring in ambulatory diabetic patients.***diabetes care*1994,**17**(5):387–396.Google Scholar - Reaven GM, Doberne L, Greenfield MS:
**Comparaison of insulin secretion and**in vivo**in-sulin action in non-obese and moderately obese individuals with non-insulin dependent diabetes mellitus.***Diabetes*1982,**31:**382–388.Google Scholar - Andres R, Swerdloff R, Pozefsky T, Coleman D:
**Manual feedback technique for the control of blood glucose concentration.**In*Automation in Analytical Chemistry*. Edited by: Skeggs Jt LT. Mediad Inc, New York; 1966:486–490.Google Scholar - Bagust A, Beale S:
**Deteriorating beta-cell function in type 2 diabetes: a long-term model.***Q J Med*2003,**96:**281–288.Google Scholar - Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC:
**Homeostasis model assessment: insulin resistance and beta cell function from fasting plasma glucose and insulin concentrations in man.***Diabetologia*1985,**28:**412–419. 10.1007/BF00280883Google Scholar - Levy JC, Matthews DR, Hermans MP:
**Correct Homeostasis Model Assessment(HOMA) evaluation uses the computer program(letter).***Diabetes Care*1998,**20:**2191–2192.Google Scholar - Bernoulli D:
**Essai d'une nouvelle analyse de la mortalité causée par la petite vérole et des avantages de 1'inoculation pour la prévenir.**In*mémoire de mathématiques et de physique*. Edited by: Académie Royale des Sciences. Paris; 1760:1–45.Google Scholar - Hethcote HW:
**A thousand and one epidemic models.**In*Frontiers in Theoretical Biology*. Edited by: Levin. 100, Springer-Verlag, Berlin; 1994:504–515. Lecture notes in BiomathGoogle Scholar - Hethcote HW:
**The Mathematics of Infectious Diseases.***SIAM REVIEW*2000,**42**(4):599–653. 10.1137/S0036144500371907MathSciNetGoogle Scholar - Boutayeb A, Kerfati A:
**Mathematical models in diabetology modelling.***Measurement and Control, C, AMSE*1994,**44:**53–63.Google Scholar - Boutayeb A, Derouich M:
**Age structured models for diabetes in East Morocco.***Mathematics and Copmuter Simulation*2002,**58:**215–229. 10.1016/S0378-4754(01)00368-8MathSciNetGoogle Scholar - Boutayeb A, Twizell EH:
**An age structured model for complications of diabetes mellitus in Morocco.***Simulation Modelling Practice and Theory*2004,**12:**77–87. 10.1016/j.simpat.2003.11.003Google Scholar - Boutayeb A, Chetouani A:
**Dynamics of a disabled population in Morocco.***Biomedical Engineering Online*2003,**2:**2. 10.1186/1475-925X-2-2Google Scholar - Boutayeb A, Twizell EH, Achouyab K, Chetouani A:
**A mathematical model for the burden of diabetes and its complications.***Biomedical Engineering Online*2004,**3:**20. 10.1186/1475-925X-3-20Google Scholar - Boutayeb A, Chetouani A, Achouyab K, Twizell EH:
**A non-linear population model of diabetes mellitus.***Journal of Applied Mathematics and computing*2006,**21:**127–139.MathSciNetGoogle Scholar - James RC, Blanchard JF, Campbell D, Clottey C, Osei W, Svenson LW, Noseworthy TW:
**A model for non-communicable disease surveillance in Canada: the prairie pilot diabetes surveillance system.***Chronic Dis Can*2004,**25**(1):7–12.Google Scholar - Paterson KR:
**Diabetes epidemiology: present and future.***Diabet Med*1999,**16**(6):446–7. 10.1046/j.1464-5491.1999.00108.xGoogle Scholar - West KM:
**Epidemiology of diabetes and its macrovascular complications.***Diabetes Care*1979,**2**(2):63–4.Google Scholar - DCCT Research Group:
**The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin dependent diabetes mellitus.***New England Journal of Medicine*1993,**329:**977–986. 10.1056/NEJM199309303291401Google Scholar - O'Neill S:
**The UK Prospective Diabetes Study(UKPDS): Questions & Answers.**[http://www.diabetes.uk] - Rohlfing CL, Wiedmeyer HM, Little RR, England JD, Tennill A, Goldstein DE:
**Defining the relationship between plasma glucose and HbAlc: analysis of glucose profiles and HbAlc in the Diabetes Control and Complications Trial.***Diabetes Care*2000,**25:**275–278.Google Scholar - Diabetes Control and Complications Trials Research Group:
**The relationship of glycemic expsure (HbA**_{ 1c }**) to the risk of development and progression of retinopathy in the diabetes control complications trial in insulin-dependent diabetes mellitus.***Diabetes*1995,**44:**968–983.Google Scholar - Diabetes Control and Complications Trials Research Group:
**The absence of a glycemic threshold for the development of long-term complication: the perspective of the diabetes control complications trial in insulin-dependent diabetes mellitus.***Diabetes*1996,**45:**1289–1298.Google Scholar - Diabetes Control and Complications Trials Research Group:
**Hypoglycemia in the Diabetes Control Complications Trial.***Diabetes*1997,**45:**271–286.Google Scholar - Stevens RJ, Kothari V, Adler AI, Stratton IM, Holman RR, (UKPDS Group):
**The UKPDS risk engine: a model for the risk of coronary heart disease in Type II diabetes (UKPDS 56).***Clinical Science*2001,**101:**671–679. 10.1042/CS20000335Google Scholar - Kothari V, Stevens RJ, Adler AI, Stratton IM, Manley SE, Neil HA, Holman RR:
**Risk of Stroke in type 2 diabetes Estimated by the UKPDS risk engine (UKPDS 60).***Stroke*2002,**33:**776–1781. 10.1161/01.STR.0000020091.07144.C7Google Scholar - Stevens RJ, Coleman RL, Adler AI, Stratton IM, Matthews DR, Holman RR:
**Risk Factors for myocardial Infarction Case Fatality and Stroke Case Fatality in Type 2 Diabete (UKPDS 66).***Diabetes Care*2003,**27**(1):201–207.Google Scholar - Stevens RJ, Stratton IM, Holman RR:
**Modeling glucose exposure as a risk factor for photocoagulation in type 2 diabetes (UKPDS 58).***Journal of Diabetes Complications*2002,**16**(6):371–376. 10.1016/S1056-8727(02)00161-7Google Scholar - Clarke PM, Gray AM, Briggs A, Farmer AJ, Fenn P, Stevens RJ, Matthews DR, Stratton IM, Holman RR, (UKPDS Group):
**A model to estimate the lifetime health outcomes of patients with Type 2 diabetes: The UKPDS Outcomes Model (UKPDS 68).***Diabetologia*2004,**47:**1747–1759. 10.1007/s00125-004-1527-zGoogle Scholar - Klein R, Klein BEK, Moss SE, Davis MD, Demets DL:
**The Wisconsin epidemiologic study of diabetic retinopathy (WESDR IX): Four year incidence and progression of diabetic retinopathy when age at diagnosis is less than 30 years.***Archives Ophtalmologic*1989,**107:**237–243.Google Scholar -
**Framingham Heart Study**[http://www.framingham.com/heart] - Diabetes Prevention Program Research Group:
**Reduction in the incidence of type 2 diabetes with lifestyle intervention or metformin.***New England Journal of Medicine*2002,**356:**393–402. 10.1056/NEJMoa012512Google Scholar - The Health Outcomes Prevention Evaluation Study investigators:
**Effects of angiotension-converting-enzyme inhibitor, ramipril on cardiovascular events in high-risk patients.***New England Journal Medicine*2000,**342:**145–153. 10.1056/NEJM200001203420301Google Scholar - The Health Outcomes Prevention Evaluation Study investigators:
**Effects of ramipril on cardiovascular and microvascular outcomes in people with diabetes mellitus: results of the HOPE study and MICRO HOPE substudy.***Lancet*2000,**355:**253–259. 10.1016/S0140-6736(99)12323-7Google Scholar -
**The Health Plan Employer Data and Information Set (HEDIS)**[http://www.ncqa.org/Programs/HEDIS] -
**Echantillon national temoin représentant des personnes diabétiques**[http://www.invs.sante.fr/publications/entred/index.html] -
**Multiple Risk Factor Intervention Trial**[http://www.gpnotebook.co.uk/cache/-1777991655.htm] -
**Heart Protection Study (HPS)**[http://www.hps.com] -
**Cholesterol and Recurrent Events (CARE)**[http://www.care.com] - Lewis EJ, Hunsicker LG, Clarke WR, Raymond P, Berl T, Rohde R, Raz I:
**The effect of angioitensin-converting-enzyme inhibition on diabetic nephropathy antagonist irbesartan in patients with nephropathy due to type 2 diabetes.***New England Journal of Medicine*1993,**329:**1456–1462. 10.1056/NEJM199311113292004Google Scholar - Irbesartin in Patients with Type-2 Diabetes and Microalbuminuria Study Group:
**The effect of Irbesartin on the development of diabetic nephropathy in patients with type 2 diabetes.***New England Journal of Medicine*2001,**345:**870–878. 10.1056/NEJMoa011489Google Scholar - Lewis EJ, Hunsicker LG, Clarke WR, Tomas P, Pohl MA, Lewis JB, Ritz E, Alkins RC, Rohde R, Raz I:
**Renoprotective effect of the angioitensin-receptor antagonist irbesartan in patients with nephropathy due to type 2 diabetes.***New England Journal of Medicine*2001,**345:**851–860. 10.1056/NEJMoa011303Google Scholar - Colhoun HM, Thomason MJ, Mackeness MI, Moton SM, Betteridge DJ, Durrington PN, Mitman GA, Neil MA, Fuller JH:
**Design of the Collaborative AtoRvastatin Diabetes Study in patients with type 2 diabetes (CARDS).***Diabetes Medicine*2002,**19**(3):201–211. 10.1046/j.1464-5491.2002.00643.xGoogle Scholar - American Diabetes Association:
**Diabetes PHD and Archimedes.**[http://www.archimedesmodel.com/pdf/archimedes-intro_selfguided_press.pdf] - Brandeau ML:
**Modeling complex Medical Decision Problems with Archimedes Model.***Annals of Internal Medicine*2005,**143:**303–304.Google Scholar - Kahn R:
**Dealing with complexity in clinical diabetes: the value of Archimedes.***Diabetes Care*2003,**26**(11):3168–3171.Google Scholar - Schlessinger L, Eddy DM:
**Archimedes: A new model simulating health care systems: the mathematical formulation.***Journal of Biomedical Informatics*2002,**35**(1):37–50. 10.1016/S1532-0464(02)00006-0Google Scholar - Eddy DM, Schlessinger L:
**Archimedes: A trial-validated model of diabetes.***Diabetes Care*2003,**26**(11):3093–3101.Google Scholar - Eddy DM, Schlessinger L:
**Validation of the Archimedes diabetes model.***Diabetes Care*2003,**26**(11):3102–3110.Google Scholar - Herman WH:
**Diabetes Modeling.***Diabetes Care*2003,**26**(11):3182–3183.MathSciNetGoogle Scholar - Bagust A, Hopkinson PK, Maier W, Currie CJ:
**An economic model of the the long-term care burden of type II diabete.***Diabetologia*2001,**44:**2140–2155. 10.1007/s001250100023Google Scholar - Eastman RC, Javitt JC, Herman WH, Dasbach EJ, Zbrozek AS, Dong F, Mannienen D, Garfield SA, Copley-Merriman C, Maier W, Eastman JF, Kotsanos J, Cowie CC, Harris M:
**Model of complications of NIDDM II: Model construction and assumptions.***Diabetes Care*1997,**20:**725–734.Google Scholar - Eastman RC, Javitt JC, Herman WH, Dasbach EJ, Copley-Merriman C, Maier W, Dong F, Mannienen D, Zbrozek AS, Kotsanos J, Garfield SA, Harris M:
**Model of complications of NIDDM. II: Analysis of the health benefits and cost-effectiveness of treating NIDDM with the goal of normoglycaemia.***Diabetes Care*1997,**20:**735–744.Google Scholar - Hernando ME, Garcia A, Perdices FJ, Torralba V, Gomez EJ, del Pozo F:
**Multi-Agent Architecture for the provision of Intelligent Telemedicine Services in Diabetes Management.**[http://www.gbt.tfo.upm.es] - Brown JB, Russel A, Chan W, Pedula K, Aickin M:
**The global diabetes model: user friendly version 3.0.***Diabetes Res Clin Pract*2000,**50**(3):S15-S46. 10.1016/S0168-8227(00)00215-1Google Scholar - Wu H:
**A case study of type 2 diabetes self-management.***Biomedical Engineering Online*2005,**4:**4. 10.1186/1475-925X-4-4Google Scholar - American Diabetes Association Consensus panel:
**Guidelines for Computer Modeling of Diabetes and its Complications. (Consensus Statement).***Diabetes Care*2004,**27**(11):2262–2265.Google Scholar - Bagust A, McEwan P:
**Guidelines for Computer Modeling of Diabetes and its Complications. (Response to American Diabetes Association Consensus Panel).***Diabetes Care*2005,**28**(11):500.Google Scholar - Khan R:
**Guidelines for Computer Modeling of Diabetes and its Complications.(Response to Bagust and McEwan).***Diabetes Care*2005,**28**(11):501.Google Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.