@article{HijaziFreitagLandwehr2023, author = {Hijazi, Saddam and Freitag, Melina A. and Landwehr, Niels}, title = {POD-Galerkin reduced order models and physics-informed neural networks for solving inverse problems for the Navier-Stokes equations}, series = {Advanced modeling and simulation in engineering sciences : AMSES}, volume = {10}, journal = {Advanced modeling and simulation in engineering sciences : AMSES}, number = {1}, publisher = {SpringerOpen}, address = {Berlin}, issn = {2213-7467}, doi = {10.1186/s40323-023-00242-2}, pages = {38}, year = {2023}, abstract = {We present a Reduced Order Model (ROM) which exploits recent developments in Physics Informed Neural Networks (PINNs) for solving inverse problems for the Navier-Stokes equations (NSE). In the proposed approach, the presence of simulated data for the fluid dynamics fields is assumed. A POD-Galerkin ROM is then constructed by applying POD on the snapshots matrices of the fluid fields and performing a Galerkin projection of the NSE (or the modified equations in case of turbulence modeling) onto the POD reduced basis. A POD-Galerkin PINN ROM is then derived by introducing deep neural networks which approximate the reduced outputs with the input being time and/or parameters of the model. The neural networks incorporate the physical equations (the POD-Galerkin reduced equations) into their structure as part of the loss function. Using this approach, the reduced model is able to approximate unknown parameters such as physical constants or the boundary conditions. A demonstration of the applicability of the proposed ROM is illustrated by three cases which are the steady flow around a backward step, the flow around a circular cylinder and the unsteady turbulent flow around a surface mounted cubic obstacle.}, language = {en} } @article{MicheletBindelliniMelinetal.2023, author = {Michelet, Robin and Bindellini, Davide and Melin, Johanna and Neumann, Uta and Blankenstein, Oliver and Huisinga, Wilhelm and Johnson, Trevor N. and Whitaker, Martin J. and Ross, Richard and Kloft, Charlotte}, title = {Insights in the maturational processes influencing hydrocortisone pharmacokinetics in congenital adrenal hyperplasia patients using a middle-out approach}, series = {Frontiers in Pharmacology}, volume = {13}, journal = {Frontiers in Pharmacology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {1663-9812}, doi = {10.3389/fphar.2022.1090554}, pages = {14}, year = {2023}, abstract = {Introduction: Hydrocortisone is the standard of care in cortisol replacement therapy for congenital adrenal hyperplasia patients. Challenges in mimicking cortisol circadian rhythm and dosing individualization can be overcome by the support of mathematical modelling. Previously, a non-linear mixed-effects (NLME) model was developed based on clinical hydrocortisone pharmacokinetic (PK) pediatric and adult data. Additionally, a physiologically-based pharmacokinetic (PBPK) model was developed for adults and a pediatric model was obtained using maturation functions for relevant processes. In this work, a middle-out approach was applied. The aim was to investigate whether PBPK-derived maturation functions could provide a better description of hydrocortisone PK inter-individual variability when implemented in the NLME framework, with the goal of providing better individual predictions towards precision dosing at the patient level. Methods: Hydrocortisone PK data from 24 adrenal insufficiency pediatric patients and 30 adult healthy volunteers were used for NLME model development, while the PBPK model and maturation functions of clearance and cortisol binding globulin (CBG) were developed based on previous studies published in the literature. Results: Clearance (CL) estimates from both approaches were similar for children older than 1 year (CL/F increasing from around 150 L/h to 500 L/h), while CBG concentrations differed across the whole age range (CBG(NLME) stable around 0.5 mu M vs. steady increase from 0.35 to 0.8 mu M for CBG (PBPK)). PBPK-derived maturation functions were subsequently included in the NLME model. After inclusion of the maturation functions, none, a part of, or all parameters were re-estimated. However, the inclusion of CL and/or CBG maturation functions in the NLME model did not result in improved model performance for the CL maturation function (\& UDelta;OFV > -15.36) and the re-estimation of parameters using the CBG maturation function most often led to unstable models or individual CL prediction bias. Discussion: Three explanations for the observed discrepancies could be postulated, i) non-considered maturation of processes such as absorption or first-pass effect, ii) lack of patients between 1 and 12 months, iii) lack of correction of PBPK CL maturation functions derived from urinary concentration ratio data for the renal function relative to adults. These should be investigated in the future to determine how NLME and PBPK methods can work towards deriving insights into pediatric hydrocortisone PK.}, language = {en} } @article{BusseSimonPetroffetal.2022, author = {Busse, David and Simon, Philipp and Petroff, David and El-Najjar, Nahed and Schmitt, Lisa and Bindellini, Davide and Dietrich, Arne and Zeitlinger, Markus and Huisinga, Wilhelm and Michelet, Robin and Wrigge, Hermann and Kloft, Charlotte}, title = {High-dosage fosfomycin results in adequate plasma and target-site exposure in morbidly obese and nonobese nonhyperfiltration patients}, series = {Antimicrobial agents and chemotherapy}, volume = {66}, journal = {Antimicrobial agents and chemotherapy}, number = {6}, publisher = {American Society for Microbiology}, address = {Washington}, issn = {0066-4804}, doi = {10.1128/aac.02302-21}, pages = {12}, year = {2022}, abstract = {The objectives of this study were the identification in (morbidly) obese and nonobese patients of (i) the most appropriate body size descriptor for fosfomycin dose adjustments and (ii) adequacy of the currently employed dosing regimens. Plasma and target site (interstitial fluid of subcutaneous adipose tissue) concentrations after fosfomycin administration (8 g) to 30 surgery patients (15 obese/15 nonobese) were obtained from a prospective clinical trial. After characterization of plasma and microdialysis-derived target site pharmacokinetics via population analysis, short-term infusions of fosfomycin 3 to 4 times daily were simulated. The adequacy of therapy was assessed by probability of pharmacokinetic/pharmacodynamic target attainment (PTA) analysis based on the unbound drug-related targets of an \%fT(>= MIC) (the fraction of time that unbound fosfomycin concentrations exceed the MIC during 24 h) of 70 and an fAUC(0-24h)/MIC (the area under the concentration-time curve from 0 to 24 h for the unbound fraction of fosfomycin relative to the MIC) of 40.8 to 83.3. Lean body weight, fat mass, and creatinine clearance calculated via adjusted body weight (ABW) (CLCRCG_ABW) of all patients (body mass index [BMI] = 20.1 to 52.0 kg/m(2)) explained a considerable proportion of between-patient pharmacokinetic variability (up to 31.0\% relative reduction). The steady-state unbound target site/plasma concentration ratio was 26.3\% lower in (morbidly) obese than nonobese patients. For infections with fosfomycin-susceptible pathogens (MIC <= 16 mg/L), intermittent "high-dosage" intravenous (i.v.) fosfomycin (8 g, three times daily) was sufficient to treat patients with a CLCRCG_ABW of,130 mL/min, irrespective of the pharmacokinetic/pharmacodynamic indices considered. For infections by Pseudomonas aeruginosa with a MIC of 32 mg/L, when the index fAUC0-24h/MIC is applied, fosfomycin might represent a promising treatment option in obese and nonobese patients, especially in combination therapy to complement beta-lactams, in which carbapenem-resistant P. aeruginosa is critical. In conclusion, fosfomycin showed excellent target site penetration in obese and nonobese patients. Dosing should be guided by renal function rather than obesity status.}, language = {en} } @article{WeineltStegemannTheloeetal.2022, author = {Weinelt, Ferdinand Anton and Stegemann, Miriam Songa and Theloe, Anja and Pf{\"a}fflin, Frieder and Achterberg, Stephan and Weber, Franz and D{\"u}bel, Lucas and Mikolajewska, Agata and Uhrig, Alexander and Kiessling, Peggy and Huisinga, Wilhelm and Michelet, Robin and Hennig, Stefanie and Kloft, Charlotte}, title = {Evaluation of a meropenem and piperacillin monitoring program in intensive care unit patients calls for the regular assessment of empirical targets and easy-to-use dosing decision tools}, series = {Antibiotics : open access journal}, volume = {11}, journal = {Antibiotics : open access journal}, number = {6}, publisher = {MDPI}, address = {Basel}, issn = {2079-6382}, doi = {10.3390/antibiotics11060758}, pages = {17}, year = {2022}, abstract = {The drug concentrations targeted in meropenem and piperacillin/tazobactam therapy also depend on the susceptibility of the pathogen. Yet, the pathogen is often unknown, and antibiotic therapy is guided by empirical targets. To reliably achieve the targeted concentrations, dosing needs to be adjusted for renal function. We aimed to evaluate a meropenem and piperacillin/tazobactam monitoring program in intensive care unit (ICU) patients by assessing (i) the adequacy of locally selected empirical targets, (ii) if dosing is adequately adjusted for renal function and individual target, and (iii) if dosing is adjusted in target attainment (TA) failure. In a prospective, observational clinical trial of drug concentrations, relevant patient characteristics and microbiological data (pathogen, minimum inhibitory concentration (MIC)) for patients receiving meropenem or piperacillin/tazobactam treatment were collected. If the MIC value was available, a target range of 1-5 x MIC was selected for minimum drug concentrations of both drugs. If the MIC value was not available, 8-40 mg/L and 16-80 mg/L were selected as empirical target ranges for meropenem and piperacillin, respectively. A total of 356 meropenem and 216 piperacillin samples were collected from 108 and 96 ICU patients, respectively. The vast majority of observed MIC values was lower than the empirical target (meropenem: 90.0\%, piperacillin: 93.9\%), suggesting empirical target value reductions. TA was found to be low (meropenem: 35.7\%, piperacillin 50.5\%) with the lowest TA for severely impaired renal function (meropenem: 13.9\%, piperacillin: 29.2\%), and observed drug concentrations did not significantly differ between patients with different targets, indicating dosing was not adequately adjusted for renal function or target. Dosing adjustments were rare for both drugs (meropenem: 6.13\%, piperacillin: 4.78\%) and for meropenem irrespective of TA, revealing that concentration monitoring alone was insufficient to guide dosing adjustment. Empirical targets should regularly be assessed and adjusted based on local susceptibility data. To improve TA, scientific knowledge should be translated into easy-to-use dosing strategies guiding antibiotic dosing.}, language = {en} } @article{PornsawadBoeckmannPanitsupakamon2022, author = {Pornsawad, Pornsarp and B{\"o}ckmann, Christine and Panitsupakamon, Wannapa}, title = {The Levenberg-Marquardt regularization for the backward heat equation with fractional derivative}, series = {Electronic transactions on numerical analysis - ETNA}, volume = {57}, journal = {Electronic transactions on numerical analysis - ETNA}, publisher = {Kent State University}, address = {Kent}, isbn = {978-3-7001-8258-0}, issn = {1068-9613}, doi = {10.1553/etna_vol57s67}, pages = {67 -- 79}, year = {2022}, abstract = {The backward heat problem with time-fractional derivative in Caputo's sense is studied. The inverse problem is severely ill-posed in the case when the fractional order is close to unity. A Levenberg-Marquardt method with a new a posteriori stopping rule is investigated. We show that optimal order can be obtained for the proposed method under a H{\"o}lder-type source condition. Numerical examples for one and two dimensions are provided.}, language = {en} } @article{StachanowNeumannBlankensteinetal.2022, author = {Stachanow, Viktoria and Neumann, Uta and Blankenstein, Oliver and Bindellini, Davide and Melin, Johanna and Ross, Richard and Whitaker, Martin J. J. and Huisinga, Wilhelm and Michelet, Robin and Kloft, Charlotte}, title = {Exploring dried blood spot cortisol concentrations as an alternative for monitoring pediatric adrenal insufficiency patients}, series = {Frontiers in pharmacology}, volume = {13}, journal = {Frontiers in pharmacology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {1663-9812}, doi = {10.3389/fphar.2022.819590}, pages = {8}, year = {2022}, abstract = {Congenital adrenal hyperplasia (CAH) is the most common form of adrenal insufficiency in childhood; it requires cortisol replacement therapy with hydrocortisone (HC, synthetic cortisol) from birth and therapy monitoring for successful treatment. In children, the less invasive dried blood spot (DBS) sampling with whole blood including red blood cells (RBCs) provides an advantageous alternative to plasma sampling. Potential differences in binding/association processes between plasma and DBS however need to be considered to correctly interpret DBS measurements for therapy monitoring. While capillary DBS samples would be used in clinical practice, venous cortisol DBS samples from children with adrenal insufficiency were analyzed due to data availability and to directly compare and thus understand potential differences between venous DBS and plasma. A previously published HC plasma pharmacokinetic (PK) model was extended by leveraging these DBS concentrations. In addition to previously characterized binding of cortisol to albumin (linear process) and corticosteroid-binding globulin (CBG; saturable process), DBS data enabled the characterization of a linear cortisol association with RBCs, and thereby providing a quantitative link between DBS and plasma cortisol concentrations. The ratio between the observed cortisol plasma and DBS concentrations varies highly from 2 to 8. Deterministic simulations of the different cortisol binding/association fractions demonstrated that with higher blood cortisol concentrations, saturation of cortisol binding to CBG was observed, leading to an increase in all other cortisol binding fractions. In conclusion, a mathematical PK model was developed which links DBS measurements to plasma exposure and thus allows for quantitative interpretation of measurements of DBS samples.}, language = {en} } @article{KellerLiuPeyerimhoff2021, author = {Keller, Matthias and Liu, Shiping and Peyerimhoff, Norbert}, title = {A note on eigenvalue bounds for non-compact manifolds}, series = {Mathematische Nachrichten}, volume = {294}, journal = {Mathematische Nachrichten}, number = {6}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {0025-584X}, doi = {10.1002/mana.201900209}, pages = {1134 -- 1139}, year = {2021}, abstract = {In this article we prove upper bounds for the Laplace eigenvalues lambda(k) below the essential spectrum for strictly negatively curved Cartan-Hadamard manifolds. Our bound is given in terms of k(2) and specific geometric data of the manifold. This applies also to the particular case of non-compact manifolds whose sectional curvature tends to -infinity, where no essential spectrum is present due to a theorem of Donnelly/Li. The result stands in clear contrast to Laplacians on graphs where such a bound fails to be true in general.}, language = {en} } @article{SaynischWagnerBaerenzungHornschildetal.2021, author = {Saynisch-Wagner, Jan and B{\"a}renzung, Julien and Hornschild, Aaron and Irrgang, Christopher and Thomas, Maik}, title = {Tide-induced magnetic signals and their errors derived from CHAMP and Swarm satellite magnetometer observations}, series = {Earth, planets and space : EPS}, volume = {73}, journal = {Earth, planets and space : EPS}, number = {1}, publisher = {Springer}, address = {Heidelberg}, issn = {1880-5981}, doi = {10.1186/s40623-021-01557-3}, pages = {11}, year = {2021}, abstract = {Satellite-measured tidal magnetic signals are of growing importance. These fields are mainly used to infer Earth's mantle conductivity, but also to derive changes in the oceanic heat content. We present a new Kalman filter-based method to derive tidal magnetic fields from satellite magnetometers: KALMAG. The method's advantage is that it allows to study a precisely estimated posterior error covariance matrix. We present the results of a simultaneous estimation of the magnetic signals of 8 major tides from 17 years of Swarm and CHAMP data. For the first time, robustly derived posterior error distributions are reported along with the reported tidal magnetic fields. The results are compared to other estimates that are either based on numerical forward models or on satellite inversions of the same data. For all comparisons, maximal differences and the corresponding globally averaged RMSE are reported. We found that the inter-product differences are comparable with the KALMAG-based errors only in a global mean sense. Here, all approaches give values of the same order, e.g., 0.09 nT-0.14 nT for M2. Locally, the KALMAG posterior errors are up to one order smaller than the inter-product differences, e.g., 0.12 nT vs. 0.96 nT for M2.}, language = {en} } @article{WormellReich2021, author = {Wormell, Caroline L. and Reich, Sebastian}, title = {Spectral convergence of diffusion maps}, series = {SIAM journal on numerical analysis / Society for Industrial and Applied Mathematics}, volume = {59}, journal = {SIAM journal on numerical analysis / Society for Industrial and Applied Mathematics}, number = {3}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {0036-1429}, doi = {10.1137/20M1344093}, pages = {1687 -- 1734}, year = {2021}, abstract = {Diffusion maps is a manifold learning algorithm widely used for dimensionality reduction. Using a sample from a distribution, it approximates the eigenvalues and eigenfunctions of associated Laplace-Beltrami operators. Theoretical bounds on the approximation error are, however, generally much weaker than the rates that are seen in practice. This paper uses new approaches to improve the error bounds in the model case where the distribution is supported on a hypertorus. For the data sampling (variance) component of the error we make spatially localized compact embedding estimates on certain Hardy spaces; we study the deterministic (bias) component as a perturbation of the Laplace-Beltrami operator's associated PDE and apply relevant spectral stability results. Using these approaches, we match long-standing pointwise error bounds for both the spectral data and the norm convergence of the operator discretization. We also introduce an alternative normalization for diffusion maps based on Sinkhorn weights. This normalization approximates a Langevin diffusion on the sample and yields a symmetric operator approximation. We prove that it has better convergence compared with the standard normalization on flat domains, and we present a highly efficient rigorous algorithm to compute the Sinkhorn weights.}, language = {en} } @article{ChangKhalilSchulze2021, author = {Chang, Der-Chen and Khalil, Sara and Schulze, Bert-Wolfgang}, title = {Analysis on regular corner spaces}, series = {The journal of geometric analysis}, volume = {31}, journal = {The journal of geometric analysis}, number = {9}, publisher = {Springer}, address = {New York}, issn = {1050-6926}, doi = {10.1007/s12220-021-00614-3}, pages = {9199 -- 9240}, year = {2021}, abstract = {We establish a new approach of treating elliptic boundary value problems (BVPs) on manifolds with boundary and regular corners, up to singularity order 2. Ellipticity and parametrices are obtained in terms of symbols taking values in algebras of BVPs on manifolds of corresponding lower singularity orders. Those refer to Boutet de Monvel's calculus of operators with the transmission property, see Boutet de Monvel (Acta Math 126:11-51, 1971) for the case of smooth boundary. On corner configuration operators act in spaces with multiple weights. We mainly study the case of upper left entries in the respective 2 x 2 operator block-matrices of such a calculus. Green operators in the sense of Boutet de Monvel (Acta Math 126:11-51, 1971) analogously appear in singular cases, and they are complemented by contributions of Mellin type. We formulate a result on ellipticity and the Fredholm property in weighted corner spaces, with parametrices of analogous kind.}, language = {en} } @article{YenvonSpechtLinetal.2022, author = {Yen, Ming-Hsuan and von Specht, Sebastian and Lin, Yen-Yu and Cotton, Fabrice and Ma, Kuo-Fong}, title = {Within- and between-event variabilities of strong-velocity pulses of moderate earthquakes within dense seismic arrays}, series = {Bulletin of the Seismological Society of America}, volume = {112}, journal = {Bulletin of the Seismological Society of America}, number = {1}, publisher = {Seismological Society of America}, address = {El Cerito, Calif.}, issn = {0037-1106}, doi = {10.1785/0120200376}, pages = {361 -- 380}, year = {2022}, abstract = {Ground motion with strong-velocity pulses can cause significant damage to buildings and structures at certain periods; hence, knowing the period and velocity amplitude of such pulses is critical for earthquake structural engineering. However, the physical factors relating the scaling of pulse periods with magnitude are poorly understood. In this study, we investigate moderate but damaging earthquakes (M-w 6-7) and characterize ground- motion pulses using the method of Shahi and Baker (2014) while considering the potential static-offset effects. We confirm that the within-event variability of the pulses is large. The identified pulses in this study are mostly from strike-slip-like earthquakes. We further perform simulations using the freq uency-wavenumber algorithm to investigate the causes of the variability of the pulse periods within and between events for moderate strike-slip earthquakes. We test the effect of fault dips, and the impact of the asperity locations and sizes. The simulations reveal that the asperity properties have a high impact on the pulse periods and amplitudes at nearby stations. Our results emphasize the importance of asperity characteristics, in addition to earthquake magnitudes for the occurrence and properties of pulses produced by the forward directivity effect. We finally quantify and discuss within- and between-event variabilities of pulse properties at short distances.}, language = {en} } @article{MolkenthinDonnerReichetal.2022, author = {Molkenthin, Christian and Donner, Christian and Reich, Sebastian and Z{\"o}ller, Gert and Hainzl, Sebastian and Holschneider, Matthias and Opper, Manfred}, title = {GP-ETAS: semiparametric Bayesian inference for the spatio-temporal epidemic type aftershock sequence model}, series = {Statistics and Computing}, volume = {32}, journal = {Statistics and Computing}, number = {2}, publisher = {Springer}, address = {Dordrecht}, issn = {0960-3174}, doi = {10.1007/s11222-022-10085-3}, pages = {25}, year = {2022}, abstract = {The spatio-temporal epidemic type aftershock sequence (ETAS) model is widely used to describe the self-exciting nature of earthquake occurrences. While traditional inference methods provide only point estimates of the model parameters, we aim at a fully Bayesian treatment of model inference, allowing naturally to incorporate prior knowledge and uncertainty quantification of the resulting estimates. Therefore, we introduce a highly flexible, non-parametric representation for the spatially varying ETAS background intensity through a Gaussian process (GP) prior. Combined with classical triggering functions this results in a new model formulation, namely the GP-ETAS model. We enable tractable and efficient Gibbs sampling by deriving an augmented form of the GP-ETAS inference problem. This novel sampling approach allows us to assess the posterior model variables conditioned on observed earthquake catalogues, i.e., the spatial background intensity and the parameters of the triggering function. Empirical results on two synthetic data sets indicate that GP-ETAS outperforms standard models and thus demonstrate the predictive power for observed earthquake catalogues including uncertainty quantification for the estimated parameters. Finally, a case study for the l'Aquila region, Italy, with the devastating event on 6 April 2009, is presented.}, language = {en} } @article{RodriguezZuluagaStolleYamazakietal.2021, author = {Rodr{\´i}guez Zuluaga, Juan and Stolle, Claudia and Yamazaki, Yosuke and Xiong, Chao and England, Scott L.}, title = {A synoptic-scale wavelike structure in the nighttime equatorial ionization anomaly}, series = {Earth and Space Science : ESS}, volume = {8}, journal = {Earth and Space Science : ESS}, number = {2}, publisher = {American Geophysical Union}, address = {Malden, Mass.}, issn = {2333-5084}, doi = {10.1029/2020EA001529}, pages = {10}, year = {2021}, abstract = {Both ground- and satellite-based airglow imaging have significantly contributed to understanding the low-latitude ionosphere, especially the morphology and dynamics of the equatorial ionization anomaly (EIA). The NASA Global-scale Observations of the Limb and Disk (GOLD) mission focuses on far-ultraviolet airglow images from a geostationary orbit at 47.5 degrees W. This region is of particular interest at low magnetic latitudes because of the high magnetic declination (i.e., about -20 degrees) and proximity of the South Atlantic magnetic anomaly. In this study, we characterize an exciting feature of the nighttime EIA using GOLD observations from October 5, 2018 to June 30, 2020. It consists of a wavelike structure of a few thousand kilometers seen as poleward and equatorward displacements of the EIA-crests. Initial analyses show that the synoptic-scale structure is symmetric about the dip equator and appears nearly stationary with time over the night. In quasi-dipole coordinates, maxima poleward displacements of the EIA-crests are seen at about +/- 12 degrees latitude and around 20 and 60 degrees longitude (i.e., in geographic longitude at the dip equator, about 53 degrees W and 14 degrees W). The wavelike structure presents typical zonal wavelengths of about 6.7 x 10(3) km and 3.3 x 10(3) km. The structure's occurrence and wavelength are highly variable on a day-to-day basis with no apparent dependence on geomagnetic activity. In addition, a cluster or quasi-periodic wave train of equatorial plasma depletions (EPDs) is often detected within the synoptic-scale structure. We further outline the difference in observing these EPDs from FUV images and in situ measurements during a GOLD and Swarm mission conjunction.}, language = {en} } @inproceedings{DemarisGrišićHuisingaetal.2020, author = {D{\´e}maris, Alise and Grišić, Ana-Marija and Huisinga, Wilhelm and Walter, Reinisch and Kloft, Charlotte}, title = {Evaluation of dosing strategies of anti-TNF alpha monoclonal antibodies using pharmacokinetic modelling and simulation}, series = {Journal of Crohn's and Colitis}, volume = {14}, booktitle = {Journal of Crohn's and Colitis}, number = {Supp. 1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {1873-9946}, doi = {10.1093/ecco-jcc/jjz203.201}, pages = {S171 -- S172}, year = {2020}, abstract = {Background: Anti-TNFα monoclonal antibodies (mAbs) are a well-established treatment for patients with Crohn's disease (CD). However, subtherapeutic concentrations of mAbs have been related to a loss of response during the first year of therapy1. Therefore, an appropriate dosing strategy is crucial to prevent the underexposure of mAbs for those patients. The aim of our study was to assess the impact of different dosing strategies (fixed dose or body size descriptor adapted) on drug exposure and the target concentration attainment for two different anti-TNFα mAbs: infliximab (IFX, body weight (BW)-based dosing) and certolizumab pegol (CZP, fixed dosing). For this purpose, a comprehensive pharmacokinetic (PK) simulation study was performed. Methods: A virtual population of 1000 clinically representative CD patients was generated based on the distribution of CD patient characteristics from an in-house clinical database (n = 116). Seven dosing regimens were investigated: fixed dose and per BW, lean BW (LBW), body surface area, height, body mass index and fat-free mass. The individual body size-adjusted doses were calculated from patient generated body size descriptor values. Then, using published PK models for IFX and CZP in CD patients2,3, for each patient, 1000 concentration-time profiles were simulated to consider the typical profile of a specific patient as well as the range of possible individual profiles due to unexplained PK variability across patients. For each dosing strategy, the variability in maximum and minimum mAb concentrations (Cmax and Cmin, respectively), area under the concentration-time curve (AUC) and the per cent of patients reaching target concentration were assessed during maintenance therapy. Results: For IFX and CZP, Cmin showed the highest variability between patients (CV ≈110\% and CV ≈80\%, respectively) with a similar extent across all dosing strategies. For IFX, the per cent of patients reaching the target (Cmin = 5 µg/ml) was similar across all dosing strategies (~15\%). For CZP, the per cent of patients reaching the target average concentration of 17 µg/ml ranged substantially (52-71\%), being the highest for LBW-adjusted dosing. Conclusion: By using a PK simulation approach, different dosing regimen of IFX and CZP revealed the highest variability for Cmin, the most commonly used PK parameter guiding treatment decisions, independent upon dosing regimen. Our results demonstrate similar target attainment with fixed dosing of IFX compared with currently recommended BW-based dosing. For CZP, the current fixed dosing strategy leads to comparable percentage of patients reaching target as the best performing body size-adjusted dosing (66\% vs. 71\%, respectively).}, language = {en} } @article{EshghiMachReichel2021, author = {Eshghi, Nasim and Mach, Thomas and Reichel, Lothar}, title = {New matrix function approximations and quadrature rules based on the Arnoldi process}, series = {Journal of computational and applied mathematics}, volume = {391}, journal = {Journal of computational and applied mathematics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0377-0427}, doi = {10.1016/j.cam.2021.113442}, pages = {12}, year = {2021}, abstract = {The Arnoldi process can be applied to inexpensively approximate matrix functions of the form f (A)v and matrix functionals of the form v*(f (A))*g(A)v, where A is a large square non-Hermitian matrix, v is a vector, and the superscript * denotes transposition and complex conjugation. Here f and g are analytic functions that are defined in suitable regions in the complex plane. This paper reviews available approximation methods and describes new ones that provide higher accuracy for essentially the same computational effort by exploiting available, but generally not used, moment information. Numerical experiments show that in some cases the modifications of the Arnoldi decompositions proposed can improve the accuracy of v*(f (A))*g(A)v about as much as performing an additional step of the Arnoldi process.}, language = {en} } @article{PohleAdamBeumer2022, author = {Pohle, Jennifer and Adam, Timo and Beumer, Larissa}, title = {Flexible estimation of the state dwell-time distribution in hidden semi-Markov models}, series = {Computational statistics \& data analysis}, volume = {172}, journal = {Computational statistics \& data analysis}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0167-9473}, doi = {10.1016/j.csda.2022.107479}, pages = {15}, year = {2022}, abstract = {Hidden semi-Markov models generalise hidden Markov models by explicitly modelling the time spent in a given state, the so-called dwell time, using some distribution defined on the natural numbers. While the (shifted) Poisson and negative binomial distribution provide natural choices for such distributions, in practice, parametric distributions can lack the flexibility to adequately model the dwell times. To overcome this problem, a penalised maximum likelihood approach is proposed that allows for a flexible and data-driven estimation of the dwell-time distributions without the need to make any distributional assumption. This approach is suitable for direct modelling purposes or as an exploratory tool to investigate the latent state dynamics. The feasibility and potential of the suggested approach is illustrated in a simulation study and by modelling muskox movements in northeast Greenland using GPS tracking data. The proposed method is implemented in the R-package PHSMM which is available on CRAN.}, language = {en} } @misc{KolbeEvans2020, author = {Kolbe, Benedikt Maximilian and Evans, Myfanwy E.}, title = {Isotopic tiling theory for hyperbolic surfaces}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1}, issn = {1866-8372}, doi = {10.25932/publishup-54428}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-544285}, pages = {30}, year = {2020}, abstract = {In this paper, we develop the mathematical tools needed to explore isotopy classes of tilings on hyperbolic surfaces of finite genus, possibly nonorientable, with boundary, and punctured. More specifically, we generalize results on Delaney-Dress combinatorial tiling theory using an extension of mapping class groups to orbifolds, in turn using this to study tilings of covering spaces of orbifolds. Moreover, we study finite subgroups of these mapping class groups. Our results can be used to extend the Delaney-Dress combinatorial encoding of a tiling to yield a finite symbol encoding the complexity of an isotopy class of tilings. The results of this paper provide the basis for a complete and unambiguous enumeration of isotopically distinct tilings of hyperbolic surfaces.}, language = {en} } @article{PathirajaLeeuwen2022, author = {Pathiraja, Sahani Darschika and Leeuwen, Peter Jan van}, title = {Multiplicative Non-Gaussian model error estimation in data assimilation}, series = {Journal of advances in modeling earth systems : JAMES}, volume = {14}, journal = {Journal of advances in modeling earth systems : JAMES}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1942-2466}, doi = {10.1029/2021MS002564}, pages = {23}, year = {2022}, abstract = {Model uncertainty quantification is an essential component of effective data assimilation. Model errors associated with sub-grid scale processes are often represented through stochastic parameterizations of the unresolved process. Many existing Stochastic Parameterization schemes are only applicable when knowledge of the true sub-grid scale process or full observations of the coarse scale process are available, which is typically not the case in real applications. We present a methodology for estimating the statistics of sub-grid scale processes for the more realistic case that only partial observations of the coarse scale process are available. Model error realizations are estimated over a training period by minimizing their conditional sum of squared deviations given some informative covariates (e.g., state of the system), constrained by available observations and assuming that the observation errors are smaller than the model errors. From these realizations a conditional probability distribution of additive model errors given these covariates is obtained, allowing for complex non-Gaussian error structures. Random draws from this density are then used in actual ensemble data assimilation experiments. We demonstrate the efficacy of the approach through numerical experiments with the multi-scale Lorenz 96 system using both small and large time scale separations between slow (coarse scale) and fast (fine scale) variables. The resulting error estimates and forecasts obtained with this new method are superior to those from two existing methods.}, language = {en} } @misc{KellerPinchoverPogorzelski2020, author = {Keller, Matthias and Pinchover, Yehuda and Pogorzelski, Felix}, title = {From hardy to rellich inequalities on graphs}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {3}, issn = {1866-8372}, doi = {10.25932/publishup-54214}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-542140}, pages = {22}, year = {2020}, abstract = {We show how to deduce Rellich inequalities from Hardy inequalities on infinite graphs. Specifically, the obtained Rellich inequality gives an upper bound on a function by the Laplacian of the function in terms of weighted norms. These weights involve the Hardy weight and a function which satisfies an eikonal inequality. The results are proven first for Laplacians and are extended to Schrodinger operators afterwards.}, language = {en} } @misc{MazzonettoSalimova2020, author = {Mazzonetto, Sara and Salimova, Diyora}, title = {Existence, uniqueness, and numerical approximations for stochastic burgers equations}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {4}, issn = {1866-8372}, doi = {10.25932/publishup-51579}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-515796}, pages = {26}, year = {2020}, abstract = {In this article, we propose an all-in-one statement which includes existence, uniqueness, regularity, and numerical approximations of mild solutions for a class of stochastic partial differential equations (SPDEs) with non-globally monotone nonlinearities. The proof of this result exploits the properties of an existing fully explicit space-time discrete approximation scheme, in particular the fact that it satisfies suitable a priori estimates. We also obtain almost sure and strong convergence of the approximation scheme to the mild solutions of the considered SPDEs. We conclude by applying the main result of the article to the stochastic Burgers equations with additive space-time white noise.}, language = {en} }