Showing 1–50 of 125 results for author: Janssen, A

NMGrad: Advancing Histopathological Bladder Cancer Grading with Weakly Supervised Deep Learning

Authors: Saul Fuster, Umay Kiraz, Trygve Eftestøl, Emiel A. M. Janssen, Kjersti Engan

Abstract: The most prevalent form of bladder cancer is urothelial carcinoma, characterized by a high recurrence rate and substantial lifetime treatment costs for patients. Grading is a prime factor for patient risk stratification, although it suffers from inconsistencies and variations among pathologists. Moreover, absence of annotations in medical imaging difficults training deep learning models. To addres… ▽ More The most prevalent form of bladder cancer is urothelial carcinoma, characterized by a high recurrence rate and substantial lifetime treatment costs for patients. Grading is a prime factor for patient risk stratification, although it suffers from inconsistencies and variations among pathologists. Moreover, absence of annotations in medical imaging difficults training deep learning models. To address these challenges, we introduce a pipeline designed for bladder cancer grading using histological slides. First, it extracts urothelium tissue tiles at different magnification levels, employing a convolutional neural network for processing for feature extraction. Then, it engages in the slide-level prediction process. It employs a nested multiple instance learning approach with attention to predict the grade. To distinguish different levels of malignancy within specific regions of the slide, we include the origins of the tiles in our analysis. The attention scores at region level is shown to correlate with verified high-grade regions, giving some explainability to the model. Clinical evaluations demonstrate that our model consistently outperforms previous state-of-the-art methods. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: https://github.com/Biomedical-Data-Analysis-Laboratory/GradeMIL

Self-Contrastive Weakly Supervised Learning Framework for Prognostic Prediction Using Whole Slide Images

Authors: Saul Fuster, Farbod Khoraminia, Julio Silva-Rodríguez, Umay Kiraz, Geert J. L. H. van Leenders, Trygve Eftestøl, Valery Naranjo, Emiel A. M. Janssen, Tahlita C. M. Zuiverloon, Kjersti Engan

Abstract: We present a pioneering investigation into the application of deep learning techniques to analyze histopathological images for addressing the substantial challenge of automated prognostic prediction. Prognostic prediction poses a unique challenge as the ground truth labels are inherently weak, and the model must anticipate future events that are not directly observable in the image. To address thi… ▽ More We present a pioneering investigation into the application of deep learning techniques to analyze histopathological images for addressing the substantial challenge of automated prognostic prediction. Prognostic prediction poses a unique challenge as the ground truth labels are inherently weak, and the model must anticipate future events that are not directly observable in the image. To address this challenge, we propose a novel three-part framework comprising of a convolutional network based tissue segmentation algorithm for region of interest delineation, a contrastive learning module for feature extraction, and a nested multiple instance learning classification module. Our study explores the significance of various regions of interest within the histopathological slides and exploits diverse learning scenarios. The pipeline is initially validated on artificially generated data and a simpler diagnostic task. Transitioning to prognostic prediction, tasks become more challenging. Employing bladder cancer as use case, our best models yield an AUC of 0.721 and 0.678 for recurrence and treatment outcome prediction respectively. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: https://github.com/Biomedical-Data-Analysis-Laboratory/HistoPrognostics

Equipping Computational Pathology Systems with Artifact Processing Pipelines: A Showcase for Computation and Performance Trade-offs

Authors: Neel Kanwal, Farbod Khoraminia, Umay Kiraz, Andres Mosquera-Zamudio, Carlos Monteagudo, Emiel A. M. Janssen, Tahlita C. M. Zuiverloon, Chunmig Rong, Kjersti Engan

Abstract: Histopathology is a gold standard for cancer diagnosis under a microscopic examination. However, histological tissue processing procedures result in artifacts, which are ultimately transferred to the digitized version of glass slides, known as whole slide images (WSIs). Artifacts are diagnostically irrelevant areas and may result in wrong deep learning (DL) algorithms predictions. Therefore, detec… ▽ More Histopathology is a gold standard for cancer diagnosis under a microscopic examination. However, histological tissue processing procedures result in artifacts, which are ultimately transferred to the digitized version of glass slides, known as whole slide images (WSIs). Artifacts are diagnostically irrelevant areas and may result in wrong deep learning (DL) algorithms predictions. Therefore, detecting and excluding artifacts in the computational pathology (CPATH) system is essential for reliable automated diagnosis. In this paper, we propose a mixture of experts (MoE) scheme for detecting five notable artifacts, including damaged tissue, blur, folded tissue, air bubbles, and histologically irrelevant blood from WSIs. First, we train independent binary DL models as experts to capture particular artifact morphology. Then, we ensemble their predictions using a fusion mechanism. We apply probabilistic thresholding over the final probability distribution to improve the sensitivity of the MoE. We developed DL pipelines using two MoEs and two multiclass models of state-of-the-art deep convolutional neural networks (DCNNs) and vision transformers (ViTs). DCNNs-based MoE and ViTs-based MoE schemes outperformed simpler multiclass models and were tested on datasets from different hospitals and cancer types, where MoE using DCNNs yielded the best results. The proposed MoE yields 86.15% F1 and 97.93% sensitivity scores on unseen data, retaining less computational cost for inference than MoE using ViTs. This best performance of MoEs comes with relatively higher computational trade-offs than multiclass models. The proposed artifact detection pipeline will not only ensure reliable CPATH predictions but may also provide quality control. △ Less

Submitted 23 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Submitted to BMC Medical Informatics and Decision Making Journal

Infinite-horizon Fuk-Nagaev inequalities

Authors: A. J. E. M. Janssen, B. Zwart

Abstract: We develop explicit bounds for the tail of the distribution for the all-time supremum of a random walk with negative drift, where the increments have a truncated heavy-tailed distribution. As an application, we consider a ruin problem in the presence of re-insurance. We develop explicit bounds for the tail of the distribution for the all-time supremum of a random walk with negative drift, where the increments have a truncated heavy-tailed distribution. As an application, we consider a ruin problem in the presence of re-insurance. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Bounding Taylor approximation errors for the exponential function in the presence of a power weight function

Authors: A. J. E. M. Janssen

Abstract: Motivated by the needs in the theory of large deviations and in the theory of Lundberg's equation with heavy-tailed distribution functions, we study for $n=0,1,...$ the maximization of $S:~\Bigl(1-e^{-s}\Bigl(1+\frac{s^1}{1!}+...+\frac{s^n}{n!}\Bigr)\Bigr)/s^δ = E_{n,δ}(s)$ over $s\geq0$, with $δ\in(0,n+1)$,… ▽ More Motivated by the needs in the theory of large deviations and in the theory of Lundberg's equation with heavy-tailed distribution functions, we study for $n=0,1,...$ the maximization of $S:~\Bigl(1-e^{-s}\Bigl(1+\frac{s^1}{1!}+...+\frac{s^n}{n!}\Bigr)\Bigr)/s^δ = E_{n,δ}(s)$ over $s\geq0$, with $δ\in(0,n+1)$, $U:~({-}1)^{n+1}\Bigl(e^{-u}-\Bigl(1-\frac{u^1}{1!}+...+({-}1)^n\,\frac{u^n}{n!} \Bigr)\Bigr)/u^δ=G_{n,δ}(u)$ over $u\geq0$ with $δ\in(n,n+1)$. We show that $E_{n,δ}(s)$ and $G_{n,δ}(u)$ have a unique maximizer $s=s_n(δ)>0$ and $u=u_n(δ)>0$ that decrease strictly from $+\infty$ at $δ=0$ and $δ=n$, respectively, to 0 at $δ=n+1$. We use Taylor's formula for truncated series with remainder in integral form to develop a criterion to decide whether a particular smooth function $S(δ)$, $δ\in(0,n+1)$, or $U(δ)$, $δ\in(n,n+1)$, respectively, is a lower/upper bound for $s_n(δ)$ and $u_n(δ)$, respectively. This criterion allows us to find lower and upper bounds for $s_n$ and $u_n$ that are reasonably tight and simple at the same time. Furthermore, as a consequence of the identities $\frac{d}{dδ}\,[{\rm ln}\,ME_{n,δ}] ={-}{\rm ln}\,s_n(δ)$ and $\frac{d}{dδ}\,[{\rm ln}\,MG_{n,δ}]={-}{\rm ln}\,u_n(δ)$, we show that $ME_{n,δ}$ and $MG_{n,δ}$ are log-convex functions of $δ\in(0,n+1)$ and $δ\in(n+1,n)$, respectively, with limiting values 1 ($δ\downarrow0$) and $1/(n+1)!$ ($δ\uparrow n+1$) for $E$, and $1/n!\,(δ\downarrow n)$ and $1/(n+1)!\,(δ\uparrow n+1)$ for $G$. The minimal values $\hat{E}_n$ and $\hat{G}_n$ of $ME_{n,δ}$ and $MG_{n,δ}$, respectively, as a function of $δ$, as well as the minimum locations $δ_{n,E}$ and $δ_{n,G}$ are determined in closed form. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Reduced-order Modeling of Modular, Position-dependent Systems with Translating Interfaces

Authors: Robert A. Egelmeers, Lars A. L. Janssen, Rob H. B. Fey, Jasper Gerritsen, Nathan van de Wouw

