Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “parameter estimation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14Linked to original sources

Estimating parameters of the family-size distribution in ascertainment sampling schemes: numerical results.

It is argued that, in any ascertainment sampling scheme using data from families of various sizes, there is never any need to assume a particular form for the (unknown) family-size distribution. There exists a simple conditional method, making no assumptions about the family-size distribution, that is always preferable to the assumption of any particular distributional form. Furthermore, the simplicity of the conditional method gives insights into properties of estimates of genetic and ascertainment parameters, which are not available when a particular form for the family-size distribution is assumed.

Biometry↗

Structure identifiability in metabolic pathways: parameter estimation in models based on the power-law formalism.

An important step in understanding a metabolic pathway is to identify its structure, in terms of the flow of material and information. In pursuing this goal, the available information for a given system is usually obtained from experiments in vitro and comes from different sources. Frequently, the final set of regulatory signals acting in the system in vivo is unclear, and some kind of test is needed on the intact system. Besides defining an appropriate experimental approach, identification of the regulatory pattern needs a theoretical framework in which the different experimental measurements can be evaluated and a final picture can be agreed on. Mathematical approaches based on sensitivity coefficients provide a useful tool for addressing this problem. Within this framework, the appropriate parameters are related to both the structure of the reaction network and the signals that regulate the target system. Thus the identification of the regulatory structure can be related to the estimation of the appropriate set of parameters. In pursuing this goal, we will show the limitations of using steady-state measurements and the usefulness of using dynamic data. We suggest a way to test the regulatory pattern in a given metabolic pathway by combining both kinds of data, and we show, by using a reference system, the potential of the method suggested.

Metabolism↗

Logistic growth curve of chickens: a comparison of techniques to estimate parameters.

Parameters of a mathematical function of growth, fit to the body weight curve of two randombred control populations of each sex of chickens from hatching through 45 weeks of age, were estimated. The logistic function was chosen from among growth formulae that express rate of gain as a function of weight at a given time and gain to be made. Two logistic parameters, growth-rate constant and age at the point of inflection, were estimated by the methods of sample quantiles and nonlinear regression from weekly mean body weights of 225 males and 281 females of the Rhode Island Red (RIR) line, and 164 males and 239 females of the White Leghorn (WL) line. Males had a larger growth-rate constant than females of the same line. The RIR line had a larger rate constant than the WL line, for each sex. Age at the point of inflection was similar for males and females in the RIR line, but smaller for males than females in the WL line. Sample quantiles yielded larger, less precise estimates of the growth-rate constant than nonlinear regression. Estimates of age at the point of inflection were usually smaller using sample quantiles.

Aging↗

Material parameter estimation with terahertz time-domain spectroscopy.

Imaging systems based on terahertz (THz) time-domain spectroscopy offer a range of unique modalities owing to the broad bandwidth, subpicosecond duration, and phase-sensitive detection of the THz pulses. Furthermore, the possibility exists for combining spectroscopic characterization or identification with imaging because the radiation is broadband in nature. To achieve this, we require novel methods for real-time analysis of THz waveforms. This paper describes a robust algorithm for extracting material parameters from measured THz waveforms. Our algorithm simultaneously obtains both the thickness and the complex refractive index of an unknown sample under certain conditions. In contrast, most spectroscopic transmission measurements require knowledge of the sample's thickness for an accurate determination of its optical parameters. Our approach relies on a model-based estimation, a gradient descent search, and the total variation measure. We explore the limits of this technique and compare the results with literature data for optical parameters of several different materials.

Algorithms↗

Endogenous model state and parameter estimation from an extensive batch experiment.

In this paper an extensive batch experiment of endogenous process behavior in an aerobic biodegradation process is presented. From these experimental data, comprising measurements of MLVSS (mixed liquor volatile suspended solids) and respiration rate, in a first step the states and unknown parameters in a four-compartmental model are reconstructed analytically. Subsequently, for a selected set of states and parameters, using the results of the previous step, a recursive state estimation procedure, in particular an Extended Kalman filter-based observer, is applied to deal with the noise properties of the data appropriately. From this it appears that the initially proposed model structure, and especially the hydrolysis term, has to be modified.

Aerobiosis↗

Stability parameter estimation at ambient temperature from studies at elevated temperatures.

