Search

Background: The minimum variance (MV) beamformer can significantly improve the image resolution in ultrasound imaging, but it has limited performance in noise reduction. We recently proposed the covariance matrix-based statistical beamforming (CMSB) for medical ultrasound imaging to reduce sidelobes and incoherent clutter.Methods: In this paper, we aim to improve the imaging performance of the MV beamformer by introducing a new pixel-based adaptive weighting approach based on CMSB, which is named as covariance matrix-based adaptive weighting (CMSAW). The proposed CMSAW estimates the mean-to-standard-deviation ratio (MSR) of a modified covariance matrix reconstructed by adaptive spatial smoothing, rotary averaging, and diagonal reducing. Moreover, adaptive diagonal reducing based on the aperture coherence is introduced in CMSAW to enhance the performance in speckle preservation.Results: The proposed CMSAW-weighted MV (CMSAW-MV) was validated through simulation, phantom experiments, and in vivo studies. The phantom experimental results show that CMSAW-MV obtains resolution improvement of 21.3% and simultaneously achieves average improvements of 96.4% and 71.8% in average contrast and generalized contrast-to-noise ratio (gCNR) for anechoic cyst, respectively, compared with MV. in vivo studies indicate that CMSAW-MV improves the noise reduction performance of MV beamformer.Conclusion: Simulation, experimental, and in vivo results all show that CMSAW-MV can improve resolution and suppress sidelobes and incoherent clutter and noise. These results demonstrate the effectiveness of CMSAW in improving the imaging performance of MV beamformer. Moreover, the proposed CMSAW with a computational complexity of [Formula: see text] has the potential to be implemented in real time using the graphics processing unit. [ABSTRACT FROM AUTHOR]

Background: With the invention of fitness trackers, it has been possible to continuously monitor a user's biometric data such as heart rates, number of footsteps taken, and amount of calories burned. This paper names the time series of these three types of biometric data, the user's "activeness", and investigates the feasibility in modeling and predicting the long-term activeness of the user.Methods: The dataset used in this study consisted of several months of biometric time-series data gathered by seven users independently. Four recurrent neural network (RNN) architectures-as well as a deep neural network and a simple regression model-were proposed to investigate the performance on predicting the activeness of the user under various length-related hyper-parameter settings. In addition, the learned model was tested to predict the time period when the user's activeness falls below a certain threshold.Results: A preliminary experimental result shows that each type of activeness data exhibited a short-term autocorrelation; and among the three types of data, the consumed calories and the number of footsteps were positively correlated, while the heart rate data showed almost no correlation with neither of them. It is probably due to this characteristic of the dataset that although the RNN models produced the best results on modeling the user's activeness, the difference was marginal; and other baseline models, especially the linear regression model, performed quite admirably as well. Further experimental results show that it is feasible to predict a user's future activeness with precision, for example, a trained RNN model could predict-with the precision of 84%-when the user would be less active within the next hour given the latest 15 min of his activeness data.Conclusions: This paper defines and investigates the notion of a user's "activeness", and shows that forecasting the long-term activeness of the user is indeed possible. Such information can be utilized by a health-related application to proactively recommend suitable events or services to the user. [ABSTRACT FROM AUTHOR]