Abstract: Many complex mechatronic systems consist of multiple interconnected dynamical subsystems, which are designed, developed, analyzed, and manufactured by multiple independent teams. To support such a design approach, a modular model framework is needed to reduce computational complexity and, at the same time, enable multiple teams to develop and analyze the subsystems in parallel. In such a modular f… ▽ More Many complex mechatronic systems consist of multiple interconnected dynamical subsystems, which are designed, developed, analyzed, and manufactured by multiple independent teams. To support such a design approach, a modular model framework is needed to reduce computational complexity and, at the same time, enable multiple teams to develop and analyze the subsystems in parallel. In such a modular framework, the subsystem models are typically interconnected by means of a static interconnection structure. However, many complex dynamical systems exhibit position-dependent behavior (e.g., induced by translating interfaces) which cannot be not captured by such static interconnection models. In this paper, a modular model framework is proposed, which allows to construct an interconnected system model, which captures the position-dependent behavior of systems with translating interfaces, such as linear guide rails, through a position-dependent interconnection structure. Additionally, this framework allows to apply model reduction on subsystem level, enabling a more effective reduction approach, tailored to the specific requirements of each subsystem. Furthermore, we show the effectiveness of this framework on an industrial wire bonder. Here, we show that including a position-dependent model of the interconnection structure 1) enables to accurately model the dynamics of a system over the operating range of the system and, 2) modular model reduction methods can be used to obtain a computationally efficient interconnected system model with guaranteed accuracy specifications. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Modular Redesign of Mechatronic Systems: Formulation of Module Specifications Guaranteeing System Dynamics Specifications

Authors: Lars A. L. Janssen, Rob H. B. Fey, Bart Besselink, Nathan van de Wouw

Abstract: Complex mechatronic systems are typically composed of interconnected modules, often developed by independent teams. This development process challenges the verification of system specifications before all modules are integrated. To address this challenge, a modular redesign framework is proposed in this paper. Herein, first, allowed changes in the dynamics (represented by frequency response functi… ▽ More Complex mechatronic systems are typically composed of interconnected modules, often developed by independent teams. This development process challenges the verification of system specifications before all modules are integrated. To address this challenge, a modular redesign framework is proposed in this paper. Herein, first, allowed changes in the dynamics (represented by frequency response functions (FRFs)) of the redesigned system are defined with respect to the original system model, which already satisfies system specifications. Second, these allowed changes in the overall system dynamics (or system redesign specifications) are automatically translated to dynamics (FRF) specifications on module level that, when satisfied, guarantee overall system dynamics (FRF) specifications. This modularity in specification management supports local analysis and verification of module design changes, enabling design teams to work in parallel without the need to iteratively rebuild the system model to check fulfilment of system FRF specifications. A modular redesign process results that shortens time-to-market and decreases redesign costs. The framework's effectiveness is demonstrated through three examples of increasing complexity, highlighting its potential to enable modular mechatronic system (re)design. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Lessons Learned from Designing an Open-Source Automated Feedback System for STEM Education

Authors: Steffen Steinert, Lars Krupp, Karina E. Avila, Anke S. Janssen, Verena Ruf, David Dzsotjan, Christian De Schryver, Jakob Karolus, Stefan Ruzika, Karen Joisten, Paul Lukowicz, Jochen Kuhn, Norbert Wehn, Stefan Küchemann

Abstract: As distance learning becomes increasingly important and artificial intelligence tools continue to advance, automated systems for individual learning have attracted significant attention. However, the scarcity of open-source online tools that are capable of providing personalized feedback has restricted the widespread implementation of research-based feedback systems. In this work, we present RATsA… ▽ More As distance learning becomes increasingly important and artificial intelligence tools continue to advance, automated systems for individual learning have attracted significant attention. However, the scarcity of open-source online tools that are capable of providing personalized feedback has restricted the widespread implementation of research-based feedback systems. In this work, we present RATsApp, an open-source automated feedback system (AFS) that incorporates research-based features such as formative feedback. The system focuses on core STEM competencies such as mathematical competence, representational competence, and data literacy. It also allows lecturers to monitor students' progress. We conducted a survey based on the technology acceptance model (TAM2) among a set of students (N=64). Our findings confirm the applicability of the TAM2 framework, revealing that factors such as the relevance of the studies, output quality, and ease of use significantly influence the perceived usefulness. We also found a linear relation between the perceived usefulness and the intention to use, which in turn is a significant predictor of the frequency of use. Moreover, the formative feedback feature of RATsApp received positive feedback, indicating its potential as an educational tool. Furthermore, as an open-source platform, RATsApp encourages public contributions to its ongoing development, fostering a collaborative approach to improve educational tools. △ Less

Submitted 19 January, 2024; originally announced January 2024.

Comments: 40 pages, 6 figures

On the theoretical framework for meniscus-guided manufacturing of large-area OPV modules

Authors: Fabian Gumpert, Annika Janßen, Robin Basu, Christoph J. Brabec, Hans-Joachim Egelhaaf, Jan Lohbreier, Andreas Distler

Abstract: For the manufacturing of thin films of solution-processable organic semiconductors, e.g. for organic photovoltaics (OPV), meniscus guided-coating techniques are the method of choice for large-scale industrial applications. However, the process requires an in-depth understanding of the respective fluid dynamics to control the resulting film thickness. In this article, we derive an analytical expres… ▽ More For the manufacturing of thin films of solution-processable organic semiconductors, e.g. for organic photovoltaics (OPV), meniscus guided-coating techniques are the method of choice for large-scale industrial applications. However, the process requires an in-depth understanding of the respective fluid dynamics to control the resulting film thickness. In this article, we derive an analytical expression to describe the layer thickness of coatings manufactured with a trapezoidal-shaped applicator as a function of various fluid and process parameters. The analytical calculations are compared with results from computational fluid dynamics (CFD) simulations and experimental data for an industrially relevant OPV active material system. The analytical calculations are compared with results from computational fluid dynamics (CFD) simulations and experimental data for an industrially relevant OPV active material system. The good agreement of all three approaches demonstrates the potential of the analytical and simulative methods to reduce time- and resource-consuming experiments to a minimum. Furthermore, our theoretical model can be used to enhance the homogeneity of large-area coatings by means of an acceleration profile of the applicator that can compensate the liquid loss during the coating process. The respective analytical expression is validated by simulated and experimentally obtained data for long-distance coatings. Finally, this approach is used to fabricate a large-area OPV module with new world record efficiency. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 23 pages, 9 figures

MSC Class: 65D99

A Dual Convolutional Neural Network Pipeline for Melanoma Diagnostics and Prognostics

Authors: Marie Bø-Sande, Edvin Benjaminsen, Neel Kanwal, Saul Fuster, Helga Hardardottir, Ingrid Lundal, Emiel A. M. Janssen, Kjersti Engan

Abstract: Melanoma is a type of cancer that begins in the cells controlling the pigment of the skin, and it is often referred to as the most dangerous skin cancer. Diagnosing melanoma can be time-consuming, and a recent increase in melanoma incidents indicates a growing demand for a more efficient diagnostic process. This paper presents a pipeline for melanoma diagnostics, leveraging two convolutional neura… ▽ More Melanoma is a type of cancer that begins in the cells controlling the pigment of the skin, and it is often referred to as the most dangerous skin cancer. Diagnosing melanoma can be time-consuming, and a recent increase in melanoma incidents indicates a growing demand for a more efficient diagnostic process. This paper presents a pipeline for melanoma diagnostics, leveraging two convolutional neural networks, a diagnosis, and a prognosis model. The diagnostic model is responsible for localizing malignant patches across whole slide images and delivering a patient-level diagnosis as malignant or benign. Further, the prognosis model utilizes the diagnostic model's output to provide a patient-level prognosis as good or bad. The full pipeline has an F1 score of 0.79 when tested on data from the same distribution as it was trained on. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: Accepted at NLDL 2024

An Ontological Model of User Preferences

Authors: Mona Abdel-Keream, Daniel Beßler, Ayden Janssen, Sascha Jongebloed, Robin Nolte, Mihai Pomarlan, Robert Porzel

Abstract: The notion of preferences plays an important role in many disciplines including service robotics which is concerned with scenarios in which robots interact with humans. These interactions can be favored by robots taking human preferences into account. This raises the issue of how preferences should be represented to support such preference-aware decision making. Several formal accounts for a notio… ▽ More The notion of preferences plays an important role in many disciplines including service robotics which is concerned with scenarios in which robots interact with humans. These interactions can be favored by robots taking human preferences into account. This raises the issue of how preferences should be represented to support such preference-aware decision making. Several formal accounts for a notion of preferences exist. However, these approaches fall short on defining the nature and structure of the options that a robot has in a given situation. In this work, we thus investigate a formal model of preferences where options are non-atomic entities that are defined by the complex situations they bring about. △ Less

Submitted 29 October, 2023; originally announced October 2023.

Mode Selection for Component Mode Synthesis with Guaranteed Assembly Accuracy

Authors: Lars A. L. Janssen, Rob H. B. Fey, Bart Besselink, Nathan van de Wouw

Abstract: In this work, a modular approach is introduced to select the most important eigenmodes for each component of a composed structural dynamics system to obtain the required accuracy of the reduced-order assembly model. To enable the use of models of complex (structural) dynamical systems in engineering practice, e.g., in a design, optimization and/or control context, the complexity of the models need… ▽ More In this work, a modular approach is introduced to select the most important eigenmodes for each component of a composed structural dynamics system to obtain the required accuracy of the reduced-order assembly model. To enable the use of models of complex (structural) dynamical systems in engineering practice, e.g., in a design, optimization and/or control context, the complexity of the models needs to be reduced. When the model consist of an assembly of multiple interconnected structural components, component mode synthesis is often the preferred model reduction method. The standard approach to component mode synthesis for such system is to select the eigenmodes of a component that are most important to accurately model the dynamic behavior of this component in a certain frequency range of interest. However, often, a more relevant goal is to obtain, in this frequency range, an accurate model of the assembly. In the proposed approach, accuracy requirements on the level of the assembly are translated to accuracy requirements on component level, by employing techniques from the field of systems and control. With these component-level requirements, the eigenmodes that are most important to accurately model the dynamic behavior of the assembly can be selected in a modular fashion. We demonstrate with two structural dynamics benchmark systems that this method based on assembly accuracy allows for a computationally efficient selection of eigenmodes that 1) guarantees satisfaction of the assembly accuracy requirements and 2) results in most cases in reduced-order models of significantly lower order with respect to the industrial standard approach in which component eigenmodes are selected using a frequency criterion. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Multivariate Regular Variation of Preferential Attachment Models

