Corrected kriging update formulae for batch-sequential data assimilation

Recently, a lot of effort has been spent in the efficient computation of kriging predictors when observations are assimilated sequentially. In particular, kriging update formulae enabling significant computational savings were derived. Taking advantage of the previous kriging mean and variance computations helps avoiding a costly matrix inversion when adding one observation to the TeX already available ones. In addition to traditional update formulae taking into account a single new observation, Emery (2009) proposed formulae for the batch-sequential case, i.e. when TeX new observations are simultaneously assimilated. However, the kriging variance and covariance formulae given in Emery (2009) for the batch-sequential case are not correct. In this paper, we fix this issue and establish correct expressions for updated kriging variances and covariances when assimilating observations in parallel. An application in sequential conditional simulation finally shows that coupling update and residual substitution approaches may enable significant speed-ups.

A sequential Cox approach for estimating the causal effect of treatment in the presence of time-dependent confounding applied to data from the Swiss HIV Cohort Study

When estimating the effect of treatment on HIV using data from observational studies, standard methods may produce biased estimates due to the presence of time-dependent confounders. Such confounding can be present when a covariate, affected by past exposure, is both a predictor of the future exposure and the outcome. One example is the CD4 cell count, being a marker for disease progression for HIV patients, but also a marker for treatment initiation and influenced by treatment. Fitting a marginal structural model (MSM) using inverse probability weights is one way to give appropriate adjustment for this type of confounding. In this paper we study a simple and intuitive approach to estimate similar treatment effects, using observational data to mimic several randomized controlled trials. Each 'trial' is constructed based on individuals starting treatment in a certain time interval. An overall effect estimate for all such trials is found using composite likelihood inference. The method offers an alternative to the use of inverse probability of treatment weights, which is unstable in certain situations. The estimated parameter is not identical to the one of an MSM, it is conditioned on covariate values at the start of each mimicked trial. This allows the study of questions that are not that easily addressed fitting an MSM. The analysis can be performed as a stratified weighted Cox analysis on the joint data set of all the constructed trials, where each trial is one stratum. The model is applied to data from the Swiss HIV cohort study.

Sequential female assessment drives complex sexual selection on bower shape in a cichlid fish

In many animals, sexual selection on male traits results from female mate choice decisions made during a sequence of courtship behaviors. We use a bower-building cichlid fish, Nyassachromis cf. microcephalus, to show how applying standard selection analysis to data on sequential female assessment provides new insights into sexual selection by mate choice. We first show that the cumulative selection differentials confirm previous results suggesting female choice favors males holding large volcano-shaped sand bowers. The sequential assessment analysis reveals these cumulative differentials are the result of selection acting on different bower dimensions during the courtship sequence; females choose to follow males courting from tall bowers, but choose to engage in premating circling with males holding bowers with large diameter platforms. The approach we present extends standard selection analysis by partitioning the variances of increasingly accurate estimates of male reproductive fitness and is applicable to systems in which sequential female assessment drives sexual selection on male traits.

KrigInv: An efficient and user-friendly implementation of batch-sequential inversion strategies based on kriging

Several strategies relying on kriging have recently been proposed for adaptively estimating contour lines and excursion sets of functions under severely limited evaluation budget. The recently released R package KrigInv 3 is presented and offers a sound implementation of various sampling criteria for those kinds of inverse problems. KrigInv is based on the DiceKriging package, and thus benefits from a number of options concerning the underlying kriging models. Six implemented sampling criteria are detailed in a tutorial and illustrated with graphical examples. Different functionalities of KrigInv are gradually explained. Additionally, two recently proposed criteria for batch-sequential inversion are presented, enabling advanced users to distribute function evaluations in parallel on clusters or clouds of machines. Finally, auxiliary problems are discussed. These include the fine tuning of numerical integration and optimization procedures used within the computation and the optimization of the considered criteria.

The effectiveness and safety of rescue treatments in 108 patients with steroid-refractory ulcerative colitis with sequential rescue therapies in a subgroup of patients

BACKGROUND: Among patients with steroid-refractory ulcerative colitis (UC) in whom a first rescue therapy has failed, a second line salvage treatment can be considered to avoid colectomy. AIM: To evaluate the efficacy and safety of second or third line rescue therapy over a one-year period. METHODS: Response to single or sequential rescue treatments with infliximab (5mg/kg intravenously (iv) at week 0, 2, 6 and then every 8weeks), ciclosporin (iv 2mg/kg/daily and then oral 5mg/kg/daily) or tacrolimus (0.05mg/kg divided in 2 doses) in steroid-refractory moderate to severe UC patients from 7 Swiss and 1 Serbian tertiary IBD centers was retrospectively studied. The primary endpoint was the one year colectomy rate. RESULTS: 60% of patients responded to the first rescue therapy, 10% went to colectomy and 30% non-responders were switched to a 2(nd) line rescue treatment. 66% of patients responded to the 2(nd) line treatment whereas 34% failed, of which 15% went to colectomy and 19% received a 3(rd) line rescue treatment. Among those, 50% patients went to colectomy. Overall colectomy rate of the whole cohort was 18%. Steroid-free remission rate was 39%. The adverse event rates were 33%, 37.5% and 30% for the first, second and third line treatment respectively. CONCLUSION: Our data show that medical intervention even with 2(nd) and 3(rd) rescue treatments decreased colectomy frequency within one year of follow up. A longer follow-up will be necessary to investigate whether sequential therapy will only postpone colectomy and what percentage of patients will remain in long-term remission.

Self-esteem development from young adulthood to old age: A cohort-sequential longitudinal study

The authors examined the development of self-esteem from young adulthood to old age. Data came from the Americans’ Changing Lives study, which includes 4 assessments across a 16-year period of a nationally representative sample of 3,617 individuals aged 25 years to 104 years. Latent growth curve analyses indicated that self-esteem follows a quadratic trajectory across the adult life span, increasing during young and middle adulthood, reaching a peak at about age 60 years, and then declining in old age. No cohort differences in the self-esteem trajectory were found. Women had lower self-esteem than did men in young adulthood, but their trajectories converged in old age. Whites and Blacks had similar trajectories in young and middle adulthood, but the self-esteem of Blacks declined more sharply in old age than did the self-esteem of Whites. More educated individuals had higher self-esteem than did less educated individuals, but their trajectories were similar. Moreover, the results suggested that changes in socioeconomic status and physical health account for the decline in self-esteem that occurs in old age.