The determination of specific kinetic constants k(i) in pH-profile studies is often undertaken at ambient temperature. However, when dealing with a drug substance that is stable at ambient temperature, the pH-profile study is conducted at a chosen elevated temperature and the kinetic parameters are given at this particular elevated temperature. But in stability studies we generally need kinetic constants at ambient or storage temperature for practical reasons (information and storage conditions of formulation). To assess this ambient kinetic information from studies at elevated temperatures, cumulative sequential steps are usually employed with very few statistical concerns on the final estimates. The statistical problems on these final estimates in cumulative procedures are highlighted in many papers. Because these stability parameters are useful for drug formulation and storage conditions, good practical decisions have to be made on the basis of statistically unbiased identified parameters. We propose in this paper a nonlinear model that allows the direct determination of specific activation energies E(ai) that are linked to the specific kinetic constants k(i). Hence, a mathematical relationship between drug concentration C, pH, temperature T, and time t is obtained. Kinetic data from acetylsalicylic acid (ASA) hydrolysis (first-order kinetics) are used to validate the model. The results show that it is possible to obtain directly, by an extrapolation procedure, the kinetic parameters (specific kinetic constants k(i), specific activation energies E(ai), and dissociation constant pK(a)) at low temperature from data gathered at elevated temperatures using more meaningful statistics.

Data Interpretation, Statistical↗

Linear and nonlinear ARMA model parameter estimation using an artificial neural network.

This paper addresses parametric system identification of linear and nonlinear dynamic systems by analysis of the input and output signals. Specifically, we investigate the relationship between estimation of the system using a feedforward neural network model and estimation of the system by use of linear and nonlinear autoregressive moving-average (ARMA) models. By utilizing a neural network model incorporating a polynomial activation function, we show the equivalence of the artificial neural network to the linear and nonlinear ARMA models. We compare the parameterization of the estimated system using the neural network and ARMA approaches by utilizing data generated by means of computer simulations. Specifically, we show that the parameters of a simulated ARMA system can be obtained from the neural network analysis of the simulated data or by conventional least squares ARMA analysis. The feasibility of applying neural networks with polynomial activation functions to the analysis of experimental data is explored by application to measurements of heart rate (HR) and instantaneous lung volume (ILV) fluctuations.

Computer Simulation↗

The amino acid requirements of growing turkeys. 1. Model construction and parameter estimation.

A model has been formulated to calculate the protein and amino acid requirements for growing turkeys. A requirement was taken as the sum of the needs for maintenance, for gain in carcass weight without feathers, and for the feathers. Parameters used in the model were obtained from a study of carcass composition along the growth period and from two balance trials with adult turkeys. The maintenance requirement was considered to be the sum of the needs for replacing skin and intestinal protein losses and for the obligatory creatinine excretion in the urine. Nitrogen retention was determined at 2.3 g/day in protein-fed adult (17 kg) male turkeys as an estimate for the replacement needs for skin protein loss. Creatinine excretion in these birds was 123 mg/day. When fed an N-free diet, 16.5 kg adult turkeys excreted a total of 3.5 g/day amino acids as a result of intestinal losses. The protein requirement for maintenance was estimated at 32 mg/day/g2/3 of body weight by summing up the different needs and dividing by .85, the coefficient of protein absorption. Protein and amino acid analyses of nude carcass and feathers and the growth curve were used to calculate the requirements for growth. These were then related to energy intake predicted by equations, which included previously determined coefficients for maintenance and growth (Hurwitz et al., 1980a), and the results of carcass lipid analysis. The calculated amino acid requirements were similar to those of the National Research Council (1977).

Amino Acids↗

Mechanical impedance of the canine diaphragm. Part 2. Theoretical model and parameter estimation.

In the paper the equation of motion of the small amplitude transverse forced vibration of a radially prestressed and circularly clamped thin membrane has been developed. The material of the membrane is considered to be homogeneous, isotropic, incompressible and viscoelastic. From the analytical solution of this equation the incremental mechanical impedance of the membrane was derived as a function of frequency, geometrical parameters and incremental viscoelastic coefficients of the material. The parameters of the model were fitted to experimental impedance data using a global optimisation procedure to obtain the incremental viscoelastic moduli of the canine diaphragm. The estimated quasi-static behaviour of the model is shown to be consistent with the results of experimental quasi-static measurements. It is concluded that the incremental viscoelastic moduli of a soft tissue and the stress dependence of these material coefficients can be determined by fitting the parameters of the model to the impedance data of that particular tissue.

Animals↗

Influence of method of model fitting and gender on estimated parameters for the relationship between starvation heat production and body size in the domestic fowl.

