Reliability of parameter estimates from models applied to respiratory impedance data

Many previous studies have fit lumped parameter models to respiratory input (Zin) and transfer (Ztr) impedance data. For frequency ranges higher than 4-32 Hz, a six-element model may be required in which an airway branch (with a resistance and inertance) is separated from a tissue branch (with a resistance, inertance, and compliance) by a shunt compliance. A sensitivity analysis is applied to predict the effects of frequency range on the accuracy of parameter estimates in this model obtained from Zin or Ztr data. Using a parameter set estimated from experimental data between 4 and 64 Hz in dogs, both Zin and Ztr were simulated from 4 to 200 Hz. Impedance sensitivity to each parameter was also calculated over this frequency range. The simulation predicted that for Zin a second resonance occurs near 80 Hz and that the impedance is considerably more sensitive to several of the parameters at frequencies surrounding this resonance than at any other frequencies. Also, unless data is obtained at very high frequencies (where the model is suspect), Zin data provides more accurate estimates than Ztr data. After adding random noise to the simulated Zin data, we attempted to extract the original parameters by using a nonlinear regression applied to three frequency ranges: 4-32, 4-64, and 4-110 Hz. Estimated parameters were substantially incorrect when using only 4- to 32-Hz or 4- to 64-Hz data, but nearly correct when fitting 4- to 110-Hz data. These results indicate that respiratory system parameters can be more accurately extracted from Zin than Ztr, and to make physiological inferences from parameter estimates based on Zin impedance data in dogs, the data must include frequencies surrounding the second resonance.     

Fabrication of polyaniline/graphene oxide composites for implementation in humidity sensing

This work reports the measurement of impedance variations under various humidity conditions at frequency ranges between 100 Hz and 5 MHz. An electrochemical polymerization process has been used in the synthesis including varying the mass ratios of graphene oxide (GO) in polyaniline. An electrochemical deposition method has been used to produce a sample film on an indium tin oxide glass slide. The percentage relative humidity (RH%) of the samples has been estimated to be 20–90%. Impedance and humidity had an inverse relationship, i.e. the impedance value decreased with an increase in humidity. In contrast with platinum capacitive humidity sensors (HS), the GO-based HS had a sensitivity of 75–99%, which was ~10-fold more than that of traditional sensors. With three different parameter weight % of GO, the frequency range have been 100 Hz to 5 MHz and RH% has been found to 20–90%. The HS showed a fast response and recovery time. Therefore, GO appears to be a useful material for building HS with high sensitivity for a comprehensive approach.     

Sensitivity analysis of respiratory parameter uncertainties: impact of criterion function form and constraints

A sensitivity analysis based on weighted least-squares regression is presented to evaluate alternative methods for fitting lumped-parameter models to respiratory impedance data. The goal is to maintain parameter accuracy simultaneously with practical experiment design. The analysis focuses on predicting parameter uncertainties using a linearized approximation for joint confidence regions. Applications are with four-element parallel and viscoelastic models for 0.125- to 4-Hz data and a six-element model with separate tissue and airway properties for input and transfer impedance data from 2-64 Hz. The criterion function form was evaluated by comparing parameter uncertainties when data are fit as magnitude and phase, dynamic resistance and compliance, or real and imaginary parts of input impedance. The proper choice of weighting can make all three criterion variables comparable. For the six-element model, parameter uncertainties were predicted when both input impedance and transfer impedance are acquired and fit simultaneously. A fit to both data sets from 4 to 64 Hz could reduce parameter estimate uncertainties considerably from those achievable by fitting either alone. For the four-element models, use of an independent, but noisy, measure of static compliance was assessed as a constraint on model parameters. This may allow acceptable parameter uncertainties for a minimum frequency of 0.275-0.375 Hz rather than 0.125 Hz. This reduces data acquisition requirements from a 16- to a 5.33- to 8-s breath holding period. These results are approximations, and the impact of using the linearized approximation for the confidence regions is discussed.     

Development of a new rapid measurement technique for fish embryo membrane permeability studies using impedance spectroscopy