Background: Gene networks in living cells can change depending on various conditions such as caused by different environments, tissue types, disease states, and development stages. Identifying the differential changes in gene networks is very important to understand molecular basis of various biological process. While existing algorithms can be used to infer two gene networks separately from gene expression data under two different conditions, and then to identify network changes, such an approach does not exploit the similarity between two gene networks, and it is thus suboptimal. A desirable approach would be clearly to infer two gene networks jointly, which can yield improved estimates of network changes. Results: In this paper, we developed a proximal gradient algorithm for differential network (ProGAdNet) inference, that jointly infers two gene networks under different conditions and then identifies changes in the network structure. Computer simulations demonstrated that our ProGAdNet outperformed existing algorithms in terms of inference accuracy, and was much faster than a similar approach for joint inference of gene networks. Gene expression data of breast tumors and normal tissues in the TCGA database were analyzed with our ProGAdNet, and revealed that 268 genes were involved in the changed network edges. Gene set enrichment analysis identified a significant number of gene sets related to breast cancer or other types of cancer that are enriched in this set of 268 genes. Network analysis of the kidney cancer data in the TCGA database with ProGAdNet also identified a set of genes involved in network changes, and the majority of the top genes identified have been reported in the literature to be implicated in kidney cancer. These results corroborated that the gene sets identified by ProGAdNet were very informative about the cancer disease status. A software package implementing the ProGAdNet, computer simulations, and real data analysis is available as Additional file 1. Conclusion: With its superior performance over existing algorithms, ProGAdNet provides a valuable tool for finding changes in gene networks, which may aid the discovery of gene-gene interactions changed under different conditions. [ABSTRACT FROM AUTHOR]

Background: Ultrasound computed tomography (USCT) is usually realized with a ring array. It can provide better imaging performance and more tissue information by emitting and receiving the ultrasound signal in different directions simultaneously. However, USCT imaging is usually applied with the synthetic aperture (SA) emission method, which leads to a long scanning time with a large number of elements on the ring array. The echo image can provide the structural information, and has a higher resolution than maps of other parameters in USCT. Hence, we proposed plane wave (PW) imaging for ring array to acquire the echo wave and reduce the scanning time considerably.Results: In this paper, an emitting and receiving process was proposed to realize plane wave imaging with a ring array. With the proposed scanning method, the number of emission events can be reduced greatly. A beamforming method based on the eigenspace-based minimum variance (ESBMV) was also combined with the scanning method. With ESBMV beamformer, the resolution and contrast ratio of reconstruction result can be maintained or even improved under a fewer-emissions condition. We validated the method using both computer simulations with Field II and phantom experiments with a ring array of 512 elements. The Verasonics® system was used to transmit and receive the ultrasound signal in the phantom experiments.Conclusions: According to the results of the experiments, the imaging results will have a better contrast ratio with a higher emitting energy. Additionally, the scanning time with the proposed method can be only one-tenth of that with the SA emission method, while the echo imaging performance still remains at a similar level or even better. [ABSTRACT FROM AUTHOR]

Background: The sample size required to power a study to a nominal level in a paired comparative diagnostic accuracy study, i.e. studies in which the diagnostic accuracy of two testing procedures is compared relative to a gold standard, depends on the conditional dependence between the two tests - the lower the dependence the greater the sample size required. A priori, we usually do not know the dependence between the two tests and thus cannot determine the exact sample size required. One option is to use the implied sample size for the maximal negative dependence, giving the largest possible sample size. However, this is potentially wasteful of resources and unnecessarily burdensome on study participants as the study is likely to be overpowered. A more accurate estimate of the sample size can be determined at a planned interim analysis point where the sample size is re-estimated.Methods: This paper discusses a sample size estimation and re-estimation method based on the maximum likelihood estimates, under an implied multinomial model, of the observed values of conditional dependence between the two tests and, if required, prevalence, at a planned interim. The method is illustrated by comparing the accuracy of two procedures for the detection of pancreatic cancer, one procedure using the standard battery of tests, and the other using the standard battery with the addition of a PET/CT scan all relative to the gold standard of a cell biopsy. Simulation of the proposed method illustrates its robustness under various conditions.Results: The results show that the type I error rate of the overall experiment is stable using our suggested method and that the type II error rate is close to or above nominal. Furthermore, the instances in which the type II error rate is above nominal are in the situations where the lowest sample size is required, meaning a lower impact on the actual number of participants recruited.Conclusion: We recommend multinomial model maximum likelihood estimation of the conditional dependence between paired diagnostic accuracy tests at an interim to reduce the number of participants required to power the study to at least the nominal level.Trial Registration: ISRCTN ISRCTN73852054 . Registered 9th of January 2015. Retrospectively registered. [ABSTRACT FROM AUTHOR]

