RESUMEN

The under-5 mortality rate (U5MR), a critical health indicator, is typically estimated from household surveys in lower and middle income countries. Spatio-temporal disaggregation of household survey data can lead to highly variable estimates of U5MR, necessitating the usage of smoothing models which borrow information across space and time. The assumptions of common smoothing models may be unrealistic when certain time periods or regions are expected to have shocks in mortality relative to their neighbors, which can lead to oversmoothing of U5MR estimates. In this paper, we develop a spatial and temporal smoothing approach based on Gaussian Markov random field models which incorporate knowledge of these expected shocks in mortality. We demonstrate the potential for these models to improve upon alternatives not incorporating knowledge of expected shocks in a simulation study. We apply these models to estimate U5MR in Rwanda at the national level from 1985 to 2019, a time period which includes the Rwandan civil war and genocide.

RESUMEN

We introduce a model that can be used for the description of the distribution of species when there is scarcity of data, based on our previous work (Ballesteros et al. J Math Biol 85(4):31, 2022). We address challenges in modeling species that are seldom observed in nature, for example species included in The International Union for Conservation of Nature's Red List of Threatened Species (IUCN 2023). We introduce a general method and test it using a case study of a near threatened species of amphibians called Plectrohyla Guatemalensis (see IUCN 2023) in a region of the UNESCO natural reserve "Tacaná Volcano", in the border between Mexico and Guatemala. Since threatened species are difficult to find in nature, collected data can be extremely reduced. This produces a mathematical problem in the sense that the usual modeling in terms of Markov random fields representing individuals associated to locations in a grid generates artificial clusters around the observations, which are unreasonable. We propose a different approach in which our random variables describe yearly averages of expectation values of the number of individuals instead of individuals (and they take values on a compact interval). Our approach takes advantage of intuitive insights from environmental properties: in nature individuals are attracted or repulsed by specific features (Ballesteros et al. J Math Biol 85(4):31, 2022). Drawing inspiration from quantum mechanics, we incorporate quantum Hamiltonians into classical statistical mechanics (i.e. Gibbs measures or Markov random fields). The equilibrium between spreading and attractive/repulsive forces governs the behavior of the species, expressed through a global control problem involving an energy operator.

Asunto(s)

RESUMEN

Gaussian graphical models are widely used to study the dependence structure among variables. When samples are obtained from multiple conditions or populations, joint analysis of multiple graphical models are desired due to their capacity to borrow strength across populations. Nonetheless, existing methods often overlook the varying levels of similarity between populations, leading to unsatisfactory results. Moreover, in many applications, learning the population-level clustering structure itself is of particular interest. In this article, we develop a novel method, called Simultaneous Clustering and Estimation of Networks via Tensor decomposition (SCENT), that simultaneously clusters and estimates graphical models from multiple populations. Precision matrices from different populations are uniquely organized as a three-way tensor array, and a low-rank sparse model is proposed for joint population clustering and network estimation. We develop a penalized likelihood method and an augmented Lagrangian algorithm for model fitting. We also establish the clustering accuracy and norm consistency of the estimated precision matrices. We demonstrate the efficacy of the proposed method with comprehensive simulation studies. The application to the Genotype-Tissue Expression multi-tissue gene expression data provides important insights into tissue clustering and gene coexpression patterns in multiple brain tissues.

RESUMEN

Network psychometrics uses graphical models to assess the network structure of psychological variables. An important task in their analysis is determining which variables are unrelated in the network, i.e., are independent given the rest of the network variables. This conditional independence structure is a gateway to understanding the causal structure underlying psychological processes. Thus, it is crucial to have an appropriate method for evaluating conditional independence and dependence hypotheses. Bayesian approaches to testing such hypotheses allow researchers to differentiate between absence of evidence and evidence of absence of connections (edges) between pairs of variables in a network. Three Bayesian approaches to assessing conditional independence have been proposed in the network psychometrics literature. We believe that their theoretical foundations are not widely known, and therefore we provide a conceptual review of the proposed methods and highlight their strengths and limitations through a simulation study. We also illustrate the methods using an empirical example with data on Dark Triad Personality. Finally, we provide recommendations on how to choose the optimal method and discuss the current gaps in the literature on this important topic.