The relationship between starvation heat production (SHP) in kJ/d and body weight (W) in kg for domestic fowl was examined with compiled calorimetric data for 294 fowls (Johnson and Farrell, 1985). Linear regression analyses using ordinary least squares methods were performed before and after transformation of the data to a logarithmic scale (base = 10). Coefficients of determination (R2) were similar for untransformed and log-transformed data (R2's = .56 and .53, respectively). The slopes for untransformed data of males and females were different (P less than 0.033). When the model parameters for: SHP = b0 Wb1 were derived from log-transformed data, there were no significant differences in exponents for males vs. females (P greater than .95). However, when the coefficients were estimated using non-linear regression techniques, a difference was detected at the level of P = 0.072. It was determined that there is no advantage in describing the relationship of SHP as a non-linear function of W, nor is it appropriate to express physiological variables as a ratio of weight raised to some mass exponent without statistical justification. Using W or W2 as a covariable for the statistical analysis of physiological data would produce equally acceptable determinations of significance.

Animals↗

A simple steady state kinetic model for alcohol dehydrogenase. Use of reaction progress curves for parameters estimation at pH 7.4.

A simple rate equation for alcohol dehydrogenase was obtained by assuming independent binding sites for ethanol and NAD+ and fully competitive inhibition by the products of the reaction, acetaldehyde and NADH. A random binding order was also assumed. The rate equation is described by six parameters: four association constants (two for the substrates and two for the products of the reaction), Vf for the forward direction, and the equilibrium constant of the reaction. The six parameters were determined at pH 7.4 by numerical analysis of progress curves of reactions started with different concentrations of ethanol and NAD+. The parameters for alcohol dehydrogenase partially purified from rat liver were: Km for ethanol = 0.746 mM, Km for NAD+ = 0.0563 mM, Km for acetaldehyde = 7.07 microM, Km for NADH = 4.77 microM and Keq = 2.36 X 10(-4). The computed values allowed a very good simulation of the experimental progress curves and little variation was observed in the kinetic parameters when the reactions were started in the presence of either NADH or acetaldehyde.

Alcohol Dehydrogenase↗

Confidence limits for the parameter estimation in the dipole localization method on the basis of spatial correlation of background EEG.

A new residual function in the inverse problem of equivalent dipole localization methods is proposed which is based on the spatial correlation of the background EEG. This residual has the advantage that it allows the calculation of confidence limits for dipole model parameters. The method was applied to VEP data, and it was studied how the localization precision depends on the recording time of the EEG. It was found that the tangential position of an equivalent dipole can be located at 99% confidence in a region of the order 7 x 7mm for a head radius of 10cm, while the 99% confidence interval of the depth estimate is approximately 1cm, with a recording time of 20 minutes. It was also observed that an EEG recording time of more than 10-15 minutes is needed to obtain stable localization precision estimates.

Brain↗

Analytical approximations of sensitivities of steady state predictions to errors in parameter estimation: II. Michaelis-Menten kinetics.

Linear sensitivity theory is used to estimate the reliability of predictions of the minimum and maximum concentrations at steady state in the Michaelis-Menten model with i.v. bolus. The dependence of the relative errors in the predictions on the errors in the pharmacokinetic parameters is derived in an analytical form. It is shown that the quality of the predictions is not equally sensitive to all errors in parameters, and that the sensitivity factors vary with the degree of saturation of the system. An example of application for a drug, such as phenytoin, is discussed. It is suggested that sensitivity analysis may be useful in design of pharmacokinetic experiments aimed at the control of steady state levels for drugs with Michaelis-Menten kinetics.

Kinetics↗

The evolution of codon preferences in Drosophila: a maximum-likelihood approach to parameter estimation and hypothesis testing.

Synonymous codon usage in related species may differ as a result of variation in mutation biases, differences in the overall strength and efficiency of selection, and shifts in codon preference-the selective hierarchy of codons within and between amino acids. We have developed a maximum-likelihood method to employ explicit population genetic models to analyze the evolution of parameters determining codon usage. The method is applied to twofold degenerate amino acids in 50 orthologous genes from D. melanogaster and D. virilis. We find that D. virilis has significantly reduced selection on codon usage for all amino acids, but the data are incompatible with a simple model in which there is a single difference in the long-term Ne, or overall strength of selection, between the two species, indicating shifts in codon preference. The strength of selection acting on codon usage in D. melanogaster is estimated to be |Nes| approximately 0.4 for most CT-ending twofold degenerate amino acids, but 1.7 times greater for cysteine and 1.4 times greater for AG-ending codons. In D. virilis, the strength of selection acting on codon usage for most amino acids is only half that acting in D. melanogaster but is considerably greater than half for cysteine, perhaps indicating the dual selection pressures of translational efficiency and accuracy. Selection coefficients in orthologues are highly correlated (rho = 0.46), but a number of genes deviate significantly from this relationship.

Amino Acids↗

Determination of a constitutive relation for passive myocardium: II. Parameter estimation.

In the first paper of this series, we proposed a new transversely isotropic pseudostrain-energy function W for describing the biomechanical behavior of excised noncontracting myocardium. The specific functional form of W was inferred directly from biaxial data to be a polynomial function of two coordinate invariant measures of the finite deformation and five material parameters. In this paper, best-fit values of the material parameters are determined from biaxial data using a nonlinear least-squares regression. These values of the parameters are shown to be well-determined, and the final constitutive relation is shown to have good predictive capabilities. Since the proposed constitutive relation describes much broader classes of in-vitro biaxial data than previously proposed relations, it may be better applicable to analyses of stress in the passive heart.

Animals↗

Tests of a linear model of visual-vestibular interaction using the technique of parameter estimation.

The goal of this study was to test whether a superposition model of smooth-pursuit and vestibuloocular reflex (VOR) eye movements could account for the stability of gaze that subjects show as they view a stationary target, during head rotations at frequencies that correspond to natural movements. Horizontal smooth-pursuit and the VOR were tested using sinusoidal stimuli with frequencies in the range 1.0-3.5 Hz. During head rotation, subjects viewed a stationary target either directly or through an optical device that required eye movements to be approximately twice the amplitude of head movements in order to maintain foveal vision of the target. The gain of compensatory eye movements during viewing through the optical device was generally greater than during direct viewing or during attempted fixation of the remembered target location in darkness. This suggests that visual factors influence the response, even at high frequencies of head rotation. During viewing through the optical device, the gain of compensatory eye movements declined as a function of the frequency of head rotation (P < 0.001) but, at any particular frequency, there was no correlation with peak head velocity 9P > 0.23), peak head acceleration (P > 0.22) or retinal slip speed (P > 0.22). The optimal values of parameters of smooth-pursuit and VOR components of a simple superposition model were estimated in the frequency domain, using the measured responses during head rotation, as each subject viewed the stationary target through the optical device. We then compared the model's prediction of smooth-pursuit gain and phase, at each frequency, with values obtained experimentally. Each subject's pursuit showed lower gain and greater phase lag than the model predicted. Smooth-pursuit performance did not improve significantly if the moving target was a 10 deg x 10 deg Amsler grid, or if sinusoidal oscillation of the target was superimposed on ramp motion. Further, subjects were still able to modulate the gain of compensatory eye movements during pseudo-random head perturbations, making improved predictor performance during visual-vestibular interactions unlikely. We conclude that the increase in gain of eye movements that compensate for head rotations when subjects view, rather than imagine, a stationary target cannot be adequately explained by superposition of VOR and smooth-pursuit signals. Instead, vision may affect VOR performance by determining the context of the behavior.

