Source author record

Michael Rosenblum

Michael Rosenblum appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Combining Covariate Adjustment with Group Sequential, Information Adaptive Designs to Improve Randomized Trial Efficiency

In clinical trials, there is potential to improve precision and reduce the required sample size by appropriately adjusting for baseline variables in the statistical analysis. This is called covariate adjustment. Despite recommendations by regulatory agencies in favor of covariate adjustment, it remains underutilized leading to inefficient trials. We address two obstacles that make it challenging to use covariate adjustment. A first obstacle is the incompatibility of many covariate adjusted estimators with commonly used boundaries in group sequential designs (GSDs). A second obstacle is the uncertainty at the design stage about how much precision gain will result from covariate adjustment. We propose a method that modifies the original estimator so that it becomes compatible with GSDs, while increasing or leaving unchanged the estimator's precision. Our approach allows the use of any asymptotically linear estimator, which covers many estimators used in randomized trials. Building on this, we propose using an information adaptive design, that is, continuing the trial until the required information level is achieved. Such a design adapts to the amount of precision gain and can lead to faster, more efficient trials, without sacrificing validity or power. We evaluate estimator performance in simulations that mimic features of a completed stroke trial.

preprint2022arXiv

Using Targeted Maximum Likelihood Estimation to Estimate Treatment Effect with Longitudinal Continuous or Binary Data: A Systematic Evaluation of 28 Diabetes Clinical Trials

The primary analysis of clinical trials in diabetes therapeutic area often involves a mixed-model repeated measure (MMRM) approach to estimate the average treatment effect for longitudinal continuous outcome, and a generalized linear mixed model (GLMM) approach for longitudinal binary outcome. In this paper, we considered another estimator of the average treatment effect, called targeted maximum likelihood estimator (TMLE). This estimator can be a one-step alternative to model either continuous or binary outcome. We compared those estimators by simulation studies and by analyzing real data from 28 diabetes clinical trials. The simulations involved different missing data scenarios, and the real data sets covered a wide range of possible distributions of the outcome and covariates in real-life clinical trials for diabetes drugs with different mechanisms of action. For all the settings, adjusted estimators tended to be more efficient than the unadjusted one. In the setting of longitudinal continuous outcome, the MMRM approach with visits and baseline variables interaction appeared to dominate the performance of the MMRM considering the main effects only for the baseline variables while showing better or comparable efficiency to the TMLE estimator in both simulations and data applications. For modeling longitudinal binary outcome, TMLE generally outperformed GLMM in terms of relative efficiency, and its avoidance of the cumbersome covariance fitting procedure from GLMM makes TMLE a more advantageous estimator.

preprint2020arXiv

High-Order Phase Reduction for Coupled Oscillators

We explore the phase reduction in networks of coupled oscillators in the higher orders of the coupling parameter. For coupled Stuart-Landau oscillators, where the phase can be introduced explicitly, we develop an analytic perturbation procedure to allow for the obtaining of the higher-order approximation explicitly. We demonstrate this by deriving the second-order phase equations for a network of three Stuart-Landau oscillators. For systems where explicit expressions of the phase are not available, we present a numerical procedure that constructs the phase dynamics equations for a small network of coupled units. We apply this approach to a network of three van der Pol oscillators and reveal components in the coupling with different scaling in the interaction strength.

preprint2020arXiv

Model-Robust Inference for Clinical Trials that Improve Precision by Stratified Randomization and Covariate Adjustment

Two commonly used methods for improving precision and power in clinical trials are stratified randomization and covariate adjustment. However, many trials do not fully capitalize on the combined precision gains from these two methods, which can lead to wasted resources in terms of sample size and trial duration. We derive consistency and asymptotic normality of model-robust estimators that combine these two methods, and show that these estimators can lead to substantial gains in precision and power. Our theorems cover a class of estimators that handle continuous, binary, and time-to-event outcomes; missing outcomes under the missing at random assumption are handled as well. For each estimator, we give a formula for a consistent variance estimator that is model-robust and that fully captures variance reductions from stratified randomization and covariate adjustment. Also, we give the first proof (to the best of our knowledge) of consistency and asymptotic normality of the Kaplan-Meier estimator under stratified randomization, and we derive its asymptotic variance. The above results also hold for the biased-coin covariate-adaptive design. We demonstrate our results using three completed, phase 3, randomized trial data sets of treatments for substance use disorder, where the variance reduction due to stratified randomization and covariate adjustment ranges from 1% to 36%.