Authors: Anja Janßen, Max Ziegenbalg

Abstract: We use the framework of multivariate regular variation to analyze the extremal behavior of preferential attachment models. To this end, we follow a directed linear preferential attachment model with a fixed number of outgoing edges per node for a random, heavy-tailed number of steps in time and treat the incoming edge count of all existing nodes as a random vector of random length. By combining ma… ▽ More We use the framework of multivariate regular variation to analyze the extremal behavior of preferential attachment models. To this end, we follow a directed linear preferential attachment model with a fixed number of outgoing edges per node for a random, heavy-tailed number of steps in time and treat the incoming edge count of all existing nodes as a random vector of random length. By combining martingale properties, moment bounds and a Breiman type theorem we show that the resulting quantity is multivariate regularly varying, both as a vector of fixed length formed by the edge counts of a finite number of oldest nodes, and also as a vector of random length viewed in sequence space. A Pólya urn representation allows us to explicitly describe the extremal dependence between the degrees with the help of Dirichlet distributions. As a by-product of our analysis we establish new results for almost sure convergence of the edge counts in sequence space as the number of nodes goes to infinity. △ Less

Submitted 4 October, 2023; originally announced October 2023.

MSC Class: 60G70; 05C80; 60F25

Balancing Privacy and Progress in Artificial Intelligence: Anonymization in Histopathology for Biomedical Research and Education

Authors: Neel Kanwal, Emiel A. M. Janssen, Kjersti Engan

Abstract: The advancement of biomedical research heavily relies on access to large amounts of medical data. In the case of histopathology, Whole Slide Images (WSI) and clinicopathological information are valuable for developing Artificial Intelligence (AI) algorithms for Digital Pathology (DP). Transferring medical data "as open as possible" enhances the usability of the data for secondary purposes but pose… ▽ More The advancement of biomedical research heavily relies on access to large amounts of medical data. In the case of histopathology, Whole Slide Images (WSI) and clinicopathological information are valuable for developing Artificial Intelligence (AI) algorithms for Digital Pathology (DP). Transferring medical data "as open as possible" enhances the usability of the data for secondary purposes but poses a risk to patient privacy. At the same time, existing regulations push towards keeping medical data "as closed as necessary" to avoid re-identification risks. Generally, these legal regulations require the removal of sensitive data but do not consider the possibility of data linkage attacks due to modern image-matching algorithms. In addition, the lack of standardization in DP makes it harder to establish a single solution for all formats of WSIs. These challenges raise problems for bio-informatics researchers in balancing privacy and progress while developing AI algorithms. This paper explores the legal regulations and terminologies for medical data-sharing. We review existing approaches and highlight challenges from the histopathological perspective. We also present a data-sharing guideline for histological data to foster multidisciplinary research and education. △ Less

Submitted 8 August, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: Accepted to FAIEMA 2023

Physical Color Calibration of Digital Pathology Scanners for Robust Artificial Intelligence Assisted Cancer Diagnosis

Authors: Xiaoyi Ji, Richard Salmon, Nita Mulliqi, Umair Khan, Yinxi Wang, Anders Blilie, Henrik Olsson, Bodil Ginnerup Pedersen, Karina Dalsgaard Sørensen, Benedicte Parm Ulhøi, Svein R Kjosavik, Emilius AM Janssen, Mattias Rantalainen, Lars Egevad, Pekka Ruusuvuori, Martin Eklund, Kimmo Kartasalo

Abstract: The potential of artificial intelligence (AI) in digital pathology is limited by technical inconsistencies in the production of whole slide images (WSIs), leading to degraded AI performance and posing a challenge for widespread clinical application as fine-tuning algorithms for each new site is impractical. Changes in the imaging workflow can also lead to compromised diagnoses and patient safety r… ▽ More The potential of artificial intelligence (AI) in digital pathology is limited by technical inconsistencies in the production of whole slide images (WSIs), leading to degraded AI performance and posing a challenge for widespread clinical application as fine-tuning algorithms for each new site is impractical. Changes in the imaging workflow can also lead to compromised diagnoses and patient safety risks. We evaluated whether physical color calibration of scanners can standardize WSI appearance and enable robust AI performance. We employed a color calibration slide in four different laboratories and evaluated its impact on the performance of an AI system for prostate cancer diagnosis on 1,161 WSIs. Color standardization resulted in consistently improved AI model calibration and significant improvements in Gleason grading performance. The study demonstrates that physical color calibration provides a potential solution to the variation introduced by different scanners, making AI-based cancer diagnostics more reliable and applicable in clinical settings. △ Less

Submitted 7 July, 2023; originally announced July 2023.

Translating Assembly Accuracy Requirements to Cut-Off Frequencies for Component Mode Synthesis

Authors: Lars A. L. Janssen, Bart Besselink, Rob H. B. Fey, Nathan van de Wouw

Abstract: One of the most popular methods for reducing the complexity of assemblies of finite element models in the field of structural dynamics is component mode synthesis. A main challenge of component mode synthesis is balancing model complexity and model accuracy, because it is difficult to predict how component reduction influences assembly model accuracy. This work introduces an approach that allows f… ▽ More One of the most popular methods for reducing the complexity of assemblies of finite element models in the field of structural dynamics is component mode synthesis. A main challenge of component mode synthesis is balancing model complexity and model accuracy, because it is difficult to predict how component reduction influences assembly model accuracy. This work introduces an approach that allows for the translation of assembly model accuracy requirements in the frequency domain to the automatic selection of the cut-off frequencies for the model-order reduction (MOR) of components. The approach is based on a mathematical approach for MOR for coupled linear systems in the field of systems and control. We show how this approach is also applicable to structural dynamics models. We demonstrate the use of this approach in the scope of component mode synthesis (CMS) methods with the aim to reduce the complexity of component models while guaranteeing accuracy requirements of the assembly model. The proposed approach is illustrated on a mechanical, three-component structural dynamics system for which reduced-order models are computed that are reduced further compared to reduction using standard methods. This results in lower simulation cost, while maintaining the required accuracy. △ Less

Submitted 2 May, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Deep Learning for Predicting Metastasis on Melanoma WSIs

Authors: Christopher Andreassen, Saul Fuster, Helga Hardardottir, Emiel A. M. Janssen, Kjersti Engan

Abstract: Northern Europe has the second highest mortality rate of melanoma globally. In 2020, the mortality rate of melanoma rose to 1.9 per 100 000 habitants. Melanoma prognosis is based on a pathologist's subjective visual analysis of the patient's tumor. This methodology is heavily time-consuming, and the prognosis variability among experts is notable, drastically jeopardizing its reproducibility. Thus,… ▽ More Northern Europe has the second highest mortality rate of melanoma globally. In 2020, the mortality rate of melanoma rose to 1.9 per 100 000 habitants. Melanoma prognosis is based on a pathologist's subjective visual analysis of the patient's tumor. This methodology is heavily time-consuming, and the prognosis variability among experts is notable, drastically jeopardizing its reproducibility. Thus, the need for faster and more reproducible methods arises. Machine learning has paved its way into digital pathology, but so far, most contributions are on localization, segmentation, and diagnostics, with little emphasis on prognostics. This paper presents a convolutional neural network (CNN) method based on VGG16 to predict melanoma prognosis as the presence of metastasis within five years. Patches are extracted from regions of interest from Whole Slide Images (WSIs) at different magnification levels used in model training and validation. Results infer that utilizing WSI patches at 20x magnification level has the best performance, with an F1 score of 0.7667 and an AUC of 0.81. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Detection and Localization of Melanoma Skin Cancer in Histopathological Whole Slide Images

Authors: Neel Kanwal, Roger Amundsen, Helga Hardardottir, Luca Tomasetti, Erling Sandoy Undersrud, Emiel A. M. Janssen, Kjersti Engan

