Minimal model S(I)=0 problem in NIDDM subjects: nonzero Bayesian estimates with credible confidence intervals.

The minimal model of glucose kinetics, in conjunction with an insulin-modified intravenous glucose tolerance test, is widely used to estimate insulin sensitivity (S(I)). Parameter estimation usually resorts to nonlinear least squares (NLS), which provides a point estimate, and its precision is expressed as a standard deviation. Applied to type 2 diabetic subjects, NLS implemented in MINMOD software often predicts S(I)=0 (the so-called "zero" S(I) problem), whereas general purpose modeling software systems, e.g., SAAM II, provide a very small S(I) but with a very large uncertainty, which produces unrealistic negative values in the confidence interval. To overcome these difficulties, in this article we resort to Bayesian parameter estimation implemented by a Markov chain Monte Carlo (MCMC) method. This approach provides in each individual the S(I) a posteriori probability density function, from which a point estimate and its confidence interval can be determined. Although NLS results are not acceptable in four out of the ten studied subjects, Bayes estimation implemented by MCMC is always able to determine a nonzero point estimate of S(I) together with a credible confidence interval. This Bayesian approach should prove useful in reanalyzing large databases of epidemiological studies.

[1]  R. DeFronzo,et al.  Pathogenesis of NIDDM: A Balanced Overview , 1992, Diabetes Care.

[2]  R N Bergman,et al.  Assessment of insulin sensitivity in vivo. , 1985, Endocrine reviews.

[3]  R. Bergman,et al.  Ability of alternative indices of insulin sensitivity to predict cardiovascular risk: comparison with the "minimal model". Insulin Resistance Atherosclerosis Study (IRAS) Investigators. , 1998, Annals of epidemiology.

[4]  Eric Walter,et al.  Identification of Parametric Models: from Experimental Data , 1997 .

[5]  C Cobelli,et al.  Minimal model SG overestimation and SI underestimation: improved accuracy by a Bayesian two-compartment model. , 1999, The American journal of physiology.

[6]  A. Festa,et al.  The relation of proinsulin and insulin to insulin sensitivity and acute insulin response in subjects with newly diagnosed Type II diabetes: the Insulin Resistance Atherosclerosis Study , 1999, Diabetologia.

[7]  G. Reaven.  Banting lecture 1988. Role of insulin resistance in human disease. , 1988, Diabetes.

[8]  R. Roman.  The mathematical modeling of metabolic and endocrine systems : E.R. Carson, C. Cobelli. and L. Finkelstein. John Wiley and Sons, New York. 1983. 394 pp., illus., index. , 1986 .

[9]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[10]  M. Rewers,et al.  Insulin sensitivity in subjects with type 2 diabetes. Relationship to cardiovascular risk factors: the Insulin Resistance Atherosclerosis Study. , 1999, Diabetes care.

[11]  D. Owens,et al.  Insulin Secretion and Sensitivity in Newly Diagnosed NIDDM Caucasians in the UK , 1996, Diabetic medicine : a journal of the British Diabetic Association.

[12]  P Vicini,et al.  The hot but not the cold minimal model allows precise assessment of insulin sensitivity in NIDDM subjects. , 1996, The American journal of physiology.

[13]  L. Tierney.  Markov Chains for Exploring Posterior Distributions , 1994 .

[14]  R. DeFronzo,et al.  Glucose clamp technique: a method for quantifying insulin secretion and resistance. , 1979, The American journal of physiology.

[15]  Claudio Cobelli,et al.  Identification of a Minimal Model of Glucose Disappearance for Estimating Insulin Sensitivity , 1979 .

[16]  C Cobelli,et al.  SAAM II: Simulation, Analysis, and Modeling Software for tracer and pharmacokinetic studies. , 1998, Metabolism: clinical and experimental.

[17]  E. Boerwinkle,et al.  The GENNID Study: A resource for mapping the genes that cause NIDDM , 1996, Diabetes Care.

[18]  Y. Z. Ider,et al.  Quantitative estimation of insulin sensitivity. , 1979, The American journal of physiology.

[19]  Peter Green,et al.  Markov chain Monte Carlo in Practice , 1996 .

[20]  G. Reaven.  Banting lecture 1988 , 1997 .

[21]  R N Bergman,et al.  Mapping Genes for NIDDM: Design of the Finland—United States Investigation of NIDDM Genetics (FUSION) Study , 1998, Diabetes Care.

[22]  R N Bergman,et al.  MINMOD: a computer program to calculate insulin sensitivity and pancreatic responsivity from the frequently sampled intravenous glucose tolerance test. , 1986, Computer methods and programs in biomedicine.

[23]  Richard M Watanabe,et al.  A Comparison Between the Minimal Model and the Glucose Clamp in the Assessment of Insulin Sensitivity Across the Spectrum of Glucose Tolerance , 1994, Diabetes.

[24]  Giovanni Sparacino,et al.  Maximum-likelihood versus maximum a posteriori parameter estimation of physiological system models: the c-peptide impulse response case study , 2000, IEEE Transactions on Biomedical Engineering.

[25]  R N Bergman,et al.  Minimal model analysis of intravenous glucose tolerance test-derived insulin sensitivity in diabetic subjects. , 1990, The Journal of clinical endocrinology and metabolism.

[26]  G. Reaven.  Role of Insulin Resistance in Human Disease , 1988, Diabetes.