Background: Bivariate random-effects models represent a widely accepted and recommended approach for meta-analysis of test accuracy studies. Standard likelihood methods routinely used for inference are prone to several drawbacks. Small sample size can give rise to unreliable inferential conclusions and convergence issues make the approach unappealing. This paper suggests a different methodology to address such difficulties.Methods: A SIMEX methodology is proposed. The method is a simulation-based technique originally developed as a correction strategy within the measurement error literature. It suits the meta-analysis framework as the diagnostic accuracy measures provided by each study are prone to measurement error. SIMEX can be straightforwardly adapted to cover different measurement error structures and to deal with covariates. The effortless implementation with standard software is an interesting feature of the method.Results: Extensive simulation studies highlight the improvement provided by SIMEX over likelihood approach in terms of empirical coverage probabilities of confidence intervals under different scenarios, independently of the sample size and the values of the correlation between sensitivity and specificity. A remarkable amelioration is obtained in case of deviations from the normality assumption for the random-effects distribution. From a computational point of view, the application of SIMEX is shown to be neither involved nor subject to the convergence issues affecting likelihood-based alternatives. Application of the method to a diagnostic review of the performance of transesophageal echocardiography for assessing ascending aorta atherosclerosis enables overcoming limitations of the likelihood procedure.Conclusions: The SIMEX methodology represents an interesting alternative to likelihood-based procedures for inference in meta-analysis of diagnostic accuracy studies. The approach can provide more accurate inferential conclusions, while avoiding convergence failure and numerical instabilities. The application of the method in the R programming language is possible through the code which is made available and illustrated using the real data example. [ABSTRACT FROM AUTHOR]

: Comparing incidental dose distributions (i.e. images) of patients with different outcomes is a straightforward way to explore dose-response hypotheses in radiotherapy. In this paper, we introduced a permutation test that compares images, such as dose distributions from radiotherapy, while tackling the multiple comparisons problem. A test statistic Tmax was proposed that summarizes the differences between the images into a single value and a permutation procedure was employed to compute the adjusted p-value. We demonstrated the method in two retrospective studies: a prostate study that relates 3D dose distributions to failure, and an esophagus study that relates 2D surface dose distributions of the esophagus to acute esophagus toxicity. As a result, we were able to identify suspicious regions that are significantly associated with failure (prostate study) or toxicity (esophagus study). Permutation testing allows direct comparison of images from different patient categories and is a useful tool for data mining in radiotherapy. [ABSTRACT FROM AUTHOR]

Background: Systems biology models tend to become large since biological systems often consist of complex networks of interacting components, and since the models usually are developed to reflect various mechanistic assumptions of those networks. Nevertheless, not all aspects of the model are equally interesting in a given setting, and normally there are parts that can be reduced without affecting the relevant model performance. There are many methods for model reduction, but few or none of them allow for a restoration of the details of the original model after the simplified model has been simulated. Results: We present a reduction method that allows for such a back-translation from the reduced to the original model. The method is based on lumping of states, and includes a general and formal algorithm for both determining appropriate lumps, and for calculating the analytical back-translation formulas. The lumping makes use of efficient methods from graph-theory and _-decomposition and is derived and exemplified on two published models for fluorescence emission in photosynthesis. The bigger of these models is reduced from 26 to 6 states, with a negligible deviation from the reduced model simulations, both when comparing simulations in the states of the reduced model and when comparing back-translated simulations in the states of the original model. The method is developed in a linear setting, but we exemplify how the same concepts and approaches can be applied to non-linear problems. Importantly, the method automatically provides a reduced model with back-translations. Also, the method is implemented as a part of the systems biology toolbox for matlab, and the matlab scripts for the examples in this paper are available in the supplementary material. Conclusions: Our novel lumping methodology allows for both automatic reduction of states using lumping, and for analytical retrieval of the original states and parameters without performing a new simulation. The two models can thus be considered as two degrees of zooming of the same model. This is a conceptually new development of model reduction approaches, which we think will stimulate much further research and will prove to be very useful in future modelling projects. [ABSTRACT FROM AUTHOR]