Information on fish embryo membrane permeability is vital in their cryopreservation. Whilst conventional volumetric measurement based assessment methods have been widely used in fish embryo membrane permeability studies, they are lengthy and reduce the capacity for multi-embryo measurement during an experimental run. A new rapid 'real-time' measurement technique is required to determine membrane permeability during cryoprotectant treatment. In this study, zebrafish (Danio rerio) embryo membrane permeability to cryoprotectants was investigated using impedance spectroscopy. An embryo holding cell, capable of holding up to 10 zebrafish embryos was built incorporating the original system electrods for measuring the impedance spectra. The holding cell was tested with deionised water and a series of KCl solutions with known conductance values to confirm the performance of the modified system. Untreated intact embryos were then tested to optimise the loading capacity and sensitivity of the system. To study the impedance changes of zebrafish embryos during cryoprotectant exposure, three, six or nine embryos at 50% epiboly stage were loaded into the holding cell in egg water, which was then removed and replaced by 0.5, 1.0, 2.0 or 3M methanol or dimethyl sulfoxide (DMSO). The impedance changes of the loaded embryos in different cryoprotectant solutions were monitored over 30 min at 22 degrees C, immediately following embryo exposure to cryoprotectants, at the frequency range of 10-10(6)Hz. The impedance changes of the embryos in egg water were used as controls. Results from this study showed that the optimum embryo loading level was six embryos per cell for each experimental run. The optimum frequency was identified at 10(3.14) or 1,380 Hz which provided good sensitivity and reproducibility. Significant impedance changes were detected after embryos were exposed to different concentrations of cryoprotectants. The results agreed well with those obtained from conventional volumetric based studies.     

Frequency domain studies of impedance characteristics of biological cells using micropipet technique. I. Erythrocyte.

This study aims at precise measurement of the membrane capacity and its frequency dependence of small biological cells using the micropipet technique. The use of AC fields as an input signal enables the magnitude and phase angle of membrane impedance to be measured at various frequencies. The micropipet technique was applied to human erythrocyte, and passive membrane capacity and conductivity were determined between 4 Hz and 10 KHz. Membrane capacity thus determined changed from 1.05 to 0.73 microF/cm2 between 4 Hz and 10 KHz. In addition to the micropipet technique, we used suspension method between 50 KHz and 10 MHz for the purpose of supplementing the new method with the one which has been in use for many years. We obtained a membrane capacity of 0.65-0.8 microF/cm2 using this technique. These values agree with the capacitance obtained with the micropipet method. Although this paper discusses only human erythrocytes, the study has been performed with lymphocytes and various forms of cancer cells. This paper is the first of the series of reports on frequency domain studies of the impedance characteristics of various biological cells.     

Analysis of lung parenchyma as a parametric porous medium.

The dynamic behavior of the lung in health and disease depends on its viscoelastic properties. To better understand these properties, several mathematical models have been utilized by many investigators. In the present work, we present a new approach that characterizes the dynamics of gas flow into a viscoelastic porous medium that models the lung structure. This problem is considered in terms of the lung input impedance on a macro level and parenchymal tissue impedance on the level of an alveolar wall. We start from a basic theoretical analysis in which macroscopic tissue deformations are represented in accordance with the linearized Navier-Stokes equations. This approach has strong theoretical underpinnings in other situations but has not been applied to analyze the impedance of the inflated lung. Our analysis provides a theoretical basis for analyzing the interaction between flow into the lungs as a biophysical diffusion process and parenchymal viscoelasticity described phenomenologically, within the frameworks of standard viscoelasticity and structural damping. This lung impedance incorporates parameters of porosity, permeability, and viscoelasticity on micro and macro levels of parenchymal tissue. The analysis shows the theoretical basis of the transformation from the impedance of alveolar walls or isolated tissue strips to that of the intact parenchyma. We also show how the loading impedance at the lung boundary may have a significant impact on the dynamic behavior of whole lung viscoelasticity. Our analysis may be useful in directing specific tests of different models and for analyzing experimental measurements of viscoelastic parameters of lung material under normal and pathological conditions.     

Application of 1D blood flow models of the human arterial network to differential pressure predictions