preprint2020arXiv

Reinforcement learning for suppression of collective activity in oscillatory ensembles

We present a use of modern data-based machine learning approaches to suppress self-sustained collective oscillations typically signaled by ensembles of degenerative neurons in the brain. The proposed hybrid model relies on two major components: an environment of oscillators and a policy-based reinforcement learning block. We report a model-agnostic synchrony control based on proximal policy optimization and two artificial neural networks in an Actor-Critic configuration. A class of physically meaningful reward functions enabling the suppression of collective oscillatory mode is proposed. The synchrony suppression is demonstrated for two models of neuronal populations -- for the ensembles of globally coupled limit-cycle Bonhoeffer-van der Pol oscillators and for the bursting Hindmarsh--Rose neurons.

preprint2016arXiv

Dynamics of weakly inhomogeneous oscillator populations: Perturbation theory on top of Watanabe-Strogatz integrability

As has been shown by Watanabe and Strogatz (WS) [Phys. Rev. Lett., 70, 2391 (1993)], a population of identical phase oscillators, sine-coupled to a common field, is a partially integrable system for any size: its dynamics reduces to equations for several collective variables. Here we develop a perturbation approach for weakly nonidentical ensembles. We calculate corrections to the WS dynamics for two types of perturbations: due to a distribution of natural frequencies and of forcing terms, and due to small white noise. We demonstrate, that in both cases the complex mean field for which the dynamical equations are written, is close up to the leading order in the perturbation to the Kuramoto order parameter. This supports validity of the dynamical reduction suggested by Ott and Antonsen [Chaos, 18, 037113 (2008)] for weakly inhomogeneous populations.

preprint2016arXiv

Improved Precision in the Analysis of Randomized Trials with Survival Outcomes, without Assuming Proportional Hazards

We present a new estimator of the restricted mean survival time in randomized trials where there is right censoring that may depend on treatment and baseline variables. The proposed estimator leverages prognostic baseline variables to obtain equal or better asymptotic precision compared to traditional estimators. Under regularity conditions and random censoring within strata of treatment and baseline variables, the proposed estimator has the following features: (i) it is interpretable under violations of the proportional hazards assumption; (ii) it is consistent and at least as precise as the Kaplan-Meier estimator under independent censoring; (iii) it remains consistent under violations of independent censoring (unlike the Kaplan-Meier estimator) when either the censoring or survival distributions are estimated consistently; and (iv) it achieves the nonparametric efficiency bound when both of these distributions are consistently estimated. We illustrate the performance of our method using simulations based on resampling data from a completed, phase 3 randomized clinical trial of a new surgical treatment for stroke; the proposed estimator achieves a 12% gain in relative efficiency compared to the Kaplan-Meier estimator. The proposed estimator has potential advantages over existing approaches for randomized trials with time-to-event outcomes, since existing methods either rely on model assumptions that are untenable in many applications, or lack some of the efficiency and consistency properties (i)-(iv). We focus on estimation of the restricted mean survival time, but our methods may be adapted to estimate any treatment effect measure defined as a smooth contrast between the survival curves for each study arm. We provide R code to implement the estimator.

preprint2015arXiv

On the equivalence of phase-oscillator and integrate-and-fire models

A quantitative comparison of various classes of oscillators (integrate-and-fire, Winfree, and Kuramoto-Daido type) is performed in the weak-coupling limit for a fully connected network of identical units. An almost perfect agreement is found, with only tiny differences among the models. We also show that the regime of self-consistent partial synchronization [SCPS] is rather general and can be observed for arbitrarily small coupling strength in any model class. As a by-product of our study, we are able to show that an integrate-and-fire model with a generic pulse shape can be always transformed into a similar model with $δ$-pulses and a suitable phase response curve.

