A critical review of mathematical models and data used in diabetology.  
Jump to Full Text  
MedLine Citation:

PMID: 16808835 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

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 costeffectiveness 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. 
Authors:

A Boutayeb; A Chetouani 
Related Documents
:

19590235  Association of secondary hyperparathyroidism with ckd progression, health care costs an... 17257715  An assessment of the external radiological impact in areas of greece with elevated natu... 17065675  Health care use and costs in the decade after identification of type 1 and type 2 diabe... 17392545  Trends in outpatient prescription drug costs in diabetic patients in germany, 19942004. 9880385  Serum glucose level and diabetes predict tissue plasminogen activatorrelated intracere... 2692795  Ontogeny of the insulin response to glucose in fetal and adult sheep. 
Publication Detail:

Type: Journal Article; Review Date: 20060629 
Journal Detail:

Title: Biomedical engineering online Volume: 5 ISSN: 1475925X ISO Abbreviation: Biomed Eng Online Publication Date: 2006 
Date Detail:

Created Date: 20060825 Completed Date: 20060918 Revised Date: 20130607 
Medline Journal Info:

Nlm Unique ID: 101147518 Medline TA: Biomed Eng Online Country: England 
Other Details:

Languages: eng Pagination: 43 Citation Subset: IM 
Affiliation:

Department of Mathematics Faculty of Sciences, Oujda, Morocco. x.boutayeb@menara.ma 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Algorithms Animals Computer Simulation Databases, Factual Diabetes Mellitus / diagnosis, drug therapy*, metabolism* Diagnosis, ComputerAssisted / methods* Drug Therapy, ComputerAssisted / methods* Glucose / metabolism* Humans Insulin / metabolism* InsulinSecreting Cells / drug effects, metabolism Models, Biological* 
Chemical  
Reg. No./Substance:

