@article{SherrattSrivastavaAinslieetal.2024, author = {Sherratt, Katharine and Srivastava, Ajitesh and Ainslie, Kylie and Singh, David E. and Cublier, Aymar and Marinescu, Maria Cristina and Carretero, Jesus and Garcia, Alberto Cascajo and Franco, Nicolas and Willem, Lander and Abrams, Steven and Faes, Christel and Beutels, Philippe and Hens, Niel and M{\"u}ller, Sebastian and Charlton, Billy and Ewert, Ricardo and Paltra, Sydney and Rakow, Christian and Rehmann, Jakob and Conrad, Tim and Sch{\"u}tte, Christof and Nagel, Kai and Abbott, Sam and Grah, Rok and Niehus, Rene and Prasse, Bastian and Sandmann, Frank and Funk, Sebastian}, title = {Characterising information gains and losses when collecting multiple epidemic model outputs}, volume = {47}, journal = {Epidemics}, publisher = {Elsevier BV}, issn = {1755-4365}, doi = {10.1016/j.epidem.2024.100765}, year = {2024}, abstract = {Collaborative comparisons and combinations of epidemic models are used as policy-relevant evidence during epidemic outbreaks. In the process of collecting multiple model projections, such collaborations may gain or lose relevant information. Typically, modellers contribute a probabilistic summary at each time-step. We compared this to directly collecting simulated trajectories. We aimed to explore information on key epidemic quantities; ensemble uncertainty; and performance against data, investigating potential to continuously gain information from a single cross-sectional collection of model results. Methods We compared July 2022 projections from the European COVID-19 Scenario Modelling Hub. Five modelling teams projected incidence in Belgium, the Netherlands, and Spain. We compared projections by incidence, peaks, and cumulative totals. We created a probabilistic ensemble drawn from all trajectories, and compared to ensembles from a median across each model's quantiles, or a linear opinion pool. We measured the predictive accuracy of individual trajectories against observations, using this in a weighted ensemble. We repeated this sequentially against increasing weeks of observed data. We evaluated these ensembles to reflect performance with varying observed data. Results. By collecting modelled trajectories, we showed policy-relevant epidemic characteristics. Trajectories contained a right-skewed distribution well represented by an ensemble of trajectories or a linear opinion pool, but not models' quantile intervals. Ensembles weighted by performance typically retained the range of plausible incidence over time, and in some cases narrowed this by excluding some epidemic shapes. Conclusions. We observed several information gains from collecting modelled trajectories rather than quantile distributions, including potential for continuously updated information from a single model collection. The value of information gains and losses may vary with each collaborative effort's aims, depending on the needs of projection users. Understanding the differing information potential of methods to collect model projections can support the accuracy, sustainability, and communication of collaborative infectious disease modelling efforts. Data availability All code and data available on Github: https://github.com/covid19-forecast-hub-europe/aggregation-info-loss}, language = {en} } @article{AmiranashviliLuedkeLietal.2024, author = {Amiranashvili, Tamaz and L{\"u}dke, David and Li, Hongwei Bran and Zachow, Stefan and Menze, Bjoern}, title = {Learning continuous shape priors from sparse data with neural implicit functions}, volume = {94}, journal = {Medical Image Analysis}, doi = {10.1016/j.media.2024.103099}, pages = {103099}, year = {2024}, abstract = {Statistical shape models are an essential tool for various tasks in medical image analysis, including shape generation, reconstruction and classification. Shape models are learned from a population of example shapes, which are typically obtained through segmentation of volumetric medical images. In clinical practice, highly anisotropic volumetric scans with large slice distances are prevalent, e.g., to reduce radiation exposure in CT or image acquisition time in MR imaging. For existing shape modeling approaches, the resolution of the emerging model is limited to the resolution of the training shapes. Therefore, any missing information between slices prohibits existing methods from learning a high-resolution shape prior. We propose a novel shape modeling approach that can be trained on sparse, binary segmentation masks with large slice distances. This is achieved through employing continuous shape representations based on neural implicit functions. After training, our model can reconstruct shapes from various sparse inputs at high target resolutions beyond the resolution of individual training examples. We successfully reconstruct high-resolution shapes from as few as three orthogonal slices. Furthermore, our shape model allows us to embed various sparse segmentation masks into a common, low-dimensional latent space — independent of the acquisition direction, resolution, spacing, and field of view. We show that the emerging latent representation discriminates between healthy and pathological shapes, even when provided with sparse segmentation masks. Lastly, we qualitatively demonstrate that the emerging latent space is smooth and captures characteristic modes of shape variation. We evaluate our shape model on two anatomical structures: the lumbar vertebra and the distal femur, both from publicly available datasets.}, language = {en} } @article{ObermeierHeimBiereetal.2022, author = {Obermeier, Patrick E and Heim, Albert and Biere, Barbara and Hage, Elias and Alchikh, Maren and Conrad, Tim and Schweiger, Brunhilde and Rath, Barbara A}, title = {Linking digital surveillance and in-depth virology to study clinical patterns of viral respiratory infections in vulnerable patient populations}, volume = {25}, journal = {iScience}, number = {5}, publisher = {Cell Press}, doi = {10.1016/j.isci.2022.104276}, year = {2022}, abstract = {To improve the identification and management of viral respiratory infections, we established a clinical and virologic surveillance program for pediatric patients fulfilling pre-defined case criteria of influenza-like illness and viral respiratory infections. The program resulted in a cohort comprising 6,073 patients (56\% male, median age 1.6 years, range 0-18.8 years), where every patient was assessed with a validated disease severity score at the point-of-care using the ViVI ScoreApp. We used machine learning and agnostic feature selection to identify characteristic clinical patterns. We tested all patients for human adenoviruses, 571 (9\%) were positive. Adenovirus infections were particularly common and mild in children ≥1 month of age but rare and potentially severe in neonates: with lower airway involvement, disseminated disease, and a 50\% mortality rate (n = 2/4). In one fatal case, we discovered a novel virus …}, language = {en} } @article{HajarolasvadiSunkaraKhavnekaretal.2022, author = {Hajarolasvadi, Noushin and Sunkara, Vikram and Khavnekar, Sagar and Beck, Florian and Brandt, Robert and Baum, Daniel}, title = {Volumetric macromolecule identification in cryo-electron tomograms using capsule networks}, volume = {23}, journal = {BMC Bioinformatics}, number = {360}, doi = {10.1186/s12859-022-04901-w}, year = {2022}, abstract = {Background: Despite recent advances in cellular cryo-electron tomography (CET), developing automated tools for macromolecule identification in submolecular resolution remains challenging due to the lack of annotated data and high structural complexities. To date, the extent of the deep learning methods constructed for this problem is limited to conventional Convolutional Neural Networks (CNNs). Identifying macromolecules of different types and sizes is a tedious and time-consuming task. In this paper, we employ a capsule-based architecture to automate the task of macro- molecule identification, that we refer to as 3D-UCaps. In particular, the architecture is composed of three components: feature extractor, capsule encoder, and CNN decoder. The feature extractor converts voxel intensities of input sub-tomograms to activities of local features. The encoder is a 3D Capsule Network (CapsNet) that takes local features to generate a low-dimensional representation of the input. Then, a 3D CNN decoder reconstructs the sub-tomograms from the given representation by upsampling. Results: We performed binary and multi-class localization and identification tasks on synthetic and experimental data. We observed that the 3D-UNet and the 3D-UCaps had an F1-score mostly above 60\% and 70\%, respectively, on the test data. In both network architectures, we observed degradation of at least 40\% in the F1-score when identifying very small particles (PDB entry 3GL1) compared to a large particle (PDB entry 4D8Q). In the multi-class identification task of experimental data, 3D-UCaps had an F1-score of 91\% on the test data in contrast to 64\% of the 3D-UNet. The better F1-score of 3D-UCaps compared to 3D-UNet is obtained by a higher precision score. We speculate this to be due to the capsule network employed in the encoder. To study the effect of the CapsNet-based encoder architecture further, we performed an ablation study and perceived that the F1-score is boosted as network depth is increased which is in contrast to the previously reported results for the 3D-UNet. To present a reproducible work, source code, trained models, data as well as visualization results are made publicly available. Conclusion: Quantitative and qualitative results show that 3D-UCaps successfully perform various downstream tasks including identification and localization of macro- molecules and can at least compete with CNN architectures for this task. Given that the capsule layers extract both the existence probability and the orientation of the molecules, this architecture has the potential to lead to representations of the data that are better interpretable than those of 3D-UNet.}, language = {en} } @article{LelievreZhang2019, author = {Leli{\`e}vre, Tony and Zhang, Wei}, title = {Pathwise estimates for effective dynamics: the case of nonlinear vectorial reaction coordinates}, journal = {Multiscale Modeling and Simulation}, number = {17}, arxiv = {http://arxiv.org/abs/1805.01928}, doi = {10.1137/18M1186034}, pages = {1019 -- 1051}, year = {2019}, abstract = {Effective dynamics using conditional expectation was proposed in [F. Legoll and T. Leli{\`e}vre, Nonlinearity, 2010] to approximate the essential dynamics of high-dimensional diffusion processes along a given reaction coordinate. The approximation error of the effective dynamics when it is used to approximate the behavior of the original dynamics has been considered in recent years. As a continuation of the previous work [F. Legoll, T. Leli{\`e}vre, and S. Olla, Stoch. Process. Appl, 2017], in this paper we obtain pathwise estimates for effective dynamics when the reaction coordinate function is either nonlinear or vector-valued.}, language = {en} } @article{HartmannSchuetteZhang2019, author = {Hartmann, Carsten and Sch{\"u}tte, Christof and Zhang, Wei}, title = {Jarzynski's equality, fluctuation theorems, and variance reduction: Mathematical analysis and numerical algorithms}, volume = {175}, journal = {Journal of Statistical Physics}, number = {6}, arxiv = {http://arxiv.org/abs/1803.09347}, doi = {10.1007/s10955-019-02286-4}, pages = {1214 -- 1261}, year = {2019}, abstract = {In this paper, we study Jarzynski's equality and fluctuation theorems for diffusion processes. While some of the results considered in the current work are known in the (mainly physics) literature, we review and generalize these nonequilibrium theorems using mathematical arguments, therefore enabling further investigations in the mathematical community. On the numerical side, variance reduction approaches such as importance sampling method are studied in order to compute free energy differences based on Jarzynski's equality.}, language = {en} } @article{ZhangKlusConradetal.2019, author = {Zhang, Wei and Klus, Stefan and Conrad, Tim and Sch{\"u}tte, Christof}, title = {Learning chemical reaction networks from trajectory data}, volume = {18}, journal = {SIAM Journal on Applied Dynamical Systems (SIADS)}, number = {4}, arxiv = {http://arxiv.org/abs/1902.04920}, doi = {10.1137/19M1265880}, pages = {2000 -- 2046}, year = {2019}, abstract = {We develop a data-driven method to learn chemical reaction networks from trajectory data. Modeling the reaction system as a continuous-time Markov chain and assuming the system is fully observed,our method learns the propensity functions of the system with predetermined basis functions by maximizing the likelihood function of the trajectory data under l^1 sparse regularization. We demonstrate our method with numerical examples using synthetic data and carry out an asymptotic analysis of the proposed learning procedure in the infinite-data limit.}, language = {en} } @article{Zhang2019, author = {Zhang, Wei}, title = {Ergodic SDEs on submanifolds and related numerical sampling schemes}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis}, arxiv = {http://arxiv.org/abs/1702.08064}, year = {2019}, abstract = {In many applications, it is often necessary to sample the mean value of certain quantity with respect to a probability measure \$\mu\$ on the level set of a smooth function ξ:R^d→R^k, 1≤k