RESUMEN

Background and objective: Pelvic bone tumors represent a harmful orthopedic condition, encompassing both benign and malignant forms. Addressing the issue of limited accuracy in current machine learning algorithms for bone tumor image segmentation, we have developed an enhanced bone tumor image segmentation algorithm. This algorithm is built upon an improved full convolutional neural network, incorporating both the fully convolutional neural network (FCNN-4s) and a conditional random field (CRF) to achieve more precise segmentation. Methodology: The enhanced fully convolutional neural network (FCNN-4s) was employed to conduct initial segmentation on preprocessed images. Following each convolutional layer, batch normalization layers were introduced to expedite network training convergence and enhance the accuracy of the trained model. Subsequently, a fully connected conditional random field (CRF) was integrated to fine-tune the segmentation results, refining the boundaries of pelvic bone tumors and achieving high-quality segmentation. Results: The experimental outcomes demonstrate a significant enhancement in segmentation accuracy and stability when compared to the conventional convolutional neural network bone tumor image segmentation algorithm. The algorithm achieves an average Dice coefficient of 93.31 %, indicating superior performance in real-time operations. Conclusion: In contrast to the conventional convolutional neural network segmentation algorithm, the algorithm presented in this paper boasts a more intricate structure, proficiently addressing issues of over-segmentation and under-segmentation in pelvic bone tumor segmentation. This segmentation model exhibits superior real-time performance, robust stability, and is capable of achieving heightened segmentation accuracy.

RESUMEN

Since 3D sensors became popular, imaged depth data are easier to obtain in the consumer sector. In applications such as defect localization on industrial objects or mass/volume estimation, precise depth data is important and, thus, benefits from the usage of multiple information sources. However, a combination of RGB images and depth images can not only improve our understanding of objects, capacitating one to gain more information about objects but also enhance data quality. Combining different camera systems using data fusion can enable higher quality data since disadvantages can be compensated. Data fusion itself consists of data preparation and data registration. A challenge in data fusion is the different resolutions of sensors. Therefore, up- and downsampling algorithms are needed. This paper compares multiple up- and downsampling methods, such as different direct interpolation methods, joint bilateral upsampling (JBU), and Markov random fields (MRFs), in terms of their potential to create RGB-D images and improve the quality of depth information. In contrast to the literature in which imaging systems are adjusted to acquire the data of the same section simultaneously, the laboratory setup in this study was based on conveyor-based optical sorting processes, and therefore, the data were acquired at different time periods and different spatial locations. Data assignment and data cropping were necessary. In order to evaluate the results, root mean square error (RMSE), signal-to-noise ratio (SNR), correlation (CORR), universal quality index (UQI), and the contour offset are monitored. With JBU outperforming the other upsampling methods, achieving a meanRMSE = 25.22, mean SNR = 32.80, mean CORR = 0.99, and mean UQI = 0.97.

RESUMEN

Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) can be used as a non-invasive method for the assessment of myocardial perfusion. The acquired images can be utilised to analyse the spatial extent and severity of myocardial ischaemia (regions with impaired microvascular blood flow). In the present paper, we propose a novel generalisable spatio-temporal hierarchical Bayesian model (GST-HBM) to automate the detection of ischaemic lesions and improve the in silico prediction accuracy by systematically integrating spatio-temporal context information. We present a computational inference procedure with an adequate trade-off between accuracy and computational efficiency, whereby model parameters are sampled from the posterior distribution with Gibbs sampling, while lower-level hyperparameters are selected using model selection strategies based on the Watanabe Akaike information criterion (WAIC). We have assessed our method on both synthetic (in silico) data with known gold-standard and 12 sets of clinical first-pass myocardial perfusion DCE-MRI datasets. We have also carried out a comparative performance evaluation with four established alternative methods: Gaussian mixture model (GMM), opening and closing operations based on Gaussian mixture model (GMMC&Omax), Markov random field constrained Gaussian mixture model (GMM-MRF) and model-based hierarchical Bayesian model (M-HBM). Our results show that the proposed GST-HBM method achieves much higher in silico prediction accuracy than the established alternative methods. Furthermore, this method appears to provide a more robust delineation of ischaemic lesions in datasets affected by spatially variant noise.

