Bayesian semiparametric analysis of short- and long- term hazard ratios with covariates
|
|
- Evelyn French
- 5 years ago
- Views:
Transcription
1 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
2 Contents Motivating example Model Extension to covariates Bayesian inference Data analysis
3 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA)
4 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA) ESA are used to relieve chemotherapy induced anemia in patients with cancer
5 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA) ESA are used to relieve chemotherapy induced anemia in patients with cancer EPO is a glycoprotein hormone that controls red blood cell production
6 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA) ESA are used to relieve chemotherapy induced anemia in patients with cancer EPO is a glycoprotein hormone that controls red blood cell production In contrast, other studies have shown that ESA treatments can compromise survival for some types of cancer
7 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA) ESA are used to relieve chemotherapy induced anemia in patients with cancer EPO is a glycoprotein hormone that controls red blood cell production In contrast, other studies have shown that ESA treatments can compromise survival for some types of cancer It was believed that effects of ESA was induced via the canonical EPO receptor (EpoR)
8 Study Survival Study: Ovarian cancer patients who have undergone a treatment with erythoprotein (EPO) stimulating agents (ESA) ESA are used to relieve chemotherapy induced anemia in patients with cancer EPO is a glycoprotein hormone that controls red blood cell production In contrast, other studies have shown that ESA treatments can compromise survival for some types of cancer It was believed that effects of ESA was induced via the canonical EPO receptor (EpoR) Objective: test whether the new ephrin-type B receptor 4 (EphB4) can better explain the effect of ESA.
9 Study n = 174 patients
10 Study n = 174 patients Covariates: X 1 = ESA treatment indicator (0 or 1). 51% received treat. X 2 = EphB4 level of expression (low or high). 39% high expsn. X 3 = EpoR level of expression (low or high) X 4 = Cytoreduction (suboptimal and optimal) X 5 = Disease s stage (I V) X 6 = Age (20 92). Average of 58 years
11 Study n = 174 patients Covariates: X 1 = ESA treatment indicator (0 or 1). 51% received treat. X 2 = EphB4 level of expression (low or high). 39% high expsn. X 3 = EpoR level of expression (low or high) X 4 = Cytoreduction (suboptimal and optimal) X 5 = Disease s stage (I V) X 6 = Age (20 92). Average of 58 years Tested proportionality assumption for each covariate. Found X 1 and X 4 with non proportional behaviour
12 Study K M survival curves Time Figure : KM curves for ESA treatment indicator.
13 Definitions In survival analysis, let
14 Definitions In survival analysis, let h T (t) hazard rate for the treatment group
15 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group
16 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group We define the hazard ratio as HR(t) = h T (t) h C (t)
17 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group We define the hazard ratio as HR(t) = h T (t) h C (t) There are three most common behaviours for HR
18 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group We define the hazard ratio as HR(t) = h T (t) h C (t) There are three most common behaviours for HR HR(t) = c proportional hazards
19 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group We define the hazard ratio as HR(t) = h T (t) h C (t) There are three most common behaviours for HR HR(t) = c proportional hazards lim t HR(t) = 1 proportional odds model
20 Definitions In survival analysis, let h T (t) hazard rate for the treatment group h C (t) hazard rate for the control group We define the hazard ratio as HR(t) = h T (t) h C (t) There are three most common behaviours for HR HR(t) = c proportional hazards lim t HR(t) = 1 proportional odds model Non-monotone HR(t). If HR(t) < 1 and HR(t) > 1, for some t, the hazards cross
21 Semiparametric model Yang and Prentice (2005) proposed a model that accounts for the previous three behaviours of a HR
22 Semiparametric model Yang and Prentice (2005) proposed a model that accounts for the previous three behaviours of a HR Their model is S C (t) = {1 + R(t)} 1 and S T (t) = {1 + λθ } θ R(t) This implies that HR(t) = λθ λ + (θ λ)/{1 + R(t)} where λ, θ > 0 and R(t) is monotone nondecreasing
23 Semiparametric model The model s name is justified by lim HR(t) = λ and lim HR(t) = θ t 0 t
24 Semiparametric model The model s name is justified by lim HR(t) = λ and lim HR(t) = θ t 0 t If λ = θ proportional hazards
25 Semiparametric model The model s name is justified by lim HR(t) = λ and lim HR(t) = θ t 0 t If λ = θ proportional hazards If θ = 1 proportional odds model
26 Semiparametric model The model s name is justified by lim HR(t) = λ and lim HR(t) = θ t 0 t If λ = θ proportional hazards If θ = 1 proportional odds model If (λ 1)(θ 1) < 0 hazard functions cross
27 Semiparametric model The model s name is justified by lim HR(t) = λ and lim HR(t) = θ t 0 t If λ = θ proportional hazards If θ = 1 proportional odds model If (λ 1)(θ 1) < 0 hazard functions cross The model is semiparametric because it includes: (λ, θ) and R(t) = 1 S C (t) S C (t) is the odds function of the control group
28 Semiparametric model S(t) Figure : Varying (λ, θ): (1, 1) (solid line), (2, 0.5) (dashed line), and (0.5, 3) (dotted-dashed line). t
29 Extensions Nieto-Barajas (2013) extended the model to include covariates with S(t x i ) = { 1 + λ } θi i R(t) θ i λ i = exp(δ x i ) and θ i = exp(γ x i ) where δ and γ do not have intercept
30 Extensions Nieto-Barajas (2013) extended the model to include covariates with S(t x i ) = { 1 + λ } θi i R(t) θ i λ i = exp(δ x i ) and θ i = exp(γ x i ) where δ and γ do not have intercept Here we can have more than two groups
31 Extensions Nieto-Barajas (2013) extended the model to include covariates with S(t x i ) = { 1 + λ } θi i R(t) θ i λ i = exp(δ x i ) and θ i = exp(γ x i ) where δ and γ do not have intercept Here we can have more than two groups HR i (t) = h(t x i )/h(t 0) has the same interpretation as in the two groups model
32 Extensions Nieto-Barajas (2013) extended the model to include covariates with S(t x i ) = { 1 + λ } θi i R(t) θ i λ i = exp(δ x i ) and θ i = exp(γ x i ) where δ and γ do not have intercept Here we can have more than two groups HR i (t) = h(t x i )/h(t 0) has the same interpretation as in the two groups model Some individuals can have prop. hazards, some others prop. odds and some others crossing hazards
33 Follow a Bayesian approach for inference on δ, γ and R(t)
34 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t):
35 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function
36 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) =
37 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.)
38 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.) There is a vast list of priors for c.h.f.:
39 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.) There is a vast list of priors for c.h.f.: beta and beta-stacy (NTR) (Hjort, 1990; Walker and Muliere, 1997)
40 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.) There is a vast list of priors for c.h.f.: beta and beta-stacy (NTR) (Hjort, 1990; Walker and Muliere, 1997) mixtures of gamma and weighted gamma (Dykstra and Laud, 1981; Lo and Weng, 1989)
41 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.) There is a vast list of priors for c.h.f.: beta and beta-stacy (NTR) (Hjort, 1990; Walker and Muliere, 1997) mixtures of gamma and weighted gamma (Dykstra and Laud, 1981; Lo and Weng, 1989) Lévy-driven proceses (Nieto-Barajas and Walker, 2004)
42 Follow a Bayesian approach for inference on δ, γ and R(t) Concentrate first on R(t): is monotone nondecreasing function lim t 0 R(t) = 0 and lim t R(t) = R(t) behaves like a cumulative hazard function (c.h.f.) There is a vast list of priors for c.h.f.: beta and beta-stacy (NTR) (Hjort, 1990; Walker and Muliere, 1997) mixtures of gamma and weighted gamma (Dykstra and Laud, 1981; Lo and Weng, 1989) Lévy-driven proceses (Nieto-Barajas and Walker, 2004) We will use the class of Lévy-driven processes
43 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds)
44 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t
45 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t lim t 0 k(t, s) = 0 for all s 0
46 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t lim t 0 k(t, s) = 0 for all s 0 This is a general class
47 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t lim t 0 k(t, s) = 0 for all s 0 This is a general class If k(t, s) = h(t, s)i (s t)β(s), and r(s) a gamma p. L&W(89)
48 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t lim t 0 k(t, s) = 0 for all s 0 This is a general class If k(t, s) = h(t, s)i (s t)β(s), and r(s) a gamma p. L&W(89) h(t, s) = (t s), and r(s) an extended gamma p. D&L(81)
49 Lévy-driven prior: mixture of kernel k(t, s) with respect to a Lévy process (or IAP) r(s) R(t) = 0 k(t, s)r(ds) k(t, s) is nondecreasing a right-continuous as a function of t lim t 0 k(t, s) = 0 for all s 0 This is a general class If k(t, s) = h(t, s)i (s t)β(s), and r(s) a gamma p. L&W(89) h(t, s) = (t s), and r(s) an extended gamma p. D&L(81) k(t, s) = { 1 e a(t s)} I (s t)), and r(s) a finite p. N-B&W(04)
50 For the kernel we will take
51 For the kernel we will take Location Weibull kernel k(t, s) = {1 e a(t s)b }I (s t), a, b > 0
52 For the kernel we will take Location Weibull kernel k(t, s) = {1 e a(t s)b }I (s t), a, b > 0 Mean Weibull kernel k(t, s) = 1 e {Γ(1+b 1 )t/s} b, b > 0
53 R(t) Figure : Kernel comparison. Loc. exp. (dotted line), loc. Weib. (solid line), and mean Weib. (dashed line). t
54 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take
55 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take A non-homogeneous version of the GG p. ρ(dv s) = Γ(1 ɛ) 1 v (1+ɛ) e β(s)v dv, where
56 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take A non-homogeneous version of the GG p. ρ(dv s) = Γ(1 ɛ) 1 v (1+ɛ) e β(s)v dv, where ɛ {(0, 1) { 1}}
57 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take A non-homogeneous version of the GG p. ρ(dv s) = Γ(1 ɛ) 1 v (1+ɛ) e β(s)v dv, where ɛ {(0, 1) { 1}} β(s) = cs d, c > 0, d { 1, 0, 1} s 0
58 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take A non-homogeneous version of the GG p. ρ(dv s) = Γ(1 ɛ) 1 v (1+ɛ) e β(s)v dv, where ɛ {(0, 1) { 1}} β(s) = cs d, c > 0, d { 1, 0, 1} s 0 Together with α(s) = ωp 0 (s)
59 For the Lévy intensity ν(dv, ds) = ρ(dv s)α(ds) we will take A non-homogeneous version of the GG p. ρ(dv s) = Γ(1 ɛ) 1 v (1+ɛ) e β(s)v dv, where ɛ {(0, 1) { 1}} β(s) = cs d, c > 0, d { 1, 0, 1} s 0 Together with α(s) = ωp 0 (s) In notation r nh-gg(ɛ, β( ), ω; P 0 )
60 E{S(t)} t E{S(t)} t E{S(t)} t E{S(t)} E{S(t)} E{S(t)} t Figure : Prior E{S(t λ, θ)}. Loc.Weib.(a, b) = (1, 2) (top row), mean Weib.b = 2 (bottom row). Decreasing β (first column), constant β (second column), and increasing β (third column). (λ, θ): (1, 1) (solid line), (2, 0.5) (dashed line), and (0.5, 3) (dotted line). t t
61 For the vector parameters δ and γ we will take
62 For the vector parameters δ and γ we will take δ s N(0, σδ 2 ), for s = 1,..., p
63 For the vector parameters δ and γ we will take δ s N(0, σδ 2 ), for s = 1,..., p γ s N(0, σγ), 2 for s = 1,..., p
64 For the vector parameters δ and γ we will take δ s N(0, σδ 2 ), for s = 1,..., p γ s N(0, σγ), 2 for s = 1,..., p δ s and γ s independent a-priori
65 Posterior characterization: via a Gibbs sampler with [r t, δ, γ], [δ t, r, γ], [γ t, r, δ]
66 Posterior characterization: via a Gibbs sampler with [r t, δ, γ], [δ t, r, γ], [γ t, r, δ] The tricky part is to find [r t, δ, γ].
67 Posterior characterization: via a Gibbs sampler with { E [r t, δ, γ], [δ t, r, γ], [γ t, r, δ] The tricky part is to find [r t, δ, γ]. This will be done by computing { E e } φ(s)r(ds) P(A Z) e φ(s)r(ds) } A = E{P(A Z)} where A = (A 1,..., A n ) is the data with A i either an exact or right censored obs.
68 To find [r t, δ, γ] we require latent variables: A pair (u i, w i ) for each exact obs., and A single u i for each right censored obs.
69 To find [r t, δ, γ] we require latent variables: A pair (u i, w i ) for each exact obs., and A single u i for each right censored obs. The u i s are dependent and the w i s are independent
70 To find [r t, δ, γ] we require latent variables: A pair (u i, w i ) for each exact obs., and A single u i for each right censored obs. s are independent [r t, w] = [r t, w, u][u t, w]d u, with [r t, w, u] an updated IAP with fixed jump locations at w i s The u i s are dependent and the w i
71 To find [r t, δ, γ] we require latent variables: A pair (u i, w i ) for each exact obs., and A single u i for each right censored obs. s are independent [r t, w] = [r t, w, u][u t, w]d u, with [r t, w, u] an updated IAP with fixed jump locations at w i s To simplify computations on [δ t, r, γ] and [γ t, r, δ] we introduce: Another latent variable v i for each obs. The u i s are dependent and the w i
72 To find [r t, δ, γ] we require latent variables: A pair (u i, w i ) for each exact obs., and A single u i for each right censored obs. s are independent [r t, w] = [r t, w, u][u t, w]d u, with [r t, w, u] an updated IAP with fixed jump locations at w i s To simplify computations on [δ t, r, γ] and [γ t, r, δ] we introduce: Another latent variable v i for each obs. These v i s are independent with gamma cond. dist. The u i s are dependent and the w i
73 Data analysis Implemented our semiparametric model with 6 covariates Same covariates were used in λ i and θ i
74 Data analysis Implemented our semiparametric model with 6 covariates Same covariates were used in λ i and θ i Table : LPML statistics for different prior settings. (ω, ɛ) (c, d) Loc. Weib. K. mean-s Weib. K. (1, 1) (1, 0.5) (1, 0) (1, 1) (1, 1) (10, 0.5) (1, 0) (1, 1) (1, 1) (1, 0) Prop. Hazards
75 Data analysis R(t) h(t 0) t t Figure : Post. estimates. R(t) (left) and baseline hazard h(t 0) (right). Using mean Weibull kernel (solid) and using location Weibull kernel (dashed).
76 Data analysis Table : Posterior estimates of covariate coefficients δ and γ. Loc. Weib. K. Parameter Mean 95% CI P(< 0) δ ( 0.160, 0.883) δ (0.885, 2.060) δ ( 0.401, 0.944) δ ( 1.219, 0.136) δ ( 0.135, 0.494) δ ( 0.617, 3.719) γ (0.558, 1.685) γ (0.250, 1.591) γ ( 0.894, 0.049) γ ( 0.242, 0.952) γ ( 0.279, 0.338) γ ( 0.834, 5.750) 0.091
77 Data analysis h(t) Figure : Posterior predictive hazard rates. Hypothetical individuals x 1 = (0, 0, 0, 1, 3, 92), x 2 = (0, 1, 0, 1, 3, 20), x 3 = (1, 0, 0, 1, 3, 92), and x 4 = (1, 1, 0, 1, 3, 20). t
78 References Yang, S. and Prentice, R. (2005). Semiparametric analysis of short-term and long-term hazard ratios with two sample survival data. Biometrika 92, Nieto-Barajas, L.E. (2013). Bayesian semiparametric analysis of short- and long-term hazard ratios with covariates. Computational Statistics and Data Analysis. In press. DOI: /j.csda
Asymptotics 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 informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Igor Prünster University of Turin, Collegio Carlo Alberto and ICER Joint work with P. Di Biasi and G. Peccati Workshop on Limit Theorems and Applications Paris, 16th January
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 informationarxiv: v2 [math.st] 27 May 2014
Full Bayesian inference with hazard mixture models Julyan Arbel a, Antonio Lijoi b,a, Bernardo Nipoti c,a, arxiv:145.6628v2 [math.st] 27 May 214 a Collegio Carlo Alberto, via Real Collegio, 3, 124 Moncalieri,
More informationTruncation error of a superposed gamma process in a decreasing order representation
Truncation error of a superposed gamma process in a decreasing order representation Julyan Arbel Inria Grenoble, Université Grenoble Alpes julyan.arbel@inria.fr Igor Prünster Bocconi University, Milan
More informationBayesian Nonparametric Inference Methods for Mean Residual Life Functions
Bayesian Nonparametric Inference Methods for Mean Residual Life Functions Valerie Poynor Department of Applied Mathematics and Statistics, University of California, Santa Cruz April 28, 212 1/3 Outline
More informationReinforced urns and the subdistribution beta-stacy process prior for competing risks analysis
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfè 1 Stefano Peluso 2 Pietro Muliere 1 1 Università Commerciale Luigi Bocconi 2 Università Cattolica
More informationIntroduction to BGPhazard
Introduction to BGPhazard José A. García-Bueno and Luis E. Nieto-Barajas February 11, 2016 Abstract We present BGPhazard, an R package which computes hazard rates from a Bayesian nonparametric view. This
More informationA nonparametric Bayesian approach to inference for non-homogeneous. Poisson processes. Athanasios Kottas 1. (REVISED VERSION August 23, 2006)
A nonparametric Bayesian approach to inference for non-homogeneous Poisson processes Athanasios Kottas 1 Department of Applied Mathematics and Statistics, Baskin School of Engineering, University of California,
More informationOn the Truncation Error of a Superposed Gamma Process
On the Truncation Error of a Superposed Gamma Process Julyan Arbel and Igor Prünster Abstract Completely random measures (CRMs) form a key ingredient of a wealth of stochastic models, in particular in
More informationOn the posterior structure of NRMI
On the posterior structure of NRMI Igor Prünster University of Turin, Collegio Carlo Alberto and ICER Joint work with L.F. James and A. Lijoi Isaac Newton Institute, BNR Programme, 8th August 2007 Outline
More informationOrder-q stochastic processes. Bayesian nonparametric applications
Order-q dependent stochastic processes in Bayesian nonparametric applications Department of Statistics, ITAM, Mexico BNP 2015, Raleigh, NC, USA 25 June, 2015 Contents Order-1 process Application in survival
More informationReinforced urns and the subdistribution beta-stacy process prior for competing risks analysis
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfè Department of Decision Sciences, Bocconi University, 20136 Milan, Italy. E-mail: andrea.arfe@phd.unibocconi.it
More informationFull Bayesian inference with hazard mixture models
ISSN 2279-9362 Full Bayesian inference with hazard mixture models Julyan Arbel Antonio Lijoi Bernardo Nipoti No. 381 December 214 www.carloalberto.org/research/working-papers 214 by Julyan Arbel, Antonio
More informationBayesian semiparametric inference for the accelerated failure time model using hierarchical mixture modeling with N-IG priors
Bayesian semiparametric inference for the accelerated failure time model using hierarchical mixture modeling with N-IG priors Raffaele Argiento 1, Alessandra Guglielmi 2, Antonio Pievatolo 1, Fabrizio
More informationNonparametric Bayesian Inference for Mean Residual. Life Functions in Survival Analysis
Nonparametric Bayesian Inference for Mean Residual Life Functions in Survival Analysis VALERIE POYNOR 1 and ATHANASIOS KOTTAS 2 1 California State University, Fullerton 2 University of California, Santa
More informationDependence structures with applications to actuarial science
with applications to actuarial science Department of Statistics, ITAM, Mexico Recent Advances in Actuarial Mathematics, OAX, MEX October 26, 2015 Contents Order-1 process Application in survival analysis
More informationChapter 2 Inference on Mean Residual Life-Overview
Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate
More informationBAYESIAN INFERENCE OF HIDDEN GAMMA WEAR PROCESS MODEL FOR SURVIVAL DATA WITH TIES
Statistica Sinica 25 (215), 1613-1635 doi:http://dx.doi.org/1.575/ss.212.351 BAYESIAN INFERENCE OF HIDDEN GAMMA WEAR PROCESS MODEL FOR SURVIVAL DATA WITH TIES Arijit Sinha 1, Zhiyi Chi 2 and Ming-Hui Chen
More informationBayesian Nonparametric Autoregressive Models via Latent Variable Representation
Bayesian Nonparametric Autoregressive Models via Latent Variable Representation Maria De Iorio Yale-NUS College Dept of Statistical Science, University College London Collaborators: Lifeng Ye (UCL, London,
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 informationModel Selection in Bayesian Survival Analysis for a Multi-country Cluster Randomized Trial
Model Selection in Bayesian Survival Analysis for a Multi-country Cluster Randomized Trial Jin Kyung Park International Vaccine Institute Min Woo Chae Seoul National University R. Leon Ochiai International
More informationMultivariate Survival Analysis
Multivariate Survival Analysis Previously we have assumed that either (X i, δ i ) or (X i, δ i, Z i ), i = 1,..., n, are i.i.d.. This may not always be the case. Multivariate survival data can arise in
More informationFlexible Regression Modeling using Bayesian Nonparametric Mixtures
Flexible Regression Modeling using Bayesian Nonparametric Mixtures Athanasios Kottas Department of Applied Mathematics and Statistics University of California, Santa Cruz Department of Statistics Brigham
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 informationA Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks
A Bayesian Nonparametric Approach to Causal Inference for Semi-competing risks Y. Xu, D. Scharfstein, P. Mueller, M. Daniels Johns Hopkins, Johns Hopkins, UT-Austin, UF JSM 2018, Vancouver 1 What are semi-competing
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 informationA Bayesian semiparametric partially PH model for clustered time-to-event data
A Bayesian semiparametric partially PH model for clustered time-to-event data Bernardo Nipoti School of Computer Science and Statistics, Trinity College, Dublin, Ireland Alejandro Jara Department of Statistics,
More informationA Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness
A Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness A. Linero and M. Daniels UF, UT-Austin SRC 2014, Galveston, TX 1 Background 2 Working model
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 informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
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 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 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 informationA Nonparametric Approach Using Dirichlet Process for Hierarchical Generalized Linear Mixed Models
Journal of Data Science 8(2010), 43-59 A Nonparametric Approach Using Dirichlet Process for Hierarchical Generalized Linear Mixed Models Jing Wang Louisiana State University Abstract: In this paper, we
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 informationLecture 8: The Metropolis-Hastings Algorithm
30.10.2008 What we have seen last time: Gibbs sampler Key idea: Generate a Markov chain by updating the component of (X 1,..., X p ) in turn by drawing from the full conditionals: X (t) j Two drawbacks:
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 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 informationarxiv: v1 [math.st] 29 Nov 2018
Reinforced urns and the subdistribution beta-stacy process prior for competing risks analysis Andrea Arfé Stefano Peluso Pietro Muliere Accepted manuscript, Scandinavian Journal of Statistics, 2018 arxiv:1811.12304v1
More informationA Bayesian Semi-parametric Survival Model with Longitudinal Markers
A Bayesian Semi-parametric Survival Model with Longitudinal Markers Song Zhang, Peter Müller and Kim-Anh Do Song.Zhang@utsouthwestern.edu Abstract We consider inference for data from a clinical trial of
More informationBayesian Statistics. Debdeep Pati Florida State University. April 3, 2017
Bayesian Statistics Debdeep Pati Florida State University April 3, 2017 Finite mixture model The finite mixture of normals can be equivalently expressed as y i N(µ Si ; τ 1 S i ), S i k π h δ h h=1 δ h
More informationSpatial Survival Analysis via Copulas
1 / 38 Spatial Survival Analysis via Copulas Tim Hanson Department of Statistics University of South Carolina, U.S.A. International Conference on Survival Analysis in Memory of John P. Klein Medical College
More informationNonparametric Bayesian modeling for dynamic ordinal regression relationships
Nonparametric Bayesian modeling for dynamic ordinal regression relationships Athanasios Kottas Department of Applied Mathematics and Statistics, University of California, Santa Cruz Joint work with Maria
More informationBayesian estimation of the discrepancy with misspecified parametric models
Bayesian estimation of the discrepancy with misspecified parametric models Pierpaolo De Blasi University of Torino & Collegio Carlo Alberto Bayesian Nonparametrics workshop ICERM, 17-21 September 2012
More informationBayesian Analysis of Survival Models with Bathtub Hazard Rates Purushottam W. Laud, Ph.D. Paul Damien, Ph.D. Stephen G. Walker, Ph.D. Technical Report
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
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 informationPackage SimSCRPiecewise
Package SimSCRPiecewise July 27, 2016 Type Package Title 'Simulates Univariate and Semi-Competing Risks Data Given Covariates and Piecewise Exponential Baseline Hazards' Version 0.1.1 Author Andrew G Chapple
More informationGaussian Processes for Survival Analysis
Gaussian Processes for Survival Analysis Tamara Fernández Department of Statistics, University of Oxford. Oxford, UK. fernandez@stats.ox.ac.uk Nicolás Rivera Department of Informatics, King s College London.
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationBayesian Nonparametric Accelerated Failure Time Models for Analyzing Heterogeneous Treatment Effects
Bayesian Nonparametric Accelerated Failure Time Models for Analyzing Heterogeneous Treatment Effects Nicholas C. Henderson Thomas A. Louis Gary Rosner Ravi Varadhan Johns Hopkins University September 28,
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 informationConstruction of Dependent Dirichlet Processes based on Poisson Processes
1 / 31 Construction of Dependent Dirichlet Processes based on Poisson Processes Dahua Lin Eric Grimson John Fisher CSAIL MIT NIPS 2010 Outstanding Student Paper Award Presented by Shouyuan Chen Outline
More informationMarginal Survival Modeling through Spatial Copulas
1 / 53 Marginal Survival Modeling through Spatial Copulas Tim Hanson Department of Statistics University of South Carolina, U.S.A. University of Michigan Department of Biostatistics March 31, 2016 2 /
More informationBayesian Nonparametrics
Bayesian Nonparametrics Peter Orbanz Columbia University PARAMETERS AND PATTERNS Parameters P(X θ) = Probability[data pattern] 3 2 1 0 1 2 3 5 0 5 Inference idea data = underlying pattern + independent
More informationThe coxvc_1-1-1 package
Appendix A The coxvc_1-1-1 package A.1 Introduction The coxvc_1-1-1 package is a set of functions for survival analysis that run under R2.1.1 [81]. This package contains a set of routines to fit Cox models
More informationConstrained estimation for binary and survival data
Constrained estimation for binary and survival data Jeremy M. G. Taylor Yong Seok Park John D. Kalbfleisch Biostatistics, University of Michigan May, 2010 () Constrained estimation May, 2010 1 / 43 Outline
More informationIn contrast, parametric techniques (fitting exponential or Weibull, for example) are more focussed, can handle general covariates, but require
Chapter 5 modelling Semi parametric We have considered parametric and nonparametric techniques for comparing survival distributions between different treatment groups. Nonparametric techniques, such as
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 informationCompound Random Measures
Compound Random Measures Jim Griffin (joint work with Fabrizio Leisen) University of Kent Introduction: Two clinical studies 3 CALGB8881 3 CALGB916 2 2 β 1 1 β 1 1 1 5 5 β 1 5 5 β Infinite mixture models
More informationHybrid Censoring; An Introduction 2
Hybrid Censoring; An Introduction 2 Debasis Kundu Department of Mathematics & Statistics Indian Institute of Technology Kanpur 23-rd November, 2010 2 This is a joint work with N. Balakrishnan Debasis Kundu
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 informationA Fully Nonparametric Modeling Approach to. BNP Binary Regression
A Fully Nonparametric Modeling Approach to Binary Regression Maria Department of Applied Mathematics and Statistics University of California, Santa Cruz SBIES, April 27-28, 2012 Outline 1 2 3 Simulation
More informationAnalysis of Cure Rate Survival Data Under Proportional Odds Model
Analysis of Cure Rate Survival Data Under Proportional Odds Model Yu Gu 1,, Debajyoti Sinha 1, and Sudipto Banerjee 2, 1 Department of Statistics, Florida State University, Tallahassee, Florida 32310 5608,
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 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 informationHeriot-Watt University
Heriot-Watt University Heriot-Watt University Research Gateway Prediction of settlement delay in critical illness insurance claims by using the generalized beta of the second kind distribution Dodd, Erengul;
More informationSTAT 331. Accelerated Failure Time Models. Previously, we have focused on multiplicative intensity models, where
STAT 331 Accelerated Failure Time Models Previously, we have focused on multiplicative intensity models, where h t z) = h 0 t) g z). These can also be expressed as H t z) = H 0 t) g z) or S t z) = e Ht
More informationBayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features. Yangxin Huang
Bayesian Inference on Joint Mixture Models for Survival-Longitudinal Data with Multiple Features Yangxin Huang Department of Epidemiology and Biostatistics, COPH, USF, Tampa, FL yhuang@health.usf.edu January
More informationFoundations of Nonparametric Bayesian Methods
1 / 27 Foundations of Nonparametric Bayesian Methods Part II: Models on the Simplex Peter Orbanz http://mlg.eng.cam.ac.uk/porbanz/npb-tutorial.html 2 / 27 Tutorial Overview Part I: Basics Part II: Models
More informationPackage ICBayes. September 24, 2017
Package ICBayes September 24, 2017 Title Bayesian Semiparametric Models for Interval-Censored Data Version 1.1 Date 2017-9-24 Author Chun Pan, Bo Cai, Lianming Wang, and Xiaoyan Lin Maintainer Chun Pan
More informationLecture 1. Introduction Statistics Statistical Methods II. Presented January 8, 2018
Introduction Statistics 211 - Statistical Methods II Presented January 8, 2018 linear models Dan Gillen Department of Statistics University of California, Irvine 1.1 Logistics and Contact Information Lectures:
More informationAn adaptive truncation method for inference in Bayesian nonparametric models
An adaptive truncation method for inference in Bayesian nonparametric models arxiv:1308.045v [stat.co] 1 May 014 J.E. Griffin School of Mathematics, Statistics and Actuarial Science, University of Kent,
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 informationLecture 3. Univariate Bayesian inference: conjugate analysis
Summary Lecture 3. Univariate Bayesian inference: conjugate analysis 1. Posterior predictive distributions 2. Conjugate analysis for proportions 3. Posterior predictions for proportions 4. Conjugate analysis
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 informationDependent Random Measures and Prediction
Dependent Random Measures and Prediction Igor Prünster University of Torino & Collegio Carlo Alberto 10th Bayesian Nonparametric Conference Raleigh, June 26, 2015 Joint wor with: Federico Camerlenghi,
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 information8. Parametric models in survival analysis General accelerated failure time models for parametric regression
8. Parametric models in survival analysis 8.1. General accelerated failure time models for parametric regression The accelerated failure time model Let T be the time to event and x be a vector of covariates.
More informationDirect likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach
Direct likelihood inference on the cause-specific cumulative incidence function: a flexible parametric regression modelling approach Sarwar I Mozumder 1, Mark J Rutherford 1, Paul C Lambert 1,2 1 Biostatistics
More informationST745: Survival Analysis: Parametric
ST745: Survival Analysis: Parametric Eric B. Laber Department of Statistics, North Carolina State University January 13, 2015 ...the statistician knows... that in nature there never was a normal distribution,
More informationBayesian nonparametric models for bipartite graphs
Bayesian nonparametric models for bipartite graphs François Caron Department of Statistics, Oxford Statistics Colloquium, Harvard University November 11, 2013 F. Caron 1 / 27 Bipartite networks Readers/Customers
More informationSurvival Analysis. Stat 526. April 13, 2018
Survival Analysis Stat 526 April 13, 2018 1 Functions of Survival Time Let T be the survival time for a subject Then P [T < 0] = 0 and T is a continuous random variable The Survival function is defined
More informationRobust and accurate inference via a mixture of Gaussian and terrors
Robust and accurate inference via a mixture of Gaussian and terrors Hyungsuk Tak ICTS/SAMSI Time Series Analysis for Synoptic Surveys and Gravitational Wave Astronomy 22 March 2017 Joint work with Justin
More informationDiscussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs
Discussion of Missing Data Methods in Longitudinal Studies: A Review by Ibrahim and Molenberghs Michael J. Daniels and Chenguang Wang Jan. 18, 2009 First, we would like to thank Joe and Geert for a carefully
More information3003 Cure. F. P. Treasure
3003 Cure F. P. reasure November 8, 2000 Peter reasure / November 8, 2000/ Cure / 3003 1 Cure A Simple Cure Model he Concept of Cure A cure model is a survival model where a fraction of the population
More informationCASE STUDY: Bayesian Incidence Analyses from Cross-Sectional Data with Multiple Markers of Disease Severity. Outline:
CASE STUDY: Bayesian Incidence Analyses from Cross-Sectional Data with Multiple Markers of Disease Severity Outline: 1. NIEHS Uterine Fibroid Study Design of Study Scientific Questions Difficulties 2.
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
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 informationSTAT 6350 Analysis of Lifetime Data. Probability Plotting
STAT 6350 Analysis of Lifetime Data Probability Plotting Purpose of Probability Plots Probability plots are an important tool for analyzing data and have been particular popular in the analysis of life
More informationRecent Advances in Bayesian Spatial Survival Modeling
1 / 47 Recent Advances in Bayesian Spatial Survival Modeling Tim Hanson Department of Statistics University of South Carolina, U.S.A. University of Alabama Department of Mathematics April 10, 2015 2 /
More informationHmms with variable dimension structures and extensions
Hmm days/enst/january 21, 2002 1 Hmms with variable dimension structures and extensions Christian P. Robert Université Paris Dauphine www.ceremade.dauphine.fr/ xian Hmm days/enst/january 21, 2002 2 1 Estimating
More informationChapter 2. Data Analysis
Chapter 2 Data Analysis 2.1. Density Estimation and Survival Analysis The most straightforward application of BNP priors for statistical inference is in density estimation problems. Consider the generic
More informationA Bayesian nonparametric dynamic AR model for multiple time series analysis
A Bayesian nonparametric dynamic AR model for multiple time series analysis Luis E. Nieto Barajas (joint with Fernando Quintana) Department of Statistics ITAM, Mexico COBAL V 8 June 2017 Luis E. Nieto
More informationA Bayesian Nonparametric Model for Predicting Disease Status Using Longitudinal Profiles
A Bayesian Nonparametric Model for Predicting Disease Status Using Longitudinal Profiles Jeremy Gaskins Department of Bioinformatics & Biostatistics University of Louisville Joint work with Claudio Fuentes
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 informationThe Log-generalized inverse Weibull Regression Model
The Log-generalized inverse Weibull Regression Model Felipe R. S. de Gusmão Universidade Federal Rural de Pernambuco Cintia M. L. Ferreira Universidade Federal Rural de Pernambuco Sílvio F. A. X. Júnior
More informationBayesian non-parametric model to longitudinally predict churn
Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics
More informationThe Gaussian process latent variable model with Cox regression
The Gaussian process latent variable model with Cox regression James Barrett Institute for Mathematical and Molecular Biomedicine (IMMB), King s College London James Barrett (IMMB, KCL) GP Latent Variable
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
More information