# DDMODEL00000227: Glucose kinetics in humans: a unique model for different tests

Short description:
This is a model for glucose kinetics in humans, including a physiology-based mechanism for glucose uptake saturation. It is able to describe the relationship between glucose utilization and glucose and insulin concentrations during different tests, such as the clamp, the oral glucose tolerance test and the mixed-meal test. By Roberto Bizzotto, Institute of Neuroscience, CNR, Padua, Italy. Copyright 2016
 Format: PharmML 0.8.x (0.8.1) Related Publication: .hiddenContent {display:none;} Glucose uptake saturation explains glucose kinetics profiles measured by different tests. Bizzotto R, Natali A, Gastaldelli A, Muscelli E, Krssak M, Brehm A, Roden M, Ferrannini E, Mari A American journal of physiology. Endocrinology and metabolism, 8/2016, Volume 311, Issue 2, pages: E346-57 Affiliation: CNR Institute of Neuroscience, Padua, Italy; roberto.bizzotto@isib.cnr.it. Abstract: It is known that for a given insulin level glucose clearance depends on glucose concentration. However, a quantitative representation of the concomitant effects of hyperinsulinemia and hyperglycemia on glucose clearance, necessary to describe heterogeneous tests such as euglycemic and hyperglycemic clamps and oral tests, is lacking. Data from five studies (123 subjects) using a glucose tracer and including all the above tests in normal and diabetic subjects were collected. A mathematical model was developed in which glucose utilization was represented as a Michaelis-Menten function of glucose with constant Km and insulin-controlled Vmax, consistently with the basic notions of glucose transport. Individual values for the model parameters were estimated using a population approach. Tracer data were accurately fitted in all tests. The estimated Km was 3.88 (2.83-5.32) mmol/l [median (interquartile range)]. Median model-derived glucose clearance at 600 pmol/l insulin was reduced from 246 to 158 ml·min(-1)·m(-2) when glucose was raised from 5 to 10 mmol/l. The model reproduced the characteristic lack of increase in glucose clearance when moderate hyperinsulinemia was accompanied by hyperglycemia. In all tests, insulin sensitivity was inversely correlated with BMI, as expected (R(2) = 0.234, P = 0.0001). In conclusion, glucose clearance in euglycemic and hyperglycemic clamps and oral tests can be described with a unifying model, consistent with the notions of glucose transport and able to reproduce the suppression of glucose clearance due to hyperglycemia observed in previous studies. The model may be important for the design of reliable glucose homeostasis simulators. Contributors: Roberto Bizzotto
 Context of model development: Mechanistic Understanding; Discrepancy between implemented model and original publication: The model in MDL/PharmML language differs from the one used in the related publication in two aspects: 1) it does not consider that the subject undergoing a paired test is unique; 2) it is used as a simulation model, as performing parameter estimation on the provided dataset ("Simulated_glucoseKinetics.csv") would be meaningless, because of the used model inputs.; Long technical model description: Long_technical_model_description_glucoseKinetics.txt; Model compliance with original publication: No; Model implementation requiring submitter’s additional knowledge: No; Modelling context description: This model was developed to be able to describe glucose kinetics during different tests, such as the clamp, the oral glucose tolerance test and the mixed-meal test, with a unique set of equations and a unique parameter distribution. The success of the development process allowed to state that the physiological mechanisms underlying glucose kinetics are similar after oral and intravenous glucose administration. Moreover, the developed model quantitatively explained, for the first time, an important phenomenon: in the presence of hyperglycaemia, hyperinsulinemia fails to increase glucose clearance as would expected if glucose clearance were unaffected by the glucose levels. This model may have future applications in areas such as the development of more accurate insulin sensitivity indices from tests other than the glucose clamp and for the study of insulin sensitivity in circumstances in which the glucose clamp would not be feasible or appropriate. Furthermore, the model may be an essential component of new accurate glucose homeostasis simulators.; Modelling task in scope: estimation; simulation; Nature of research: Clinical research & Therapeutic use; Therapeutic/disease area: Endocrinology;
 Validation Status: Annotations are correct. Certification Comment: This model is not certified.