Abstract: Melanoma diagnosed and treated in its early stages can increase the survival rate. A projected increase in skin cancer incidents and a dearth of dermatopathologists have emphasized the need for computational pathology (CPATH) systems. CPATH systems with deep learning (DL) models have the potential to identify the presence of melanoma by exploiting underlying morphological and cellular features. Th… ▽ More Melanoma diagnosed and treated in its early stages can increase the survival rate. A projected increase in skin cancer incidents and a dearth of dermatopathologists have emphasized the need for computational pathology (CPATH) systems. CPATH systems with deep learning (DL) models have the potential to identify the presence of melanoma by exploiting underlying morphological and cellular features. This paper proposes a DL method to detect melanoma and distinguish between normal skin and benign/malignant melanocytic lesions in Whole Slide Images (WSI). Our method detects lesions with high accuracy and localizes them on a WSI to identify potential regions of interest for pathologists. Interestingly, our DL method relies on using a single CNN network to create localization maps first and use them to perform slide-level predictions to determine patients who have melanoma. Our best model provides favorable patch-wise classification results with a 0.992 F1 score and 0.99 sensitivity on unseen data. The source code is https://github.com/RogerAmundsen/Melanoma-Diagnosis-and-Localization-from-Whole-Slide-Images-using-Convolutional-Neural-Networks. △ Less

Submitted 4 November, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted at EUSIPCO 23

Journal ref: 10.23919/EUSIPCO58844.2023.10290087

Modular Model Reduction of Interconnected Systems: A Top-Down Approach

Authors: Lars A. L. Janssen, Bart Besselink, Rob H. B. Fey, Nathan van de Wouw

Abstract: Models of complex systems often consist of multiple interconnected subsystem/component models that are developed by multi-disciplinary teams of engineers or scientists. To ensure that such interconnected models can be applied for the purpose of simulation and/or control, a reduced-order model for the interconnected dynamics is needed. In the scope of this paper, we pursue this goal by subsystem re… ▽ More Models of complex systems often consist of multiple interconnected subsystem/component models that are developed by multi-disciplinary teams of engineers or scientists. To ensure that such interconnected models can be applied for the purpose of simulation and/or control, a reduced-order model for the interconnected dynamics is needed. In the scope of this paper, we pursue this goal by subsystem reduction to warrant modularity of the reduction approach. Clearly, by reducing the complexity of the subsystem models, not only the accuracy of the subsystem models is affected, but, consequently, also the accuracy of the interconnected model. It is typically difficult to predict a priori how the interconnected model accuracy is affected precisely by the subsystem reduction. In this work, we address this challenge by introducing a top-down approach which enables the translation of given accuracy requirements of the interconnected system model to accuracy requirements at subsystem model level, by using mathematical tools from robust performance analysis. This allows for the independent reduction of subsystem models while guaranteeing the desired accuracy of the interconnected model. In addition, we show how this top-down approach can be used to significantly reduce the interconnected model in an illustrative structural dynamics case study. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Tail-dependence, exceedance sets, and metric embeddings

Authors: Anja Janßen, Sebastian Neblung, Stilian Stoev

Abstract: There are many ways of measuring and modeling tail-dependence in random vectors: from the general framework of multivariate regular variation and the flexible class of max-stable vectors down to simple and concise summary measures like the matrix of bivariate tail-dependence coefficients. This paper starts by providing a review of existing results from a unifying perspective, which highlights conn… ▽ More There are many ways of measuring and modeling tail-dependence in random vectors: from the general framework of multivariate regular variation and the flexible class of max-stable vectors down to simple and concise summary measures like the matrix of bivariate tail-dependence coefficients. This paper starts by providing a review of existing results from a unifying perspective, which highlights connections between extreme value theory and the theory of cuts and metrics. Our approach leads to some new findings in both areas with some applications to current topics in risk management. We begin by using the framework of multivariate regular variation to show that extremal coefficients, or equivalently, the higher-order tail-dependence coefficients of a random vector can simply be understood in terms of random exceedance sets, which allows us to extend the notion of Bernoulli compatibility. In the special but important case of bi-variate tail-dependence, we establish a correspondence between tail-dependence matrices and $L^1$- and $\ell_1$-embeddable finite metric spaces via the spectral distance, which is a metric on the space of jointly $1$-Fréchet random variables. Namely, the coefficients of the cut-decomposition of the spectral distance and of the Tawn-Molchanov max-stable model realizing the corresponding bi-variate extremal dependence coincide. We show that line metrics are rigid and if the spectral distance corresponds to a line metric, the higher order tail-dependence is determined by the bi-variate tail-dependence matrix. Finally, the correspondence between $\ell_1$-embeddable metric spaces and tail-dependence matrices allows us to revisit the realizability problem, i.e. checking whether a given matrix is a valid tail-dependence matrix. We confirm a conjecture of Shyamalkumar & Tao (2020) that this problem is NP-complete. △ Less

Submitted 2 December, 2022; originally announced December 2022.

MSC Class: 60G70 (Primary) 51K05; 60E05; 68R12; 68Q25 (Secondary)

Modular Model Reduction of Interconnected Systems: A Robust Performance Analysis Perspective

Authors: Lars A. L. Janssen, Bart Besselink, Rob H. B. Fey, Nathan van de Wouw

Abstract: Many complex engineering systems consist of multiple subsystems that are developed by different teams of engineers. To analyse, simulate and control such complex systems, accurate yet computationally efficient models are required. Modular model reduction, in which the subsystem models are reduced individually, is a practical and an efficient method to obtain accurate reduced-order models of such c… ▽ More Many complex engineering systems consist of multiple subsystems that are developed by different teams of engineers. To analyse, simulate and control such complex systems, accurate yet computationally efficient models are required. Modular model reduction, in which the subsystem models are reduced individually, is a practical and an efficient method to obtain accurate reduced-order models of such complex systems. However, when subsystems are reduced individually, without taking their interconnections into account, the effect on stability and accuracy of the resulting reduced-order interconnected system is difficult to predict. In this work, a mathematical relation between the accuracy of reduced-order linear-time invariant subsystem models and (stability and accuracy of) resulting reduced-order interconnected linear time-invariant model is introduced. This result can subsequently be used in two ways. Firstly, it can be used to translate accuracy characteristics of the reduced-order subsystem models directly to accuracy properties of the interconnected reduced-order model. Secondly, it can also be used to translate specifications on the interconnected system model accuracy to accuracy requirements on subsystem models that can be used for fit-for-purpose reduction of the subsystem models. These applications of the proposed analysis framework for modular model reduction are demonstrated on an illustrative structural dynamics example. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Invariance properties of limiting point processes and applications to clusters of extremes

Authors: Anja Janßen, Johan Segers

Abstract: Motivated by examples from extreme value theory we introduce the general notion of a cluster process as a limiting point process of returns of a certain event in a time series. We explore general invariance properties of cluster processes which are implied by stationarity of the underlying time series under minimal assumptions. Of particular interest are the cluster size distributions, where we in… ▽ More Motivated by examples from extreme value theory we introduce the general notion of a cluster process as a limiting point process of returns of a certain event in a time series. We explore general invariance properties of cluster processes which are implied by stationarity of the underlying time series under minimal assumptions. Of particular interest are the cluster size distributions, where we introduce the two notions of inspected and typical cluster sizes and derive general properties of and connections between them. While the extremal index commonly used in extreme value theory is often interpreted as the inverse of a "mean cluster size", we point out that this only holds true for the expected value of the typical cluster size, caused by an effect very similar to the inspection paradox in renewal theory. △ Less

Submitted 2 November, 2023; v1 submitted 22 July, 2022; originally announced July 2022.

MSC Class: 60G10; 60G70

The MICADO Atmospheric Dispersion Corrector: Optomechanical design, expected performance and calibration techniques

Authors: J. A. van den Born, R. Romp, A. W. Janssen, R. Navarro, W. Jellema, E. Tolstoy, B. Jayawardhana, M. Hartl

Abstract: The differential refraction of light passing through the atmosphere can have a severe impact on image quality if no atmospheric dispersion corrector (ADC) is used. For the Extremely Large Telescope (ELT) this holds true well into the infrared. MICADO, the near-infrared imaging camera for the ELT, will employ a cryogenic ADC consisting of two counter-rotating Amici prisms with diameters of 125 mm.… ▽ More The differential refraction of light passing through the atmosphere can have a severe impact on image quality if no atmospheric dispersion corrector (ADC) is used. For the Extremely Large Telescope (ELT) this holds true well into the infrared. MICADO, the near-infrared imaging camera for the ELT, will employ a cryogenic ADC consisting of two counter-rotating Amici prisms with diameters of 125 mm. The mechanism will reduce the atmospheric dispersion to below 2.5 milli arcseconds (mas), with a set goal of 1 mas. In this report, we provide an overview of the current status of the ADC in development for MICADO. We summarise the optomechanical design and discuss how the cryogenic environment impacts the performance. We will also discuss our plan to use a diffraction mask in the cold pupil to calibrate and validate the performance once the instrument is fully integrated. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 16 pages, 12 figures, submitted to Proceedings of SPIE Astronomical Telescopes & Instrumentation 2022

Asymptotic analysis of Emden-Fowler type equation with an application to power flow models

Authors: M. H. M. Christianen, A. J. E. M. Janssen, M. Vlasiou, B. Zwart

Abstract: Emden-Fowler type equations are nonlinear differential equations that appear in many fields such as mathematical physics, astrophysics and chemistry. In this paper, we perform an asymptotic analysis of a specific Emden-Fowler type equation that emerges in a queuing theory context as an approximation of voltages under a well-known power flow model. Thus, we place Emden-Fowler type equations in the… ▽ More Emden-Fowler type equations are nonlinear differential equations that appear in many fields such as mathematical physics, astrophysics and chemistry. In this paper, we perform an asymptotic analysis of a specific Emden-Fowler type equation that emerges in a queuing theory context as an approximation of voltages under a well-known power flow model. Thus, we place Emden-Fowler type equations in the context of electrical engineering. We derive properties of the continuous solution of this specific Emden-Fowler type equation and study the asymptotic behavior of its discrete analog. We conclude that the discrete analog has the same asymptotic behavior as the classical continuous Emden-Fowler type equation that we consider. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 33 pages, 3 figures