preprint2014arXiv

Chimera-like states in an ensemble of globally coupled oscillators

We demonstrate emergence of a complex state in a homogeneous ensemble of globally coupled identical oscillators, reminiscent of chimera states in locally coupled oscillator lattices. In this regime some part of the ensemble forms a regularly evolving cluster, while all other units irregularly oscillate and remain asynchronous. We argue that chimera emerges because of effective bistability which dynamically appears in the originally monostable system due to internal delayed feedback in individual units. Additionally, we present two examples of chimeras in bistable systems with frequency-dependent phase shift in the global coupling.

preprint2014arXiv

Effective phase connectivity from observations

We present a novel approach for recovery of the directional connectivity of a small oscillator network by means of the phase dynamics reconstruction from multivariate time series data. The main idea is to use a triplet analysis instead of the traditional pair-wise one. Our technique reveals effective phase connectivity which is generally not equivalent to structural one. We demonstrate that by comparing the coupling functions from all possible triplets of oscillators, we are able to achieve in the reconstruction a good separation between existing and non-existing connections, and thus reliably reproduce the network structure.

preprint2014arXiv

interAdapt -- An Interactive Tool for Designing and Evaluating Randomized Trials with Adaptive Enrollment Criteria

The interAdapt R package is designed to be used by statisticians and clinical investigators to plan randomized trials. It can be used to determine if certain adaptive designs offer tangible benefits compared to standard designs, in the context of investigators' specific trial goals and constraints. Specifically, interAdapt compares the performance of trial designs with adaptive enrollment criteria versus standard (non-adaptive) group sequential trial designs. Performance is compared in terms of power, expected trial duration, and expected sample size. Users can either work directly in the R console, or with a user-friendly shiny application that requires no programming experience. Several added features are available when using the shiny application. For example, the application allows users to immediately download the results of the performance comparison as a csv-table, or as a printable, html-based report.

preprint2014arXiv

Solitary state at the edge of synchrony in ensembles with attractive and repulsive interaction

We discuss the desynchronization transition in networks of globally coupled identical oscillators with attractive and repulsive interactions. We show that, if attractive and repulsive groups act in antiphase or close to that, a solitary state emerges with a single repulsive oscillator split up from the others fully synchronized. With further increase of the repulsing strength, the synchronized cluster becomes fuzzy and the dynamics is given by a variety of stationary states with zero common forcing. Intriguingly, solitary states represent the natural link between coherence and incoherence. The phenomenon is described analytically for phase oscillators with sine coupling and demonstrated numerically for more general amplitude models.

preprint2014arXiv

Targeted Maximum Likelihood Estimation using Exponential Families

Targeted maximum likelihood estimation (TMLE) is a general method for estimating parameters in semiparametric and nonparametric models. Each iteration of TMLE involves fitting a parametric submodel that targets the parameter of interest. We investigate the use of exponential families to define the parametric submodel. This implementation of TMLE gives a general approach for estimating any smooth parameter in the nonparametric model. A computational advantage of this approach is that each iteration of TMLE involves estimation of a parameter in an exponential family, which is a convex optimization problem for which software implementing reliable and computationally efficient methods exists. We illustrate the method in three estimation problems, involving the mean of an outcome missing at random, the parameter of a median regression model, and the causal effect of a continuous exposure, respectively. We conduct a simulation study comparing different choices for the parametric submodel, focusing on the first of these problems. To the best of our knowledge, this is the first study investigating robustness of TMLE to different specifications of the parametric submodel. We find that the choice of submodel can have an important impact on the behavior of the estimator in finite samples.

preprint2013arXiv

Optimal Tests of Treatment Effects for the Overall Population and Two Subpopulations in Randomized Trials, using Sparse Linear Programming

