-
Study Duration Prediction for Clinical Trials with Time-to-Event Endpoints Using Mixture Distributions Accounting for Heterogeneous Population
Authors:
Hong Zhang,
Jie Pu,
Shibing Deng,
Satrajit Roychoudhury,
Haitao Chu,
Douglas Robinson
Abstract:
In the era of precision medicine, more and more clinical trials are now driven or guided by biomarkers, which are patient characteristics objectively measured and evaluated as indicators of normal biological processes, pathogenic processes, or pharmacologic responses to therapeutic interventions. With the overarching objective to optimize and personalize disease management, biomarker-guided clinic…
▽ More
In the era of precision medicine, more and more clinical trials are now driven or guided by biomarkers, which are patient characteristics objectively measured and evaluated as indicators of normal biological processes, pathogenic processes, or pharmacologic responses to therapeutic interventions. With the overarching objective to optimize and personalize disease management, biomarker-guided clinical trials increase the efficiency by appropriately utilizing prognostic or predictive biomarkers in the design. However, the efficiency gain is often not quantitatively compared to the traditional all-comers design, in which a faster enrollment rate is expected (e.g. due to no restriction to biomarker positive patients) potentially leading to a shorter duration. To accurately predict biomarker-guided trial duration, we propose a general framework using mixture distributions accounting for heterogeneous population. Extensive simulations are performed to evaluate the impact of heterogeneous population and the dynamics of biomarker characteristics and disease on the study duration. Several influential parameters including median survival time, enrollment rate, biomarker prevalence and effect size are identitied. Re-assessments of two publicly available trials are conducted to empirically validate the prediction accuracy and to demonstrate the practical utility. The R package \emph{detest} is developed to implement the proposed method and is publicly available on CRAN.
△ Less
Submitted 31 December, 2023;
originally announced January 2024.
-
A Portfolio Rebalancing Approach for the Indian Stock Market
Authors:
Jaydip Sen,
Arup Dasgupta,
Subhasis Dasgupta,
Sayantani Roychoudhury
Abstract:
This chapter presents a calendar rebalancing approach to portfolios of stocks in the Indian stock market. Ten important sectors of the Indian economy are first selected. For each of these sectors, the top ten stocks are identified based on their free-float market capitalization values. Using the ten stocks in each sector, a sector-specific portfolio is designed. In this study, the historical stock…
▽ More
This chapter presents a calendar rebalancing approach to portfolios of stocks in the Indian stock market. Ten important sectors of the Indian economy are first selected. For each of these sectors, the top ten stocks are identified based on their free-float market capitalization values. Using the ten stocks in each sector, a sector-specific portfolio is designed. In this study, the historical stock prices are used from January 4, 2021, to September 20, 2023 (NSE Website). The portfolios are designed based on the training data from January 4, 2021 to June 30, 2022. The performances of the portfolios are tested over the period from July 1, 2022, to September 20, 2023. The calendar rebalancing approach presented in the chapter is based on a yearly rebalancing method. However, the method presented is perfectly flexible and can be adapted for weekly or monthly rebalancing. The rebalanced portfolios for the ten sectors are analyzed in detail for their performances. The performance results are not only indicative of the relative performances of the sectors over the training (i.e., in-sample) data and test (out-of-sample) data, but they also reflect the overall effectiveness of the proposed portfolio rebalancing approach.
△ Less
Submitted 15 October, 2023;
originally announced October 2023.
-
Application of BadNets in Spam Filters
Authors:
Swagnik Roychoudhury,
Akshaj Kumar Veldanda
Abstract:
Spam filters are a crucial component of modern email systems, as they help to protect users from unwanted and potentially harmful emails. However, the effectiveness of these filters is dependent on the quality of the machine learning models that power them. In this paper, we design backdoor attacks in the domain of spam filtering. By demonstrating the potential vulnerabilities in the machine learn…
▽ More
Spam filters are a crucial component of modern email systems, as they help to protect users from unwanted and potentially harmful emails. However, the effectiveness of these filters is dependent on the quality of the machine learning models that power them. In this paper, we design backdoor attacks in the domain of spam filtering. By demonstrating the potential vulnerabilities in the machine learning model supply chain, we highlight the need for careful consideration and evaluation of the models used in spam filters. Our results show that the backdoor attacks can be effectively used to identify vulnerabilities in spam filters and suggest the need for ongoing monitoring and improvement in this area.
△ Less
Submitted 18 July, 2023;
originally announced July 2023.
-
On fixed and uncertain mixture prior weights
Authors:
Beat Neuenschwander,
Simon Wandel,
Satrajit Roychoudhury,
Heinz Schmidli
Abstract:
This paper focuses on the specification of the weights for the components of mixture priors.
This paper focuses on the specification of the weights for the components of mixture priors.
△ Less
Submitted 30 June, 2023; v1 submitted 27 June, 2023;
originally announced June 2023.
-
Duration of and time to response in oncology clinical trials from the perspective of the estimand framework
Authors:
Hans-Jochen Weber,
Stephen Corson,
Jiang Li,
Francois Mercier,
Satrajit Roychoudhury,
Martin Oliver Sailer,
Stephen Sun,
Alexander Todd,
Godwin Yung
Abstract:
Duration of response (DOR) and time to response (TTR) are typically evaluated as secondary endpoints in early-stage clinical studies in oncology when efficacy is assessed by the best overall response (BOR) and presented as the overall response rate (ORR). Despite common use of DOR and TTR in particular in single-arm studies, the definition of these endpoints and the questions they are intended to…
▽ More
Duration of response (DOR) and time to response (TTR) are typically evaluated as secondary endpoints in early-stage clinical studies in oncology when efficacy is assessed by the best overall response (BOR) and presented as the overall response rate (ORR). Despite common use of DOR and TTR in particular in single-arm studies, the definition of these endpoints and the questions they are intended to answer remain unclear. Motivated by the estimand framework, we present relevant scientific questions of interest for DOR and TTR and propose corresponding estimand definitions. We elaborate on how to deal with relevant intercurrent events which should follow the same considerations as implemented for the primary response estimand. A case study in mantle cell lymphoma illustrates the implementation of relevant estimands of DOR and TTR. We close the paper with practical recommendations to implement DOR and TTR in clinical study protocols.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Investigating the Electronic Structure of Prospective Water-splitting Oxide BaCe$_{0.25}$Mn$_{0.75}$O$_{3-δ}$ Before and After Thermal Reduction
Authors:
Subhayan Roychoudhury,
Sarah Shulda,
Anuj Goyal,
Robert Bell,
Sami Sainio,
Nicholas Strange,
James Eujin Park,
Eric N. Coker,
Stephan Lany,
David Ginley,
David Prendergast
Abstract:
BaCe$_{0.25}$Mn$_{0.75}$O$_{3-δ}$ (BCM), a non-stoichiometric oxide closely resembling a perovskite crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar energy, oxygen-vacancies can be created in BCM and the reduced crystal so obtained can, in turn, produce H2 by stripping oxygen f…
▽ More
BaCe$_{0.25}$Mn$_{0.75}$O$_{3-δ}$ (BCM), a non-stoichiometric oxide closely resembling a perovskite crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar energy, oxygen-vacancies can be created in BCM and the reduced crystal so obtained can, in turn, produce H2 by stripping oxygen from H2O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM, is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge x-ray absorption spectroscopy (XAS). The computed projected density-of-states (PDOS) and orbital-plots are used to propose a simplified model for orbital-mixing between the oxygen and the ligand atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and to categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that, as a consequence of dielectric screening, the change in electron-density resulting from the reduction is strongly localized around the oxygen vacancy. Our experimental studies reveal a marked lowering of the first O K-edge peak in the reduced crystal which is shown to result from a diminished O-2p contribution to the frontier unoccupied orbitals, in accordance with the tight-binding scheme. Our study paves the way for investigation of the working-mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
Excitonic effects in X-ray absorption spectra of fluoride salts and their surfaces
Authors:
Ana Sanz-Matias,
Subhayan Roychoudhury,
Xuefei Feng,
Feipeng Yang,
Li Cheng Kao,
Kevin R. Zavadil,
Jinghua Guo,
David Prendergast
Abstract:
Operando X-ray absorption spectroscopy (XAS) can probe the electrode-electrolyte interface with single-digit nanometer depth resolution and offers a wealth of insight into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure aboutthe excited at…
▽ More
Operando X-ray absorption spectroscopy (XAS) can probe the electrode-electrolyte interface with single-digit nanometer depth resolution and offers a wealth of insight into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure aboutthe excited atoms. To this end, we explore fluorine K-edge XAS of mono- (Li, Na, K) and di-valent (Mg, Ca, Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials, despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the delta-SCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced, despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.
△ Less
Submitted 23 September, 2022;
originally announced September 2022.
-
Efficient core-excited state orbital perspective on calculating X-ray absorption transitions in determinant framework
Authors:
Subhayan Roychoudhury,
David Prendergast
Abstract:
X-ray absorption spectroscopy (XAS) is an explicit probe of the unoccupied electronic structure of materials and an invaluable tool for fingerprinting various electronic properties and phenomena. Computational methods capable of simulating and analysing such spectra are therefore in high demand for complementing the experimental results and for extracting valuable insights therefrom. In particular…
▽ More
X-ray absorption spectroscopy (XAS) is an explicit probe of the unoccupied electronic structure of materials and an invaluable tool for fingerprinting various electronic properties and phenomena. Computational methods capable of simulating and analysing such spectra are therefore in high demand for complementing the experimental results and for extracting valuable insights therefrom. In particular, a recently proposed first-principles approach titled Many-Body XAS (MBXAS), which approximates the final (initial) state as a Slater determinant constructed from Kohn-Sham (KS) orbitals optimized in absence (presence) of the relevant core-hole has shown promising prospects in evaluating the transition amplitudes. In this article, we show that the MBXAS approach can be rederived using a transition operator expressed entirely in the basis of core-excited state KS orbitals and that this reformulation offers substantial practical and conceptual advantages. In addition to circumventing previous issues of convergence with respect to the number of unoccupied ground-state orbitals, the aforementioned representation reduces the computational expense by rendering the calculation of such orbitals unnecessary altogether. The reformulated approach also provides a direct pathway for comparing the many-body approximation with the so-called single-particle treatment and indicates the relative importance in observed XAS intensity of the relaxation of the valence occupied subspace induced by the core excitation. Finally, using the core-excited state basis, we define auxiliary orbitals for x-ray absorption and demonstrate their utility in explaining the spectral intensity by contrasting them with single-particle approximations to the excited state.
△ Less
Submitted 23 August, 2022;
originally announced August 2022.
-
CleaRIXS: A Fast and Accurate First-Principles Method for Simulation and Analysis of Resonant Inelastic X-ray Scattering
Authors:
Subhayan Roychoudhury,
David Prendergast
Abstract:
Resonant Inelastic X-ray Scattering (RIXS), which probes the occupied and unoccupied electronic subspaces in an interrelated fashion, is one of the most detailed, complex and information-rich experimental techniques employed in the investigation of electronic structure across physics, chemistry and materials science. We introduce an ab-initio, accurate and efficient computational framework for sim…
▽ More
Resonant Inelastic X-ray Scattering (RIXS), which probes the occupied and unoccupied electronic subspaces in an interrelated fashion, is one of the most detailed, complex and information-rich experimental techniques employed in the investigation of electronic structure across physics, chemistry and materials science. We introduce an ab-initio, accurate and efficient computational framework for simulation and analysis of RIXS spectra, which combines two diverse and established approaches to modeling electronic excited states. The Core-Hole Linear-Response RIXS (CleaRIXS) method not only ensures accurate incorporation of the interaction of electrons with core and valence holes, but also automatically maps the salient RIXS features to the relevant electronic excitations and de-excitations. Through direct comparison with previous methanol C K-edge RIXS measurements [J. Phys. Chem. A 120, 2260 (2016)], we show the efficacy of the formalism in modeling different regions of the RIXS spectrum and in gaining physical insight regarding their origins. The importance of including the valence electron-hole interactions is explored, in addition to the connection between CleaRIXS and determinant-based approaches for simulating X-ray absorption and non-resonant X-ray emission. A qualitative discussion of additional RIXS features arising from possible vibrationally-mediated processes is also presented. CleaRIXS provides a robust and extendable framework for prediction and interpretation of RIXS processes and for the simulation of complex electronic excited states in general.
△ Less
Submitted 18 July, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
The Predictive Individual Effect for Survival Data
Authors:
Beat Neuenschwander,
Satrajit Roychoudhury,
Simon Wandel,
Kannan Natarajan,
Emmanuel Zuber
Abstract:
The call for patient-focused drug development is loud and clear, as expressed in the 21st Century Cures Act and in recent guidelines and initiatives of regulatory agencies. Among the factors contributing to modernized drug development and improved health-care activities are easily interpretable measures of clinical benefit. In addition, special care is needed for cancer trials with time-to-event e…
▽ More
The call for patient-focused drug development is loud and clear, as expressed in the 21st Century Cures Act and in recent guidelines and initiatives of regulatory agencies. Among the factors contributing to modernized drug development and improved health-care activities are easily interpretable measures of clinical benefit. In addition, special care is needed for cancer trials with time-to-event endpoints if the treatment effect is not constant over time. We propose the predictive individual effect which is a patient-centric and tangible measure of clinical benefit under a wide variety of scenarios. It can be obtained by standard predictive calculations under a rank preservation assumption that has been used previously in trials with treatment switching. We discuss four recent Oncology trials that cover situations with proportional as well as non-proportional hazards (delayed treatment effect or crossing of survival curves). It is shown that the predictive individual effect offers valuable insights beyond p-values, estimates of hazard ratios or differences in median survival. Compared to standard statistical measures, the predictive individual effect is a direct, easily interpretable measure of clinical benefit. It facilitates communication among clinicians, patients, and other parties and should therefore be considered in addition to standard statistical results.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
Changes in polarization dictate necessary approximations for modeling electronic de-excitation intensity: an application to X-ray emission
Authors:
Subhayan Roychoudhury,
Leonardo A. Cunha,
Martin Head-Gordon,
David Prendergast
Abstract:
We systematically investigate the underlying relations among different levels of approximation for simulating electronic de-excitations, with a focus on modeling X-ray emission spectroscopy (XES). Using Fermi's golden rule and explicit modeling of the initial, core-excited state and the final, valence-hole state, we show that XES can be accurately modeled by using orbital optimization for the vari…
▽ More
We systematically investigate the underlying relations among different levels of approximation for simulating electronic de-excitations, with a focus on modeling X-ray emission spectroscopy (XES). Using Fermi's golden rule and explicit modeling of the initial, core-excited state and the final, valence-hole state, we show that XES can be accurately modeled by using orbital optimization for the various final states within a Slater-determinant framework. However, in this paper, we introduce a much cheaper approach reliant only on a single self-consistent field for all the final states, and show that it is typically sufficient. Further approximations reveal that these fundamentally many-body transitions can be reasonably approximated by projections of ground state orbitals, but that the ground state alone is insufficient. Furthermore, except in cases where the core-ionization induces negligible changes in polarization, linear-response approaches within the adiabatic approximation will have difficulty in accurately modeling de-excitation to the core level. Therefore, change in the net dipole moment of the valence electrons can serve as a metric for the validity of the linear-response approximation.
△ Less
Submitted 18 December, 2021;
originally announced December 2021.
-
Estimands in Hematologic Oncology Trials
Authors:
Steven Sun,
Hans-Jochen Weber,
Emily Butler,
Kaspar Rufibach,
Satrajit Roychoudhury
Abstract:
The estimand framework included in the addendum to the ICH E9 guideline facilitates discussions to ensure alignment between the key question of interest, the analysis, and interpretation. Therapeutic knowledge and drug mechanism play a crucial role in determining the strategy and defining the estimand for clinical trial designs. Clinical trials in patients with hematological malignancies often pre…
▽ More
The estimand framework included in the addendum to the ICH E9 guideline facilitates discussions to ensure alignment between the key question of interest, the analysis, and interpretation. Therapeutic knowledge and drug mechanism play a crucial role in determining the strategy and defining the estimand for clinical trial designs. Clinical trials in patients with hematological malignancies often present unique challenges for trial design due to complexity of treatment options and existence of potential curative but highly risky procedures, e.g. stem cell transplant or treatment sequence across different phases (induction, consolidation, maintenance). Here, we illustrate how to apply the estimand framework in hematological clinical trials and how the estimand framework can address potential difficulties in trial result interpretation.
This paper is a result of a cross-industry collaboration to connect the International Conference on Harmonisation (ICH) E9 addendum concepts to applications. Three randomized phase 3 trials will be used to consider common challenges including intercurrent events in hematologic oncology trials to illustrate different scientific questions and the consequences of the estimand choice for trial design, data collection, analysis, and interpretation. Template language for describing estimand in both study protocols and statistical analysis plans is suggested for statisticians' reference.
△ Less
Submitted 1 October, 2020;
originally announced October 2020.
-
Principal Stratum Strategy: Potential Role in Drug Development
Authors:
Björn Bornkamp,
Kaspar Rufibach,
Jianchang Lin,
Yi Liu,
Devan V. Mehrotra,
Satrajit Roychoudhury,
Heinz Schmidli,
Yue Shentu,
Marcel Wolbers
Abstract:
A randomized trial allows estimation of the causal effect of an intervention compared to a control in the overall population and in subpopulations defined by baseline characteristics. Often, however, clinical questions also arise regarding the treatment effect in subpopulations of patients, which would experience clinical or disease related events post-randomization. Events that occur after treatm…
▽ More
A randomized trial allows estimation of the causal effect of an intervention compared to a control in the overall population and in subpopulations defined by baseline characteristics. Often, however, clinical questions also arise regarding the treatment effect in subpopulations of patients, which would experience clinical or disease related events post-randomization. Events that occur after treatment initiation and potentially affect the interpretation or the existence of the measurements are called {\it intercurrent events} in the ICH E9(R1) guideline. If the intercurrent event is a consequence of treatment, randomization alone is no longer sufficient to meaningfully estimate the treatment effect. Analyses comparing the subgroups of patients without the intercurrent events for intervention and control will not estimate a causal effect. This is well known, but post-hoc analyses of this kind are commonly performed in drug development. An alternative approach is the principal stratum strategy, which classifies subjects according to their potential occurrence of an intercurrent event on both study arms. We illustrate with examples that questions formulated through principal strata occur naturally in drug development and argue that approaching these questions with the ICH E9(R1) estimand framework has the potential to lead to more transparent assumptions as well as more adequate analyses and conclusions. In addition, we provide an overview of assumptions required for estimation of effects in principal strata. Most of these assumptions are unverifiable and should hence be based on solid scientific understanding. Sensitivity analyses are needed to assess robustness of conclusions.
△ Less
Submitted 8 February, 2021; v1 submitted 12 August, 2020;
originally announced August 2020.
-
Alternative Analysis Methods for Time to Event Endpoints under Non-proportional Hazards: A Comparative Analysis
Authors:
Ray S. Lin,
Ji Lin,
Satrajit Roychoudhury,
Keaven M. Anderson,
Tianle Hu,
Bo Huang,
Larry F Leon,
Jason JZ Liao,
Rong Liu,
Xiaodong Luo,
Pralay Mukhopadhyay,
Rui Qin,
Kay Tatsuoka,
Xuejing Wang,
Yang Wang,
Jian Zhu,
Tai-Tsang Chen,
Renee Iacona,
Cross-Pharma Non-proportional Hazards Working Group
Abstract:
The log-rank test is most powerful under proportional hazards (PH). In practice, non-PH patterns are often observed in clinical trials, such as in immuno-oncology; therefore, alternative methods are needed to restore the efficiency of statistical testing. Three categories of testing methods were evaluated, including weighted log-rank tests, Kaplan-Meier curve-based tests (including weighted Kaplan…
▽ More
The log-rank test is most powerful under proportional hazards (PH). In practice, non-PH patterns are often observed in clinical trials, such as in immuno-oncology; therefore, alternative methods are needed to restore the efficiency of statistical testing. Three categories of testing methods were evaluated, including weighted log-rank tests, Kaplan-Meier curve-based tests (including weighted Kaplan-Meier and Restricted Mean Survival Time, RMST), and combination tests (including Breslow test, Lee's combo test, and MaxCombo test). Nine scenarios representing the PH and various non-PH patterns were simulated. The power, type I error, and effect estimates of each method were compared. In general, all tests control type I error well. There is not a single most powerful test across all scenarios. In the absence of prior knowledge regarding the PH or non-PH patterns, the MaxCombo test is relatively robust across patterns. Since the treatment effect changes overtime under non-PH, the overall profile of the treatment effect may not be represented comprehensively based on a single measure. Thus, multiple measures of the treatment effect should be pre-specified as sensitivity analyses to evaluate the totality of the data.
△ Less
Submitted 20 September, 2019;
originally announced September 2019.
-
Beyond p-values: a phase II dual-criterion design with statistical significance and clinical relevance
Authors:
Satrajit Roychoudhury,
Nicolas Scheuer,
Beat Neuenschwander
Abstract:
Background: Well-designed phase II trials must have acceptable error rates relative to a pre-specified success criterion, usually a statistically significant p-value. Such standard designs may not always suffice from a clinical perspective because clinical relevance may call for more. For example, proof-of-concept in phase II often requires not only statistical significance but also a sufficiently…
▽ More
Background: Well-designed phase II trials must have acceptable error rates relative to a pre-specified success criterion, usually a statistically significant p-value. Such standard designs may not always suffice from a clinical perspective because clinical relevance may call for more. For example, proof-of-concept in phase II often requires not only statistical significance but also a sufficiently large effect estimate.
Purpose: We propose dual-criterion designs to complement statistical significance with clinical relevance, discuss their methodology, and illustrate their implementation in phase II.
Methods: Clinical relevance requires the effect estimate to pass a clinically motivated threshold (the decision value). In contrast to standard designs, the required effect estimate is an explicit design input whereas study power is implicit. The sample size for a dual-criterion design needs careful considerations of the study's operating characteristics (type-I error, power).
Results: Dual-criterion designs are discussed for a randomized controlled and a single-arm phase II trial, including decision criteria, sample size calculations, decisions under various data scenarios, and operating characteristics. The designs facilitate GO/NO-GO decisions due to their complementary statistical-clinical criterion.
Conclusion: To improve evidence-based decision-making, a formal yet transparent quantitative framework is important. Dual-criterion designs offer an appealing statistical-clinical compromise, which may be preferable to standard designs if evidence against the null hypothesis alone does not suffice for an efficacy claim.
△ Less
Submitted 21 August, 2019;
originally announced August 2019.
-
Bayesian leveraging of historical control data for a clinical trial with time-to-event endpoint
Authors:
Satrajit Roychoudhury,
Beat Neuenschwander
Abstract:
The recent 21st Century Cures Act propagates innovations to accelerate the discovery, development, and delivery of 21st century cures. It includes the broader application of Bayesian statistics and the use of evidence from clinical expertise. An example of the latter is the use of trial-external (or historical) data, which promises more efficient or ethical trial designs. We propose a Bayesian met…
▽ More
The recent 21st Century Cures Act propagates innovations to accelerate the discovery, development, and delivery of 21st century cures. It includes the broader application of Bayesian statistics and the use of evidence from clinical expertise. An example of the latter is the use of trial-external (or historical) data, which promises more efficient or ethical trial designs. We propose a Bayesian meta-analytic approach to leveraging historical data for time-to-event endpoints, which are common in oncology and cardiovascular diseases. The approach is based on a robust hierarchical model for piecewise exponential data. It allows for various degrees of between trial-heterogeneity and for leveraging individual as well as aggregate data. An ovarian carcinoma trial and a non-small-cell cancer trial illustrate methodological and practical aspects of leveraging historical data for the analysis and design of time-to-event trials.
△ Less
Submitted 6 February, 2020; v1 submitted 20 August, 2019;
originally announced August 2019.
-
Robust Design and Analysis of Clinical Trials With Non-proportional Hazards: A Straw Man Guidance from a Cross-pharma Working Group
Authors:
Satrajit Roychoudhury,
Keaven M Anderson,
Jiabu Ye,
Pralay Mukhopadhyay
Abstract:
Loss of power and clear description of treatment differences are key issues in designing and analyzing a clinical trial where non-proportional hazard is a possibility. A log-rank test may be very inefficient and interpretation of the hazard ratio estimated using Cox regression is potentially problematic. In this case, the current ICH E9 (R1) addendum would suggest designing a trial with a clinical…
▽ More
Loss of power and clear description of treatment differences are key issues in designing and analyzing a clinical trial where non-proportional hazard is a possibility. A log-rank test may be very inefficient and interpretation of the hazard ratio estimated using Cox regression is potentially problematic. In this case, the current ICH E9 (R1) addendum would suggest designing a trial with a clinically relevant estimand, e.g., expected life gain. This approach considers appropriate analysis methods for supporting the chosen estimand. However, such an approach is case specific and may suffer lack of power for important choices of the underlying alternate hypothesis distribution. On the other hand, there may be a desire to have robust power under different deviations from proportional hazards. Also, we would contend that no single number adequately describes treatment effect under non-proportional hazards scenarios. The cross-pharma working group has proposed a combination test to provide robust power under a variety of alternative hypotheses. These can be specified for primary analysis at the design stage and methods appropriately accounting for combination test correlations are efficient for a variety of scenarios. We have provided design and analysis considerations based on a combination test under different non-proportional hazard types and present a straw man proposal for practitioners. The proposals are illustrated with real life example and simulation.
△ Less
Submitted 12 January, 2021; v1 submitted 19 August, 2019;
originally announced August 2019.
-
Spin-Phonon coupling parameters from maximally localized Wannier functions and first principles electronic structure: the case of durene single crystal
Authors:
Subhayan Roychoudhury,
Stefano Sanvito
Abstract:
Spin-orbit interaction is an important vehicle for spin relaxation. At finite temperature lattice vibrations modulate the spin-orbit interaction and thus generate a mechanism for spin-phonon coupling, which needs to be incorporated in any quantitative analysis of spin transport. Starting from a density functional theory \textit{ab initio} electronic structure, we calculate spin-phonon matrix eleme…
▽ More
Spin-orbit interaction is an important vehicle for spin relaxation. At finite temperature lattice vibrations modulate the spin-orbit interaction and thus generate a mechanism for spin-phonon coupling, which needs to be incorporated in any quantitative analysis of spin transport. Starting from a density functional theory \textit{ab initio} electronic structure, we calculate spin-phonon matrix elements over the basis of maximally localized Wannier functions. Such coupling terms form an effective Hamiltonian to be used to extract thermodynamic quantities, within a multiscale approach particularly suitable for organic crystals. The symmetry of the various matrix elements are analyzed by using the $Γ$-point phonon modes of a one-dimensional chain of Pb atoms. Then the method is employed to extract the spin-phonon coupling of solid durene, a high-mobility crystal organic semiconducting. Owing to the small masses of carbon and hydrogen spin-orbit is weak in durene and so is the spin-phonon coupling. Most importantly we demonstrate that the largest contribution to the spin-phonon interaction originates from Holstein-like phonons, namely from internal molecular vibrations.
△ Less
Submitted 12 July, 2018;
originally announced July 2018.
-
XDFT: an efficient first-principles method for neutral excitations in molecules
Authors:
Subhayan Roychoudhury,
Stefano Sanvito,
David D. O'Regan
Abstract:
State-of-the-art methods for calculating neutral excitation energies are typically demanding and limited to single electron-hole pairs and their composite plasmons. Here we introduce excitonic density-functional theory (XDFT) a computationally light, generally applicable, first-principles technique for calculating neutral excitations based on generalized constrained DFT. In order to simulate an M-…
▽ More
State-of-the-art methods for calculating neutral excitation energies are typically demanding and limited to single electron-hole pairs and their composite plasmons. Here we introduce excitonic density-functional theory (XDFT) a computationally light, generally applicable, first-principles technique for calculating neutral excitations based on generalized constrained DFT. In order to simulate an M-particle excited state of an N-electron system, XDFT automatically optimizes a constraining potential to confine N-M electrons within the ground-state Kohn-Sham valence subspace. We demonstrate the efficacy of XDFT by calculating the lowest single-particle singlet and triplet excitation energies of the well-known Thiel molecular test set, with results which are in excellent agreement with time-dependent DFT. Furthermore, going beyond the capability of adiabatic time-dependent DFT, we show that XDFT can successfully capture double excitations. Overall our method makes optical gaps, excition bindings and oscillator strengths readily accessible at a computational cost comparable to that of standard DFT. As such, XDFT appears as an ideal candidate to work within high-throughput discovery frameworks and within linear-scaling methods for large systems.
△ Less
Submitted 4 March, 2018;
originally announced March 2018.
-
Wannier-function-based constrained DFT with nonorthogonality-correcting Pulay forces in application to the reorganization effects in graphene-adsorbed pentacene
Authors:
Subhayan Roychoudhury,
David D. O'Regan,
Stefano Sanvito
Abstract:
Pulay terms arise in the Hellman-Feynman forces in electronic structure calculations when one employs a basis set made of localized orbitals that move with their host atoms. If the total energy of the system depends on a subspace population defined in terms of the localized orbitals across multiple atoms, then unconventional Pulay terms will emerge due to the variation of the orbital nonorthogonal…
▽ More
Pulay terms arise in the Hellman-Feynman forces in electronic structure calculations when one employs a basis set made of localized orbitals that move with their host atoms. If the total energy of the system depends on a subspace population defined in terms of the localized orbitals across multiple atoms, then unconventional Pulay terms will emerge due to the variation of the orbital nonorthogonality with ionic translation. Here, we derive the required exact expressions for such terms, which cannot be eliminated by orbital orthonormalization. We have implemented these corrected ionic forces within the linear-scaling density functional theory (DFT) package ONETEP, and have used constrained DFT to calculate the reorganization energy of a pentacene molecule adsorbed on a graphene flake. The calculations are performed by including ensemble DFT, corrections for periodic boundary conditions, and empirical Van der Waals interactions. For this system we find that tensorially invariant population analysis yields an adsorbate subspace population that is very close to integer-valued when based upon nonorthogonal Wannier functions, and also but less precisely when using pseudoatomic functions. Thus, orbitals can provide a very effective population analysis for constrained DFT. Our calculations show that the reorganization energy of the adsorbed pentacene is typically lower than that of pentacene in the gas phase. We attribute this effect to steric hindrance.
△ Less
Submitted 8 May, 2018; v1 submitted 5 February, 2018;
originally announced February 2018.
-
Spin-orbit Hamiltonian for organic crystals from first principles electronic structure and Wannier functions
Authors:
Subhayan Roychoudhury,
Stefano Sanvito
Abstract:
Spin-orbit coupling in organic crystals is responsible for many spin-relaxation phenomena, going from spin diffusion to intersystem crossing. With the goal of constructing effective spin-orbit Hamiltonians to be used in multiscale approaches to the thermodynamical properties of organic crystals, we present a method that combines density functional theory with the construction of Wannier functions.…
▽ More
Spin-orbit coupling in organic crystals is responsible for many spin-relaxation phenomena, going from spin diffusion to intersystem crossing. With the goal of constructing effective spin-orbit Hamiltonians to be used in multiscale approaches to the thermodynamical properties of organic crystals, we present a method that combines density functional theory with the construction of Wannier functions. In particular we show that the spin-orbit Hamiltonian constructed over maximally localised Wannier functions can be computed by direct evaluation of the spin-orbit matrix elements over the Wannier functions constructed in absence of spin-orbit interaction. This eliminates the prob- lem of computing the Wannier functions for almost degenerate bands, a problem always present with the spin-orbit-split bands of organic crystals. Examples of the method are presented for isolated molecules, for mono-dimensional chains of Pb and C atoms and for triarylamine-based one-dimansional single crystals.
△ Less
Submitted 9 February, 2017; v1 submitted 25 November, 2016;
originally announced November 2016.
-
Charge transfer energies of benzene physisorbed on a graphene sheet from constrained density functional theory
Authors:
Subhayan Roychoudhury,
Carlo Motta,
Stefano Sanvito
Abstract:
Constrained density functional theory (CDFT) is used to evaluate the energy level alignment of a benzene molecule as it approaches a graphene sheet. Within CDFT the problem is conveniently mapped onto evaluating total energy differences between different charge-separated states, and it does not consist in determining a quasi-particle spectrum. We demonstrate that the simple local density approxima…
▽ More
Constrained density functional theory (CDFT) is used to evaluate the energy level alignment of a benzene molecule as it approaches a graphene sheet. Within CDFT the problem is conveniently mapped onto evaluating total energy differences between different charge-separated states, and it does not consist in determining a quasi-particle spectrum. We demonstrate that the simple local density approximation provides a good description of the level aligmnent along the entire binding curve, with excellent agreement to experiments at an infinite separation and to GW calculations close to the bonding distance. The method also allows us to explore the effects due to the presence of graphene structural defects and of multiple molecules. In general all our results can be reproduced by a classical image charge model taking into account the finite dielectric constant of graphene.
△ Less
Submitted 27 October, 2015;
originally announced October 2015.
-
On the structure of a family of probability generating functions induced by shock models
Authors:
Satrajit Roychoudhury,
Manish C. Bhattacharjee
Abstract:
We explore conditions for a class of functions defined via an integral representation to be a probability generating function of some positive integer valued random variable. Interest in and research on this question is motivated by an apparently surprising connection between a family of classic shock models due to Esary et. al. (1973) and the negatively aging nonparametric notion of ``strongly…
▽ More
We explore conditions for a class of functions defined via an integral representation to be a probability generating function of some positive integer valued random variable. Interest in and research on this question is motivated by an apparently surprising connection between a family of classic shock models due to Esary et. al. (1973) and the negatively aging nonparametric notion of ``strongly decreasing failure rate'' (SDFR) introduced by Bhattacharjee (2005). A counterexample shows that there exist probability generating functions with our integral representation which are not discrete SDFR, but when used as shock resistance probabilities can give rise to a SDFR survival distribution in continuous time.
△ Less
Submitted 15 May, 2008;
originally announced May 2008.