• Model owner: Roberto Bizzotto
• Submitted: Oct 14, 2016 5:54:08 PM
##### Revisions
• Version: 15
• Submitted on: Jun 30, 2017 8:02:50 AM
• Submitted by: Roberto Bizzotto
• With comment: Edited model metadata online.
• Version: 13
• Submitted on: Oct 14, 2016 6:12:56 PM
• Submitted by: Roberto Bizzotto
• With comment: Edited model metadata online.
• Version: 11
• Submitted on: Oct 14, 2016 6:05:42 PM
• Submitted by: Roberto Bizzotto
• With comment: Edited model metadata online.
• Version: 9
• Submitted on: Oct 14, 2016 5:54:08 PM
• Submitted by: Roberto Bizzotto
• With comment: Edited model metadata online.

### Name

Glucose kinetics in humans: a unique model for different tests

### Description

This is a model for glucose kinetics in humans, including a physiology-based mechanism for glucose uptake saturation. It is able to describe the relationship between glucose utilization and glucose and insulin concentrations during different tests, such as the clamp, the oral glucose tolerance test and the mixed-meal test. By Roberto Bizzotto, Institute of Neuroscience, CNR, Padua, Italy. Copyright 2016

 T

### Function Definitions

 $\mathrm{additiveError}:\mathrm{real}\left(\mathrm{additive}:\mathrm{real}\right)=\mathrm{additive}$

### Parameter Model: $\mathrm{pm}$

#### Random Variables

${\mathrm{eta_KmG}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_KmG}\right)$
${\mathrm{eta_Vmax0}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_Vmax0}\right)$
${\mathrm{eta_Emax}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_Emax}\right)$
${\mathrm{eta_gamma}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_gamma}\right)$
${\mathrm{eta_KmI}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_KmI}\right)$
${\mathrm{eta_t12I}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_t12I}\right)$
${\mathrm{eta_t12G}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_t12G}\right)$
${\mathrm{eta_V}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_V}\right)$
${\mathrm{eta_flambda3}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_flambda3}\right)$
${\mathrm{eta_flambda2}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_flambda2}\right)$
${\mathrm{eta_w1}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_w1}\right)$
${\mathrm{eta_fw2}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_fw2}\right)$
${\mathrm{eta_F}}_{\mathrm{vm_mdl.ID}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.var_F}\right)$
${\mathrm{epsilon}}_{\mathrm{vm_err.DV}}~\mathrm{Normal2}\left(\mathrm{mean}=0,\mathrm{var}=\mathrm{pm.sigma}\right)$

#### Population Parameters

$\mathrm{typ_KmG}$
$\mathrm{typ_Vmax0}$
$\mathrm{typ_Emax}$
$\mathrm{typ_gamma}$
$\mathrm{typ_KmI}$
$\mathrm{typ_t12I}$
$\mathrm{typ_t12G}$
$\mathrm{typ_V}$
$\mathrm{typ_flambda3}$
$\mathrm{typ_flambda2}$
$\mathrm{typ_w1}$
$\mathrm{typ_fw2}$
$\mathrm{typ_F}$
$\mathrm{var_KmG}$
$\mathrm{var_Vmax0}$
$\mathrm{var_Emax}$
$\mathrm{var_gamma}$
$\mathrm{var_KmI}$
$\mathrm{var_t12I}$
$\mathrm{var_t12G}$
$\mathrm{var_V}$
$\mathrm{var_flambda3}$
$\mathrm{var_flambda2}$
$\mathrm{var_w1}$
$\mathrm{var_fw2}$
$\mathrm{var_F}$
$\mathrm{corr_gamma_KmI}$
$\mathrm{alpha}$
$\mathrm{sigma}$

#### Individual Parameters