A new application of 1D models of the human arterial network is proposed. We take advantage of the sensitivity of the models predictions for the pressure profiles within the main aorta to key model parameter values. We propose to use the patterns in the predicted differences from a base case as a way to infer to the most probable changes in the parameter values. We demonstrate this application using an impedance model that we have recently developed (Johnson, 2010). The input model parameters are all physiologically related, such as the geometric dimensions of large arteries, various blood properties, vessel elasticity, etc. and can therefore be patient specific. As a base case, nominal values from the literature are used. The necessary information to characterize the smaller arteries, arterioles, and capillaries is taken from a physical scaling model (West, 1999). Model predictions for the effective impedance of the human arterial system closely agree with experimental data available in the literature. The predictions for the pressure wave development along the main arteries are also found in qualitative agreement with previous published results. The model has been further validated against our own measured pressure data in the carotid and radial arteries, obtained from healthy individuals. Upon changes in the value of key model parameters, we show that the differences seen in the pressure profiles correspond to qualitatively different patterns for different parameters. This suggests the possibility of using the model in interpreting multiple pressure data of healthy/diseased individuals.     

Optimal parametric sensitivity control for the estimation of kinetic parameters in bioreactors

In this paper the well-known problem of optimal input design is considered. In particular, the focus is on input design for the estimation of kinetic parameters in bioreactors. The problem is formulated as follows: given the model structure (f,g), which is assumed to be affine in the input, and the specific parameter of interest theta;(k) find a feedback law that maximizes the sensitivity of the model output to the parameter under different flow conditions in the bioreactor and, possibly, minimize the input or state costs. Analytical solutions to these problems are presented. As an example a bioreactor with a biomass that grows according to the well-known Monod kinetics is considered.     

Real-time measurement of PMA-induced cellular alterations by microelectrode array-based impedance spectroscopy

For a feasible and cost-effective impedance measurement of cellular alterations in real-time, we combined commercially available microelectrode arrays (MEAs), consisting of 60 microelectrodes, with a conventional impedance analyzer. For proof of principle, a breast carcinoma cell line (MCF-7) was cultured on MEAs, and cellular alterations were measured by impedance spectroscopy at a frequency ranging from 10 Hz to 1 MHz. Cells were stimulated with phorbol 12-myristate 13-acetate (PMA) at different concentrations to activate protein kinase C (PKC)-mediated extra- and intracellular changes. By addition of 0.03 microM PMA, an increase of the relative impedance (Z(rel)) was observed after 10 min with a maximum at 1 kHz. Moreover a gradual elevation of the impedance was measured 60 min after stimulation with PMA. If 0.3 microM PMA was applied, the maximal amplitude of the relative impedance after 60 min shifted from 1 kHz (0.03 microM PMA) to 150 Hz. Subsequently, the impedance was further increased up to 90 min after PMA application, after which the impedance reduced after 240 min. Since we could use MEAs for at least 10 times without affecting the sensitivity, our study revealed that commercially available MEAs comprising nanocolumnar titanium nitrite electrodes are suitable microstructures for a highly reproducible and cost-effective multisite measurement of intracellular processes by impedance spectroscopy.     

Muscular performance modeling of the upper limb in static postures

The purpose of the present paper is to describe and evaluate the polynomial models for predicting the muscular work capacity of the upper limb during sustained holding tasks. This research was concerned with the relationship between indicators of performance, i.e., specific posture or specific level of maximum voluntary contraction (MVC), and then modeling the functional data based on experimental results to estimate factors that may have an effect on task performance. To this end, we designed an experiment using 10 subjects in which each subject performed sustained isometric shoulder and elbow flexion endurance exercise under 27 conditions [3 shoulder angles (SA)x3 elbow angles (EA)x3 levels of %MVC]. Throughout all experiments, subjective perception of effort was assessed using the Borg scale, every 60, 30, and 10 s during the 20%, 40%, and 60% MVC tests, respectively. Proposal models were represented by three approaches: model A: estimation of endurance time (ET), with input variables such as SA, EA, and %MVC; model B: estimation of recommendation time (RT, the time during which the operator was able to maintain a position under the desired condition), with input variables such as SA, EA, %MVC, and required rate on the Borg scale; and model C: estimation of limit strength or %MVC, with input variables such as SA, EA, request limit time for work (LT), and required rate on the Borg scale. Statistical analysis indicated that the three proposal estimation models based on polynomial regression functions showed high significance (p<0.0001). The proposal models suggested and recommended the possibility of finding the best positions entailing the reduction and minimization of total muscular strain from manual material handling tasks in different work situations, with the consequent increase in work efficiency.     