Heavy-traffic single-server queues and the transform method

Authors: M. A. A. Boon, A. J. E. M. Janssen, J. S. H. van Leeuwaarden

Abstract: Heavy-traffic limit theory deals with queues that operate close to criticality and face severe queueing times. Let $W$ denote the steady-state waiting time in the ${\rm GI}/{\rm G}/1$ queue. Kingman (1961) showed that $W$, when appropriately scaled, converges in distribution to an exponential random variable as the system's load approaches 1. The original proof of this famous result uses the trans… ▽ More Heavy-traffic limit theory deals with queues that operate close to criticality and face severe queueing times. Let $W$ denote the steady-state waiting time in the ${\rm GI}/{\rm G}/1$ queue. Kingman (1961) showed that $W$, when appropriately scaled, converges in distribution to an exponential random variable as the system's load approaches 1. The original proof of this famous result uses the transform method. Starting from the Laplace transform of the pdf of $W$ (Pollaczek's contour integral representation), Kingman showed convergence of transforms and hence weak convergence of the involved random variables. We apply and extend this transform method to obtain convergence of moments with error assessment. We also demonstrate how the transform method can be applied to so-called nearly deterministic queues in a Kingman-type and a Gaussian heavy-traffic regime. We demonstrate numerically the accuracy of the various heavy-traffic approximations. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Innovation Begets Innovation and Concentration: The Case of Upstream Oil & Gas in the North Sea

Authors: Michele Fioretti, Alessandro Iaria, Aljoscha Janssen, Clément Mazet-Sonilhac, Robert K. Perrons

Abstract: We investigate the effect of technology adoption on competition by leveraging a unique dataset on production, costs, and asset characteristics for North Sea upstream oil & gas companies. Relying on heterogeneity in the geological suitability of fields and a landmark decision of the Norwegian Supreme Court that increased the returns of capital investment in Norway relative to the UK, we show that t… ▽ More We investigate the effect of technology adoption on competition by leveraging a unique dataset on production, costs, and asset characteristics for North Sea upstream oil & gas companies. Relying on heterogeneity in the geological suitability of fields and a landmark decision of the Norwegian Supreme Court that increased the returns of capital investment in Norway relative to the UK, we show that technology adoption increases market concentration. Firms with prior technology-specific know-how specialize more in fields suitable for the same technology but also invest more in high-risk-high-return fields (e.g., ultra-deep recovery), diversifying their technology portfolio and ultimately gaining larger shares of the North Sea market. Our analyses illustrate how technology adoption can lead to market concentration both directly through specialization and indirectly via experimentation. △ Less

Submitted 27 May, 2022; v1 submitted 26 May, 2022; originally announced May 2022.

Comparison of stability regions for a line distribution network with stochastic load demands

Authors: M. H. M. Christianen, J. Cruise, A. J. E. M. Janssen, S. Shneer, M. Vlasiou, B. Zwart

Abstract: We compare stability regions for different power flow models in the process of charging electric vehicles (EVs) by considering their random arrivals, their stochastic demand for energy at charging stations, and the characteristics of the electricity distribution network. We assume the distribution network is a line with charging stations located on it. We consider the Distflow and the Linearized D… ▽ More We compare stability regions for different power flow models in the process of charging electric vehicles (EVs) by considering their random arrivals, their stochastic demand for energy at charging stations, and the characteristics of the electricity distribution network. We assume the distribution network is a line with charging stations located on it. We consider the Distflow and the Linearized Distflow models and we assume that EVs have an exponential charging requirement, that voltage drops on the distribution network stay under control and that the number of charging stations $N$ goes to infinity. We investigate the stability of utility-optimizing power allocations in large distribution networks for both power flow models by controlling the arrival rate of EVs to charging stations. For both power flow models, we show that to obtain stability, the maximum feasible arrival rate, i.e. stability region of vehicles is decaying as $1/N^2$, and the difference between those arrival rates is up to constants, which we compare explicitly. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: 40 pages, 2 figures

Cluster based inference for extremes of time series

Authors: Holger Drees, Anja Janßen, Sebastian Neblung

Abstract: We introduce a new type of estimator for the spectral tail process of a regularly varying time series. The approach is based on a characterizing invariance property of the spectral tail process, which is incorporated into the new estimator via a projection technique. We show uniform asymptotic normality of this estimator, both in the case of known and of unknown index of regular variation. In a si… ▽ More We introduce a new type of estimator for the spectral tail process of a regularly varying time series. The approach is based on a characterizing invariance property of the spectral tail process, which is incorporated into the new estimator via a projection technique. We show uniform asymptotic normality of this estimator, both in the case of known and of unknown index of regular variation. In a simulation study the new procedure shows a more stable performance than previously proposed estimators. △ Less

Submitted 15 March, 2021; originally announced March 2021.

MSC Class: 62G32; 62M10; 62G05; 60F17

MOSAIC: the high-multiplex and multi-IFU spectrograph for the ELT

Authors: Rubén Sánchez-Janssen, Francois Hammer, Simon Morris, Jean-Gabriel Cuby, Lex Kaper, Matthias Steinmetz, Jose Afonso, Beatriz Barbuy, Edwin Bergin, Alexis Finoguenov, Jesús Gallego, Susan Kassin, Christopher Miller, Goran Östlin, Laura Pentericci, Daniel Schaerer, Bodo Ziegler, Fanny Chemla, Gavin Dalton, Fatima De Frondat, Chris Evans, David Le Mignant, Mathieu Puech, Myriam Rodrigues, Sylvestre Taburet , et al. (18 additional authors not shown)

Abstract: MOSAIC is the planned multi-object spectrograph for the 39m Extremely Large Telescope (ELT). Conceived as a multi-purpose instrument, it offers both high multiplex and multi-IFU capabilities at a range of intermediate to high spectral resolving powers in the visible and the near-infrared. MOSAIC will enable unique spectroscopic surveys of the faintest sources, from the oldest stars in the Galaxy a… ▽ More MOSAIC is the planned multi-object spectrograph for the 39m Extremely Large Telescope (ELT). Conceived as a multi-purpose instrument, it offers both high multiplex and multi-IFU capabilities at a range of intermediate to high spectral resolving powers in the visible and the near-infrared. MOSAIC will enable unique spectroscopic surveys of the faintest sources, from the oldest stars in the Galaxy and beyond to the first populations of galaxies that completed the reionisation of the Universe--while simultaneously opening up a wide discovery space. In this contribution we present the status of the instrument ahead of Phase B, showcasing the key science cases as well as introducing the updated set of top level requirements and the adopted architecture. The high readiness level will allow MOSAIC to soon enter the construction phase, with the goal to provide the ELT community with a world-class MOS capability as soon as possible after the telescope first light. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: Accepted to SPIE

A Fourier optics approach to evaluate the astrometric performance of MICADO

Authors: J. A. van den Born, W. Jellema, R. Navarro, E. Tolstoy, B. Jayawardhana, A. W. Janssen

Abstract: We present our investigation into the impact of wavefront errors on high accuracy astrometry using Fourier Optics. MICADO, the upcoming near-IR imaging instrument for the Extremely Large Telescope, will offer capabilities for relative astrometry with an accuracy of 50 micro arcseconds (μas). Due to the large size of the point spread function (PSF) compared to the astrometric requirement, the detai… ▽ More We present our investigation into the impact of wavefront errors on high accuracy astrometry using Fourier Optics. MICADO, the upcoming near-IR imaging instrument for the Extremely Large Telescope, will offer capabilities for relative astrometry with an accuracy of 50 micro arcseconds (μas). Due to the large size of the point spread function (PSF) compared to the astrometric requirement, the detailed shape and position of the PSF on the detector must be well understood. Furthermore, because the atmospheric dispersion corrector of MICADO is a moving component within an otherwise mostly static instrument, it might not be sufficient to perform a simple pre-observation calibration. Therefore, we have built a Fourier Optics framework, allowing us to evaluate the small changes in the centroid position of the PSF as a function of wavefront error. For a complete evaluation, we model both the low order surface form errors, using Zernike polynomials, and the mid- and high-spatial frequencies, using Power Spectral Density analysis. The described work will then make it possible, performing full diffractive beam propagation, to assess the expected astrometric performance of MICADO. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 13 pages, 13 figures, to be submitted to the SPIE Astronomical Telescopes & Instrumentation 2020 conference

MOSAIC on the ELT: high-multiplex spectroscopy to unravel the physics of stars and galaxies from the dark ages to the present-day

Authors: F. Hammer, S. Morris, J. G. Cuby, L. Kaper, M. Steinmetz, J. Afonso, B. Barbuy, E. Bergin, A. Finogenov, J. Gallego, S. Kassin, C. Miller, G. Ostlin, L. Penterricci, D. Schaerer, B. Ziegler, F. Chemla, G. Dalton, F. De Frondat, C. Evans, D. Le Mignant, M. Puech, M. Rodrigues, R. Sanchez-Janssen, S. Taburet , et al. (14 additional authors not shown)

Abstract: The powerful combination of the cutting-edge multi-object spectrograph MOSAIC with the world largest telescope, the ELT, will allow us to probe deeper into the Universe than was possible. MOSAIC is an extremely efficient instrument in providing spectra for the numerous faint sources in the Universe, including the very first galaxies and sources of cosmic reionization. MOSAIC has a high multiplex i… ▽ More The powerful combination of the cutting-edge multi-object spectrograph MOSAIC with the world largest telescope, the ELT, will allow us to probe deeper into the Universe than was possible. MOSAIC is an extremely efficient instrument in providing spectra for the numerous faint sources in the Universe, including the very first galaxies and sources of cosmic reionization. MOSAIC has a high multiplex in the NIR and in the VIS, in addition to multi-Integral Field Units (Multi-IFUs) in NIR. As such it is perfectly suited to carry out an inventory of dark matter (from rotation curves) and baryons in the cool-warm gas phases in galactic haloes at z=3-4. MOSAIC will enable detailed maps of the intergalactic medium at z=3, the evolutionary history of dwarf galaxies during a Hubble time, the chemistry directly measured from stars up to several Mpc. Finally, it will measure all faint features seen in cluster gravitational lenses or in streams surrounding nearby galactic halos, providing MOSAIC to be a powerful instrument with an extremely large space of discoveries. The preliminary design of MOSAIC is expected to begin next year, and its level of readiness is already high, given the instrumental studies made by the team. △ Less

Submitted 12 February, 2021; v1 submitted 6 November, 2020; originally announced November 2020.

Comments: 8 pages, 4 Figures, Pre-edited version, to appear in the ESO Messenger No.182 - Quarter 1 2021 - Version identical to the Edited one

Permutation inference in factorial survival designs with the CASANOVA

Authors: Marc Ditzhaus, Arnold Janssen, Markus Pauly

Abstract: We propose inference procedures for general nonparametric factorial survival designs with possibly right-censored data. Similar to additive Aalen models, null hypotheses are formulated in terms of cumulative hazards. Thereby, deviations are measured in terms of quadratic forms in Nelson-Aalen-type integrals. Different to existing approaches this allows to work without restrictive model assumptions… ▽ More We propose inference procedures for general nonparametric factorial survival designs with possibly right-censored data. Similar to additive Aalen models, null hypotheses are formulated in terms of cumulative hazards. Thereby, deviations are measured in terms of quadratic forms in Nelson-Aalen-type integrals. Different to existing approaches this allows to work without restrictive model assumptions as proportional hazards. In particular, crossing survival or hazard curves can be detected without a significant loss of power. For a distribution-free application of the method, a permutation strategy is suggested. The resulting procedures' asymptotic validity as well as their consistency are proven and their small sample performances are analyzed in extensive simulations. Their applicability is finally illustrated by analyzing an oncology data set. △ Less

Submitted 24 June, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

AGN feedback in a galaxy merger: Multi-phase, galaxy-scale outflows including a fast molecular gas blob ~6 kpc away from IRAS F08572+3915

Authors: R. Herrera-Camus, A. Janssen, E. Sturm, D. Lutz, S. Veilleux, R. Davies, T. Shimizu, E. González-Alfonso, D. S. N. Rupke, L. Tacconi, R. Genzel, C. Cicone, R. Maiolino, A. Contursi, J. Graciá-Carpio

Abstract: To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, findin… ▽ More To understand the role that AGN feedback plays in galaxy evolution we need in-depth studies of the multi-phase structure and energetics of galaxy-wide outflows. In this work we present new, deep ($\sim$50 hr) NOEMA CO(1-0) line observations of the molecular gas in the powerful outflow driven by the AGN in the ultra-luminous infrared galaxy IRAS F08572+3915. We spatially resolve the outflow, finding that its most likely configuration is a wide-angle bicone aligned with the kinematic major axis of the rotation disk. The molecular gas in the wind reaches velocities up to approximately $\pm$1200 km s$^{-1}$ and transports nearly 20% of the molecular gas mass in the system. We detect a second outflow component located $\sim$6 kpc north-west from the galaxy moving away at $\sim$900 km s$^{-1}$, which could be the result of a previous episode of AGN activity. The total mass and energetics of the outflow, which includes contributions from the ionized, neutral, warm and cold molecular gas phases is strongly dominated by the cold molecular gas. In fact, the molecular mass outflow rate is higher than the star formation rate, even if we only consider the gas in the outflow that is fast enough to escape the galaxy, which accounts for about $\sim$40% of the total mass of the outflow. This results in an outflow depletion time for the molecular gas in the central $\sim$1.5 kpc region of only $\sim3$ Myr, a factor of $\sim2$ shorter than the depletion time by star formation activity. △ Less

Submitted 14 November, 2019; originally announced November 2019.

Comments: 17 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 635, A47 (2020)

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

Authors: D. Lutz, E. Sturm, A. Janssen, S. Veilleux, S. Aalto, C. Cicone, A. Contursi, R. I. Davies, C. Feruglio, J. Fischer, S. Garcia-Burillo, R. Genzel, E. González-Alfonso, J. Gracía-Carpio R. Herrera-Camus, R. Maiolino, A. Schruba, T. Shimizu, A. Sternberg, L. J. Tacconi, A. Weiß

Abstract: We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transit… ▽ More We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged) △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: 37 pages. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A134 (2020)

