Spearman simultaneous estimation for a compartmental model Technical report no. 9
Compartmental model, represented as system of linear combination of exponentials with common exponential parameters, for tracer experiments
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Compartmental model, represented as system of linear combination of exponentials with common exponential parameters, for tracer experiments
Simultaneous estimation procedure for parameters in multiple equation regression model
A new two-compartment model for the passage of particles through the gastro-intestinal tract of ruminants is proposed. In this model, a gamma distribution of lifetimes is introduced in the first compartment; thereby, passage from that compartment becomes time-dependent. This modification is strongly suggested by the physical alteration which certain substances, e.g. hay particles, undergo in the digestive process. The proposed model is applied to experimental data.
Biomathematical model and least square estimation from time series data of discrete particle population in stochastic compartmental system
Electrical impedance spectrometry involves measurement of the complex resistance of a load at multiple frequencies. With this information in the form of impedance magnitude and phase, or resistance and reactance, basic structure or function of the load can be estimated. The "load" targeted for measurement and estimation in this study consisted of the water-bearing tissues of the human calf. It was proposed and verified that by measuring the electrical impedance of the human calf and fitting this data to a model of fluid compartments, the lumped-model volume of intracellular and extracellular spaces could be estimated, By performing this estimation over time, the volume dynamics during application of stimuli which affect the direction of gravity can be viewed. The resulting data can form a basis for further modeling and verification of cardiovascular and compartmental modeling of fluid reactions to microgravity as well as countermeasures to the headward shift of fluid during head-down tilt or spaceflight.
A substantial foundation of tools for model-based reasoning with incomplete knowledge was developed: QSIM (a qualitative simulation program) and its extensions for qualitative simulation; Q2, Q3 and their successors for quantitative reasoning on a qualitative framework; and the CC (component-connection) and QPC (Qualitative Process Theory) model compilers for building QSIM QDE (qualitative differential equation) models starting from different ontological assumptions. Other model-compilers for QDE's, e.g., using bond graphs or compartmental models, have been developed elsewhere. These model-building tools will support automatic construction of qualitative models from physical specifications, and further research into selection of appropriate modeling viewpoints. For monitoring and diagnosis, plausible hypotheses are unified against observations to strengthen or refute the predicted behaviors. In MIMIC (Model Integration via Mesh Interpolation Coefficients), multiple hypothesized models of the system are tracked in parallel in order to reduce the 'missing model' problem. Each model begins as a qualitative model, and is unified with a priori quantitative knowledge and with the stream of incoming observational data. When the model/data unification yields a contradiction, the model is refuted. When there is no contradiction, the predictions of the model are progressively strengthened, for use in procedure planning and differential diagnosis. Only under a qualitative level of description can a finite set of models guarantee the complete coverage necessary for this performance. The results of this research are presented in several publications. Abstracts of these published papers are presented along with abtracts of papers representing work that was synergistic with the NASA grant but funded otherwise. These 28 papers include but are not limited to: 'Combined qualitative and numerical simulation with Q3'; 'Comparative analysis and qualitative integral representations'; 'Model-based monitoring of dynamic systems'; 'Numerical behavior envelopes for qualitative models'; 'Higher-order derivative constraints in qualitative simulation'; and 'Non-intersection of trajectories in qualitative phase space: a global constraint for qualitative simulation.'
An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials with an Investigative New Drug (IND) protocol. The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trials with INSCOP. Methods: Twelve healthy human subjects were administered three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min and 24 h after dosing and scopolamine concentrations were measured by using a validated LC‐MS‐MS assay. Pharmacokinetic Compartmental models, using actual dosing and sampling times, were built using Phoenix (version 1.2). Model selection was based on the likelihood ratio test on the difference of criteria (‐2LL) and comparison of the quality of fit plots. Results: The best structural model for INSCOP (minimal ‐2LL= 502.8) was established. It consisted of one compartment each for plasma, saliva and urine, respectively, which were connected with linear transport processes except the nonlinear PK process from plasma to saliva compartment. The best‐fit estimates of PK parameters from individual PK compartmental analysis and Population PK model analysis were shown in Tables 1 and 2, respectively. Conclusion: A population PK model that could predict population and individual PK of scopolamine in plasma, saliva and urine after dosing was developed and validated. Incorporating a non‐linear transfer from plasma to saliva compartments resulted in a significantly improved model fitting. The model could be used to predict scopolamine plasma concentrations from salivary and urinary drug levels, allowing non‐invasive therapeutic monitoring of scopolamine in space and other remote environments.
An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials for an Investigative New Drug (IND). The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial with INSCOP. METHODS: Twelve healthy human subjects were administered three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min to 24 h after dosing and scopolamine concentrations measured by using a validated LC-MS-MS assay. Pharmacokinetic Compartmental models, using actual dosing and sampling times, were built using Phoenix (version 1.2). Model discrimination was performed, by minimizing the Akaike Information Criteria (AIC), maximizing the coefficient of determination (r²) and by comparison of the quality of fit plots. RESULTS: The best structural model to describe scopolamine disposition after INSCOP administration (minimal AIC =907.2) consisted of one compartment for plasma, saliva and urine respectively that were inter-connected with different rate constants. The estimated values of PK parameters were compiled in Table 1. The model fitting exercises revealed a nonlinear PK for scopolamine between plasma and saliva compartments for K21, Vmax and Km. CONCLUSION: PK model for INSCOP was developed and for the first time it satisfactorily predicted the PK of scopolamine in plasma, saliva and urine after INSCOP administration. Using non-linear PK yielded the best structural model to describe scopolamine disposition between plasma and saliva compartments, and inclusion of non-linear PK resulted in a significant improved model fitting. The model can be utilized to predict scopolamine plasma concentration using saliva and/or urine data that allows non-invasive assessment of pharmacotherapeutics of scopolamine in space and other remote environments without requiring blood sampling.
An intranasal gel dosage formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness (SMS). The bioavailability and pharmacokinetics (PK) were evaluated under IND (Investigational New Drug) guidelines. The aim of the project was to develop a PK model that can predict the relationships among plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial protocol with INSCOP. Twelve healthy human subjects were administered at three dose levels (0.1, 0.2 and 0.4 mg) of INSCOP. Serial blood, saliva and urine samples were collected between 5 min to 24 h after dosing and scopolamine concentrations were measured by using a validated LC‐MS‐MS assay. PK compartmental models, using actual dosing and sampling time, were established using Phoenix (version 1.2). Model selection was based on a likelihood ratio test on the difference of criteria (‐2LL (i.e. log-likelihood ratio test)) and comparison of the quality of fit plots. The results: Predictable correlations among scopolamine concentrations in compartments of plasma, saliva and urine were established, and for the first time the model satisfactorily predicted the population and individual PK of INSCOP in plasma, saliva and urine. The model can be utilized to predict the INSCOP plasma concentration by saliva and urine data, and it will be useful for monitoring the PK of scopolamine in space and other remote environments using non‐invasive sampling of saliva and/or urine.
Calcium loss from bones during space flight creates a risk for astronauts who travel into space, and may prohibit space flights to other planets. The problem of calcium loss during space flight has been studied using animal models, bed rest (as a ground-based model), and humans in-flight. In-flight studies have typically documented bone loss by comparing bone mass before and after flight. To identify changes in metabolism leading to bone loss, we have performed kinetic studies using stable isotopes of calcium. Oral (Ca-43) and intravenous (Ca-46) tracers were administered to subjects (n=3), three-times before flight, once in-flight (after 110 days), and three times post-flight (on landing day, and 9 days and 3 months after flight). Samples of blood, saliva, urine, and feces were collected for up to 5 days after isotope administration, and were analyzed for tracer enrichment. Tracer data in tissues were analyzed using a compartmental model for calcium metabolism and the WinSAAM software. The model was used to: account for carryover of tracer between studies, fit data for all studies using the minimal number of changes between studies, and calculate calcium absorption, excretion, bone calcium deposition and bone calcium resorption. Results showed that fractional absorption decreased by 50% during flight and that bone resorption and urinary excretion increased by 50%. Results were supported by changes in biochemical markers of bone metabolism. Inflight bone loss of approximately 250 mg Ca/d resulted from decreased calcium absorption combined with increased bone resorption and excretion. Further studies will assess the time course of these changes during flight, and the effectiveness of countermeasures to mitigate flight-induced bone loss. The overall goal is to enable human travel beyond low-Earth orbit, and to allow for better understanding and treatment of bone diseases on Earth.
A program has been written for the identification of parameters in certain linear systems. These systems appear in biomedical problems, particularly in compartmental models of pharmacokinetics. The method presented here assumes that some of the state variables are regularly modified by jump conditions. This simulates administration of drugs following some prescribed drug regime. Parameters are identified by a least-square fit of the linear differential system to a set of experimental observations. The method is especially suited when the interval of observation of the system is very long.
Realistic, biophysically-based compartmental models were constructed of several primary sensory interneurons in the cricket cercal sensory system. A dynamic atlas of the afferent input to these cells was used to set spatio-temporal parameters for the simulated stimulus-dependent synaptic inputs. We examined the roles of dendritic morphology, passive membrane properties, and active conductances on the frequency tuning of the neurons. The sensitivity of narrow-band low pass interneurons could be explained entirely by the electronic structure of the dendritic arbors and the dynamic sensitivity of the SIZ. The dynamic characteristics of interneurons with higher frequency sensitivity required models with voltage-dependent dendritic conductances.
Studies of calcium kinetics require administration of tracer doses of calcium and subsequent repeated sampling of biological fluids. This study was designed to develop techniques that would allow estimation of calcium kinetics by using small (micrograms) doses of isotopes instead of the more common large (mg) doses to minimize tracer perturbation of the system and reduce cost, and to explore the use of saliva sampling as an alternative to blood sampling. Subjects received an oral dose (133 micrograms) of 43Ca and an i.v. dose (7.7 micrograms) of 46Ca. Isotopic enrichment in blood, urine, saliva and feces was well above thermal ionization mass spectrometry measurement precision up to 170 h after dosing. Fractional calcium absorptions determined from isotopic ratios in blood, urine and saliva were similar. Compartmental modeling revealed that kinetic parameters determined from serum or saliva data were similar, decreasing the necessity for blood samples. It is concluded from these results that calcium kinetics can be assessed with micrograms doses of stable isotopes, thereby reducing tracer costs and with saliva samples, thereby reducing the amount of blood needed.
We cut serial sections through the medial part of the rat vestibular macula for transmission electron microscopic (TEM) examination, computer-assisted 3-D reconstruction, and compartmental modeling. The ultrastructural research showed that many primary vestibular neurons have an unmyelinated segment, often branched, that extends between the heminode (putative site of the spike initiation zone) and the expanded terminal(s) (calyx, calyces). These segments, termed the neuron branches, and the calyces frequently have spine-like processes of various dimensions with bouton endings that morphologically are afferent, efferent, or reciprocal to other macular neural elements. The major questions posed by this study were whether small details of morphology, such as the size and location of neuronal processes or synapses, could influence the output of a vestibular afferent, and whether a knowledge of morphological details could guide the selection of values for simulation parameters. The conclusions from our simulations are (1) values of 5.0 k omega cm2 for membrane resistivity and 1.0 nS for synaptic conductance yield simulations that best match published physiological results; (2) process morphology has little effect on orthodromic spread of depolarization from the head (bouton) to the spike initiation zone (SIZ); (3) process morphology has no effect on antidromic spread of depolarization to the process head; (4) synapses do not sum linearly; (5) synapses are electrically close to the SIZ; and (6) all whole-cell simulations should be run with an active SIZ.
An improved understanding of spaceflight-induced ocular pathology, including the loss of visual acuity, globe flattening, optic disk edema and distension of the optic nerve and optic nerve sheath, is of keen interest to space medicine. Cephalad fluid shift causes a profoundly altered distribution of fluid within the compartments of the head and body, and may indirectly generate phenomena that are biomechanically relevant to visual function, such as choroidal engorgement, compromised drainage of blood and cerebrospinal fluid (CSF), and altered translaminar pressure gradient posterior to the eye. The experimental body of evidence with respect to the consequences of fluid shift has not yet been able to provide a definitive picture of the sequence of events. On earth, elevated intracranial pressure (ICP) is associated with idiopathic intracranial hypertension (IIH), which can produce ocular pathologies that look similar to those seen in some astronauts returning from long-duration flight. However, the clinically observable features of the Visual Impairment and Intracranial Pressure (VIIP) syndrome in space and IIH on earth are not entirely consistent. Moreover, there are at present no experimental measurements of ICP in microgravity. By its very nature, physiological measurements in spaceflight are sparse, and the space environment does not lend itself to well-controlled experiments. In the absence of such data, numerical modeling can play a role in the investigation of biomechanical causal pathways that are suspected of involvement in VIIP. In this work, we describe the conceptual framework for modeling the altered compartmental fluid distribution that represents an equilibrium fluid distribution resulting from the loss of hydrostatic pressure gradient.
The integration of two types of physiological system simulations is presented. The long term model is a circulatory system model which simulates long term blood flow variations and compartmental fluid shifts. The short term models simulate transient phenomena of the respiratory, thermoregulatory, and pulsatile cardiovascular systems as they respond to stimuli such as LBNP, exercise, and environmental gaseous variations. An overview of the interfacing approach is described. Descriptions of the variable interface for long term to short term and between the three short term models are given.
Compact analytic solution to distribution theory of nonsteady state model illustrated with biological data - age dependency in compartmental analysis
Mechanical model for propellant sloshing in irregular compartmentation of S-4B Workshop configuration during roll and lateral excitations