Multilevel statistical models for psychosomatic research.
Explore the source record for details and available documents.
SEARCH · Search PubMed
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.
Explore the source record for details and available documents.
We study an optimal nonparametric regression model for a threshold detector exposed to a noisy, subthreshold signal. The problem of recovering the signal is similar to that faced by neurons in nervous systems, although our model is intended to be normative rather than realistic. In our approach, the time-integrating activity of the neuron is modeled by kernel regression. Several aspects of the performance of the model are studied, including the existence of an optimal amount of noise (stochastic resonance). We construct a sequential, data-driven procedure for estimating the subthreshold signal. The performance of our model for threshold data is compared with kernel estimation for fully observed data. Finally, we discuss differences between our estimator and the best estimator for a constant signal.
We present a Bayesian dynamical inference method for characterizing cardiorespiratory (CR) dynamics in humans by inverse modeling from blood pressure time-series data. The technique is applicable to a broad range of stochastic dynamical models and can be implemented without severe computational demands. A simple nonlinear dynamical model is found that describes a measured blood pressure time series in the primary frequency band of the CR dynamics. The accuracy of the method is investigated using model-generated data with parameters close to the parameters inferred in the experiment. The connection of the inferred model to a well-known beat-to-beat model of the baroreflex is discussed.
Explore the source record for details and available documents.
This paper presents a framework for weighted fusion of several Active Shape and Active Appearance Models. The approach is based on the eigenspace fusion method proposed by Hall et al., which has been extended to fuse more than two weighted eigenspaces using unbiased mean and covariance matrix estimates. To evaluate the performance of fusion, a comparative assessment on segmentation precision as well as facial verification tests are performed using the AR, EQUINOX, and XM2VTS databases. Based on the results, it is concluded that the fusion is useful when the model needs to be updated online or when the original observations are absent.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Longitudinal data of 50-yr-old men, followed for 10 yr, were analyzed to identify variables which may be associated with healthy tooth surface loss. Three types of longitudinal models were used, arriving at similar conclusions. Predictors of sound tooth surface loss with aging included pocket depth, mobility, and recession.
Explore the source record for details and available documents.
In comprehensive functional genomics projects, systematic analysis of phenotypes is vital. However, conventional phenotypic screening is done mainly by imprecise visual observation of qualitative traits, and, therefore, in silico screening techniques for quantitative traits are required. In this report, we propose in silico phenotypic screening method that utilizes a Gaussian mixture model for the trait distribution in the offspring of a mutagenized line and the likelihood ratio test between the estimated Gaussian mixture model and the wild-type single Gaussian model. In order to evaluate the proposed method, we performed a screening experiment using real trait data of Arabidopsis. In this experiment, the proposed screening method properly distinguished the mutant line from the wild-type line. Furthermore, we conducted power analysis of the proposed method and two conventional methods under various simulated conditions of sample size and distribution of trait frequency. The result of the power analysis confirmed the effectiveness of the proposed method compared to the conventional methods.
The diurnal pattern of growth hormone (GH) serum levels depends on the frequency and amplitude of GH secretory events, the kinetics of GH infusion into and clearance from the circulation, and the feedback of GH on its secretion. We present a two-dimensional linear differential equation model based on these physiological principles to describe GH diurnal patterns. The model characterizes the onset times of the secretory events, the secretory event amplitudes, as well as the infusion, clearance, and feedback half-lives of GH. We illustrate the model by using maximum likelihood methods to fit it to GH measurements collected in 12 normal, healthy women during 8 h of scheduled sleep and a 16-h circadian constant-routine protocol. We assess the importance of the model components by using parameter standard error estimates and Akaike's Information Criterion. During sleep, both the median infusion and clearance half-life estimates were 13.8 min, and the median number of secretory events was 2. During the constant routine, the median infusion half-life estimate was 12.6 min, the median clearance half-life estimate was 11.7 min, and the median number of secretory events was 5. The infusion and clearance half-life estimates and the number of secretory events are consistent with current published reports. Our model gave an excellent fit to each GH data series. Our analysis paradigm suggests an approach to decomposing GH diurnal patterns that can be used to characterize the physiological properties of this hormone under normal and pathological conditions.
Neurohormone secretion is viewed here as a variable (unknown) admixture of basal and pulsatile release mechanisms, convolved with individually fitted biexponential elimination kinetics. This construct allows maximum-likelihood estimates of both (regulated and constitutive) components of hormone secretion. Thereby, we infer that a prolonged slow-component half-life of gonadotropin removal and amplified pulsatile (and total) daily luteinizing hormone (LH) secretion rates jointly explicate the postmenopausal elevation in serum LH concentrations without a necessary rise in basal LH secretion rates. This biomathematical formulation should be useful in exploring other neuroregulatory mechanisms that underlie single or dual alterations in the basal versus pulsatile modes of hormone secretion.
Studies assessing changes in maximal aerobic capacity (VO2 max) associated with aging have traditionally employed the ratio of VO2 max to body weight. Log-linear, ordinary least-squares, and weighted least-squares models may avoid some of the inherent weaknesses associated with the use of ratios. In this study we used four different methods to examine the age-associated decline in VO2 max in a cross-sectional sample of 276 healthy men, aged 45-80 yr. Sixty-one of the men were aerobically trained athletes, and the remainder were sedentary. The model that accounted for the largest proportion of variance was a weighted least-squares model that included age, fat-free mass, and an indicator variable denoting exercise training status. The model accounted for 66% of the variance in VO2 max and satisfied all the important general linear model assumptions. The other approaches failed to satisfy one or more of these assumptions. The results indicated that VO2 max declines at the same rate in athletic and sedentary men (0.24 l/min or 9%/decade) and that 35% of this decline (0.08 l . min-1 . decade-1) is due to the age-associated loss of fat-free mass.
When neuronal excitability is increased in area CA3 of the hippocampus in vitro, the pyramidal cells generate periodic bursts of action potentials that are synchronized across the network. We have previously provided evidence that synaptic depression at the excitatory recurrent collateral synapses in the CA3 network terminates each population burst so that the next burst cannot begin until these synapses have recovered. These findings raise the possibility that burst timing can be described in terms of the probability of recovery of this population of synapses. Here we demonstrate that when neuronal excitability is changed in the CA3 network, the mean and variance of the interburst interval change in a manner that is consistent with a timing mechanism comprised of a pool of exponentially relaxing pacemakers. The relaxation time constant of these pacemakers is the same as the time constant describing the recovery from activity-dependent depression of recurrent collateral synapses. Recovery was estimated from the rate of spontaneous transmitter release versus time elapsed since the last CA3 burst. Pharmacological and long-term alterations of synaptic strength and network excitability affected CA3 burst timing as predicted by the cumulative binomial distribution if the burst pace-maker consists of a pool of recovering recurrent synapses. These findings indicate that the recovery of a pool of synapses from burst-induced depression is a sufficient explanation for burst timing in the in vitro CA3 neuronal network. These findings also demonstrate how information regarding the nature of a pacemaker can be derived from the temporal pattern of synchronous network activity. This information could also be extracted from less accessible networks such as those generating interictal epileptiform discharges in vivo.
PURPOSE: To analyze the sensitivity of various sensory tests adjusted for glaucomatous optic disc damage. METHODS: In a cross-sectional study, the results of testing of 196 control subjects (age range, 18-69 years) and 308 patients with chronic open-angle glaucoma (age range, 18-70 years) were included. The perimetric mean defect (MD), a temporal contrast sensitivity test (TCS), a spatiotemporal contrast sensitivity test (STCS), the peak time of a blue-on-yellow visual evoked potential (BYVEP), and the amplitude of a pattern-reversal electroretinogram (PERG) were evaluated by a specific logistic regression model. This model included glaucomatous damage, quantified by neuroretinal rim area corrected for disc size, as a covariate of sensitivity. RESULTS: Sensitivity of diagnostic tests increased for all procedures with increasing loss of neuroretinal rim area. With progressing optic disc damage, MD and STCS showed higher sensitivity than did TCS. BYVEP showed a higher sensitivity than PERG in all disease stages. In general, the psychophysical tests were more sensitive than the electrophysiological ones. CONCLUSIONS: The specific model used in this study was an appropriate tool to analyze the sensitivity of several sensory glaucoma tests in relation to disease stage. Moreover, tests that were more sensitive in early disease stages (TCS) and others that were more sensitive in more advanced stages (MD, STCS) were identified.
In order to investigate the relationship between the pneumatization of temporal bones and the occurrence of otitis media in Greenlandic Inuit, 36 Greenlandic Inuit were examined by radiography of the temporal bones. The pneumatized cell area was measured planimetrically. All subjects answered a questionnaire on infectious middle ear disease (IMED), and an objective otologic examination was performed. Nine persons of 34 (26%) reported IMED in childhood, and there was IMED reported in all pneumatized areas below 400 mm2. Based on bilateral area measures, a polychotomous logistic regression model was applied. The occurrence of IMED was shown to be associated with smaller areas, and unilateral IMED was associated with pronounced asymmetry. The model has enabled risk assessment, and 8 persons were designated by the model as having had IMED with a sensitivity of 67% (confidence interval 30% to 93%) and a specificity of 92% (confidence interval 74% to 99%). When the model was applied to a historical anthropological Inuit material from 1700 to 1800 AD, 6 of 56 crania were designated as having had IMED. This method has improved the accuracy of estimating the occurrence of IMED in ancient populations.
Experiments offer more reliable evidence on causation than observational studies, which is not to gainsay the contribution to knowledge from observation. Experiments should be analyzed as experiments, not as observational studies. A simple comparison of rates might be just the right tool, with little value added by "sophisticated" models. This article discusses current models for causation, as applied to experimental and observational data. The intention-to-treat principle and the effect of treatment on the treated will also be discussed. Flaws in per-protocol and treatment-received estimates will be demonstrated.
BACKGROUND: Marine turtles deposit their eggs on tropical or subtropical beaches during discrete nesting seasons that span several months. The number and distribution of nests laid during a nesting season provide vital information on various aspects of marine turtle ecology and conservation. RESULTS: In the case of leatherback sea turtles nesting in French Guiana, we developed a mathematical model to explore the phenology of their nesting season, derived from an incomplete nest count dataset. We detected 3 primary components in the nest distribution of leatherbacks: an overall shape that corresponds to the arrival and departure of leatherback females in the Guianas region, a sinusoidal pattern with a period of approximately 10 days that is related to physiological constraints of nesting female leatherbacks, and a sinusoidal pattern with a period of approximately 15 days that likely reflects the influence of spring high tides on nesting female turtles. CONCLUSION: The model proposed here offers a variety of uses for both marine turtles and also other taxa when individuals are observed in a particular location for only part of the year.