Biometrical Journal最新文献

In laboratory medicine, due to the lack of sample availability and resources, measurements of many quantities of interest are commonly collected over a few samples, making statistical inference particularly challenging. In this context, several hypotheses can be tested, and studies are not often powered accordingly. We present a semiparametric Bayesian approach to effectively test multiple hypotheses applied to an experiment that aims to identify cytokines involved in Crohn's disease (CD) infection that may be ongoing in multiple tissues. We assume that the positive correlation commonly observed between cytokines is caused by latent groups of effects, which in turn result from a common cause. These clusters are effectively modeled through a Dirichlet Process (DP) that is one of the most popular choices as nonparametric prior in Bayesian statistics and has been proven to be a powerful tool for model-based clustering. We use a spike–slab distribution as the base measure of the DP. The nonparametric part has been included in an additive model whose parametric component is a Bayesian hierarchical model. We include simulations that empirically demonstrate the effectiveness of the proposed testing procedure in settings that mimic our application's sample size and data structure. Our CD data analysis shows strong evidence of a cytokine gradient in the external intestinal tissue.

Bayesian inference and the use of posterior or posterior predictive probabilities for decision making have become increasingly popular in clinical trials. The current practice in Bayesian clinical trials relies on a hybrid Bayesian-frequentist approach where the design and decision criteria are assessed with respect to frequentist operating characteristics such as power and type I error rate conditioning on a given set of parameters. These operating characteristics are commonly obtained via simulation studies. The utility of Bayesian measures, such as “assurance,” that incorporate uncertainty about model parameters in estimating the probabilities of various decisions in trials has been demonstrated. However, the computational burden remains an obstacle toward wider use of such criteria. In this article, we propose methodology which utilizes large sample theory of the posterior distribution to define parametric models for the sampling distribution of the posterior summaries used for decision making. The parameters of these models are estimated using a small number of simulation scenarios, thereby refining these models to capture the sampling distribution for small to moderate sample size. The proposed approach toward the assessment of conditional and marginal operating characteristics and sample size determination can be considered as simulation-assisted rather than simulation-based. It enables formal incorporation of uncertainty about the trial assumptions via a design prior and significantly reduces the computational burden for the design of Bayesian trials in general.

Landmarking is an alternative to complex multistate models when the aim is to calculate dynamic predictions. We develop the concept of landmarking for the case of left truncation and competing risks from the application background of drug safety assessment in pregnancy. The method is illustrated with a cohort study of the German Embryotox Pharmacovigilance Institute in Berlin to assess if the risk or the cumulative incidence of adverse pregnancy outcomes, like spontaneous abortions (SABs), is increased in fluoroquinolone-exposed women. Furthermore, we conduct an extensive simulation study to compare the dynamic predictions and coefficient estimates obtained by landmarking to those from nonparametric multistate models and classical time-dependent covariate Cox regression. The results from the simulation study indicate that attenuation of the effects is present in the landmark estimates, also in the complex setting considered here, but the estimates are still close to those from the multistate models. Regarding the Berlin fluoroquinolone data, the fluoroquinolone exposure of a pregnant woman in the first trimester seems to increase her cumulative incidence of elective termination of pregnancy over women never exposed before, but there is no evidence of a significantly increased risk or cumulative incidence in exposed women for SABs. This supports previous results on the same data, which were driven from an analysis without landmarking methods.

The modified Poisson and least-squares regression analyses for binary outcomes have been widely used as effective multivariable analysis methods to provide risk ratio and risk difference estimates in clinical and epidemiological studies. However, there is no certain evidence that assessed their operating characteristics under small and sparse data settings and no effective methods have been proposed for these regression analyses to address this issue. In this article, we show that the modified Poisson regression provides seriously biased estimates under small and sparse data settings. In addition, the modified least-squares regression provides unbiased estimates under these settings. We further show that the ordinary robust variance estimators for both of the methods have certain biases under situations that involve small or moderate sample sizes. To address these issues, we propose the Firth-type penalized methods for the modified Poisson and least-squares regressions. The adjustment methods lead to a more accurate and stable risk ratio estimator under small and sparse data settings, although the risk difference estimator is not invariant. In addition, to improve the inferences of the effect measures, we provide an improved robust variance estimator for these regression analyses. We conducted extensive simulation studies to assess the performances of the proposed methods under real-world conditions and found that the accuracies of the point and interval estimations were markedly improved by the proposed methods. We illustrate the effectiveness of these methods by applying them to a clinical study of epilepsy.

The analysis of multiple time-to-event outcomes in a randomized controlled clinical trial can be accomplished with existing methods. However, depending on the characteristics of the disease under investigation and the circumstances in which the study is planned, it may be of interest to conduct interim analyses and adapt the study design if necessary. Due to the expected dependency of the endpoints, the full available information on the involved endpoints may not be used for this purpose. We suggest a solution to this problem by embedding the endpoints in a multistate model. If this model is Markovian, it is possible to take the disease history of the patients into account and allow for data-dependent design adaptations. To this end, we introduce a flexible test procedure for a variety of applications, but are particularly concerned with the simultaneous consideration of progression-free survival (PFS) and overall survival (OS). This setting is of key interest in oncological trials. We conduct simulation studies to determine the properties for small sample sizes and demonstrate an application based on data from the NB2004-HR study.

The understanding of species interactions and ecosystem dynamics hinges upon the study of ecological niches. Quantifying the overlap of Hutchinsonian-niches has garnered significant attention, with many recent publications addressing the issue. Prior work on estimating niche overlap often did not provide confidence intervals or assumed multivariate normality, seriously limiting applications in ecology, and biodiversity research. This paper extends a nonparametric approach, previously applied to the two-species case, to multiple species. For estimation, a consistent plug-in estimator based on rank sums is proposed and its asymptotic distribution is derived under weak conditions. The novel methodology is then applied to a study comparing the ecological niches of the Eurasian eagle owl, common buzzard, and red kite. These species share a habitat in Central Europe but exhibit distinct population trends. The analysis explores their breeding habitat preferences, considering the intricate competition dynamics and utilizing the nonparametric approach to niche overlap estimation. Our proposed method provides a valuable inferential tool for the quantitative evaluation of differences and overlap between niches.

In this work, a method to regularize Cox frailty models is proposed that accommodates time-varying covariates and time-varying coefficients and is based on the full likelihood instead of the partial likelihood. A particular advantage of this framework is that the baseline hazard can be explicitly modeled in a smooth, semiparametric way, for example, via P-splines. Regularization for variable selection is performed via a lasso penalty and via group lasso for categorical variables while a second penalty regularizes wiggliness of smooth estimates of time-varying coefficients and the baseline hazard. Additionally, adaptive weights are included to stabilize the estimation. The method is implemented in the R function coxlasso, which is now integrated into the package PenCoxFrail, and will be compared to other packages for regularized Cox regression.