$\mathrm{ln}\left(\mathrm{KmG}\right)=\mathrm{ln}\left(\mathrm{pm.typ_KmG}\right)+\mathrm{pm.eta_KmG}$
$\mathrm{ln}\left(\mathrm{Vmax0}\right)=\mathrm{ln}\left(\mathrm{pm.typ_Vmax0}\right)+\mathrm{pm.eta_Vmax0}$
$\mathrm{ln}\left(\mathrm{Emax}\right)=\mathrm{ln}\left(\mathrm{pm.typ_Emax}\right)+\mathrm{pm.eta_Emax}$
$\mathrm{ln}\left(\mathrm{gamma}\right)=\mathrm{ln}\left(\mathrm{pm.typ_gamma}\right)+\mathrm{pm.eta_gamma}$
$\mathrm{ln}\left(\mathrm{KmI}\right)=\mathrm{ln}\left(\mathrm{pm.typ_KmI}\right)+\mathrm{pm.eta_KmI}$
$\mathrm{ln}\left(\mathrm{t12I}\right)=\mathrm{ln}\left(\mathrm{pm.typ_t12I}\right)+\mathrm{pm.eta_t12I}$
$\mathrm{ln}\left(\mathrm{t12G}\right)=\mathrm{ln}\left(\mathrm{pm.typ_t12G}\right)+\mathrm{pm.eta_t12G}$
$\mathrm{ln}\left(V\right)=\mathrm{ln}\left(\mathrm{pm.typ_V}\right)+\mathrm{pm.eta_V}$
$\mathrm{logit}\left(\mathrm{flambda3}\right)=\mathrm{logit}\left(\mathrm{pm.typ_flambda3}\right)+\mathrm{pm.eta_flambda3}$
$\mathrm{logit}\left(\mathrm{flambda2}\right)=\mathrm{logit}\left(\mathrm{pm.typ_flambda2}\right)+\mathrm{pm.eta_flambda2}$
$\mathrm{logit}\left(\mathrm{w1}\right)=\mathrm{logit}\left(\mathrm{pm.typ_w1}\right)+\mathrm{pm.eta_w1}$
$\mathrm{logit}\left(\mathrm{fw2}\right)=\mathrm{logit}\left(\mathrm{pm.typ_fw2}\right)+\mathrm{pm.eta_fw2}$
$\mathrm{ln}\left(F\right)=\mathrm{ln}\left(\mathrm{pm.typ_F}\right)+\mathrm{pm.eta_F}$

#### Random Variable Correlation

$\mathrm{corr}\left(\mathrm{eta_gamma},\mathrm{eta_KmI}\right)=\mathrm{pm.corr_gamma_KmI}$

### Structural Model: $\mathrm{sm}$

#### Variables

