Goodness-of-fit criteria for survival data

The definition of an appropriate measure for goodness-of-fit in case of survival data comparable to R^2 in linear regression is difficult due to censored observations. In this paper, a variety of answers based on different residuals and variance of survival curves are presented together with a newly introduced criterion. In univariate simulation studies, the presented criteria are examined with respect to their dependence on the value of the coefficient associated with the covariate; underlying covariate distribution and censoring percentage in the data. Investigation of the relations between the values of the different criteria indicates strong dependencies, although the absolute values show high discrepancies and the criteria building processes differ substantially

Symmetric Diblock Copolymers in Thin Films (I): Phase stability in Self-Consistent Field Calculations and Monte Carlo Simulations

We investigate the phase behavior of symmetric AB diblock copolymers confined into a thin film. The film boundaries are parallel, impenetrable and attract the A component of the diblock copolymer. Using a self-consistent field technique [M.W. Matsen, J.Chem.Phys. {\bf 106}, 7781 (1997)], we study the ordered phases as a function of incompatibility $\chi$ and film thickness in the framework of the Gaussian chain model. For large film thickness and small incompatibility, we find first order transitions between phases with different number of lamellae which are parallel oriented to the film boundaries. At high incompatibility or small film thickness, transitions between parallel oriented and perpendicular oriented lamellae occur. We compare the self-consistent field calculations to Monte Carlo simulations of the bond fluctuation model for chain length N=32. In the simulations we quench several systems from $\chi N=0$ to $\chi N=30$ and monitor the morphology into which the diblock copolymers assemble. Three film thicknesses are investigated, corresponding to parallel oriented lamellae with 2 and 4 interfaces and a perpendicular oriented morphology. Good agreement between self-consistent field calculations and Monte Carlo simulations is found.Comment: to appear in J.Chem.Phy

Monte Carlo simulations of copolymers at homopolymer interfaces: Interfacial structure as a function of the copolymer density

By means of extensive Monte Carlo simulations of the bond fluctuation model, we study the effect of adding AB diblock copolymers on the properties of an interface between demixed homopolymer phases. The parameters are chosen such that the homopolymers are strongly segregated, and the whole range of copolymer concentrations in the two phase coexistence region is scanned. We compare the ``mushroom'' regime, in which copolymers are diluted and do not interact with each other, with the ``wet brush'' regime, where copolymers overlap and stretch, but are still swollen by the homopolymers. A ``dry brush'' regime is never entered for our choice of chain lengths. ``Intrinsic'' profiles are calculated using a block analysis method introduced by us in earlier work. We discuss density profiles, orientational profiles and contact number profiles. In general, the features of the profiles are similar at all copolymer concentrations, however, the profiles in the concentrated regime are much broader than in the dilute regime. The results compare well with self-consistent field calculations.Comment: to appear in J. Chem. Phy

Implementation of complex interactions in a Cox regression framework

The standard Cox proportional hazards model has been extended by functionally describable interaction terms. The first of which are related to neural networks by adopting the idea of transforming sums of weighted covariables by means of a logistic function. A class of reasonable weight combinations within the logistic transformation is described. Apart from the standard covariable product interaction, a product of logistically transformed covariables has also been included in the analysis of performance of the new terms. An algorithm combining likelihood ratio tests and AIC criterion has been defined for model choice. The critical values of the likelihood ratio test statistics had to be corrected in order to guarantee a maximum type I error of 5% for each interaction term. The new class of interaction terms allows interpretation of functional relationships between covariables with more flexibility and can easily be implemented in standard software packages

Dephasing-assisted Gain and Loss in Mesoscopic Quantum Systems

Motivated by recent experiments, we analyse the phonon-assisted steady-state gain of a microwave field driving a double quantum-dot in a resonator. We apply the results of our companion paper, which derives the complete set of fourth-order Lindblad dissipators using Keldysh methods, to show that resonator gain and loss are substantially affected by dephasing-assisted dissipative processes in the quantum-dot system. These additional processes, which go beyond recently proposed polaronic theories, are in good quantitative agreement with experimental observationsComment: 5 pages, 3 Figures, published together with arXiv:1608.0416

Equation of state and critical behavior of polymer models: A quantitative comparison between Wertheim's thermodynamic perturbation theory and computer simulations

We present an application of Wertheim's Thermodynamic Perturbation Theory (TPT1) to a simple coarse grained model made of flexibly bonded Lennard-Jones monomers. We use both the Reference Hyper-Netted-Chain (RHNC) and Mean Spherical approximation (MSA) integral equation theories to describe the properties of the reference fluid. The equation of state, the density dependence of the excess chemical potential, and the critical points of the liquid--vapor transition are compared with simulation results and good agreement is found. The RHNC version is somewhat more accurate, while the MSA version has the advantage of being almost analytic. We analyze the scaling behavior of the critical point of chain fluids according to TPT1 and find it to reproduce the mean field exponents: The critical monomer density is predicted to vanish as $n^{-1/2}$ upon increasing the chain length $n$ while the critical temperature is predicted to reach an asymptotic finite temperature that is attained as $n^{-1/2}$. The predicted asymptotic finite critical temperature obtained from the RHNC and MSA versions of TPT1 is found to be in good agreement with the $\Theta$ point of our polymer model as obtained from the temperature dependence of the single chain conformations.Comment: to appear in J.Chem.Phy

ReDecay: A novel approach to speed up the simulation at LHCb

With the steady increase in the precision of flavour physics measurements collected during LHC Run 2, the LHCb experiment requires simulated data samples of larger and larger sizes to study the detector response in detail. The simulation of the detector response is the main contribution to the time needed to simulate full events. This time scales linearly with the particle multiplicity. Of the dozens of particles present in the simulation only the few participating in the signal decay under study are of interest, while all remaining particles mainly affect the resolutions and efficiencies of the detector. This paper presents a novel development for the LHCb simulation software which re-uses the rest of the event from previously simulated events. This approach achieves an order of magnitude increase in speed and the same quality compared to the nominal simulation