Asunto(s)

RESUMEN

Let Z(t)=exp2BH(t)-t2H,t∈R with BH(t),t∈R a standard fractional Brownian motion (fBm) with Hurst parameter H∈(0,1] and define for x non-negative the Berman function BZ(x)=EI{ϵ0(RZ)>x}ϵ0(RZ)∈(0,∞),where the random variable R independent of Z has survival function 1/x,x⩾1 and ϵ0(RZ)=∫RIRZ(t)>1dt.In this paper we consider a general random field (rf) Z that is a spectral rf of some stationary max-stable rf X and derive the properties of the corresponding Berman functions. In particular, we show that Berman functions can be approximated by the corresponding discrete ones and derive interesting representations of those functions which are of interest for Monte Carlo simulations presented in this article.

RESUMEN

BACKGROUND: In the healthcare domain today, despite the substantial adoption of electronic health information systems, a significant proportion of medical reports still exist in paper-based formats. As a result, there is a significant demand for the digitization of information from these paper-based reports. However, the digitization of paper-based laboratory reports into a structured data format can be challenging due to their non-standard layouts, which includes various data types such as text, numeric values, reference ranges, and units. Therefore, it is crucial to develop a highly scalable and lightweight technique that can effectively identify and extract information from laboratory test reports and convert them into a structured data format for downstream tasks. METHODS: We developed an end-to-end Natural Language Processing (NLP)-based pipeline for extracting information from paper-based laboratory test reports. Our pipeline consists of two main modules: an optical character recognition (OCR) module and an information extraction (IE) module. The OCR module is applied to locate and identify text from scanned laboratory test reports using state-of-the-art OCR algorithms. The IE module is then used to extract meaningful information from the OCR results to form digitalized tables of the test reports. The IE module consists of five sub-modules, which are time detection, headline position, line normalization, Named Entity Recognition (NER) with a Conditional Random Fields (CRF)-based method, and step detection for multi-column. Finally, we evaluated the performance of the proposed pipeline on 153 laboratory test reports collected from Peking University First Hospital (PKU1). RESULTS: In the OCR module, we evaluate the accuracy of text detection and recognition results at three different levels and achieved an averaged accuracy of 0.93. In the IE module, we extracted four laboratory test entities, including test item name, test result, test unit, and reference value range. The overall F1 score is 0.86 on the 153 laboratory test reports collected from PKU1. With a single CPU, the average inference time of each report is only 0.78 s. CONCLUSION: In this study, we developed a practical lightweight pipeline to digitalize and extract information from paper-based laboratory test reports in diverse types and with different layouts that can be adopted in real clinical environments with the lowest possible computing resources requirements. The high evaluation performance on the real-world hospital dataset validated the feasibility of the proposed pipeline.

Asunto(s)

RESUMEN

Many applications in science and engineering involve data defined at specific geospatial locations, which are often modeled as random fields. The modeling of a proper correlation function is essential for the probabilistic calibration of the random fields, but traditional methods were developed with the assumption to have observations with evenly spaced data. Available methods dealing with irregularly spaced data generally require either interpolation or computationally expensive solutions. Instead, we propose a simple approach based on least square regression to estimate the autocorrelation function. We first tested our methodology on an artificially produced dataset to assess the performance of our method. The accuracy of the method and its robustness to the level of noise in the data indicate that it is suitable for use in realistic problems. In addition, the methodology was used on a major application, the modeling of animal species connected with zoonotic diseases. Understanding the population dynamics of reservoirs of zoonotic diseases, such as bats, is a crucial first step to predict and prevent potential spillover of deadly viruses like Ebola. Due to the limited data on bats across Africa, their density and migrations can only be studied with probabilistic numerical models based on samples of the ecological bare carrying capacity (K0). For this purpose, the bare carrying capacity was modeled as a random field and its statistics calibrated with the available data. The bare carrying capacity of bats was found to be denser in central Africa. This is because climatic and environmental conditions are more suitable for the survival of bats. The proposed methodology for random field calibration was shown to be a promising approach, which can cope with large gaps in data and with complex applications involving large geographical areas and high resolution.