Sharp bound on the radial derivatives of the Zernike circle polynomials (disk polynomials)

Authors: A. J. E. M. Janssen

Abstract: We sharpen the bound $n^{2k}$ on the maximum modulus of the $k^{\rm th}$ radial derivative of the Zernike circle polynomials (disk polynomials) of degree $n$ to $n^2(n^2-1^2)\cdot ... \cdot(n^2-(k-1)^2)/2^k(1/2)_k$. This bound is obtained from a result of Koornwinder on the non-negativity of connection coefficients of the radial parts of the circle polynomials when expanded into a series of Chebys… ▽ More We sharpen the bound $n^{2k}$ on the maximum modulus of the $k^{\rm th}$ radial derivative of the Zernike circle polynomials (disk polynomials) of degree $n$ to $n^2(n^2-1^2)\cdot ... \cdot(n^2-(k-1)^2)/2^k(1/2)_k$. This bound is obtained from a result of Koornwinder on the non-negativity of connection coefficients of the radial parts of the circle polynomials when expanded into a series of Chebyshev polynomials of the first kind. The new bound is shown to be sharp for, for instance, Zernike circle polynomials of degree $n$ and azimuthal order $m$ when $m=O(\sqrt{n})$ by using an explicit expression for the connection coefficients in terms of squares of Jacobi polynomials evaluated at 0. Keywords: Zernike circle polynomial, disk polynomial, radial derivative, Chebyshev polynomial, connection coefficient, Gegenbauer polynomial. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Dependence correction of multiple tests with applications to sparsity

Authors: Marc Ditzhaus, Arnold Janssen

Abstract: The present paper establishes new multiple procedures for simultaneous testing of a large number of hypotheses under dependence. Special attention is devoted to experiments with rare false hypotheses. This sparsity assumption is typically for various genome studies when a portion of remarkable genes should be detected. The aim is to derive tests which control the false discovery rate (FDR) always… ▽ More The present paper establishes new multiple procedures for simultaneous testing of a large number of hypotheses under dependence. Special attention is devoted to experiments with rare false hypotheses. This sparsity assumption is typically for various genome studies when a portion of remarkable genes should be detected. The aim is to derive tests which control the false discovery rate (FDR) always at finite sample size. The procedures are compared for the set up of dependent and independent $p$-values. It turns out that the FDR bounds differ by a dependency factor which can be used as a correction quantity. We offer sparsity modifications and improved dependence tests which generalize the Benjamini-Yekutieli test and adaptive tests in the sense of Storey. As a byproduct, an early stopped test is presented in order to bound the number of rejections. The new procedures perform well for real genome data examples. △ Less

Submitted 30 September, 2019; originally announced September 2019.

MSC Class: 62G10

Multiscale Deep Neural Networks for Multiclass Tissue Classification of Histological Whole-Slide Images

Authors: Rune Wetteland, Kjersti Engan, Trygve Eftestøl, Vebjørn Kvikstad, Emilius A. M. Janssen

Abstract: Correct treatment of urothelial carcinoma patients is dependent on accurate grading and staging of the cancer tumour. This is determined manually by a pathologist by examining the histological whole-slide images (WSI). The large size of these images makes this a time-consuming and challenging task. The WSI contain a variety of tissue types, and a method for defining diagnostic relevant regions wou… ▽ More Correct treatment of urothelial carcinoma patients is dependent on accurate grading and staging of the cancer tumour. This is determined manually by a pathologist by examining the histological whole-slide images (WSI). The large size of these images makes this a time-consuming and challenging task. The WSI contain a variety of tissue types, and a method for defining diagnostic relevant regions would have several advantages for visualization as well as further input to automated diagnosis systems. We propose an automatic multiscale method for classification of tiles from WSI of urothelial carcinoma patients into six classes. Three architectures based on convolutional neural network (CNN) were tested: MONO-CNN (400x), DI-CNN (100x/400x) and TRI-CNN (25x/100x/400x). The preliminary results show that the two multiscale models performed significantly better than the mono-scale model, achieving an F1-score of 0.986, substantiating that utilising multiple scales in the model aids the classification accuracy. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Comments: International Conference on Medical Imaging with Deep Learning (MIDL), 2019 [arXiv:1907.08612]

Report number: MIDL/2019/ExtendedAbstract/SyebijnnYN

$k$-means clustering of extremes