Background: Assistive technologies aim to increase quality of life, reduce dependence on care giver and on the long term care system. Several studies have demonstrated the effectiveness in the use of assistive technology for environment control and communication systems. The progress of brain-computer interfaces (BCI) research together with exoskeleton enable a person with motor impairment to interact with new elements in the environment. This paper aims to evaluate the environment control interface (ECI) developed under the AIDE project conditions, a multimodal interface able to analyze and extract relevant information from the environments as well as from the identification of residual abilities, behaviors, and intentions of the user.Methods: This study evaluated the ECI in a simulated scenario using a two screen layout: one with the ECI and the other with a simulated home environment, developed for this purpose. The sensorimotor rhythms and the horizontal oculoversion, acquired through BCI2000, a multipurpose standard BCI platform, were used to online control the ECI after the user training and system calibration. Eight subjects with different neurological diseases and spinal cord injury participated in this study. The subjects performed simulated activities of daily living (ADLs), i.e. actions in the simulated environment as drink, switch on a lamp or raise the bed head, during ten minutes in two different modes, AIDE mode, using a prediction model, to recognize the user intention facilitating the scan, and Manual mode, without a prediction model.Results: The results show that the mean task time spent in the AIDE mode was less than in the Manual, i.e the users were able to perform more tasks in the AIDE mode during the same time. The results showed a statistically significant differences with p<0.001. Regarding the steps, i.e the number of abstraction levels crossed in the ECI to perform an ADL, the users performed one step in the 90% of the tasks using the AIDE mode and three steps, at least, were necessary in the Manual mode. The user's intention prediction was performed through conditional random fields (CRF), with a global accuracy about 87%.Conclusions: The environment analysis and the identification of the user's behaviors can be used to predict the user intention opening a new paradigm in the design of the ECIs. Although the developed ECI was tested only in a simulated home environment, it can be easily adapted to a real environment increasing the user independence at home. [ABSTRACT FROM AUTHOR]

Background: Atherosclerotic plaque progression and rupture are believed to be associated with mechanical stress conditions. In this paper, patient-specific in vivo intravascular ultrasound (IVUS) coronary plaque image data were used to construct computational models with fluid-structure interaction (FSI) and cyclic bending to investigate correlations between plaque wall thickness and both flow shear stress and plaque wall stress conditions.Methods: IVUS data were acquired from 10 patients after voluntary informed consent. The X-ray angiogram was obtained prior to the pullback of the IVUS catheter to determine the location of the coronary artery stenosis, vessel curvature and cardiac motion. Cyclic bending was specified in the model representing the effect by heart contraction. 3D anisotropic FSI models were constructed and solved to obtain flow shear stress (FSS) and plaque wall stress (PWS) values. FSS and PWS values were obtained for statistical analysis. Correlations with p < 0.05 were deemed significant.Results: Nine out of the 10 patients showed positive correlation between wall thickness and flow shear stress. The mean Pearson correlation r-value was 0.278 ± 0.181. Similarly, 9 out of the 10 patients showed negative correlation between wall thickness and plaque wall stress. The mean Pearson correlation r-value was -0.530 ± 0.210.Conclusion: Our results showed that plaque vessel wall thickness correlated positively with FSS and negatively with PWS. The patient-specific IVUS-based modeling approach has the potential to be used to investigate and identify possible mechanisms governing plaque progression and rupture and assist in diagnosis and intervention procedures. This represents a new direction of research. Further investigations using more patient follow-up data are warranted. [ABSTRACT FROM AUTHOR]

A conference paper about CaliBayes is presented. CaliBayes is an integration of GRID based simulation and data resources for Bayesian calibration of biological models. The goal is to find a combination of parameters that give the best match to available experimental data, together with some idea of the uncertainty in the parameter value. The calibration algorithm selects parameter combinations at which the simulator will be run.