We propose new, optimal methods for analyzing randomized trials, when it is suspected that treatment effects may differ in two predefined subpopulations. Such sub-populations could be defined by a biomarker or risk factor measured at baseline. The goal is to simultaneously learn which subpopulations benefit from an experimental treatment, while providing strong control of the familywise Type I error rate. We formalize this as a multiple testing problem and show it is computationally infeasible to solve using existing techniques. Our solution involves a novel approach, in which we first transform the original multiple testing problem into a large, sparse linear program. We then solve this problem using advanced optimization techniques. This general method can solve a variety of multiple testing problems and decision theory problems related to optimal trial design, for which no solution was previously available. In particular, we construct new multiple testing procedures that satisfy minimax and Bayes optimality criteria. For a given optimality criterion, our new approach yields the optimal tradeoff? between power to detect an effect in the overall population versus power to detect effects in subpopulations. We demonstrate our approach in examples motivated by two randomized trials of new treatments for HIV.

preprint2013arXiv

Synchrony suppression in ensembles of coupled oscillators via adaptive vanishing feedback

Synchronization and emergence of a collective mode is a general phenomenon, frequently observed in ensembles of coupled self-sustained oscillators of various natures. In several circumstances, in particular in cases of neurological pathologies, this state of the active medium is undesirable. Destruction of this state by a specially designed stimulation is a challenge of high clinical relevance. Typically, the precise effect of an external action on the ensemble is unknown, since the microscopic description of the oscillators and their interactions are not available. We show, that desynchronization in case of a large degree of uncertainty about important features of the system is nevertheless possible; it can by achieved by virtue of a feedback loop with an additional adaptation of parameters. The adaptation also ensures desynchronization of ensembles with non-stationary, time-varying parameters. We perform the stability analysis of the feedback-controlled system and demonstrate efficient destruction of synchrony for several models, including those of spiking and bursting neurons.

preprint2011arXiv

Optimal Phase Description of Chaotic Oscillators

We introduce an optimal phase description of chaotic oscillations by generalizing the concept of isochrones. On chaotic attractors possessing a general phase description, we define the optimal isophases as Poincaré surfaces showing return times as constant as possible. The dynamics of the resultant optimal phase is maximally decoupled of the amplitude dynamics, and provides a proper description of phase resetting of chaotic oscillations. The method is illustrated with the Rössler and Lorenz systems.

preprint2011arXiv

Reconstructing phase dynamics of oscillator networks

We generalize our recent approach to reconstruction of phase dynamics of coupled oscillators from data [B. Kralemann et al., Phys. Rev. E, 77, 066205 (2008)] to cover the case of small networks of coupled periodic units. Starting from the multivariate time series, we first reconstruct genuine phases and then obtain the coupling functions in terms of these phases. The partial norms of these coupling functions quantify directed coupling between oscillators. We illustrate the method by different network motifs for three coupled oscillators and for random networks of five and nine units. We also discuss nonlinear effects in coupling.

preprint2010arXiv

Collective Phase Chaos in the Dynamics of Interacting Oscillator Ensembles

We study chaotic behavior of order parameters in two coupled ensembles of self-sustained oscillators. Coupling within each of these ensembles is switched on and off alternately, while the mutual interaction between these two subsystems is arranged through quadratic nonlinear coupling. We show numerically that in the course of alternating Kuramoto transitions to synchrony and back to asynchrony, the exchange of excitations between two subpopulations proceeds in such a way that their collective phases are governed by an expanding circle map similar to the Bernoulli map. We perform the Lyapunov analysis of the dynamics and discuss finite-size effects.

preprint2010arXiv

Partially integrable dynamics of ensembles of nonidentical oscillators

We consider ensembles of sine-coupled phase oscillators consisting of subpopulations of identical units, with a general heterogeneous coupling between subpopulations. Using the Watanabe-Strogatz ansatz we reduce the dynamics of the ensemble to a relatively small number of dynamical variables plus microscopic constants of motion. This reduction is independent of the sizes of subpopulations and remains valid in the thermodynamic limits, where these sizes or/and the number of subpopulations are infinite. We demonstrate that the approach to the dynamics of such systems, recently proposed by Ott and Antonsen, corresponds to a particular choice of microscopic constants of motion. The theory is applied to the standard Kuramoto model and to the description of two interacting subpopulations, exhibiting a chimera state. Furthermore, we analyze the dynamics of the extension of the Kuramoto model for the case of nonlinear coupling and demonstrate the multistability of synchronous states.