Parametric Signal Estimation Using the Cumulative Distribution Transform.
Ontology highlight
ABSTRACT:
SUBMITTER: Rubaiyat AHM
PROVIDER: S-EPMC7392180 | biostudies-literature | 2020
REPOSITORIES: biostudies-literature
Publications
Parametric Signal Estimation Using the Cumulative Distribution Transform.
Rubaiyat Abu Hasnat Mohammad AHM Hallam Kyla M KM Nichols Jonathan M JM Hutchinson Meredith N MN Li Shiying S Rohde Gustavo K GK
IEEE transactions on signal processing : a publication of the IEEE Signal Processing Society 20200525
We present a new method for estimating signal model parameters using the Cumulative Distribution Transform (CDT). Our approach minimizes the Wasserstein distance between measured and model signals. We derive some useful properties of the CDT and show that the resulting estimation problem, while nonlinear in the original signal domain, becomes a linear least squares problem in the transform domain. Furthermore, we discuss the properties of the estimator in the presence of noise and present a nove ...[more]
Similar Datasets
Project description:We present a new supervised image classification method applicable to a broad class of image deformation models. The method makes use of the previously described Radon Cumulative Distribution Transform (R-CDT) for image data, whose mathematical properties are exploited to express the image data in a form that is more suitable for machine learning. While certain operations such as translation, scaling, and higher-order transformations are challenging to model in native image space, we show the R-CDT can capture some of these variations and thus render the associated image classification problems easier to solve. The method - utilizing a nearest-subspace algorithm in the R-CDT space - is simple to implement, non-iterative, has no hyper-parameters to tune, is computationally efficient, label efficient, and provides competitive accuracies to state-of-the-art neural networks for many types of classification problems. In addition to the test accuracy performances, we show improvements (with respect to neural network-based methods) in terms of computational efficiency (it can be implemented without the use of GPUs), number of training samples needed for training, as well as out-of-distribution generalization. The Python code for reproducing our results is available at [1].
| S-EPMC9032314 | biostudies-literature
Project description:Separating the decay signal from diffusion-weighted scans into two or more components can be challenging. The phasor technique is well established in the field of optical microscopy for visualization and separation of fluorescent dyes with different lifetimes. The use of the phasor technique for separation of diffusion-weighted decay signals was recently proposed. In this study, we investigate the added value of this technique for fitting decay models and visualization of decay rates. Phasor visualization was performed in five glioblastoma patients. Using simulations, the influence of incorrect diffusivity values and of the number of b-values on fitting a three-component model with fixed diffusivities (dubbed "unmixing") was investigated for both a phasor-based fit and a linear least squares (LLS) fit. Phasor-based intravoxel incoherent motion (IVIM) fitting was compared with nonlinear least squares (NLLS) and segmented fitting (SF) methods in terms of accuracy and precision. The distributions of the parameter estimates of simulated data were compared with those obtained in a healthy volunteer. In the phasor visualizations of two glioblastoma patients, a cluster of points was observed that was not seen in healthy volunteers. The identified cluster roughly corresponded to the enhanced edge region of the tumor of two glioblastoma patients visible on fluid-attenuated inversion recovery (FLAIR) images. For fitting decay models the usefulness of the phasor transform is less pronounced, but the additional knowledge gained from the geometrical configuration of phasor space can aid fitting routines. This has led to slightly improved fitting results for the IVIM model: phasor-based fitting yielded parameter maps with higher precision than the NLLS and SF methods for parameters f and D (interquartile range [IQR] for f: NLLS 27, SF 12, phasor 5.7%; IQR for D: NLLS 0.28, SF 0.18, phasor 0.10 μm2 /s). For unmixing, LLS fitting slightly but consistently outperformed phasor-based fitting in all of the tested scenarios.
| S-EPMC7685171 | biostudies-literature
Project description:Non-parametric and semi-parametric resampling procedures are widely used to perform support estimation in computational biology and bioinformatics. Among the most widely used methods in this class is the standard bootstrap method, which consists of random sampling with replacement. While not requiring assumptions about any particular parametric model for resampling purposes, the bootstrap and related techniques assume that sites are independent and identically distributed (i.i.d.). The i.i.d. assumption can be an over-simplification for many problems in computational biology and bioinformatics. In particular, sequential dependence within biomolecular sequences is often an essential biological feature due to biochemical function, evolutionary processes such as recombination, and other factors. To relax the simplifying i.i.d. assumption, we propose a new non-parametric/semi-parametric sequential resampling technique that generalizes "Heads-or-Tails" mirrored inputs, a simple but clever technique due to Landan and Graur. The generalized procedure takes the form of random walks along either aligned or unaligned biomolecular sequences. We refer to our new method as the SERES (or "SEquential RESampling") method. To demonstrate the performance of the new technique, we apply SERES to estimate support for the multiple sequence alignment problem. Using simulated and empirical data, we show that SERES-based support estimation yields comparable or typically better performance compared to state-of-the-art methods.
| S-EPMC7164268 | biostudies-literature
Project description:Given the latent stratum membership, principal stratification models with continuous outcomes naturally fit in the parametric estimation framework of Gaussian mixtures. However, with models that are not nonparametrically identified, relying on parametric mixture modeling has been mostly discouraged as a way of identifying principal effects. This study revisits this rather deserted use of parametric mixture modeling, which may open up various possibilities in principal stratification modeling. The main problem with using the parametric mixture modeling approach is that it is hard to assess the quality of principal effect estimates given its reliance on parametric conditions. As a way of assessing the estimation quality in this situation, this study proposes that we use parametric mixture modeling in two different ways, with and without the assurance of nonparametric identification. The key identifying assumption employed in this study is the moving exclusion restriction, a flexible version of the standard exclusion restriction assumption. This assumption is used as a temporary vehicle to help assess the quality of principal effect estimates obtained relying on parametric mixture modeling. The study presents promising results, showing the possibility of using parametric mixture modeling as an accessible tool for causal inference.
| S-EPMC9232942 | biostudies-literature
Project description:Mounting evidence for the role of sleep spindles in neuroplasticity has led to an increased interest in these non-rapid eye movement (NREM) sleep oscillations. It has been hypothesized that fast and slow spindles might play a different role in memory processing. Here, we present a new sleep spindle detection algorithm utilizing a continuous wavelet transform (CWT) and individual adjustment of slow and fast spindle frequency ranges. Eighteen nap recordings of ten subjects were used for algorithm validation. Our method was compared with both a human scorer and a commercially available SIESTA spindle detector. For the validation set, mean agreement between our detector and human scorer measured during sleep stage 2 using kappa coefficient was 0.45, whereas mean agreement between our detector and SIESTA algorithm was 0.62. Our algorithm was also applied to sleep-related memory consolidation data previously analyzed with a SIESTA detector and confirmed previous findings of significant correlation between spindle density and declarative memory consolidation. We then applied our method to a study in monozygotic (MZ) and dizygotic (DZ) twins, examining the genetic component of slow and fast sleep spindle parameters. Our analysis revealed strong genetic influence on variance of all slow spindle parameters, weaker genetic effect on fast spindles, and no effects on fast spindle density and number during stage 2 sleep.
| S-EPMC4652604 | biostudies-literature
Project description:BackgroundT-wave alternans (TWA) has been associated with increased vulnerability to ventricular tachyarrhythmias and sudden cardiac death. However, both random (white) noise and (patho)physiologic processes (ie, premature ventricular contractions and heart and respiration rates) may hamper TWA estimation and therefore lessen its clinical utility for risk stratification.ObjectiveTo investigate the effect of random noise and certain (patho)physiologic processes on the estimation of TWA by using the fast Fourier transform method and to develop methods to overcome these potential sources of error.MethodsWe used a combination of human electrocardiogram data and computer simulations to assess the effects of a premature ventricular contraction and random and colored noise on the accuracy of TWA estimation.ResultsWe quantitatively demonstrate that replacing a "bad" beat with an odd/even median beat is a more accurate approach than replacing it with the overall average or the overall median beat. We also show that phase resetting may have a significant effect on alternans estimation and that estimation of alternans by using frequencies >0.4922 cycles/beat in a 128-point fast Fourier transform provides the most accurate approach for estimating the alternans when phase resetting is likely to occur. In addition, our data demonstrate that the number of indeterminate TWA tests due to high levels of noise can be reduced when the alternans voltage exceeds a new higher threshold. Furthermore, the amplitude of random noise has a significant effect on alternans estimation and the alternans voltage threshold should be adjusted for noise levels >1.8 μV. Finally, we quantitatively demonstrate that colored noise may lead to a false-positive or a false-negative result. We propose methods to estimate the effect of these (patho)physiologic processes on the alternans estimation in order to determine whether a TWA test is likely to be a true positive or a true negative.ConclusionThis study introduces novel methods to overcome potential sources of error in the estimation of TWA. These methods may improve the utility of TWA either for ambulatory monitoring or for clinical risk stratification for ventricular arrhythmias and sudden cardiac death.
| S-EPMC3288719 | biostudies-literature
Project description:Double electron-electron resonance (DEER) spectroscopy measures distance distributions between spin labels in proteins, yielding important structural and energetic information about conformational landscapes. Analysis of an experimental DEER signal in terms of a distance distribution is a nontrivial task due to the ill-posed nature of the underlying mathematical inversion problem. This work introduces a Bayesian probabilistic inference approach to analyze DEER data, using a multi-Gauss mixture model for the distance distribution. The method uses Markov chain Monte Carlo (MCMC) sampling to determine a posterior probability distribution over model parameter space. This distribution contains all the information available from the data, including a full quantification of the uncertainty about the parameters. The corresponding uncertainty about the distance distribution is captured via an ensemble of posterior predictive distributions. Several synthetic examples illustrate the method. An experimental example shows the importance of model checking and comparison using residual analysis and Bayes factors. Overall, the Bayesian approach allows for more robust inference about protein conformations from DEER spectroscopy.
| S-EPMC7846514 | biostudies-literature
Project description:We suggest a procedure for estimating Nth degree polynomial approximations to unknown (or known) probability density functions (PDFs) based on N statistical moments from each distribution. The procedure is based on the method of moments and is setup algorithmically to aid applicability and to ensure rigor in use. In order to show applicability, polynomial PDF approximations are obtained for the distribution families Normal, Log-Normal, Weibull as well as for a bimodal Weibull distribution and a data set of anonymized household electricity use. The results are compared with results for traditional PDF series expansion methods of Gram-Charlier type. It is concluded that this procedure is a comparatively simple procedure that could be used when traditional distribution families are not applicable or when polynomial expansions of probability distributions might be considered useful approximations. In particular this approach is practical for calculating convolutions of distributions, since such operations become integrals of polynomial expressions. Finally, in order to show an advanced applicability of the method, it is shown to be useful for approximating solutions to the Smoluchowski equation.
| S-EPMC5386244 | biostudies-literature
Project description:MotivationAlthough the majority of gene histories found in a clade of organisms are expected to be generated by a common process (e.g. the coalescent process), it is well known that numerous other coexisting processes (e.g. horizontal gene transfers, gene duplication and subsequent neofunctionalization) will cause some genes to exhibit a history distinct from those of the majority of genes. Such 'outlying' gene trees are considered to be biologically interesting, and identifying these genes has become an important problem in phylogenetics.ResultsWe propose and implement kdetrees, a non-parametric method for estimating distributions of phylogenetic trees, with the goal of identifying trees that are significantly different from the rest of the trees in the sample. Our method compares favorably with a similar recently published method, featuring an improvement of one polynomial order of computational complexity (to quadratic in the number of trees analyzed), with simulation studies suggesting only a small penalty to classification accuracy. Application of kdetrees to a set of Apicomplexa genes identified several unreliable sequence alignments that had escaped previous detection, as well as a gene independently reported as a possible case of horizontal gene transfer. We also analyze a set of Epichloë genes, fungi symbiotic with grasses, successfully identifying a contrived instance of paralogy.Availability and implementationOur method for estimating tree distributions and identifying outlying trees is implemented as the R package kdetrees and is available for download from CRAN.
| S-EPMC4176058 | biostudies-literature
Project description:It is generally difficult to define reasonable parameters and interpret their values in mathematical models of social phenomena. Rather than directly fitting abstract parameters against empirical data, we should define some concrete parameters to denote the sociocultural factors relevant for particular phenomena, and compute the values of these parameters based upon the corresponding empirical data. Taking the example of modeling studies of language competition, we propose a language diffusion principle and two language inheritance principles to compute two critical parameters, namely the impacts and inheritance rates of competing languages, in our language competition model derived from the Lotka-Volterra competition model in evolutionary biology. These principles assign explicit sociolinguistic meanings to those parameters and calculate their values from the relevant data of population censuses and language surveys. Using four examples of language competition, we illustrate that our language competition model with thus-estimated parameter values can reliably replicate and predict the dynamics of language competition, and it is especially useful in cases lacking direct competition data.
| S-EPMC3683775 | biostudies-literature