RESUMEN

Background and objective: Bone tumor is a kind of harmful orthopedic disease, there are benign and malignant points. Aiming at the problem that the accuracy of the existing machine learning algorithm for bone tumor image segmentation is not high, a bone tumor image segmentation algorithm based on improved full convolutional neural network which consists fully convolutional neural network (FCNN-4s) and conditional random field (CRF). Methodology: The improved fully convolutional neural network (FCNN-4s) was used to perform coarse segmentation on preprocessed images. Batch normalization layers were added after each convolutional layer to accelerate the convergence speed of network training and improve the accuracy of the trained model. Then, a fully connected conditional random field (CRF) was fused to refine the bone tumor boundary in the coarse segmentation results, achieving the fine segmentation effect. Results: The experimental results show that compared with the traditional convolutional neural network bone tumor image segmentation algorithm, the algorithm has a great improvement in segmentation accuracy and stability, the average Dice can reach 91.56%, the real-time performance is better. Conclusion: Compared with the traditional convolutional neural network segmentation algorithm, the algorithm in this paper has a more refined structure, which can effectively solve the problem of over-segmentation and under-segmentation of bone tumors. The segmentation prediction has better real-time performance, strong stability, and can achieve higher segmentation accuracy.

RESUMEN

To date, several POS taggers have been introduced to facilitate the success of semantic analysis for different languages. However, the task of POS tagging becomes a bit intricate in morphologically complex languages, like Amharic. In this paper, we evaluated different models such as bidirectional long short term memory, convolutional neural network in combination with bidirectional long short term memory, and conditional random field for Amharic POS tagging. Various features, both language-dependent and -independent, have been explored in a conditional random field model. Besides, word-level and character-level features are analyzed in deep neural network models. A convolutional neural network is utilized for encoding features at the word and character level. Each model's performance has evaluated on the dataset that contained 321 K tokens and manually tagged with 31 POS tags. Lastly, the best performance obtained by an end-to-end deep neural network model, convolutional neural network in combination with bidirectional long term short memory and conditional random field, is 97.23% accuracy. This is the highest accuracy for Amharic POS tagging task and is competent with contemporary taggers currently existing in different languages.

RESUMEN

Texture segmentation plays a crucial role in the domain of image analysis and its recognition. Noise is inextricably linked to images, just like it is with every signal received by sensing, which has an impact on how well the segmentation process performs in general. Recent literature reveals that the research community has started recognizing the domain of noisy texture segmentation for its work towards solutions for the automated quality inspection of objects, decision support for biomedical images, facial expressions identification, retrieving image data from a huge dataset and many others. Motivated by the latest work on noisy textures, during our work being presented here, Brodatz and Prague texture images are contaminated with Gaussian and salt-n-pepper noise. A three-phase approach is developed for the segmentation of textures contaminated by noise. In the first phase, these contaminated images are restored using techniques with excellent performance as per the recent literature. In the remaining two phases, segmentation of the restored textures is carried out by a novel technique developed using Markov Random Fields (MRF) and objective customization of the Median Filter based on segmentation performance metrics. When the proposed approach is evaluated on Brodatz textures, an improvement of up to 16% segmentation accuracy for salt-n-pepper noise with 70% noise density and 15.1% accuracy for Gaussian noise (with a variance of 50) has been made in comparison with the benchmark approaches. On Prague textures, accuracy is improved by 4.08% for Gaussian noise (with variance 10) and by 2.47% for salt-n-pepper noise with 20% noise density. The approach in the present study can be applied to a diversified class of image analysis applications spanning a wide spectrum such as satellite images, medical images, industrial inspection, geo-informatics, etc.