$\mathrm{INS}$
$\mathrm{GLU}$
$\mathrm{TOBS}$
$\mathrm{T1}$
$\mathrm{GLU1}$
$\mathrm{INS1}$
$\mathrm{VHL}=700$
$\mathrm{deltaHL}=15$
$\mathrm{delta}=10$
$\mathrm{w2}=\left(1-\mathrm{pm.w1}\right)\cdot \mathrm{pm.fw2}$
$\mathrm{w3}=1-\mathrm{pm.w1}-\mathrm{sm.w2}$
$\mathrm{lambda1}=\frac{\frac{\left(\mathrm{pm.w1}\cdot \mathrm{pm.flambda2}\cdot \mathrm{pm.flambda3}+\mathrm{sm.w2}\cdot \mathrm{pm.flambda3}+\mathrm{sm.w3}\right)}{\mathrm{pm.flambda2}\cdot \mathrm{pm.flambda3}}\cdot \mathrm{sm.delta}\cdot \mathrm{pm.F}}{\left(\mathrm{sm.delta}\cdot \left(\mathrm{pm.V}-\mathrm{sm.VHL}\right)-\mathrm{pm.F}\right)}$
$\mathrm{lambda2}=\mathrm{sm.lambda1}\cdot \mathrm{pm.flambda2}$
$\mathrm{lambda3}=\mathrm{sm.lambda2}\cdot \mathrm{pm.flambda3}$
$\mathrm{c1}=\frac{\mathrm{sm.deltaHL}\cdot \mathrm{pm.F}}{\left(\mathrm{sm.deltaHL}\cdot \mathrm{sm.VHL}-2\cdot \mathrm{pm.F}\right)}$
$\mathrm{c2}=\frac{-\mathrm{sm.deltaHL}\cdot \mathrm{pm.F}}{\left(\mathrm{sm.deltaHL}\cdot \mathrm{sm.VHL}-\mathrm{pm.F}\right)}$
$I=\frac{\left(T-\mathrm{sm.T1}\right)}{\left(\mathrm{sm.TOBS}-\mathrm{sm.T1}\right)}\cdot \left(\mathrm{sm.INS}-\mathrm{sm.INS1}\right)+\mathrm{sm.INS1}$
$\mathrm{GL}=\frac{\left(T-\mathrm{sm.T1}\right)}{\left(\mathrm{sm.TOBS}-\mathrm{sm.T1}\right)}\cdot \left(\mathrm{sm.GLU}-\mathrm{sm.GLU1}\right)+\mathrm{sm.GLU1}$
$\mathrm{t0}=0$
$\begin{array}{c}\frac{d}{dT}\mathrm{X1}=\frac{\left(\mathrm{sm.GL}-\mathrm{sm.X1}\right)\cdot \mathrm{ln}\left(2\right)}{\mathrm{pm.t12G}}\\ \mathrm{X1}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}X=\frac{\left(\mathrm{sm.X1}-\mathrm{sm.X}\right)\cdot \mathrm{ln}\left(2\right)}{\mathrm{pm.t12G}}\\ X\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{Z1}=\frac{\left(\mathrm{sm.I}-\mathrm{sm.Z1}\right)\cdot \mathrm{ln}\left(2\right)}{\mathrm{pm.t12I}}\\ \mathrm{Z1}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}Z=\frac{\left(\mathrm{sm.Z1}-\mathrm{sm.Z}\right)\cdot \mathrm{ln}\left(2\right)}{\mathrm{pm.t12I}}\\ Z\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\mathrm{Vmax}=\mathrm{pm.Vmax0}+\frac{\mathrm{pm.Emax}\cdot {\mathrm{sm.Z}}^{\mathrm{pm.gamma}}}{\left({\mathrm{pm.KmI}}^{\mathrm{pm.gamma}}+{\mathrm{sm.Z}}^{\mathrm{pm.gamma}}\right)}$
$\mathrm{cl}=\frac{\mathrm{sm.Vmax}}{\left(\mathrm{pm.KmG}+\mathrm{sm.X}\right)}$
$E=\frac{\mathrm{sm.cl}}{\mathrm{pm.F}}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xHL1}=\mathrm{sm.c2}\cdot \mathrm{sm.xHL1}+\mathrm{sm.Gv}\\ \mathrm{xHL1}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xHL2}=-\mathrm{sm.deltaHL}\cdot \mathrm{sm.xHL2}+\mathrm{sm.Gv}\\ \mathrm{xHL2}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$G=\mathrm{sm.c1}\cdot \mathrm{sm.xHL1}-\mathrm{sm.c1}\cdot \mathrm{sm.xHL2}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xPER1}=-\mathrm{sm.lambda1}\cdot \mathrm{sm.xPER1}+\mathrm{pm.w1}\cdot \left(1-\mathrm{sm.E}\right)\cdot \mathrm{sm.G}\\ \mathrm{xPER1}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xPER2}=-\mathrm{sm.lambda2}\cdot \mathrm{sm.xPER2}+\mathrm{sm.w2}\cdot \left(1-\mathrm{sm.E}\right)\cdot \mathrm{sm.G}\\ \mathrm{xPER2}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xPER3}=-\mathrm{sm.lambda3}\cdot \mathrm{sm.xPER3}+\mathrm{sm.w3}\cdot \left(1-\mathrm{sm.E}\right)\cdot \mathrm{sm.G}\\ \mathrm{xPER3}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\begin{array}{c}\frac{d}{dT}\mathrm{xPER4}=\mathrm{sm.lambda1}\cdot \mathrm{sm.xPER1}+\mathrm{sm.lambda2}\cdot \mathrm{sm.xPER2}+\mathrm{sm.lambda3}\cdot \mathrm{sm.xPER3}-\mathrm{sm.delta}\cdot \mathrm{sm.xPER4}\\ \mathrm{xPER4}\left(T=\mathrm{sm.t0}\right)=0\end{array}$
$\mathrm{Gv}=\mathrm{sm.delta}\cdot \mathrm{sm.xPER4}$

#### Variables

$\mathrm{depot}\left(\mathrm{ADM}=1,\mathrm{TARGET}=\mathrm{sm.X1}\right)$
$\mathrm{depot}\left(\mathrm{ADM}=2,\mathrm{TARGET}=\mathrm{sm.X}\right)$
$\mathrm{depot}\left(\mathrm{ADM}=3,\mathrm{TARGET}=\mathrm{sm.Z1}\right)$
$\mathrm{depot}\left(\mathrm{ADM}=4,\mathrm{TARGET}=\mathrm{sm.Z}\right)$
$\mathrm{depot}\left(\mathrm{ADM}=5,\mathrm{TARGET}=\mathrm{sm.xHL1},P=\frac{1}{\mathrm{pm.F}}\right)$
$\mathrm{depot}\left(\mathrm{ADM}=6,\mathrm{TARGET}=\mathrm{sm.xHL2},P=\frac{1}{\mathrm{pm.F}}\right)$

### Observation Model: $\mathrm{om1}$

#### Continuous Observation

