Computational Methods for Multiplicative Intensity Models Using Weighted Gamma Processes: Proportional Hazards, Marked Point Processes, and Panel Count Data

By Ishwaran, Hemant; James, Lancelot F. | Journal of the American Statistical Association, March 2004 | Go to article overview

Computational Methods for Multiplicative Intensity Models Using Weighted Gamma Processes: Proportional Hazards, Marked Point Processes, and Panel Count Data


Ishwaran, Hemant, James, Lancelot F., Journal of the American Statistical Association


1. INTRODUCTION

Aalen (1975, 1978) developed a unified theory for nonparametric inference in multiplicative intensity models from a frequentist perspective. This treatment included, for example, the life-testing model, the multiple decrement model, birth and death processes, and branching processes. A Bayesian treatment for the real line was given by Lo and Weng (1989), who modeled hazard rates in the multiplicative intensity model as mixtures of a known kernel k with a finite measure [mu] modeled as a weighted gamma measure on the real line. That is, a hazard r is modeled as

r(x|[mu]) = [[integral].sub.R] k(x, v) [mu](dv). (1)

Lo and Weng (1989) showed that using arbitrary kernels provides the user with a great deal of flexibility; for instance, the choice of the kernel k(x, v) = I{v [less than or equal to] x} gives monotoneincreasing hazards as considered by Dykstra and Laud (1981), whereas kernels k(x, v) = I{|x - a| [greater than or equal to] v} and k(x, v) = I{|x - a| [less than or equal to] v} give U-shaped hazards (with minimum and maximum at a) similar to those of Glaser (1980), and normal density kernels k(x, v) = exp(-.5(x - v)[.sup.2]/[[tau].sup.2])/square root of (2[pi][[tau].sup.2]) can be used to estimate hazards without shape restriction.

In this article we develop a general approach to Bayesian inference for hazard (intensity) rates in nonparametric and semiparametric multiplicative intensity models by incorporating kernel mixtures of spatial weighted gamma process priors. This approach extends the work of Lo and Weng (1989) and Dykstra and Laud (1981) from a nonparametric setting on the real line to the nonparametric and semiparametric settings over general spaces and applies to the nonparametric multiplicative intensity models considered by Aalen (1975, 1978) and their semiparametric extensions developed by Andersen, Borgan, Gill, and Keiding (1993, chap. III). Models that fall within this framework that have been considered using gamma and weighted gamma processes from a Bayesian perspective include Markov models used in survival analysis subject to certain types of censoring, filtering, and truncation (Arjas and Gasbarra 1994; Laud, Smith, and Damien 1996; Gasbarra and Karia 2000), as well as Poisson point process models used in reliability (Kuo and Ghosh 1997), forest ecology (Wolpert and Ickstadt 1998a), and health exposure analysis (Best, Ickstadt, and Wolpert 2000). Another related application was given by Ibrahim, Chen, and MacEachern (1999) who used a weighted gamma process to select variables in Cox proportional hazards models.

A major contribution of this article is to develop a unified computational treatment of these problems from a Bayesian perspective. As was shown by Lo and Weng (1989) (see also Lo, Brunner, and Chan 1996) the posterior for multiplicative intensity models under weighted gamma processes share common structural features with posterior distributions for models subject to the Dirichlet process (i.e., Dirichlet process mixture models). (See Lo 1984 for background and posterior descriptions of Dirichlet process mixture models.) Recently, James (2003) extended these results to an abstract semiparametric setting (see Sec. 3), thus providing explicit calculus for relating posteriors for spatial semiparametric intensity models to posteriors for Dirichlet process mixture models. This equivalence, in combination with our use of a weighted gamma process approximation (Sec. 3), allows us to use efficient Dirichlet process computational procedures and to approximate the laws for general functionals of interest. An important aspect is that we avoid ad hoc methods used to approximate likelihoods. For example, in survival analysis problems we do not discretize time, as is often done to simplify computations. Another nice benefit of using Dirichlet process methods is that they are well understood and have a rich literature that can be drawn on. …

The rest of this article is only available to active members of Questia

Already a member? Log in now.

Notes for this article

Add a new note
If you are trying to select text to create highlights or citations, remember that you must now click or tap on the first word, and then click or tap on the last word.
One moment ...
Default project is now your active project.
Project items

Items saved from this article

This article has been saved
Highlights (0)
Some of your highlights are legacy items.

Highlights saved before July 30, 2012 will not be displayed on their respective source pages.

You can easily re-create the highlights by opening the book page or article, selecting the text, and clicking “Highlight.”

Citations (0)
Some of your citations are legacy items.

Any citation created before July 30, 2012 will labeled as a “Cited page.” New citations will be saved as cited passages, pages or articles.

We also added the ability to view new citations from your projects or the book or article where you created them.

Notes (0)
Bookmarks (0)

You have no saved items from this article

Project items include:
  • Saved book/article
  • Highlights
  • Quotes/citations
  • Notes
  • Bookmarks
Notes
Cite this article

Cited article

Style
Citations are available only to our active members.
Buy instant access to cite pages or passages in MLA, APA and Chicago citation styles.

(Einhorn, 1992, p. 25)

(Einhorn 25)

1. Lois J. Einhorn, Abraham Lincoln, the Orator: Penetrating the Lincoln Legend (Westport, CT: Greenwood Press, 1992), 25, http://www.questia.com/read/27419298.

Cited article

Computational Methods for Multiplicative Intensity Models Using Weighted Gamma Processes: Proportional Hazards, Marked Point Processes, and Panel Count Data
Settings

Settings

Typeface
Text size Smaller Larger Reset View mode
Search within

Search within this article

Look up

Look up a word

  • Dictionary
  • Thesaurus
Please submit a word or phrase above.
Print this page

Print this page

Why can't I print more than one page at a time?

Help
Full screen

matching results for page

    Questia reader help

    How to highlight and cite specific passages

    1. Click or tap the first word you want to select.
    2. Click or tap the last word you want to select, and you’ll see everything in between get selected.
    3. You’ll then get a menu of options like creating a highlight or a citation from that passage of text.

    OK, got it!

    Cited passage

    Style
    Citations are available only to our active members.
    Buy instant access to cite pages or passages in MLA, APA and Chicago citation styles.

    "Portraying himself as an honest, ordinary person helped Lincoln identify with his audiences." (Einhorn, 1992, p. 25).

    "Portraying himself as an honest, ordinary person helped Lincoln identify with his audiences." (Einhorn 25)

    "Portraying himself as an honest, ordinary person helped Lincoln identify with his audiences."1

    1. Lois J. Einhorn, Abraham Lincoln, the Orator: Penetrating the Lincoln Legend (Westport, CT: Greenwood Press, 1992), 25, http://www.questia.com/read/27419298.

    Cited passage

    Thanks for trying Questia!

    Please continue trying out our research tools, but please note, full functionality is available only to our active members.

    Your work will be lost once you leave this Web page.

    Buy instant access to save your work.

    Already a member? Log in now.

    Author Advanced search

    Oops!

    An unknown error has occurred. Please click the button below to reload the page. If the problem persists, please try again in a little while.