Asunto(s)

RESUMEN

Coronary artery disease is one of the most common types of cardiovascular disease. Death from coronary heart disease is influenced by genetic factors in both women and men. In this article, we propose a novel Bayesian variable selection framework for the identification of important genetic variants associated with coronary artery disease disease status. Instead of treating each feature independently as in conventional Bayesian variable selection methods, we propose an innovative prior for the inclusion probabilities of genetic variants that accounts for their ordering structure. We assume that neighboring variants are more likely to be selected together as they tend to be highly correlated and have similar biological functions. Additionally, we propose to group participating subjects based on underlying population structure and fit separate regressions, so that the regression coefficients can better reflect different disease risks in different population groups. Our approach borrows strength across regression models through an innovative prior inspired by the Markov random fields. The proposed framework can improve variable selection and prediction performances as demonstrated in the simulation studies. We also apply the proposed framework to the CATHeterization GENetics data with binary Coronary artery disease disease status.

Asunto(s)

RESUMEN

Viruses are a vastly underestimated component of biodiversity that occur as diverse communities across hierarchical scales from the landscape level to individual hosts. The integration of community ecology with disease biology is a powerful, novel approach that can yield unprecedented insights into the abiotic and biotic drivers of pathogen community assembly. Here, we sampled wild plant populations to characterize and analyze the diversity and co-occurrence structure of within-host virus communities and their predictors. Our results show that these virus communities are characterized by diverse, non-random coinfections. Using a novel graphical network modeling framework, we demonstrate how environmental heterogeneity influences the network of virus taxa and how the virus co-occurrence patterns can be attributed to non-random, direct statistical virus-virus associations. Moreover, we show that environmental heterogeneity changed virus association networks, especially through their indirect effects. Our results highlight a previously underestimated mechanism of how environmental variability can influence disease risks by changing associations between viruses that are conditional on their environment.

Asunto(s)

RESUMEN

With advances in modern worlds technology, huge datasets that show dependencies in space as well as in time occur frequently in practice. As an example, several monitoring stations at different geographical locations track hourly concentration measurements of a number of air pollutants for several years. Such a dataset contains thousands of multivariate observations, thus, proper statistical analysis needs to account for dependencies in space and time between and among the different monitored variables. To simplify the consequent multivariate spatio-temporal statistical analysis it might be of interest to detect linear transformations of the original observations that result in straightforward interpretative, spatio-temporally uncorrelated processes that are also highly likely to have a real physical meaning. Blind source separation (BSS) represents a statistical methodology which has the aim to recover so-called latent processes, that exactly meet the former requirements. BSS was already successfully used in sole temporal and sole spatial applications with great success, but, it was not yet introduced for the spatio-temporal case. In this contribution, a reasonable and innovative generalization of BSS for multivariate space-time random fields (stBSS), under second-order stationarity, is proposed, together with two space-time extensions of the well-known algorithms for multiple unknown signals extraction (stAMUSE) and the second-order blind identification (stSOBI) which solve the formulated problem. Furthermore, symmetry and separability properties of the model are elaborated and connections to the space-time linear model of coregionalization and to the classical principal component analysis are drawn. Finally, the usefulness of the new methods is shown in a thorough simulation study and on a real environmental application.

RESUMEN

This paper introduces a new modeling framework for the statistical analysis of point patterns on a manifold Md, defined by a connected and compact two-point homogeneous space, including the special case of the sphere. The presented approach is based on temporal Cox processes driven by a L2(Md)-valued log-intensity. Different aggregation schemes on the manifold of the spatiotemporal point-referenced data are implemented in terms of the time-varying discrete Jacobi polynomial transform of the log-risk process. The n-dimensional microscale point pattern evolution in time at different manifold spatial scales is then characterized from such a transform. The simulation study undertaken illustrates the construction of spherical point process models displaying aggregation at low Legendre polynomial transform frequencies (large scale), while regularity is observed at high frequencies (small scale). K-function analysis supports these results under temporal short, intermediate and long range dependence of the log-risk process.