Microcirculation impedance analysis in cat lung

An analysis of pulsatile microcirculation in cat lung, with special attention to the pulmonary microvascular impedance, is presented. A theoretical calculation is made on the basis of a complete set of experimental data on the morphology and elasticity of cat's pulmonary capillary sheets. The transfer matrix of the pulmonary microvascular impedance is obtained. The input impedance at the capillary entrance and exit are determined. The input impedance at the pulmonary arterial trunk is compared under various physiological conditions. It is shown that although the impact of pulmonary microcirculation on the relationship between the steady mean flow and pressure in the pulmonary arteries and veins is decisively large, the influence of the alveolar microcirculation on the input impedance at the pulmonary arterial trunk is small.     

The controversy about spinal neuronal nitric oxide synthase: under which conditions is it up- or downregulated?

In recent years, the regulation of the synthesis of nitric oxide (NO) in the central nervous system has attracted much interest because it has been shown that NO is involved in a wide variety of functions such as neuroprotection, neurotoxicity, neurotransmission, and neuroplas- ticity under physiological and pathophysiological conditions. However, the use of different detection techniques for neuronal nitric oxide synthase (nNOS), different animal species, and different experimental lesions has led to contradictory results concerning the direction of changes in spinal nNOS expression. This paper summarizes the available data on the expression on nNOS in the spinal cord under physiological and pathological conditions and tries to extract some of the basic mechanisms that underlie neuronal up- or downregulation of this enzyme. Wherever possible, results obtained with the NADPH-dependent diaphorase reaction are also included for reasons of comparison. The main conclusion is that changes in spinal nNOS expression critically depend on the type of afferent fibres activated by a specific lesion as well as the intensity and duration of input to the spinal cord. This input may be further modified by supraspinal influences. Thus the exact composition of these factors, which is undoubtfully highly variable between different experimental models, appears to determine whether the spinal NO system responds with an up- or downregulation of nNOS expression or in a bidirectional way. With regard to the diaphorase reaction it is becoming increasingly clear that under pathological conditions data obtained with this reaction differ markedly from those obtained with immunohistochemical visualization of nNOS.     

Estimation of hurdle clearance parameters using a monocular human motion tracking method

This paper presents a method of monocular human motion tracking for estimation of hurdle clearance kinematic parameters. The analysis involved 10 image sequences of five hurdlers at various training levels. Recording of the sequences was carried out under simulated starting conditions of a 110 m hurdle race. The parameters were estimated using the particle swarm optimization algorithm and they are based on analysis of the images recorded with a 100 Hz camera. The proposed method does not involve using any special clothes, markers, inertial sensors, etc. As the quality criteria, the mean absolute error and mean relative error were used. The level of computed errors justifies the use of this method to estimate hurdle clearance parameters.     

Metallic oxide CdIn2O4 films for the label free electrochemical detection of DNA hybridization

DNA functionalised semiconductor metallic oxide electrodes have been developed for the direct electrochemical detection of DNA hybridization, without labelling or the introduction of a redox couple. Conductive CdIn(2)O(4) thin films with controlled properties were deposited on glass substrates using an aerosol pyrolysis technique. The films exhibit a polycrystalline microstructure with a surface roughness of 1.5 nm (r.m.s.) and an electrical resistivity ranging between 1 and 3 x 10(-3) Omega cm. These electrodes were functionalised using hydroxylation and silanisation steps, to allow the binding of DNA probe sequences (20 bases). The electrical detection of DNA hybridization with complementary sequences has been performed using electrochemical impedance spectrometry (EIS) measuring the variation of impedance before and after hybridization. Two set-ups were used, a standard set-up including three electrodes and a set-up including two symmetrical electrodes. In both configurations, a significant increase of the impedance modulus, more particularly of the real part of the impedance (160-225% according to the electrochemical cell used) has been obtained over a frequency range of 10-10(5)Hz. DNA hybridization has also been systematically confirmed using the fluorescence spectrometry. This study emphasizes the high sensitivity of the CdIn(2)O(4) as a working electrode for the detection of biological events occurring at the electrode surface.     

Dual assessment of airway area profile and respiratory input impedance from a single transient wave.