Authors: Anja Janßen, Phyllis Wan

Abstract: The $k$-means clustering algorithm and its variant, the spherical $k$-means clustering, are among the most important and popular methods in unsupervised learning and pattern detection. In this paper, we explore how the spherical $k$-means algorithm can be applied in the analysis of only the extremal observations from a data set. By making use of multivariate extreme value analysis we show how it c… ▽ More The $k$-means clustering algorithm and its variant, the spherical $k$-means clustering, are among the most important and popular methods in unsupervised learning and pattern detection. In this paper, we explore how the spherical $k$-means algorithm can be applied in the analysis of only the extremal observations from a data set. By making use of multivariate extreme value analysis we show how it can be adopted to find "prototypes" of extremal dependence and we derive a consistency result for our suggested estimator. In the special case of max-linear models we show furthermore that our procedure provides an alternative way of statistical inference for this class of models. Finally, we provide data examples which show that our method is able to find relevant patterns in extremal observations and allows us to classify extremal events. △ Less

Submitted 25 May, 2019; v1 submitted 5 April, 2019; originally announced April 2019.

MSC Class: 62G32; 62H30; 60G70

Counting cliques and cycles in scale-free inhomogeneous random graphs

Authors: A. J. E. M. Janssen, Johan S. H. van Leeuwaarden, Seva Shneer

Abstract: Scale-free networks contain many small cliques and cycles. We model such networks as inhomogeneous random graphs with regularly varying infinite-variance weights. For these models, the number of cliques and cycles have exact integral expressions amenable to asymptotic analysis. We obtain various asymptotic descriptions for how the average number of cliques and cycles, of any size, grow with the ne… ▽ More Scale-free networks contain many small cliques and cycles. We model such networks as inhomogeneous random graphs with regularly varying infinite-variance weights. For these models, the number of cliques and cycles have exact integral expressions amenable to asymptotic analysis. We obtain various asymptotic descriptions for how the average number of cliques and cycles, of any size, grow with the network size. For the cycle asymptotics we invoke the theory of circulant matrices. △ Less

Submitted 11 December, 2018; originally announced December 2018.

On a minimum distance procedure for threshold selection in tail analysis

Authors: Holger Drees, Anja Janßen, Sidney I. Resnick, Tiandong Wang

Abstract: Power-law distributions have been widely observed in different areas of scientific research. Practical estimation issues include how to select a threshold above which observations follow a power-law distribution and then how to estimate the power-law tail index. A minimum distance selection procedure (MDSP) is proposed in Clauset et al. (2009) and has been widely adopted in practice, especially in… ▽ More Power-law distributions have been widely observed in different areas of scientific research. Practical estimation issues include how to select a threshold above which observations follow a power-law distribution and then how to estimate the power-law tail index. A minimum distance selection procedure (MDSP) is proposed in Clauset et al. (2009) and has been widely adopted in practice, especially in the analyses of social networks. However, theoretical justifications for this selection procedure remain scant. In this paper, we study the asymptotic behavior of the selected threshold and the corresponding power-law index given by the MDSP. We find that the MDSP tends to choose too high a threshold level and leads to Hill estimates with large variances and root mean squared errors for simulated data with Pareto-like tails. △ Less

Submitted 12 February, 2020; v1 submitted 15 November, 2018; originally announced November 2018.

MSC Class: 60G70; 62E20; 60G15; 62G30; 05C80

Journal ref: SIAM Journal on Mathematics of Data Science 2020 2:1, 75-102

The MICADO first light imager for the ELT: overview, operation, simulation

Authors: Richard Davies, João Alves, Yann Clénet, Florian Lang-Bardl, Harald Nicklas, Jörg-Uwe Pott, Roberto Ragazzoni, Eline Tolstoy, Paola Amico, Heiko Anwand-Heerwart, Santiago Barboza, Lothar Barl, Pierre Baudoz, Ralf Bender, Naidu Bezawada, Peter Bizenberger, Wilfried Boland, Piercarlo Bonifacio, Bruno Borgo, Tristan Buey, Frédéric Chapron, Fanny Chemla, Mathieu Cohen, Oliver Czoske, Vincent Deo , et al. (76 additional authors not shown)

Abstract: MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and m… ▽ More MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve. △ Less

Submitted 26 July, 2018; originally announced July 2018.

Comments: Proc SPIE 10702. SPIE's copyright notice: "Copyright 2018 Society of PhotoOptical Instrumentation Engineers. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modification of the content of the paper are prohibited."

Implementation and benchmarking of a crosstalk-free method for wavefront Zernike coefficients reconstruction using Shack-Hartmann sensor data

Authors: R. S. Beisheuvel, A. J. E. M. Janssen, P. Pozzi, S. F. Pereira

Abstract: In wavefront characterization, often the combination of a Shack-Hartmann sensor and a reconstruction method utilizing the Cartesian derivatives of Zernike circle polynomials (the least-squares method, to be called here Method A) is used, which is known to introduce crosstalk. In \citep{janssen2014zernike} a crosstalk-free analytic expression of the LMS estimator of the wavefront Zernike coefficien… ▽ More In wavefront characterization, often the combination of a Shack-Hartmann sensor and a reconstruction method utilizing the Cartesian derivatives of Zernike circle polynomials (the least-squares method, to be called here Method A) is used, which is known to introduce crosstalk. In \citep{janssen2014zernike} a crosstalk-free analytic expression of the LMS estimator of the wavefront Zernike coefficients is given in terms of wavefront partial derivatives (leading to what we call Method B). Here, we show an implementation of this analytic result where the derivative data are obtained using the Shack-Hartmann sensor and compare it with the conventional least-squares method. △ Less

Submitted 9 July, 2018; originally announced July 2018.

Comments: The first author (R. S. Biesheuvel) has a new affiliation since this paper has been written: DEMCON Focal BV, Institutenweg 25A, 7521 PH Enschede, The Netherlands

The ELT-MOS (MOSAIC): towards the construction phase

Authors: Simon Morris, François Hammer, Pascal Jagourel, Christopher J. Evans, Mathieu Puech, Gavin B. Dalton, Myriam Rodrigues, Ruben Sanchez-Janssen, Ewan Fitzsimons, Beatriz Barbuy, Jean-Gabriel Cuby, Lex Kaper, Martin Roth, Gérard Rousset, Richard Myers, Olivier Le Fèvre, Alexis Finogenov, Jari Kotilainen, Bruno Castilho, Goran Ostlin, Sofia Feltzing, Andreas Korn, Jesus Gallego, África Castillo Morales, Jorge Iglesias-Páramo , et al. (28 additional authors not shown)

Abstract: When combined with the huge collecting area of the ELT, MOSAIC will be the most effective and flexible Multi-Object Spectrograph (MOS) facility in the world, having both a high multiplex and a multi-Integral Field Unit (Multi-IFU) capability. It will be the fastest way to spectroscopically follow-up the faintest sources, probing the reionisation epoch, as well as evaluating the evolution of the dw… ▽ More When combined with the huge collecting area of the ELT, MOSAIC will be the most effective and flexible Multi-Object Spectrograph (MOS) facility in the world, having both a high multiplex and a multi-Integral Field Unit (Multi-IFU) capability. It will be the fastest way to spectroscopically follow-up the faintest sources, probing the reionisation epoch, as well as evaluating the evolution of the dwarf mass function over most of the age of the Universe. MOSAIC will be world-leading in generating an inventory of both the dark matter (from realistic rotation curves with MOAO fed NIR IFUs) and the cool to warm-hot gas phases in z=3.5 galactic haloes (with visible wavelenth IFUs). Galactic archaeology and the first massive black holes are additional targets for which MOSAIC will also be revolutionary. MOAO and accurate sky subtraction with fibres have now been demonstrated on sky, removing all low Technical Readiness Level (TRL) items from the instrument. A prompt implementation of MOSAIC is feasible, and indeed could increase the robustness and reduce risk on the ELT, since it does not require diffraction limited adaptive optics performance. Science programmes and survey strategies are currently being investigated by the Consortium, which is also hoping to welcome a few new partners in the next two years. △ Less

Submitted 9 July, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 6 pages, 3 figures, Submitted to SPIE Proceedings, AS18

On the consistency of adaptive multiple tests

Authors: Marc Ditzhaus, Arnold Janssen

Abstract: Much effort has been done to control the "false discovery rate" (FDR) when $m$ hypotheses are tested simultaneously. The FDR is the expectation of the "false discovery proportion" $\text{FDP}=V/R$ given by the ratio of the number of false rejections $V$ and all rejections $R$. In this paper, we have a closer look at the FDP for adaptive linear step-up multiple tests. These tests extend the well kn… ▽ More Much effort has been done to control the "false discovery rate" (FDR) when $m$ hypotheses are tested simultaneously. The FDR is the expectation of the "false discovery proportion" $\text{FDP}=V/R$ given by the ratio of the number of false rejections $V$ and all rejections $R$. In this paper, we have a closer look at the FDP for adaptive linear step-up multiple tests. These tests extend the well known Benjamini and Hochberg test by estimating the unknown amount $m_0$ of the true null hypotheses. We give exact finite sample formulas for higher moments of the FDP and, in particular, for its variance. Using these allows us a precise discussion about the consistency of adaptive step-up tests. We present sufficient and necessary conditions for consistency on the estimators $\widehat m_0$ and the underlying probability regime. We apply our results to convex combinations of generalized Storey type estimators with various tuning parameters and (possibly) data-driven weights. The corresponding step-up tests allow a flexible adaptation. Moreover, these tests control the FDR at finite sample size. We compare these tests to the classical Benjamini and Hochberg test and discuss the advantages of it. △ Less