$Y=\mathrm{sm.G}+\mathrm{additiveError}\left(\mathrm{additive}=\mathrm{pm.alpha}\right)+\mathrm{pm.epsilon}$

## External Dataset

 OID $\mathrm{nm_ds}$ Tool Format Monolix

### File Specification

 Format $\mathrm{csv}$ Delimiter comma File Location Simulated_glucoseKinetics.csv

### Column Definitions

Column ID Position Column Type Value Type
$\mathrm{TIME}$
$1$
$\mathrm{idv}$
$\mathrm{real}$
$\mathrm{DV}$
$2$
$\mathrm{dv}$
$\mathrm{real}$
$\mathrm{MDV}$
$3$
$\mathrm{mdv}$
$\mathrm{int}$
$\mathrm{AMT}$
$4$
$\mathrm{dose}$
$\mathrm{real}$
$\mathrm{RATE}$
$5$
$\mathrm{rate}$
$\mathrm{real}$
$\mathrm{ID}$
$6$
$\mathrm{id}$
$\mathrm{int}$
$\mathrm{INS}$
$7$
$\mathrm{reg}$
$\mathrm{real}$
$\mathrm{GLU}$
$8$
$\mathrm{reg}$
$\mathrm{real}$
$\mathrm{CMT}$
$9$
$\mathrm{adm}$
$\mathrm{int}$
$\mathrm{TOBS}$
$\mathrm{10}$
$\mathrm{reg}$
$\mathrm{real}$
$\mathrm{T1}$
$\mathrm{11}$
$\mathrm{reg}$
$\mathrm{real}$
$\mathrm{GLU1}$
$\mathrm{12}$
$\mathrm{reg}$
$\mathrm{real}$
$\mathrm{INS1}$
$\mathrm{13}$
$\mathrm{reg}$
$\mathrm{real}$

### Column Mappings

Column Ref Modelling Mapping
$TIME$
$T$
$DV$
$\mathrm{om1.Y}$
$AMT$
$\left\{\begin{array}{lll}\text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}1& \text{if}& \mathrm{AMT}=\text{1}\\ \text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}2& \text{if}& \mathrm{AMT}=\text{2}\\ \text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}3& \text{if}& \mathrm{AMT}=\text{3}\\ \text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}4& \text{if}& \mathrm{AMT}=\text{4}\\ \text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}5& \text{if}& \mathrm{AMT}=\text{5}\\ \text{blk:}\phantom{\rule{1mm}{0ex}}\mathrm{sm},\phantom{\rule{2mm}{0ex}}\text{adm:}\phantom{\rule{1mm}{0ex}}6& \text{if}& \mathrm{AMT}=\text{6}\end{array}$
$ID$
$\mathrm{vm_mdl.ID}$
$INS$
$\mathrm{sm.INS}$
$GLU$
$\mathrm{sm.GLU}$
$TOBS$
$\mathrm{sm.TOBS}$
$T1$
$\mathrm{sm.T1}$
$GLU1$
$\mathrm{sm.GLU1}$
$INS1$
$\mathrm{sm.INS1}$

## Simulation Step

 OID $\mathrm{simulStep_1}$

### Variable Assignments

Variable Value
pm.typ_KmG
$3.88$
pm.typ_Vmax0
$338$
pm.typ_Emax
$4812$
pm.typ_gamma
$1.62$
pm.typ_KmI
$784$
pm.typ_t12I
$15.9$
pm.typ_t12G
$0.7$
pm.typ_V
$12648$
pm.typ_flambda3
$0.0582$
pm.typ_flambda2
$0.154$
pm.typ_w1
$0.609$
pm.typ_fw2
$0.901$
pm.typ_F
$2688$
pm.var_KmG
$0.219$
pm.var_Vmax0
$0$
pm.var_Emax
$0.112$
pm.var_gamma
$0.111$
pm.var_KmI
$0.263$
pm.var_t12I
$0.151$
pm.var_t12G
$0$
pm.var_V
$0.0557$
pm.var_flambda3
$0.179$
pm.var_flambda2
$0$
pm.var_w1
$0.773$
pm.var_fw2
$0$
pm.var_F
$0$
pm.corr_gamma_KmI
$-0.44$
pm.alpha
$0.014$
pm.sigma
$1$

### Operations

#### Operation: $1$

 Op Type generic
Name Value

## Step Dependencies

Step OID Preceding Steps
$\mathrm{simulStep_1}$