RESUMEN

The functionality enhancement of ferroelectrics by local polar clusters called polar nanoregions (PNRs) is one of the current interests in materials science. KTa1-xNbxO3 (KTN) with perovskite structure is a well-known electro-optic crystal with a large Kerr effect. The existence of PNRs in relaxor-like ferroelectric Nb-rich KTN with homovalent B-site cations is controversial. This paper reviews recent progress in understanding precursor dynamics in Nb-rich KTN crystals studied using Brillouin scattering. The intense central peak (CP) and significant softening of sound velocity are observed above the Curie temperature (TC) due to the polarization fluctuations in PNRs. The effects of Li-doping, defects, and electric fields on the growth and/or creation of PNRs are found using changes in acoustic properties. The electric-field-induced TC, which is shifted to higher values with increases in applied voltage, including critical endpoint (CEP) and field gradient by trapped electrons, are discussed as well. This new knowledge may give new insight into advanced functionality in perovskite ferroelectrics.

RESUMEN

Mosquito vectors of eastern equine encephalitis virus (EEEV) and West Nile virus (WNV) in the USA reside within broad multi-species assemblages that vary in spatial and temporal composition, relative abundances and vector competence. These variations impact the risk of pathogen transmission and the operational management of these species by local public health vector control districts. However, most models of mosquito vector dynamics focus on single species and do not account for co-occurrence probabilities between mosquito species pairs across environmental gradients. In this investigation, we use for the first time conditional Markov Random Fields (CRF) to evaluate spatial co-occurrence patterns between host-seeking mosquito vectors of EEEV and WNV around sampling sites in Manatee County, Florida. Specifically, we aimed to: (i) quantify correlations between mosquito vector species and other mosquito species; (ii) quantify correlations between mosquito vectors and landscape and climate variables; and (iii) investigate whether the strength of correlations between species pairs are conditional on landscape or climate variables. We hypothesized that either mosquito species pairs co-occur in patterns driven by the landscape and/or climate variables, or these vector species pairs are unconditionally dependent on each other regardless of the environmental variables. Our results indicated that landscape and bioclimatic covariates did not substantially improve the overall model performance and that the log abundances of the majority of WNV and EEEV vector species were positively dependent on other vector and non-vector mosquito species, unconditionally. Only five individual mosquito vectors were weakly dependent on environmental variables with one exception, Culiseta melanura, the primary vector for EEEV, which showed a strong correlation with woody wetland, precipitation seasonality and average temperature of driest quarter. Our analyses showed that majority of the studied mosquito species' abundance and distribution are insignificantly better predicted by the biotic correlations than by environmental variables. Additionally, these mosquito vector species may be habitat generalists, as indicated by the unconditional correlation matrices between species pairs, which could have confounded our analysis, but also indicated that the approach could be operationalized to leverage species co-occurrences as indicators of vector abundances in unsampled areas, or under scenarios where environmental variables are not informative.

Asunto(s)

RESUMEN

We revisit several conditionally formulated Gaussian Markov random fields, known as the intrinsic conditional autoregressive model, the proper conditional autoregressive model, and the Leroux et al. conditional autoregressive model, as well as convolution models such as the well known Besag, York and Mollie model, its (adaptive) re-parameterization, and its scaled alternatives, for their roles of modelling underlying spatial risks in Bayesian disease mapping. Analytic and simulation studies, with graphic visualizations, and disease mapping case studies, present insights and critique on these models for their nature and capacities in characterizing spatial dependencies, local influences, and spatial covariance and correlation functions, and in facilitating stabilized and efficient posterior risk prediction and inference. It is illustrated that these models are Gaussian (Markov) random fields of different spatial dependence, local influence, and (covariance) correlation functions and can play different and complementary roles in Bayesian disease mapping applications.

Asunto(s)