Submitted 8 January, 2018; originally announced January 2018.

MSC Class: 62G10; 62G20

Dynamic adaptive procedures that control the false discovery rate

Authors: Peter MacDonald, Kun Liang, Arnold Janssen

Abstract: In the multiple testing problem with independent tests, the classical linear step-up procedure controls the false discovery rate (FDR) at level $π_0α$, where $π_0$ is the proportion of true null hypotheses and $α$ is the target FDR level. Adaptive procedures can improve power by incorporating estimates of $π_0$, which typically rely on a tuning parameter. Fixed adaptive procedures set their tuning… ▽ More In the multiple testing problem with independent tests, the classical linear step-up procedure controls the false discovery rate (FDR) at level $π_0α$, where $π_0$ is the proportion of true null hypotheses and $α$ is the target FDR level. Adaptive procedures can improve power by incorporating estimates of $π_0$, which typically rely on a tuning parameter. Fixed adaptive procedures set their tuning parameters before seeing the data and can be shown to control the FDR in finite samples. We develop theoretical results for dynamic adaptive procedures whose tuning parameters are determined by the data. We show that, if the tuning parameter is chosen according to a left-to-right stopping time rule, the corresponding dynamic adaptive procedure controls the FDR in finite samples. Examples include the recently proposed right-boundary procedure and the widely used lowest-slope procedure, among others. Simulation results show that the right-boundary procedure is more powerful than other dynamic adaptive procedures under independence and mild dependence conditions. △ Less

Submitted 28 August, 2019; v1 submitted 6 December, 2017; originally announced December 2017.

Comments: To appear in Electronic Journal of Statistics

Spitzer's identity for discrete random walks

Authors: A. J. E. M. Janssen, Johan S. H. van Leeuwaarden

Abstract: Spitzer's identity describes the position of a reflected random walk over time in terms of a bivariate transform. Among its many applications in probability theory are congestion levels in queues and random walkers in physics. We present a new derivation of Spitzer's identity under the assumption that the increments of the random walk have bounded jumps to the left. This mild assumption facilitate… ▽ More Spitzer's identity describes the position of a reflected random walk over time in terms of a bivariate transform. Among its many applications in probability theory are congestion levels in queues and random walkers in physics. We present a new derivation of Spitzer's identity under the assumption that the increments of the random walk have bounded jumps to the left. This mild assumption facilitates a proof of Spitzer's identity that only uses basic properties of analytic functions and contour integration. The main novelty, believed to be of broader interest, is a reversed approach that recognizes a factored polynomial expression as the outcome of Cauchy's formula. △ Less

Submitted 26 October, 2017; originally announced October 2017.

Detectability of nonparametric signals: higher criticism versus likelihood ratio

Authors: Marc Ditzhaus, Arnold Janssen

Abstract: We study the signal detection problem in high dimensional noise data (possibly) containing rare and weak signals. Log-likelihood ratio (LLR) tests depend on unknown parameters, but they are needed to judge the quality of detection tests since they determine the detection regions. The popular Tukey's higher criticism (HC) test was shown to achieve the same completely detectable region as the LLR te… ▽ More We study the signal detection problem in high dimensional noise data (possibly) containing rare and weak signals. Log-likelihood ratio (LLR) tests depend on unknown parameters, but they are needed to judge the quality of detection tests since they determine the detection regions. The popular Tukey's higher criticism (HC) test was shown to achieve the same completely detectable region as the LLR test does for different (mainly) parametric models. We present a novel technique to prove this result for very general signal models, including even nonparametric $p$-value models. Moreover, we address the following questions which are still pending since the initial paper of Donoho and Jin: What happens on the border of the completely detectable region, the so-called detection boundary? Does HC keep its optimality there? In particular, we give a complete answer for the heteroscedastic normal mixture model. As a byproduct, we give some new insights about the LLR test's behavior on the detection boundary by discussing, among others, Pitmans's asymptotic efficiency as an application of Le Cam's theory. △ Less

Submitted 7 August, 2018; v1 submitted 21 September, 2017; originally announced September 2017.

MSC Class: Primary 62G10; 62G20; secondary 62G32

Design and Processing of Invertible Orientation Scores of 3D Images for Enhancement of Complex Vasculature

Authors: M. H. J. Janssen, A. J. E. M. Janssen, E. J. Bekkers, J. Olivan Bescos, R. Duits

Abstract: The enhancement and detection of elongated structures in noisy image data is relevant for many biomedical imaging applications. To handle complex crossing structures in 2D images, 2D orientation scores $U: \mathbb{R} ^ 2\times S ^ 1 \rightarrow \mathbb{C}$ were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores… ▽ More The enhancement and detection of elongated structures in noisy image data is relevant for many biomedical imaging applications. To handle complex crossing structures in 2D images, 2D orientation scores $U: \mathbb{R} ^ 2\times S ^ 1 \rightarrow \mathbb{C}$ were introduced, which already showed their use in a variety of applications. Here we extend this work to 3D orientation scores $U: \mathbb{R} ^ 3 \times S ^ 2\rightarrow \mathbb{C}$. First, we construct the orientation score from a given dataset, which is achieved by an invertible coherent state type of transform. For this transformation we introduce 3D versions of the 2D cake-wavelets, which are complex wavelets that can simultaneously detect oriented structures and oriented edges. Here we introduce two types of cake-wavelets, the first uses a discrete Fourier transform, the second is designed in the 3D generalized Zernike basis, allowing us to calculate analytical expressions for the spatial filters. Finally, we show two applications of the orientation score transformation. In the first application we propose an extension of crossing-preserving coherence enhancing diffusion via our invertible orientation scores of 3D images which we apply to real medical image data. In the second one we develop a new tubularity measure using 3D orientation scores and apply the tubularity measure to both artificial and real medical data. △ Less

Submitted 27 November, 2017; v1 submitted 7 July, 2017; originally announced July 2017.

MSC Class: 62H35; 65T60; 58J65; 37L05

ISM conditions in z~0.2 Lyman-Break Analogs

Authors: A. Contursi, A. J. Baker, S. Berta, B. Magnelli, D. Lutz, J. Fischer, A. Verma, M. Nielbock, J. Grácia Carpio, S. Veilleux, E. Sturm, R. Davies, R. Genzel, S. Hailey-Dunsheath, R. Herrera-Camus, A. Janssen, A. Poglitsch, A. Sternberg, L. J. Tacconi

Abstract: We present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar mass… ▽ More We present an analysis of far--infrared (FIR) [CII] and [OI] fine structure line and continuum observations obtained with $Herschel$/PACS, and CO(1-0) observations obtained with the IRAM Plateau de Bure Interferometer, of Lyman Break Analogs (LBAs) at $z\sim 0.2$. The principal aim of this work is to determine the typical ISM properties of $z\sim 1-2$ Main Sequence (MS) galaxies, with stellar masses between $10^{9.5}$ and $10^{11}$ $M_{\odot}$, which are currently not easily detectable in all these lines even with ALMA and NOEMA. We perform PDR modeling and apply different IR diagnostics to derive the main physical parameters of the FIR emitting gas and dust and we compare the derived ISM properties to those of galaxies on and above the MS at different redshifts. We find that the ISM properties of LBAs are quite extreme (low gas temperature, high density and thermal pressure) with respect to those found in local normal spirals and more active local galaxies. LBAs have no [CII] deficit despite having the high specific star formation rates (sSFRs) typical of starbursts. Although LBAs lie above the local MS, we show that their ISM properties are more similar to those of high-redshift MS galaxies than of local galaxies above the main sequence. This data set represents an important reference for planning future ALMA [CII] observations of relatively low-mass MS galaxies at the epoch of the peak of the cosmic star formation. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 19 pages, 12 Figures,8 Tables, Accepted for publication in A&A

Clustering Spectrum of scale-free networks

Authors: Clara Stegehuis, Remco van der Hofstad, Johan S. H. van Leeuwaarden, A. J. E. M Janssen

Abstract: Real-world networks often have power-law degrees and scale-free properties such as ultra-small distances and ultra-fast information spreading. In this paper, we study a third universal property: three-point correlations that suppress the creation of triangles and signal the presence of hierarchy. We quantify this property in terms of $\bar c(k)$, the probability that two neighbors of a degree-$k$… ▽ More Real-world networks often have power-law degrees and scale-free properties such as ultra-small distances and ultra-fast information spreading. In this paper, we study a third universal property: three-point correlations that suppress the creation of triangles and signal the presence of hierarchy. We quantify this property in terms of $\bar c(k)$, the probability that two neighbors of a degree-$k$ node are neighbors themselves. We investigate how the clustering spectrum $k\mapsto\bar c(k)$ scales with $k$ in the hidden variable model and show that $c(k)$ follows a {\it universal curve} that consists of three $k$-ranges where $\bar c(k)$ remains flat, starts declining, and eventually settles on a power law $\bar c(k)\sim k^{-α}$ with $α$ depending on the power law of the degree distribution. We test these results against ten contemporary real-world networks and explain analytically why the universal curve properties only reveal themselves in large networks. △ Less

Submitted 5 October, 2017; v1 submitted 6 June, 2017; originally announced June 2017.

Journal ref: Phys. Rev. E 96, 042309 (2017)