Adult↗

Discrimination of the isozymes of human placental hexosaminidase by kinetic parameter estimation.

Human placental beta-N-acetylhexosaminidase (EC 3.2.1.52) (HEX) is a lysosomal glycosyl hydrolase with an acidic pH optimum. Four isozymes (HEX B, HEX I1, HEX I2, and HEX A) have been isolated from human placenta. HEX BA derived from the subunit rearrangement of HEX A was also prepared. To determine if the isozymes of hexosaminidase differ in their kinetic parameters, the conditions for 4-methylumbelliferyl-beta-D-N-acetylglucosaminide hydrolysis were optimized for each isozyme. The Km values were essentially similar and varied from 0.64 +/- 0.06 for HEX BA to 0.85 +/- 0.13 for HEX I1. The Vmax values were similar only for HEX I1 (3.90 +/- 0.28 kat kg-1) and HEX I2 (4.40 +/- 0.17). Vmax values varied significantly between HEX A (9.68 +/- 0.52), HEX B (8.00 +/- 0.75), HEX BA (4.81 +/- 0.17), and the HEX I values.

Electrophoresis, Polyacrylamide Gel↗

Dynamic evaluation of human vestibulo-ocular function using white noise rotation stimulus and linear system parameter estimation techniques.

White noise acceleration inputs were used to determine the human VOR transfer function both for normal subjects and for patients falling into two pilot categories: unilateral labyrinthectomy and reduced bilateral responses. The systematic patterns shown in the transfer function of the pilot abnormal categories as compared to the normal data suggests one method of classifying test results (table I). Frequency domain linear systems parameter fits were also made using the same data. The changes in these fit parameters, when pilot abnormal data is compared to normal data, suggests the use of the parameter fits themselves as a second classification scheme (fig. 1). The second scheme is not appropriate in cases where the response is unrelated to the stimulus.

Eye Movements↗