Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

R: st: Population attributable fractions (PAFs) in discrete-time survival analysis. -punaf-

From   "Carlo Lazzaro" <>
To   <>
Subject   R: st: Population attributable fractions (PAFs) in discrete-time survival analysis. -punaf-
Date   Mon, 1 Jul 2013 13:21:58 +0200

I suppose that Angelo refers to the following reference (access to the full
text conditional on subscription to Stat Med):

Samuelsen SO, Eide GE. Attributable fractions with survival data. Stat Med.
2008 Apr 30;27(9):1447-67.

Kind regards,
-----Messaggio originale-----
[] Per conto di Roger B. Newson
Inviato: lunedì 1 luglio 2013 12:57
Oggetto: Re: st: Population attributable fractions (PAFs) in discrete-time
survival analysis. -punaf-

Yes, you can use -punaf- after a generalized linear model (GLM) with a
complementary log-log link and a binomial error function. Or after any other
GLM that gives positive-valued conditional expectations (which includes
proportions and also Gamma and inverse-Gaussian means).

For proportional-hazard models (and also for case-control data), there is a
package -punafcc-, which you can also download from SSC, and which estimates
population attributable hazard factions (after proportional-hazard
regressions), or population attributable fractions (after logit regressions
on case-control data).

Angelo has not given the Samuelsen & Eide (2008) reference on PAHFs in full.
However, I would guess that the PAHFs of that reference would be either the
same as, or similar to, those produced by -punafcc-. I would very much like
to know the full reference, so I can read it and find out more.

I hope this helps.

Best wishes


Roger B Newson BSc MSc DPhil
Lecturer in Medical Statistics
Respiratory Epidemiology and Public Health Group National Heart and Lung
Institute Imperial College London Royal Brompton Campus Room 33, Emmanuel
Kaye Building 1B Manresa Road London SW3 6LR UNITED KINGDOM
Tel: +44 (0)20 7352 8121 ext 3381
Fax: +44 (0)20 7351 8322
Web page:
Departmental Web page:

Opinions expressed are those of the author, not of the institution.

On 01/07/2013 00:13, Angelo Belardi wrote:
> Dear All,
> I am working on discrete-time proportional hazard models with a 
> non-parametric baseline hazard function, using -cloglog- in 
> person-period formatted datasets.
> I would like to additionally calculate population attributable 
> fractions (PAFs) in these models.
> However, I have never worked with PAFs in survival analyses before and 
> therefore don't know which functions to use and how to correctly 
> interpret the results.
> Previously, I calculated PAFs in STATA with the -punaf- package from 
> Roger Newson, e.g.
> for logistic regressions.
> Can I use -punaf- here as well, just after calculating the estimates 
> over -cloglog-?
> Or is there another function/package for this situation?
> Or would it be better to calculate population attributable hazard 
> fractions (PAHFs) as described in Samuelsen & Eide (2008)?
> Thanks for any help or advice on the subject.
> Regards,
> Angelo
> Ref:
> S. O. Samuelsen, G. E. Eide, Statist. Med. 27, 1447 (2008).
> *
> *   For searches and help try:
> *
> *
> *
*   For searches and help try:

*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index