Bayesian Analysis of Survival Models with Bathtub Hazard Rates Purushottam W. Laud, Ph.D. Paul Damien, Ph.D. Stephen G. Walker, Ph.D. Technical Report
|
|
- Warren Holmes
- 5 years ago
- Views:
Transcription
1 Bayesian Analysis of Survival Models with Bathtub Hazard Rates Purushottam W. Laud, Ph.D. Paul Damien, Ph.D. Stephen G. Walker, Ph.D. Technical Report 31 June 1998 Division of Biostatistics Medical College of Wisconsin 871 Watertown Plank Road Milwaukee WI Phone: (414)
2 Bayesian Analysis of Survival Models with Bathtub Hazard Rates Purushottam W. Laud 1 Paul Damien 2 Stephen G. Walker 3 1 Division of Biostatistics, Medical College of Wisconsin, Milwaukee, USA. 2 University of Michigan Business School, Ann Arbor, USA. 3 Department of Mathematics, Imperial College, London, UK. SUMMARY Nonparametric Bayesian inference using bathtub hazard rates in survival analysis is developed. A new simulation method, Auxiliary Random Functions, is introduced. When used within a Gibbs sampler, this method enables a unied treatment of exact, right-censored, left-censored, left-trucated and interval censored data, with and without covariates. The models and methods are exemplied via illustrative analyses. Key Words: Extended gamma process, Censored data, Gibbs sampling, Auxiliary functions, Covariates, Bioassay. 1 Introduction In survival analysis, often times, hazard specic prior information and/or hazard specic inference may be available and/or desirable. Modelling increasing or decreasing hazard rates is relatively straightforward; see, for example, Dykstra & Laud (1981). In reality, however, hazard rates are typically both decreasing and increasing, possibly with substantial uncertainty regarding the point where the monotonicity changes. Amman (1984) suggests a model for such hazard rates, appropriately termed bathtubs, using a combined extended gamma process. In this context, a unied nonparametric Bayesian inference solution for bathtub hazard rates given dierent types of censored and/or truncated data with or without covariates has not been available. This is partly due to the complexity of the modelling structure and partly due to the one-o, often dicult, computational approaches needed to sample from the appropriate posterior distributions. To the best of our knowledge, there are only three approaches available 1
3 for hazard rates. The other two, which do not model the shape, are : (i) the conditional gamma martingale of Arjas and Gasbarra (1994) and (ii) the discretized correlated constant hazards of Leonard (1978) and Gamerman (1991). For an excellent review of the general literature on Bayesian survival analysis, see Sinha and Dey (1997). This paper develops a non/semi-parametric Bayesian analysis of bathtub shaped hazard rates given dierent types of censored data. Computations are rendered straightforward based on a new stochastic substitution method. This method, which we term Auxiliary Random Functions (ARF), is particularly suited to simulate from the posterior conditional distributions in a Gibbs sampler (Smith and Roberts, 1993) in the present context. 2 Bathub Processes Here we collect, briey, the key properties of the extended gamma process that will be of use in the rest of the discussion; see, Dykstra & Laud (1981) and Amman (1984) for details. To begin, we consider increasing hazard rates. Let G(; ) denote the gamma distribution with shape parameter and scale parameter 1= >. For =, we dene this distribution to be degenerate at zero. For >, its density is ( x?1 exp(?x) g(xj; ) = ; x > ;?() ; x. Let (t), t, be a non-decreasing left continuous real valued function such that () = and (t); t, be a positive right continuous real valued function with left hand limits existing and bounded away from and 1. Take Z(t), t, to be a gamma process with parameter function (). That is, Z() =, Z(t) has independent increments and for t > s, Z(t)? Z(s) is G((t)? (s); 1). Considering a version of this process such that all its sample paths are non-decreasing and left continuous, let Z r(t) = [(s)]?1 dz(s): [;t) This process fr(t); t g is called the Extended Gamma (EG) (or the Increasing Extended, IEG,)process and we denote it by r(t) I?((); ()): In a manner similar to the description above one can dene an EG process for decreasing hazard rates as in Laud (1977). Take (), () and Z() as above. In addition, let Z(1) G((1); 1) be such that Z(1)?lim t!1 Z(t) is independent of the rest of the process, (1) > lim t!1 (t) and (1) >. Now dene Z Z r(t) = [(s)]?1 dz(s) = [t;1] [t;1) [(s)]?1 dz(s) + [(1)]?1 [Z(1)? lim 2 t!1 Z(t)] :
4 The main purpose of the values at innity is to preclude lim t!1 r(t) = so that R [;1) r(t)dt is innite and the survival time described by r() does not have positive probability at innity. It is clear that such a process generates nonincreasing hazard rates. We call it the Decreasing Extended Gamma (DEG) and denote it by r(t) D?((); ()): Amman (1984) uses the above development and denes a combined process by r(t) = r I (t) + r D (t); where the superscript, inherited also by the parameter functions, indicates the direction of monotonicity, and the two component processes are independent. Amman termed this a U- shaped process whereas elsewhere in the reliability and survival analysis literature one also encounters the usage \bathtub-shaped hazard rates". It should be noted that such a combined process does not necessarily generate hazard rates with this shape. It does, however, add a great deal of exibility in modelling the shape. For example, the expected hazard rate can be arranged to decrease initially and to increase beyond a certain time point by a suitable choice of the parameter functions. A parametric choice for the expectation can be a Weibull distribution with shape parameter less than one for the decreasing component and shape parameter exceeding two for the increasing component. In general, one can require with Since we get the conditions and d 2 d2 E(r(t)) = dt2 dt 2 E(rI (t)) + d2 dt 2 E(rD (t)) Z E(r(t)) = d E(r(t)) = for some t > : dt [;t) Z [ I (s)]?1 d I (s) + [t;1] [ D (s)]?1 d D (s); I (t) I (2) (t)? I (1) (t)i (1) (t) [ I (t)] 2? D (t) D (2) (t)? D (1) (t)d (1) (t) [ D (t)] 2 I (t) (1) I (t)? D (t) (1) = for some t > D (t) where the parenthesized subscript i on a function denotes its i th derivative. Such a prior, while maintaining a bathtub expectation, allows the hazard rates to have more general shapes. The data can then either reinforce the bathtub shape or indicate otherwise via the posterior expectation. 3
5 Several aspects of the posterior distribution are given by Dykstra and Laud (1981), Laud (1977) and Amman (1984). These can be summarized as follows : Fact 2.1 : The EG (or DEG) process is conjugate with respect to right censored data. In particular, given an observation X > x, writing R x r(t)dt as R 1 (x? R t)+ dr(t) (or 1 min(x; t)dr(t)) it follows that the posterior process is again EG (or DEG) with I ()(or D ()) unchanged and the parameter modied as ~ I (t) = I (t) + (x? t) + ~ D (t) = D (t) + min(x; t) : Here the superscript + on a quantity denotes its positive part. For a combined process, the posterior is simply the combined process made with these functions. Fact 2.2 : The posterior with respect to exact data is a mixture of EG(DEG) processes. The mixing measure is n-dimensional, where n is the number of exact observations, and has both continuous and jump components. For a combined process, the posterior is again a mixture of the corresponding combined processes. Fact 2.3 : Although analytic expressions for the posterior mean and variance functions are given by the above mentioned authors, their numerical evaluation is virtually impossible, especially as the sample size increases to more than around 2. They do not even consider the possibility of including covariates as a consequence. Our goal in the rest of this paper is to provide a full Bayesian analysis of the combined hazard rate model for various types of censored data as well as exact observations, both with and without covariates. We note here that the rest of the discussion applies also to data arising in a bioassay because a bioassay can be treated as a censored data problem. In addition, these methods also work for estimating the intensity rate function of a nonhomogeneous Poisson process because the form of the likelihood there is the same as that for survival data (see, for example, Lawless (1982)). Notation: In the following, [Aj:] and [A] denote the conditional and marginal distributions of A, respectively; note that A could be a process as well. Also, where the context is clear, [A] will imply the prior distribution of A. I() refers to the indicator function. 4
6 3 Posterior Computations Via ARFs Damien et al. (1998) develop a general method to sample from posterior distributions using auxiliary variables as follows. Suppose the required conditional distribution for a random variable X is denoted f. The basic idea is to introduce an auxiliary variable U, construct the joint density of U and X, with marginal density for X given by f, and then extend the Gibbs sampler to include the extra full conditional for U. In a Bayesian context, consider the posterior density given by f(x) / l(x)(x) and suppose it is not possible to sample directly from f. The general idea proposed by Damien et al. is to introduce an auxiliary variable Y, dened on the interval (; 1) or more strictly the interval (; l(^x)), where ^x maximises l(:), and dene the joint density with X by f(x; y) / I(y < l(x))(x): The full conditional for Y is U(; l(x)), a uniform variable on the interval (; l(x)), and the full conditional for X is, restricted to the set A y = fx : l(x) > yg. In this paper, we develop a simulation method that is particularly suited to the context of sampling monotone hazard rate processes in a manner similar in spirit to that of Damien et al. (1998). Whereas Damien et al. use latent variables to put a function in place which depends on a nite number of parameters, the latent variables used in this paper put an entire random function into place. Hence we refer to it as the auxiliary random function (ARF) method. To introduce the method, we rst consider only the increasing hazard rate case without covariates, dropping the superscript I. We then implement it for combined processes and various survival data scenarios. Exact observations The contribution to the likelihood resulting from one exact observation at x is given by: Z x f(x) = r(x) F (x) = r(x) exp? r(s)ds ; where F is the survival function. Based on Fact 2.1, it is clear that F (x), being the contribution of a right-censored observation, updates the prior parameters of the EG process. In the above equation, from an ARF perspective, the following question arises: how do we substitute the random function r() so that sampling from the posterior is accomplished readily? Take the cumulative distribution of a random variable Y given the hazard rate process r() as: P (Y yjr()) = r(y) I( y x) + I(y x): r(x) The distribution of Y is dened in terms of a random function, hence the phrase \auxiliary random function". It is in this sense that the ARF algorithm is dierent from Damien et al.'s 5
7 algorithm. To sample Y jr(), let Y = r?1 (Ur(x)) where U U(; 1) and r?1 (w) = infftjr(t) wg. To address sampling r()jy, write the density of Y jr() as [yjr()] = dr(y) I( y x): r(x) Combining this with the prior and the likelihood we arrive at Z x [r()jy] / [r()]dr(y) exp? r(s)ds ; where dr(y) denotes the increment of r() at y. It is clear that this term combines naturally with the EG process, adding a unit jump at y to the shape parameter (). Thus, as a result of the ARF substitution and Fact 2.1, the posterior full conditional for r() is once again an EG process with the parameters updated by: ~(s) = (s) + I(s y) and ~ (s) = (s) + (x? s) + : The case of the decreasing hazard rates is similar with the following change : the auxiliary variable Y has the survival function P (Y > yjr()) = r(y) I(y > x) r(x) and possibly takes the value \1". This results in a DEG process for the full conditional for r() with parameters ~(s) = (s) + I(s y) and ~ (s) = (s) + min(x; s): Turning to the case of the combined process, we write the posterior distribution as [r D (); r I ()jx = x] / [r D ()][r I ()]fr D (x) + r I (x)g exp(? Now the stochastic substitution takes the form Z x fr D (s) + r I (s)gds): [yjx = x; r D (); r I ()] = dri (y)i( y x) + f?dr D (y)gi(y > x) : r D (x) + r I (x) Straightforward algebra results in where [r D (); r I ()jy; X = x] / EG(~ I ; ~ I )EG(~ D ; ~ D ); ~ D (s) = D (s) + I(y > x)i(s y); ~ I (s) = I (s) + I( y x)i(s y) and the two scale parameter functions as given in Fact
8 In each of the cases considered above, the full conditional of the hazard rate after the ARF substitution reduces to EG, DEG or a combination. The increments of the EG(DEG) process can be simulated by using the Bondesson (1981) method as described in Laud, Smith and Damien (1996) or, as is often adequate, simply approximated as gamma variates. Left censored data The contribution to the likelihood from a left-censored observaton, X x, is Z x F (x) = 1? exp? r(s)ds : First, using the auxiliary variate W jr() with hazard rate r() restricted to [; x], we have and [wjr(); X x] = r(w) exp(? R w r(s)ds) 1? exp(? R x I( w x); r(s)ds) Z w [r()jw; X x] / [r()]r(w) exp? r(s)ds : Now, in a manner similar to the exact data scenario, we introduce an ARF Y and proceed as described above. Interval censored data The contribution of an observation, a < X b to the likelihood is Z a Z! b F (a)? F (b) = exp? r(s)ds? exp? r(s)ds : Employing W jr() with hazard rate r() restricted to (a; b] yields the same full conditional for r() as in the left censored data scenario above. Left truncated data Here one observes X (or a censored version of it) only if it exceeds an observed left truncation time L. Thus the contribution of X = x; L = l to the likelihood is Z f(x) x? F (l) I(l x) = r(x) exp r(t)dt I(l x): 7 l
9 Once again, the ARF Y as dened above suces to substitute for r(x). It can be shown easily that the eect of the term e?r x r(t)dt l is to modify the scale functions given in Fact 2.1 as ~ I (t) = I (t) + (x? max(t; l)) + and ~ D (t) = D (t) + (min(x; t)? l) + : The unique conjugacy property of the extended gamma process for the latent model is highlighted in the Appendix. 4 Cox Model Computations Via ARFs Consider now the extension of the nonparametric analysis to include covariates. Employing the Cox model, Cox (1972), and initially tackling a single right-censored observation, let and z denote, respectively, the vectors of regression parameters and covariates. The likelihood function for a right censored observation is given by L(r(); ) = P (X > xjr(); ) = exp?e z Z x r(s)ds : Conditioning on and following the development leading to Fact 2.1, we easily arrive at [r()j; X > x] EG (); () ~ where ~() = () + e z (x? s) + : A similar result obtains for decreasing hazards, yielding ~ () = () + e z min(x; s). In the case of exact data, [r()j; X = x] = [r()]e z r(x) exp?e z Z x r(s)ds : Clearly, the problem now yields to the same ARF substitution employed in the case without covariates. In addition, the decreasing hazards follow suit, leading to the same solution for the combined process as in the previous section. To examine the full conditional distribution of, let denote the usual \noncensoring" indicator and write the likelihood as Z x L(r(); ; x; ) / e z exp?e z r(s)ds : Simulating from [j] is easily accomplished using the auxiliary variable method discussed in Damien et al. Omitting straightforward algebraic details, the simulation of j is given by the following full conditional in the Gibbs sampler : [vjall else] / U(; e z j j ); 8
10 [wjall else] / E(1)I(H(x)e z < w); [ j jall else] / [ j ]fi(a < j < b; z j > ) + I(b < j < a; z j < ) + I(z j = )g; where U and E(1) denote the uniform and standard exponential distributions, respectively and a = ln(v) z j, b = ln(w?ln(h(x))?p z k6=j k k z j ; H(x) = R x r(s)ds. 5 Illustrative Analyses We describe a simulated example without covariates and a case study involving covariates. Example 1 We simulated samples of size 5, 2 and 5 from a mixture of two gamma densities with common scale.25, and shape parameters.5 and 5, with mixing weights.7 and.3 respectively. The corresponding bathtub hazard rates are depicted in Figure 1 (see, Glaser(198) for conditions under which mixtures have bathtub hazards). We consider two dierent prior means: Prior1: r(t) = 4exp(?2t + :8t; Prior2:r(t) = exp(?2t) + :4t; these are depicted in Figure 1. Based on the model described in Section 3, the posterior mean and the sample-based two-standard-deviation limits were calculated; these are also shown in Figure 1. The shape of the posterior mean is seen to respond to the prior modeling via the combined process. Example 2 We consider here data for 548 chronic myelogeneous leukemia (CML) patients between the ages of 15 and 55 years who received an HLA-identical sibling bone marrow transplant. These transplants were reported to the International Bone Marrow Transplant Registry, a working group of over 3 transplant centers worldwide that contribute transplant and follow-up data to a statistical center at the Medical College of Wisconsin. These 548 patients were diagnosed between 1983 and 1991 as suering from CML. Approximately half of the observations are right censored corresponding to patients that are alive or, in a few cases, lost to follow-up. More detailed description of the data are available in Zhang and Klein (1998) and the references therein. For the purposes of this illustration, we consider the survival time from the transplant date as possibly related to the following covariates: (i) Time from diagnosis to transplant. (ii) Gender. (iii) Whether or not the CML diagnosis was made after (It is well recognized in the related medical literature that substantial technical improvements became available in 1988.) (iv) Whether or not the patient's spleen size exceeded 1cm. (v) Whether or not the patient's age exceeded 34. Employing the proportional hazards regression model of Cox (1972), it is thought desirable to model the baseline hazard rate as initially high, reecting complications related to the 9
11 transplant, decreasing thereafter and possibly increasing as patients relapse into the disease. Implementing the analysis outlined in the previous sections, the posterior distributions of the regression coecients indicated that only gender and diagnosis date aected survival. The results presented in Figures 2 through 4 were obtained using these two covariates and their interaction. With the intention of leaving the likelihood nearly unaected by the prior, uniform prior distributions for the regression coecients were chosen, each on the interval (?5; 5). In 5 samples from the Gibbs chain after discarding the initial 1, the values for the interaction were all positive, and those for the rst two coecients all negative, providing substantial evidence that all three terms inuence the survival time. It is customary in survival analysis of this type to address one's inferences to the more readily interpretable relative risk (or risk ratio) as opposed to the original regression coecient. For a binary covariate, it is the ratio of the hazards for the two groups and is constant over time in Cox's model. Figure 2 presents posterior densities of four such ratios to illuminate the nature of the interaction. From the two solid line curves one notes that the gender gap in risk prior to 1988 diminished considerably in the later period. Looked at another way, via the dashed curves, risk for males lowered substantially as shown by the curve to the left lying entirely below.85. For females, however, the procedural improvements of 1988 and later do not appear to have resulted in decreased risk. Albeit by a small amount, the posterior probability of an increased risk indeed exceeds.5. Figure 3 compares post-transplant survival in the four groups. At 25 months, there is much posterior uncertainty in dierentiating among three of the four groups; only males with pre diagnoses can be condently declared to have a lower survival probability than the rest. The posterior mean survival probabilities over eight years are presented in part (b). It should be noted that the posterior uncertainties for these curves are readily available via the posterior samples. However, their clear depiction is somewhat challenging. Finally, Figure 4 addresses the baseline hazard rate and survival. It presents the posterior mean curves and two-sigma limits. The prior means are also shown. The prior two-sigma limits are omitted, however, since the proper but vague choice of prior parameter functions rendered these limits o-scale. 6 Discussion In this paper we have provided a complete Bayesian solution to modelling bathtub hazard rates nonparametrically. A new stochastic simulation method was introduced to sample the posterior processes. We also provided a solution for the semi-parametric model with covariates. Natural extensions are to consider time-dependent covariates and to detect when data contradict the proportional hazards assumption. From the perspective of the physicians and scientists wishing 1
12 to use these methods, based on the stored posterior samples, graphical implementations that can respond to queries regarding substantively interesting functionals would be attractive indeed. Acknowledgements This research was funded in part by an NSF grant, an EPSRC Realising Our Potential Award, an EPSRC Fellowship, and a Faculty Research Award by the University of Michigan Business School. We thank the IBMTR for providing the data in Example 2. Appendix We are working with the joint data model with likelihood component Z f(s; tjr())dtds = dr(s) exp (t? u) + dr(u) I( s t)dt; where the function r() is the parameter and s is latent; i.e. unobserved, whereas t is observed. We show that if r() has independent increments then we have a desirable conjugacy property if these increments are gamma variables. This can be seen if we write the likelihood for one observation as Y dr(u) I(u=s) exp??(t? u) + dr(u) : u Therefore, if dr() G((); ()), the updates given in Section 3 are immediate. The expression above demonstrates the unique conjugacy property of the extended gamma process. 11
13 References Amman, L. (1984). Bayesian nonparametric inference for quantal response data. Annals of Statistics, 12, Arjas, E. and Gasbarra, D. (1994). Nonparametric Bayesian inference from right censored survival data using the Gibbs sampler. Statistica Sinica, 4, Bondesson, L. (1982). On simulation from innitely divisible distributions. Adv. App. Prob. 14, Cox, D.R. (1972). Regression models and life tables (with discussion). J. Roy. Statist. Soc. B, 34, Damien, P., Wakeeld, J.C. and Walker, S.G. (1998). Gibbs sampling for Bayesian nonconjugate models using auxiliary variables. To appear in the J. Roy. Statist. Soc. B. Dykstra, R.L. and Laud, P.W. (1981). A Bayesian nonparametric approach to reliability. Annals of Statistics, 9, Gamerman, D. (1991). Dynamic Bayesian models for survival data. Applied Statistics, 4, Glaser, R.E. (198). Bathtub and related failure rate characterizations. J. Am. Statist. Assoc., 75, Laud, P.W. (1977). Bayesian nonparametric inference in reliability. Ph.D. Dissertation, University of Missouri, Columbia, MO. Laud, P.W., Smith, A.F.M. and Damien, P. (1996). Monte Carlo methods for approximating a posterior hazard process. Statistics and Computing, 6, Lawless, J.F. (1982). New York. Statistical Models & Methods For Lifetime Data, Wiley Publications, Leonard, T. (1978). Density estimation, stochastic processes and prior information. J. Roy. 12
14 Statist. Soc. B, 4, Sinha, D. and Dey, D.K. (1997). Semiparametric Bayesian analysis of survival data. J. Am. Statist. Assoc., 92, Smith, A.F.M. and Roberts, G.O. (1993). Bayesian computations via the Gibbs sampler and related Markov chain Monte Carlo methods. J. Roy. Statist. Soc. B 55, Zhang, M.J. and Klein, J.P. (1998). Condence bands for the dierence of two survival curves under proportional hazards model. Technical Report 29, Division of Biostatistics, Medical College of Wisconsin, Milwaukee, Wisconsin. 13
15 14
16 15
17 16
18 17
Statistical Inference for Stochastic Epidemic Models
Statistical Inference for Stochastic Epidemic Models George Streftaris 1 and Gavin J. Gibson 1 1 Department of Actuarial Mathematics & Statistics, Heriot-Watt University, Riccarton, Edinburgh EH14 4AS,
More informationBAYESIAN NON-PARAMETRIC SIMULATION OF HAZARD FUNCTIONS
Proceedings of the 2009 Winter Simulation Conference. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. BAYESIAN NON-PARAETRIC SIULATION OF HAZARD FUNCTIONS DITRIY BELYI Operations
More informationMultistate models in survival and event history analysis
Multistate models in survival and event history analysis Dorota M. Dabrowska UCLA November 8, 2011 Research supported by the grant R01 AI067943 from NIAID. The content is solely the responsibility of the
More informationDynamic Prediction of Disease Progression Using Longitudinal Biomarker Data
Dynamic Prediction of Disease Progression Using Longitudinal Biomarker Data Xuelin Huang Department of Biostatistics M. D. Anderson Cancer Center The University of Texas Joint Work with Jing Ning, Sangbum
More informationDuration Analysis. Joan Llull
Duration Analysis Joan Llull Panel Data and Duration Models Barcelona GSE joan.llull [at] movebarcelona [dot] eu Introduction Duration Analysis 2 Duration analysis Duration data: how long has an individual
More informationSurvival Analysis Math 434 Fall 2011
Survival Analysis Math 434 Fall 2011 Part IV: Chap. 8,9.2,9.3,11: Semiparametric Proportional Hazards Regression Jimin Ding Math Dept. www.math.wustl.edu/ jmding/math434/fall09/index.html Basic Model Setup
More informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationThe iterative convex minorant algorithm for nonparametric estimation
The iterative convex minorant algorithm for nonparametric estimation Report 95-05 Geurt Jongbloed Technische Universiteit Delft Delft University of Technology Faculteit der Technische Wiskunde en Informatica
More informationRonald Christensen. University of New Mexico. Albuquerque, New Mexico. Wesley Johnson. University of California, Irvine. Irvine, California
Texts in Statistical Science Bayesian Ideas and Data Analysis An Introduction for Scientists and Statisticians Ronald Christensen University of New Mexico Albuquerque, New Mexico Wesley Johnson University
More informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationCTDL-Positive Stable Frailty Model
CTDL-Positive Stable Frailty Model M. Blagojevic 1, G. MacKenzie 2 1 Department of Mathematics, Keele University, Staffordshire ST5 5BG,UK and 2 Centre of Biostatistics, University of Limerick, Ireland
More informationMulti-state Models: An Overview
Multi-state Models: An Overview Andrew Titman Lancaster University 14 April 2016 Overview Introduction to multi-state modelling Examples of applications Continuously observed processes Intermittently observed
More informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
More informationGROUPED SURVIVAL DATA. Florida State University and Medical College of Wisconsin
FITTING COX'S PROPORTIONAL HAZARDS MODEL USING GROUPED SURVIVAL DATA Ian W. McKeague and Mei-Jie Zhang Florida State University and Medical College of Wisconsin Cox's proportional hazard model is often
More informationBayesian Nonparametric Regression for Diabetes Deaths
Bayesian Nonparametric Regression for Diabetes Deaths Brian M. Hartman PhD Student, 2010 Texas A&M University College Station, TX, USA David B. Dahl Assistant Professor Texas A&M University College Station,
More informationInference on the Univariate Frailty Model: An Approach Bayesian Reference Analysis
Inference on the Univariate Frailty Model: An Approach Bayesian Reference Analysis Vera Lucia D. Tomazella and Camila Bertini Martins Universidade Federal de São Carlos, SP-Brasil (vera@ufscar.br; cacamila-bertini@yahoo.com.br)
More informationOther Survival Models. (1) Non-PH models. We briefly discussed the non-proportional hazards (non-ph) model
Other Survival Models (1) Non-PH models We briefly discussed the non-proportional hazards (non-ph) model λ(t Z) = λ 0 (t) exp{β(t) Z}, where β(t) can be estimated by: piecewise constants (recall how);
More informationGoodness-of-Fit Tests With Right-Censored Data by Edsel A. Pe~na Department of Statistics University of South Carolina Colloquium Talk August 31, 2 Research supported by an NIH Grant 1 1. Practical Problem
More informationThe simple slice sampler is a specialised type of MCMC auxiliary variable method (Swendsen and Wang, 1987; Edwards and Sokal, 1988; Besag and Green, 1
Recent progress on computable bounds and the simple slice sampler by Gareth O. Roberts* and Jerey S. Rosenthal** (May, 1999.) This paper discusses general quantitative bounds on the convergence rates of
More informationYou know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?
You know I m not goin diss you on the internet Cause my mama taught me better than that I m a survivor (What?) I m not goin give up (What?) I m not goin stop (What?) I m goin work harder (What?) Sir David
More informationFrailty Models and Copulas: Similarities and Differences
Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationPhilosophy and Features of the mstate package
Introduction Mathematical theory Practice Discussion Philosophy and Features of the mstate package Liesbeth de Wreede, Hein Putter Department of Medical Statistics and Bioinformatics Leiden University
More informationA general mixed model approach for spatio-temporal regression data
A general mixed model approach for spatio-temporal regression data Thomas Kneib, Ludwig Fahrmeir & Stefan Lang Department of Statistics, Ludwig-Maximilians-University Munich 1. Spatio-temporal regression
More informationMultistate models and recurrent event models
Multistate models Multistate models and recurrent event models Patrick Breheny December 10 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/22 Introduction Multistate models In this final lecture,
More informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and
More informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationRank Regression with Normal Residuals using the Gibbs Sampler
Rank Regression with Normal Residuals using the Gibbs Sampler Stephen P Smith email: hucklebird@aol.com, 2018 Abstract Yu (2000) described the use of the Gibbs sampler to estimate regression parameters
More informationOn the Fisher Bingham Distribution
On the Fisher Bingham Distribution BY A. Kume and S.G Walker Institute of Mathematics, Statistics and Actuarial Science, University of Kent Canterbury, CT2 7NF,UK A.Kume@kent.ac.uk and S.G.Walker@kent.ac.uk
More informationAnalysis of competing risks data and simulation of data following predened subdistribution hazards
Analysis of competing risks data and simulation of data following predened subdistribution hazards Bernhard Haller Institut für Medizinische Statistik und Epidemiologie Technische Universität München 27.05.2013
More informationContents. Part I: Fundamentals of Bayesian Inference 1
Contents Preface xiii Part I: Fundamentals of Bayesian Inference 1 1 Probability and inference 3 1.1 The three steps of Bayesian data analysis 3 1.2 General notation for statistical inference 4 1.3 Bayesian
More informationTowards a Bayesian model for Cyber Security
Towards a Bayesian model for Cyber Security Mark Briers (mbriers@turing.ac.uk) Joint work with Henry Clausen and Prof. Niall Adams (Imperial College London) 27 September 2017 The Alan Turing Institute
More information5. Parametric Regression Model
5. Parametric Regression Model The Accelerated Failure Time (AFT) Model Denote by S (t) and S 2 (t) the survival functions of two populations. The AFT model says that there is a constant c > 0 such that
More informationBayesian nonparametric inference for random distributions and related functions
J. R. Statist. Soc. B(1999) 61, Part 3, pp. 485^527 Bayesian nonparametric inference for random distributions and related functions Stephen G. Walker, Imperial College of Science, Technology and Medicine,
More informationFahrmeir: Discrete failure time models
Fahrmeir: Discrete failure time models Sonderforschungsbereich 386, Paper 91 (1997) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner Discrete failure time models Ludwig Fahrmeir, Universitat
More informationMultistate models and recurrent event models
and recurrent event models Patrick Breheny December 6 Patrick Breheny University of Iowa Survival Data Analysis (BIOS:7210) 1 / 22 Introduction In this final lecture, we will briefly look at two other
More informationCURE FRACTION MODELS USING MIXTURE AND NON-MIXTURE MODELS. 1. Introduction
ØÑ Å Ø Ñ Ø Ð ÈÙ Ð Ø ÓÒ DOI: 10.2478/v10127-012-0001-4 Tatra Mt. Math. Publ. 51 (2012), 1 9 CURE FRACTION MODELS USING MIXTURE AND NON-MIXTURE MODELS Jorge A. Achcar Emílio A. Coelho-Barros Josmar Mazucheli
More informationOnline Appendix to: Marijuana on Main Street? Estimating Demand in Markets with Limited Access
Online Appendix to: Marijuana on Main Street? Estating Demand in Markets with Lited Access By Liana Jacobi and Michelle Sovinsky This appendix provides details on the estation methodology for various speci
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Estimation of the primary hazard ratio in the presence of a secondary covariate with non-proportional hazards An undergraduate honors thesis submitted to the Department
More informationApproximation of Survival Function by Taylor Series for General Partly Interval Censored Data
Malaysian Journal of Mathematical Sciences 11(3): 33 315 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Approximation of Survival Function by Taylor
More informationON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS. Abstract. We introduce a wide class of asymmetric loss functions and show how to obtain
ON STATISTICAL INFERENCE UNDER ASYMMETRIC LOSS FUNCTIONS Michael Baron Received: Abstract We introduce a wide class of asymmetric loss functions and show how to obtain asymmetric-type optimal decision
More informationPROD. TYPE: COM ARTICLE IN PRESS. Computational Statistics & Data Analysis ( )
COMSTA 28 pp: -2 (col.fig.: nil) PROD. TYPE: COM ED: JS PAGN: Usha.N -- SCAN: Bindu Computational Statistics & Data Analysis ( ) www.elsevier.com/locate/csda Transformation approaches for the construction
More informationPower Calculations for Preclinical Studies Using a K-Sample Rank Test and the Lehmann Alternative Hypothesis
Power Calculations for Preclinical Studies Using a K-Sample Rank Test and the Lehmann Alternative Hypothesis Glenn Heller Department of Epidemiology and Biostatistics, Memorial Sloan-Kettering Cancer Center,
More informationBayesian semiparametric analysis of short- and long- term hazard ratios with covariates
Bayesian semiparametric analysis of short- and long- term hazard ratios with covariates Department of Statistics, ITAM, Mexico 9th Conference on Bayesian nonparametrics Amsterdam, June 10-14, 2013 Contents
More informationOptimum Hybrid Censoring Scheme using Cost Function Approach
Optimum Hybrid Censoring Scheme using Cost Function Approach Ritwik Bhattacharya 1, Biswabrata Pradhan 1, Anup Dewanji 2 1 SQC and OR Unit, Indian Statistical Institute, 203, B. T. Road, Kolkata, PIN-
More informationINVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION
Pak. J. Statist. 2017 Vol. 33(1), 37-61 INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION A. M. Abd AL-Fattah, A.A. EL-Helbawy G.R. AL-Dayian Statistics Department, Faculty of Commerce, AL-Azhar
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationPackage threg. August 10, 2015
Package threg August 10, 2015 Title Threshold Regression Version 1.0.3 Date 2015-08-10 Author Tao Xiao Maintainer Tao Xiao Depends R (>= 2.10), survival, Formula Fit a threshold regression
More informationBayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units
Bayesian nonparametric estimation of finite population quantities in absence of design information on nonsampled units Sahar Z Zangeneh Robert W. Keener Roderick J.A. Little Abstract In Probability proportional
More informationFrailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Mela. P.
Frailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Melanie M. Wall, Bradley P. Carlin November 24, 2014 Outlines of the talk
More informationStochastic dominance with imprecise information
Stochastic dominance with imprecise information Ignacio Montes, Enrique Miranda, Susana Montes University of Oviedo, Dep. of Statistics and Operations Research. Abstract Stochastic dominance, which is
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods By Oleg Makhnin 1 Introduction a b c M = d e f g h i 0 f(x)dx 1.1 Motivation 1.1.1 Just here Supresses numbering 1.1.2 After this 1.2 Literature 2 Method 2.1 New math As
More informationPENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA
PENALIZED LIKELIHOOD PARAMETER ESTIMATION FOR ADDITIVE HAZARD MODELS WITH INTERVAL CENSORED DATA Kasun Rathnayake ; A/Prof Jun Ma Department of Statistics Faculty of Science and Engineering Macquarie University
More informationNonparametric Bayesian Survival Analysis using Mixtures of Weibull Distributions
Nonparametric Bayesian Survival Analysis using Mixtures of Weibull Distributions ATHANASIOS KOTTAS University of California, Santa Cruz ABSTRACT. Bayesian nonparametric methods have been applied to survival
More information2 Ishwaran and James: Gibbs sampling stick-breaking priors Our second method, the blocked Gibbs sampler, works in greater generality in that it can be
Hemant ISHWARAN 1 and Lancelot F. JAMES 2 Gibbs Sampling Methods for Stick-Breaking Priors A rich and exible class of random probability measures, which we call stick-breaking priors, can be constructed
More informationResolving GRB Light Curves
Resolving GRB Light Curves Robert L. Wolpert Duke University wolpert@stat.duke.edu w/mary E Broadbent (Duke) & Tom Loredo (Cornell) 2014 August 03 17:15{17:45 Robert L. Wolpert Resolving GRB Light Curves
More informationA COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky
A COMPARISON OF POISSON AND BINOMIAL EMPIRICAL LIKELIHOOD Mai Zhou and Hui Fang University of Kentucky Empirical likelihood with right censored data were studied by Thomas and Grunkmier (1975), Li (1995),
More informationBayesian Predictions for Exponentially Distributed Failure Times With One Change-Point
c Heldermann Verlag Economic Quality Control ISSN 0940-5151 Vol 18 2003, No. 2, 195 207 Bayesian Predictions for Exponentially Distributed Failure Times With One Change-Point Y. Abdel-Aty, J. Franz, M.A.W.
More informationREGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520
REGRESSION ANALYSIS FOR TIME-TO-EVENT DATA THE PROPORTIONAL HAZARDS (COX) MODEL ST520 Department of Statistics North Carolina State University Presented by: Butch Tsiatis, Department of Statistics, NCSU
More informationLikelihood Ratio Tests and Intersection-Union Tests. Roger L. Berger. Department of Statistics, North Carolina State University
Likelihood Ratio Tests and Intersection-Union Tests by Roger L. Berger Department of Statistics, North Carolina State University Raleigh, NC 27695-8203 Institute of Statistics Mimeo Series Number 2288
More informationMultivariate spatial modeling
Multivariate spatial modeling Point-referenced spatial data often come as multivariate measurements at each location Chapter 7: Multivariate Spatial Modeling p. 1/21 Multivariate spatial modeling Point-referenced
More information1. INTRODUCTION Propp and Wilson (1996,1998) described a protocol called \coupling from the past" (CFTP) for exact sampling from a distribution using
Ecient Use of Exact Samples by Duncan J. Murdoch* and Jerey S. Rosenthal** Abstract Propp and Wilson (1996,1998) described a protocol called coupling from the past (CFTP) for exact sampling from the steady-state
More informationAnalysis of incomplete data in presence of competing risks
Journal of Statistical Planning and Inference 87 (2000) 221 239 www.elsevier.com/locate/jspi Analysis of incomplete data in presence of competing risks Debasis Kundu a;, Sankarshan Basu b a Department
More informationMixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals. John W. Mac McDonald & Alessandro Rosina
Mixture modelling of recurrent event times with long-term survivors: Analysis of Hutterite birth intervals John W. Mac McDonald & Alessandro Rosina Quantitative Methods in the Social Sciences Seminar -
More informationPractical Bayesian Quantile Regression. Keming Yu University of Plymouth, UK
Practical Bayesian Quantile Regression Keming Yu University of Plymouth, UK (kyu@plymouth.ac.uk) A brief summary of some recent work of us (Keming Yu, Rana Moyeed and Julian Stander). Summary We develops
More informationEvaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions
Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions C. Xing, R. Caspeele, L. Taerwe Ghent University, Department
More informationPhysician Performance Assessment / Spatial Inference of Pollutant Concentrations
Physician Performance Assessment / Spatial Inference of Pollutant Concentrations Dawn Woodard Operations Research & Information Engineering Cornell University Johns Hopkins Dept. of Biostatistics, April
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
More informationSurvival Distributions, Hazard Functions, Cumulative Hazards
BIO 244: Unit 1 Survival Distributions, Hazard Functions, Cumulative Hazards 1.1 Definitions: The goals of this unit are to introduce notation, discuss ways of probabilistically describing the distribution
More informationCOMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH *
Jordan Journal of Mathematics and Statistics JJMS 4,, pp. 6-78 COMPARISON OF RELATIVE RISK FUNCTIONS OF THE RAYLEIGH DISTRIBUTION UNDER TYPE-II CENSORED SAMPLES: BAYESIAN APPROACH * Sanku Dey ABSTRACT:
More informationChapter 4 Fall Notations: t 1 < t 2 < < t D, D unique death times. d j = # deaths at t j = n. Y j = # at risk /alive at t j = n
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 4 Fall 2012 4.2 Estimators of the survival and cumulative hazard functions for RC data Suppose X is a continuous random failure time with
More informationTextbook: Survivial Analysis Techniques for Censored and Truncated Data 2nd edition, by Klein and Moeschberger
Lecturer: James Degnan Office: SMLC 342 Office hours: MW 12:00 1:00 or by appointment E-mail: jamdeg@unm.edu Please include STAT474 or STAT574 in the subject line of the email to make sure I don t overlook
More informationRene Tabanera y Palacios 4. Danish Epidemiology Science Center. Novo Nordisk A/S Gentofte. September 1, 1995
Estimation of variance in Cox's regression model with gamma frailties. Per Kragh Andersen 2 John P. Klein 3 Kim M. Knudsen 2 Rene Tabanera y Palacios 4 Department of Biostatistics, University of Copenhagen,
More informationGibbs Sampling in Latent Variable Models #1
Gibbs Sampling in Latent Variable Models #1 Econ 690 Purdue University Outline 1 Data augmentation 2 Probit Model Probit Application A Panel Probit Panel Probit 3 The Tobit Model Example: Female Labor
More informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationAnalysis of Gamma and Weibull Lifetime Data under a General Censoring Scheme and in the presence of Covariates
Communications in Statistics - Theory and Methods ISSN: 0361-0926 (Print) 1532-415X (Online) Journal homepage: http://www.tandfonline.com/loi/lsta20 Analysis of Gamma and Weibull Lifetime Data under a
More informationIntegrated Likelihood Estimation in Semiparametric Regression Models. Thomas A. Severini Department of Statistics Northwestern University
Integrated Likelihood Estimation in Semiparametric Regression Models Thomas A. Severini Department of Statistics Northwestern University Joint work with Heping He, University of York Introduction Let Y
More informationYounshik Chung and Hyungsoon Kim 968). Sharples(990) showed how variance ination can be incorporated easily into general hierarchical models, retainin
Bayesian Outlier Detection in Regression Model Younshik Chung and Hyungsoon Kim Abstract The problem of 'outliers', observations which look suspicious in some way, has long been one of the most concern
More informationCensoring mechanisms
Censoring mechanisms Patrick Breheny September 3 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Fixed vs. random censoring In the previous lecture, we derived the contribution to the likelihood
More informationParametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored Data
Columbia International Publishing Journal of Advanced Computing (2013) 1: 43-58 doi:107726/jac20131004 Research Article Parametric Maximum Likelihood Estimation of Cure Fraction Using Interval-Censored
More informationMarkov Chain Monte Carlo in Practice
Markov Chain Monte Carlo in Practice Edited by W.R. Gilks Medical Research Council Biostatistics Unit Cambridge UK S. Richardson French National Institute for Health and Medical Research Vilejuif France
More informationOr How to select variables Using Bayesian LASSO
Or How to select variables Using Bayesian LASSO x 1 x 2 x 3 x 4 Or How to select variables Using Bayesian LASSO x 1 x 2 x 3 x 4 Or How to select variables Using Bayesian LASSO On Bayesian Variable Selection
More informationStructural Econometrics: Dynamic Discrete Choice. Jean-Marc Robin
Structural Econometrics: Dynamic Discrete Choice Jean-Marc Robin 1. Dynamic discrete choice models 2. Application: college and career choice Plan 1 Dynamic discrete choice models See for example the presentation
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 24 Paper 153 A Note on Empirical Likelihood Inference of Residual Life Regression Ying Qing Chen Yichuan
More informationPrerequisite: STATS 7 or STATS 8 or AP90 or (STATS 120A and STATS 120B and STATS 120C). AP90 with a minimum score of 3
University of California, Irvine 2017-2018 1 Statistics (STATS) Courses STATS 5. Seminar in Data Science. 1 Unit. An introduction to the field of Data Science; intended for entering freshman and transfers.
More informationDepartment of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim
Tests for trend in more than one repairable system. Jan Terje Kvaly Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim ABSTRACT: If failure time data from several
More informationSemi-Competing Risks on A Trivariate Weibull Survival Model
Semi-Competing Risks on A Trivariate Weibull Survival Model Cheng K. Lee Department of Targeting Modeling Insight & Innovation Marketing Division Wachovia Corporation Charlotte NC 28244 Jenq-Daw Lee Graduate
More informationRidge analysis of mixture response surfaces
Statistics & Probability Letters 48 (2000) 3 40 Ridge analysis of mixture response surfaces Norman R. Draper a;, Friedrich Pukelsheim b a Department of Statistics, University of Wisconsin, 20 West Dayton
More informationJoint Modeling of Longitudinal Item Response Data and Survival
Joint Modeling of Longitudinal Item Response Data and Survival Jean-Paul Fox University of Twente Department of Research Methodology, Measurement and Data Analysis Faculty of Behavioural Sciences Enschede,
More informationAugustin: Some Basic Results on the Extension of Quasi-Likelihood Based Measurement Error Correction to Multivariate and Flexible Structural Models
Augustin: Some Basic Results on the Extension of Quasi-Likelihood Based Measurement Error Correction to Multivariate and Flexible Structural Models Sonderforschungsbereich 386, Paper 196 (2000) Online
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
More informationInstrumental variables estimation in the Cox Proportional Hazard regression model
Instrumental variables estimation in the Cox Proportional Hazard regression model James O Malley, Ph.D. Department of Biomedical Data Science The Dartmouth Institute for Health Policy and Clinical Practice
More informationBayesian Analysis. Bayesian Analysis: Bayesian methods concern one s belief about θ. [Current Belief (Posterior)] (Prior Belief) x (Data) Outline
Bayesian Analysis DuBois Bowman, Ph.D. Gordana Derado, M. S. Shuo Chen, M. S. Department of Biostatistics and Bioinformatics Center for Biomedical Imaging Statistics Emory University Outline I. Introduction
More informationvar D (B) = var(b? E D (B)) = var(b)? cov(b; D)(var(D))?1 cov(d; B) (2) Stone [14], and Hartigan [9] are among the rst to discuss the role of such ass
BAYES LINEAR ANALYSIS [This article appears in the Encyclopaedia of Statistical Sciences, Update volume 3, 1998, Wiley.] The Bayes linear approach is concerned with problems in which we want to combine
More informationBIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca Lenz-Tönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
More informationBayesian Graphical Models
Graphical Models and Inference, Lecture 16, Michaelmas Term 2009 December 4, 2009 Parameter θ, data X = x, likelihood L(θ x) p(x θ). Express knowledge about θ through prior distribution π on θ. Inference
More informationThe Poisson transform for unnormalised statistical models. Nicolas Chopin (ENSAE) joint work with Simon Barthelmé (CNRS, Gipsa-LAB)
The Poisson transform for unnormalised statistical models Nicolas Chopin (ENSAE) joint work with Simon Barthelmé (CNRS, Gipsa-LAB) Part I Unnormalised statistical models Unnormalised statistical models
More informationST5212: Survival Analysis
ST51: Survival Analysis 8/9: Semester II Tutorial 1. A model for lifetimes, with a bathtub-shaped hazard rate, is the exponential power distribution with survival fumction S(x) =exp{1 exp[(λx) α ]}. (a)
More information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More information