This report concerns the inference of geometric and mechanical airway characteristics based on information derived from a single transient planar wave recorded at the airway opening. We describe a new method to simultaneously measure upper airway area and respiratory input impedance by performing dual analysis of a single pressure wave. The algorithms required to reconstruct airway dimensions and mechanical characteristics were developed, implemented, and tested with reference to known physical models. Our method appears suitable to estimate, even under severe intensive care unit conditions, the respiratory system frequency response (above 10 Hz) in intubated patients and the patency of the endotracheal tube used to connect the patients to the ventilator.     

Anaerobic digestion model No. 1-based distributed parameter model of an anaerobic reactor: II. Model validation

In this study, an ADM1-based distributed parameter model was validated using experimental results obtained in a laboratory-scale 10 L UASB reactor. Sensitivity analysis of the model parameters was used to select four parameters for estimation by a numerical procedure while other parameters were accepted from ADM1 benchmark simulations. The parameter estimation procedure used measurements of liquid phase components obtained at different sampling points in the reactor and under different operating conditions. Model verification used real time fluorescence-based measurements of chemical oxygen demand and volatile fatty acids at four sampling locations in the reactor. Overall, the distributed parameter model was able to describe the distribution of liquid phase components in the reactor and adequately simulated the effect of external recirculation on degradation efficiency. The model can be used in the design, analysis and optimization of UASB reactors.     

A comparison of different methods for estimating single-trial P300 latencies

Inferences from comparative analyses of reaction time and P300 latency are stronger when the various aspects of the distribution across trials are treated in the same way for both variables. To this end, a number of studies have resorted to estimation of P300 latency at the single-trial level. This report presents a comparative evaluation of two common methods for such single-trial analysis, i.e., peak-picking and template-matching. Both methods were applied to a representative set of real data, comprising different task conditions and two age groups. Relevant scoring parameters were varied: low-pass filter settings (down to 0.94 Hz) for peak-picking, template duration (250–970 msec) and use of covariance vs. correlation for template-matching, and use of a noise-range criterion for both methods. It is concluded that peak-picking with a 3.4 Hz filter, and template-matching using covariance and template duration between 600 and 800 msec, are best in terms of sensitivity and reliability, with peak-picking surpassing template-matching. Also, the marked increase in the number of rejected trials when the noise-range criterion was applied resulted in unwanted modulation of behavioral effects of task conditions and age groups.     

Low-frequency respiratory system resistance in the normal dog during mechanical ventilation

The low-frequency resistances of the respiratory system, lung, and chest wall were investigated in four anesthetized paralyzed dogs mechanically ventilated at various frequencies between 0.08 and 0.83 Hz. The resistances were calculated by three different methods: 1) as the real part of the complex impedance obtained from regular ventilation data, 2) as the effective resistance of a two-compartment model fitted to the same data, and 3) as the resistance of a single-compartment model fitted to data obtained during sinusoidal ventilation at various frequencies. Alveolar pressures were measured by a closed-chest alveolar capsule technique and afforded a direct measure of airways resistance. All three resistance estimates were very similar and decreased markedly with frequency between 0 and 1 Hz. The real part of lung impedance at the higher frequencies investigated (around 5 Hz) closely approximated airways resistance, as predicted by the eight-parameter viscoelastic model of respiratory mechanics proposed by Bates et al. (J. Appl. Physiol. 67:2276-2285, 1989).     

A nonlinear viscoelastic model of lung tissue mechanics.

There have been a number of attempts recently to use linear models to describe the low-frequency (0-2 Hz) dependence of lung tissue resistance (Rti) and elastance (Eti). Only a few attempts, however, have been made to account for the volume dependence of these quantities, all of which require the tissues to be plastoelastic. In this paper we specifically avoid invoking plastoelasticity and develop a nonlinear viscoelastic model that is also capable of accounting for the nonlinear and frequency-dependent features of lung tissue mechanics. The model parameters were identified by fitting the model to data obtained in a previous study from dogs during sinusoidal ventilation. The model was then used to simulate pressure and flow data by use of various types of ventilation patterns similar to those that have been employed experimentally. Rti and Eti were estimated from the simulated data by use of four different estimation techniques commonly applied in respiratory mechanics studies. We found that the estimated volume dependence of Rti and Eti is sensitive to both the ventilation pattern and the estimation technique, being in error by as much as 217 and 22%, respectively.