0/Insulin; 50997/Glucose 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): Biomed Eng Online ISSN: 1475925X Publisher: BioMed Central, London 
Article Information Download PDF Copyright ? 2006 Boutayeb and Chetouani; licensee BioMed Central Ltd. openaccess: This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Received Day: 8 Month: 5 Year: 2006 Accepted Day: 29 Month: 6 Year: 2006 collection publication date: Year: 2006 Electronic publication date: Day: 29 Month: 6 Year: 2006 Volume: 5First Page: 43 Last Page: 43 ID: 1553453 Publisher Id: 1475925X543 PubMed Id: 16808835 DOI: 10.1186/1475925X543 
A critical review of mathematical models and data used in diabetology  
A Boutayeb1  Email: x.boutayeb@menara.ma 
A Chetouani1  Email: abd.chetouani@menara.ma 
1Department of Mathematics Faculty of Sciences, Oujda, Morocco 
It is now commonly admitted that diabetes is sweeping the globe as a silent epidemic largely contributing to the growing burden of noncommunicable 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: glucoseinsulin dynamics, epidemiology of diabetes and its complications, cost of diabetes and costeffectiveness 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 glucoseinsulin 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 sofar 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 nonexhaustive overview including model structure, computer algorithms, data used as input and model validation.
The majority of mathematical models proposed in the literature were devoted to the dynamics of glucoseinsulin, 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 insulinglucose 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 glucoseinsulin dynamics is thought to begin with the socalled 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 insulinexcitable 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 }= p3/p2, 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}].
However, some authors [^{19},^{36}^{}^{38}] indicated that while the minimal model has minimal number of constants (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 nonobservable 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 glucoseinsulin 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 preinjection 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 glucoseinsulin 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 glucoseinsulin 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 (nx1) state vector, u(t) is the (mx1) control vector, y(t) is the (qx1) 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.
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 devicesmicrosystems 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}].
An example of models with closedloop 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.
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 betacell [^{52}^{}^{54}].
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 nonlinear population models of diabetes mellitus, using ordinary differential equations and numerical implementation.
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 insulindependent 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 IRMA2 trial [^{89}], Irbesartan Diabetic Nephropathy Trial (IDNT) [^{90}], the Collaborative AtoRvastatin Diabetes Study (CARDS) [^{91}].
According to [^{92}], until recently, there have been four main kinds of mathematical models in health care:
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 personbyperson, objectbyobject 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 ). 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 subpopulation of the HOPE Trial (MicroHOPE) [^{82}], the Heart Protection Study (HPS) [^{86}], Cholesterol and Recurrent Events (CARE) [^{87}], the ACE Inhibitors and Diabetic Nephropathy Trial (Lewis) [^{88}], the IRMA2 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 macrovascular 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 and computer algorithms were devoted to the burden, cost and costeffectiveness of diabetes [^{99}^{}^{101}], telemedicine and home management of diabetes [^{102}^{}^{104}]. In a case study paper [^{104}], Wu proposed the following model for selfmanagement 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}].
Mathematical models constitute interesting tools for the understanding of diseases. They provide insights, improve intuitions, clarify assumptions for formal theory, allow for planning studies, estimating parameters, determining sensitivities, assessing conjectures, simulating simple and complex phenomena and providing future predictions. In the case of diabetes, simple and comprehensive models dealing with different aspects of the disease, have been used during the last three decades. In general, simple models are so simple as to be inadequate but they have the advantage of using a small number of identifiable parameters. Comprehensive models on the other hand are models which try to represent the system (biological, clinical, economic, etc...) by taking into account all interactions. This makes them very complex and generally not identifiable. In the present paper, our main objective was an overview of models and studies dealing with different aspects of diagnosis, care and management of diabetes and its complications. We presented a non exhaustive list of published models with their theoretical and applied aspects, indicating what were the hypotheses that lead authors to propose new, modified, generalized or alternative models. But we must say that we did not intend to compare all models or classify them according to whatever criterion of performance. Even when available models are simple, they are not necessary comparable. For instance, Bolie's model [^{24}] is one of the simplest models proposed to estimate the linear glucose disappearance and insulinglucose dynamics. However, since the author is one of the pioneers in this field, his work remains an unavoidable reference. Another very simple model was recently proposed in a case study by Wu [^{104}] but the purpose was a selfmanagement of type 2 diabetes. The model being a case study based on a single type 2 diabetes person, the results yielded need to be considered with caution. The discrete matrix model considered by Boutayeb and Derouich [^{59}] is also simple but the approach is completely different since it deals with the control of complications Bergman et al. [^{7}] discussed seven models before selecting the " best one" which became from then the well known minimal model (between simple inadequate and comprehensive not identifiable models), they based their selection on the following criteria:
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 socalled "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.
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 nonexhaustive panorama of the papers which have used mathematical modeling for different aspects of diabetes, including glucoseinsulin dynamics, betacell 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.
The paper is dedicated to the students : Wiam, Sara, Rajae, Oumnia, Mariam, Hanae and FatimaZohra 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.
References
The world health report 2003Today's challenges. Geneva, World Health Organization;  
International Diabetes FederationIFD report 2003  
Boutayeb A,Boutayeb S. The burden of non communicable diseases in developing countries 2005;4International Journal of Equity in Health; :2. [pmid: 15651987]  
Ratner RE. Type 2 diabetes mellitus: the grand overviewDiabet Med 1998;14:54–57.  
Jiwa F. Diabetes in the 1990s an overviewStat Bull Metrop Co 1997;78:2–8.  
The Expert Committee on the Diagnosis and Classification of Diabetes MellitusReport of the Expert Committee on the diagnosis and Classification of Diabetes MellitusDiabetes Care 1997;20:1183–1197. [pmid: 9203460]  
Bergman RN,Finegood DT,Ader M. Assessment of Insulin Sensitivity in VivoEndicrine Reviews 1985;6:45–86.  
Bergman RN. Bergman RN, Lovejoy JCThe minimal model: yesterday, today and tomorrowThe minmal model Approach and Determination of Glucose Tolerance 1997;7Boston : Louisiana state university Press; :3–50.  
Bergman RN. Novotny, Green, BostonThe minimal model of glucose regulation: a biographyMathematical Modeling in Nutrition and Health. 2001Kluwer Academic/Plenum;  
Nucci G,Cobelli C. Models of subcatuneous insulin kinetics: a critical reviewComputer Methods and Programs in Biomedicine 2000;62:249–257. [pmid: 10837910] [doi: 10.1016/S01692607(00)000717]  
Sundell J,Knuuti J. Insulin and myocardial blood flowCardiovascular Research 2003;57:312–319. [pmid: 12566104] [doi: 10.1016/S00086363(02)007186]  
Finegood DT. Bergman RN, Lovejoy JCApplication of the minimal model of glucose kineticsThe minmal model Approach and Determination of Glucose Tolerance 1997;7Boston : Louisiana state university Press; :51–122.  
Mari A. Mathematical modelling in glucose metabolism and insulin secretionCurrent Opinion Clinical Nutrition Metabolism Care 2002;5:495–501. [doi: 10.1097/0007519720020900000007]  
Lehmann ED,Deutsch T. Application of computers in diabetes care: a review. II, Computers for decision support and educationMed Inform 1995;20:303–329.  
Bellazzi R,Nucci G,Cobelli C. The Subcutaneous Route to Insulin Dependent Diabetes Therapy: ClosedLoop and Partially ClosedLoop Control Strategies for insulin Delivery and Measuring Glucose ConcentrationIEEE Engrg Medicine Biol 2001;20:54–64. [doi: 10.1109/51.897828]  
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 PatientsIEEE Engineering in Medicine and Biologie 2001;20:65–73. [doi: 10.1109/51.897829]  
Koschinsky T,Heinemann . Sensors for glucose monitoring: technical and clinical aspectsDiabetes/Metabolism Research and Reviews 2001;17:113–123. [pmid: 11307176] [doi: 10.1002/dmrr.188]  
Palerm CCR. Drug Infusion Control: An Extended Direct Model Reference Adaptive Control StrategyPhD Thesis 2003;chap 3?4Rensselear Polytechnic Institute, Troy, New York; :45–60.  
Makroglou A,Li J,Kuang Y. Mathematical models and software tools for the glucoseinsulin regulatory system and diabetes: an overviewApplied Numerical Mathematics 2006;56:559–573. [doi: 10.1016/j.apnum.2005.04.023]  
Kalergis M,De Grandpre E,Andersons C. The Role of the Glycemic Index in the Prevention and Management of Diabetes: A Review and DiscussionCanadian Journal of Diabetes 2005;29:27–38.  
Atun R,GurolUrganci I. Health expenditure: an 'investment' rather than a cost ?. 2005International Economics Programme; IEP XP 05/01.  
Boutayeb A. The double burden of communicable and non communicable diseases in developing countries : A reviewTransactions of the Royal Society of Tropical Medicine and Hygiene 2006;100:191–199. [pmid: 16274715] [doi: 10.1016/j.trstmh.2005.07.021]  
Himsworth HP,Ker RB. Insulinsensitive and insulin insensitive types of diabetes millitusCli Sci 1939;4:119–122.  
Bolie VW. Coefficients of normal blood glucose regulationJ Appl Physiol 1961;16:783–788. [pmid: 13870789]  
Ackerman E,Gatewood LC,Rosevear JW,Molnar GD. Model studies of blood glucose regulationBull Math Biophys 1965;27:21–24. [pmid: 5884133]  
Della C,Romano MR,Voehhelin MR,Seriam E. On a mathematical model for the analysis of the glucose tolerance curveDiabetes 1970;19:145–148. [pmid: 4908764]  
Serge G,Turcogl M,Varcellone G. Modelling blood glucose and insulin kinetics in normal diabetic and obese subjectsDiabetes 1973;22:94–97. [pmid: 4687648]  
Srinivasan R,Kadish AH,Sridhar R. A mathematical model for the control mechanism of freefatty acid and glucose metabolism in normal humansComp Biomed Res 1970;3:146–149. [doi: 10.1016/00104809(70)900212]  
Bergman RN,Ider YZ,Bowden CR,Cobelli C. Quantitative Estimation of Insulin SensitivityAm J Physiol 1979;23:E667–E677. [pmid: 443421]  
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 DogDiabetes 1980;29:979–990. [pmid: 7002673]  
Bergman RN,Bowden CR,Cobelli C. Cobelli, BergmanThe Minimal Model approach to quantification of factors controlling glucose disposal in manCarbohydrate Metabolism 1981;chap 13John Wiley & Sons Ltd; :269–293.  
Cobelli C,Mari A. Validation of mathematical models complex endocrinemetabomism systems. A case study on a model of glucose regulationMed & Biot Eng & Comput 1983;21:390–399.  
Cobelli C,Nucci G,del Prato S. A physiological simulation model of the glucoseinsulin system in type I diabetesDiabetes, Nutrition & Metabolism 1998;11:78–80.  
Derouich M,Boutayeb A. The effect of physical exercise on the dynamics of glucose and insulinJournal of Biomechanics 2002;35:911–917. [pmid: 12052393] [doi: 10.1016/S00219290(02)000556]  
Bergman RN. Pathogenesis and prediction of diabetes mellitus: lessons from integrative physiologyMount Sinai J Medicine 2002;60:280–290.  
De Gaetano A,Arino O. Mathematical Modelling of the Intravenous Glucose Tolerance TestJ Math Biol 2000;40:136–168. [pmid: 10743599] [doi: 10.1007/s002850050007]  
De Gaetano A,Arino O. A statistical approach to the determination of stability for dynamical systems modelling physiological processesMath Comput Modelling 2000;31:41–51. [doi: 10.1016/S08957177(00)000200]  
Li J,Kuang Y,Li B. Analysis of IVGTT GlucoseInsulin Interaction Models with time delayDiscrete and Continous Dynamical Systems Series B 2000;1:103–124.  
Mukhopadhyay A,De Gaetano A,Arino O. Modelling the intravenous glucose tolerance test: A global study for singledistributeddelay modelDiscrete and Continous Dynamical Systems Series B 2004;4:407–417.  
Cobelli C,Thomaseth K. Optimal input design for identification of compartimental models : theory and applications to a model of glucose kineticsMathematical Biosciences 1985;77:267–270. [doi: 10.1016/00255564(85)901014]  
Cobelli C,Thomaseth K. The minimal model of glucose disappearence : optimal input studiesMath Biosciences 1987;83:127–130. [doi: 10.1016/00255564(87)901076]  
Cobelli C,Thomaseth K. An optimalality of the impulse input for linear system identificationMath Biosciences 1988;89:127–129. [doi: 10.1016/00255564(88)900946]  
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 examplesMath Biosciences 89:135–137. [doi: 10.1016/00255564(88)900958]  
Lam ZH,Hwang KS,Lee JY,Chase JG,Walker GC. Active insulin infusion using optimal and derivative weighted controlMedical engineering physics 2002;24:663–672. [pmid: 12460725] [doi: 10.1016/S13504533(02)001479]  
Salam J,Charles MA. Devices for insulin administrationDiabetes care 1999;13:955–979.  
Lehman ED,Deutsch T. A physiological model of glucose insulin interaction in type 1 diabetes mellitusJ Biomed Eng 1992;14:235–242. [pmid: 1588781]  
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 1993;18Med Inform London; :83–101.  
Shimoda S,Nishida K,Sakakida M,Konno Y,ishinose K,Uehara M,Nowak T,Shichiri M. Closedloop subcutaneous isulin infusion algorithm with a short acting insulin analog for longterm clinical application of a wearable artificial endocrine pancreasFrontiers Med Biol Eng 1997;8:197–211.  
Hashiguchi Y,Sakakida M,Nishida K,Uemura T,Kajiwara K,Shichiri M. Development of a miniaturized glucose monitoring system by combining a needletype glucose sensor with microdialysis sampling method. Longterm subcutaneous tissue glucose monitoring in ambulatory diabetic patientsdiabetes care 1994;17:387–396. [pmid: 8062605]  
Reaven GM,Doberne L,Greenfield MS. Comparaison of insulin secretion and in vivo insulin action in nonobese and moderately obese individuals with noninsulin dependent diabetes mellitusDiabetes 1982;31:382–388. [pmid: 6759253]  
Andres R,Swerdloff R,Pozefsky T,Coleman D. Skeggs Jt LTManual feedback technique for the control of blood glucose concentrationAutomation in Analytical Chemistry 1966Mediad Inc, New York; :486–490.  
Bagust A,Beale S. Deteriorating betacell function in type 2 diabetes: a longterm modelQ J Med 2003;96:281–288.  
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 manDiabetologia 1985;28:412–419. [pmid: 3899825] [doi: 10.1007/BF00280883]  
Levy JC,Matthews DR,Hermans MP. Correct Homeostasis Model Assessment(HOMA) evaluation uses the computer program(letter)Diabetes Care 1998;20:2191–2192. [pmid: 9839117]  
Bernoulli D. Acad?mie Royale des SciencesEssai d'une nouvelle analyse de la mortalit? caus?e par la petite v?role et des avantages de 1'inoculation pour la pr?venirm?moire de math?matiques et de physique 1760Paris; :1–45.  
Hethcote HW. LevinA thousand and one epidemic modelsFrontiers in Theoretical Biology 1994100, SpringerVerlag, Berlin; :504–515. Lecture notes in Biomath  
Hethcote HW. The Mathematics of Infectious DiseasesSIAM REVIEW 2000;42:599–653. [doi: 10.1137/S0036144500371907]  
Boutayeb A,Kerfati A. Mathematical models in diabetology modellingMeasurement and Control, C, AMSE 1994;44:53–63.  
Boutayeb A,Derouich M. Age structured models for diabetes in East MoroccoMathematics and Copmuter Simulation 2002;58:215–229. [doi: 10.1016/S03784754(01)003688]  
Boutayeb A,Twizell EH. An age structured model for complications of diabetes mellitus in MoroccoSimulation Modelling Practice and Theory 2004;12:77–87. [doi: 10.1016/j.simpat.2003.11.003]  
Boutayeb A,Chetouani A. Dynamics of a disabled population in MoroccoBiomedical Engineering Online 2003;2:2. [pmid: 12625838] [doi: 10.1186/1475925X22]  
Boutayeb A,Twizell EH,Achouyab K,Chetouani A. A mathematical model for the burden of diabetes and its complicationsBiomedical Engineering Online 2004;3:20. [pmid: 15222886] [doi: 10.1186/1475925X320]  
Boutayeb A,Chetouani A,Achouyab K,Twizell EH. A nonlinear population model of diabetes mellitusJournal of Applied Mathematics and computing 2006;21:127–139.  
James RC,Blanchard JF,Campbell D,Clottey C,Osei W,Svenson LW,Noseworthy TW. A model for noncommunicable disease surveillance in Canada: the prairie pilot diabetes surveillance systemChronic Dis Can 2004;25:7–12. [pmid: 15298483]  
Paterson KR. Diabetes epidemiology: present and futureDiabet Med 1999;16:446–7. [pmid: 10391390] [doi: 10.1046/j.14645491.1999.00108.x]  
West KM. Epidemiology of diabetes and its macrovascular complicationsDiabetes Care 1979;2:63–4. [pmid: 520131]  
DCCT Research GroupThe effect of intensive treatment of diabetes on the development and progression of longterm complications in insulin dependent diabetes mellitusNew England Journal of Medicine 1993;329:977–986. [pmid: 8366922] [doi: 10.1056/NEJM199309303291401]  
O'Neill S. The UK Prospective Diabetes Study(UKPDS): Questions & Answers  
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 TrialDiabetes Care 2000;25:275–278. [pmid: 11815495]  
Diabetes Control and Complications Trials Research GroupThe relationship of glycemic expsure (HbA_{1c}) to the risk of development and progression of retinopathy in the diabetes control complications trial in insulindependent diabetes mellitusDiabetes 1995;44:968–983. [pmid: 7622004]  
Diabetes Control and Complications Trials Research GroupThe absence of a glycemic threshold for the development of longterm complication: the perspective of the diabetes control complications trial in insulindependent diabetes mellitusDiabetes 1996;45:1289–1298. [pmid: 8826962]  
Diabetes Control and Complications Trials Research GroupHypoglycemia in the Diabetes Control Complications TrialDiabetes 1997;45:271–286.  
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. [pmid: 11724655] [doi: 10.1042/CS20000335]  
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. [doi: 10.1161/01.STR.0000020091.07144.C7]  
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:201–207. [pmid: 14693990]  
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:371–376. [doi: 10.1016/S10568727(02)001617]  
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. [pmid: 15517152] [doi: 10.1007/s001250041527z]  
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 yearsArchives Ophtalmologic 1989;107:237–243.  
Framingham Heart Study  
Diabetes Prevention Program Research GroupReduction in the incidence of type 2 diabetes with lifestyle intervention or metforminNew England Journal of Medicine 2002;356:393–402. [doi: 10.1056/NEJMoa012512]  
The Health Outcomes Prevention Evaluation Study investigatorsEffects of angiotensionconvertingenzyme inhibitor, ramipril on cardiovascular events in highrisk patientsNew England Journal Medicine 2000;342:145–153. [doi: 10.1056/NEJM200001203420301]  
The Health Outcomes Prevention Evaluation Study investigatorsEffects of ramipril on cardiovascular and microvascular outcomes in people with diabetes mellitus: results of the HOPE study and MICRO HOPE substudyLancet 2000;355:253–259. [pmid: 10675071] [doi: 10.1016/S01406736(99)123237]  
The Health Plan Employer Data and Information Set (HEDIS)  
Echantillon national temoin repr?sentant des personnes diab?tiques  
Multiple Risk Factor Intervention Trial  
Heart Protection Study (HPS)  
Cholesterol and Recurrent Events (CARE)  
Lewis EJ,Hunsicker LG,Clarke WR,Raymond P,Berl T,Rohde R,Raz I. The effect of angioitensinconvertingenzyme inhibition on diabetic nephropathy antagonist irbesartan in patients with nephropathy due to type 2 diabetesNew England Journal of Medicine 1993;329:1456–1462. [pmid: 8413456] [doi: 10.1056/NEJM199311113292004]  
Irbesartin in Patients with Type2 Diabetes and Microalbuminuria Study GroupThe effect of Irbesartin on the development of diabetic nephropathy in patients with type 2 diabetesNew England Journal of Medicine 2001;345:870–878. [pmid: 11565519] [doi: 10.1056/NEJMoa011489]  
Lewis EJ,Hunsicker LG,Clarke WR,Tomas P,Pohl MA,Lewis JB,Ritz E,Alkins RC,Rohde R,Raz I. Renoprotective effect of the angioitensinreceptor antagonist irbesartan in patients with nephropathy due to type 2 diabetesNew England Journal of Medicine 2001;345:851–860. [pmid: 11565517] [doi: 10.1056/NEJMoa011303]  
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:201–211. [doi: 10.1046/j.14645491.2002.00643.x]  
American Diabetes AssociationDiabetes PHD and Archimedes  
Brandeau ML. Modeling complex Medical Decision Problems with Archimedes ModelAnnals of Internal Medicine 2005;143:303–304. [pmid: 16103475]  
Kahn R. Dealing with complexity in clinical diabetes: the value of ArchimedesDiabetes Care 2003;26:3168–3171. [pmid: 14578256]  
Schlessinger L,Eddy DM. Archimedes: A new model simulating health care systems: the mathematical formulationJournal of Biomedical Informatics 2002;35:37–50. [pmid: 12415725] [doi: 10.1016/S15320464(02)000060]  
Eddy DM,Schlessinger L. Archimedes: A trialvalidated model of diabetesDiabetes Care 2003;26:3093–3101. [pmid: 14578245]  
Eddy DM,Schlessinger L. Validation of the Archimedes diabetes modelDiabetes Care 2003;26:3102–3110. [pmid: 14578246]  
Herman WH. Diabetes ModelingDiabetes Care 2003;26:3182–3183. [pmid: 14578260]  
Bagust A,Hopkinson PK,Maier W,Currie CJ. An economic model of the the longterm care burden of type II diabeteDiabetologia 2001;44:2140–2155. [pmid: 11793015] [doi: 10.1007/s001250100023]  
Eastman RC,Javitt JC,Herman WH,Dasbach EJ,Zbrozek AS,Dong F,Mannienen D,Garfield SA,CopleyMerriman C,Maier W,Eastman JF,Kotsanos J,Cowie CC,Harris M. Model of complications of NIDDM II: Model construction and assumptionsDiabetes Care 1997;20:725–734. [pmid: 9135934]  
Eastman RC,Javitt JC,Herman WH,Dasbach EJ,CopleyMerriman 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 costeffectiveness of treating NIDDM with the goal of normoglycaemiaDiabetes Care 1997;20:735–744. [pmid: 9135935]  
Hernando ME,Garcia A,Perdices FJ,Torralba V,Gomez EJ,del Pozo F. MultiAgent Architecture for the provision of Intelligent Telemedicine Services in Diabetes Management  
Brown JB,Russel A,Chan W,Pedula K,Aickin M. The global diabetes model: user friendly version 3.0Diabetes Res Clin Pract 2000;50:S15–S46. [pmid: 11080561] [doi: 10.1016/S01688227(00)002151]  
Wu H. A case study of type 2 diabetes selfmanagementBiomedical Engineering Online 2005;4:4. [pmid: 15644138] [doi: 10.1186/1475925X44]  
American Diabetes Association Consensus panelGuidelines for Computer Modeling of Diabetes and its Complications. (Consensus Statement)Diabetes Care 2004;27:2262–2265. [pmid: 15333499]  
Bagust A,McEwan P. Guidelines for Computer Modeling of Diabetes and its Complications. (Response to American Diabetes Association Consensus Panel)Diabetes Care 2005;28:500. [pmid: 15677831]  
Khan R. Guidelines for Computer Modeling of Diabetes and its Complications.(Response to Bagust and McEwan)Diabetes Care 2005;28:501. 
Article Categories:

Previous Document: Pvalue based visualization of codon usage data.
Next Document: In vitro hydrodynamic properties of the Miethke ProGAV hydrocephalus shunt.