Parametric frailty models for clustered data with arbitrary censoring: application to effect of male circumcision on HPV clearance.
BACKGROUND
In epidemiological studies, subjects are often followed for a period during which study outcomes are measured at selected time points, such as by diagnostic testing performed on biological samples collected at each visit. Although test results may indicate the presence or absence of a disease or condition, they cannot provide information on when exactly it occurred. Such study designs generate arbitrarily censored time-to-event data, which can include left, interval and right censoring. Adding to this complexity, the data may be clustered such that observations within the same cluster are not independent, such as time to recovery of an infectious disease of family or community members. This data structure is observed when evaluating circumcision's effect on clearance of penile high risk human papillomavirus (HR-HPV) infections using data collected from the male circumcision(MC) trial conducted in Rakai, Uganda, where the multiple infections within individual and HPV testings performed at trial follow-up visits gave rise to the clustered data with arbitrary censoring.


METHODS
We describe the use of parametric proportional hazards frailty models and accelerated failure time frailty models to examine the relationship between explanatory variables and the survival outcomes that are subject to arbitrary censoring, while accounting for the correlation within clusters. Standard software such as SAS can be used for parameter estimation.


RESULTS
Circumcision's effect on HPV infection was a secondary end point in the Rakai MC trial, and HPV genotyping was conducted for penile samples of a subset of trial participants collected at enrollment, 6, 12 and 24-month follow up visits. At enrollment, 36.7% intervention arm men (immediate circumcision) and 36.6% control arm men (delayed circumcision at 2 years) were infected with HR-HPV, with the number of infections per man being 1-5. The proposed models were used to examine whether MC facilitated clearance of the prevalent infections. Results show that clearance of multiple infections within each man is highly correlated, and clearance was 60% faster if a man was circumcised.


CONCLUSIONS
Parametric frailty models provide viable ways to study the relationship between exposure variables and clustered survival outcome that is subject to arbitrary censoring, as is often observed in HPV epidemiology studies.
Kong
Xiangrong
X
Department of Population, Family and Reproductive Health, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD, USA. xikong@jhsph.edu
Archer
Kellie J
KJ
Moulton
Lawrence H
LH
Gray
Ronald H
RH
Wang
Mei-Cheng
MC
D002944:Circumcision, Male; D016000:Cluster Analysis; D005838:Genotype; D006801:Humans; D008297:Male; D015233:Models, Statistical; D027383:Papillomaviridae / Q000235:genetics / Q000502:physiology*; D030361:Papillomavirus Infections / Q000517:prevention & control / Q000821:virology*; D016016:Proportional Hazards Models; D014454:Uganda
Int J Cancer. 2006 Oct 1;119(7):1623-9
16646070
Biometrics. 2000 Dec;56(4):1016-22
11129456
Stat Med. 2008 Jul 20;27(16):3105-24
18038445
Stat Med. 2004 Oct 30;23(20):3177-92
15449337
N Engl J Med. 2009 Mar 26;360(13):1298-309
19321868
J Infect Dis. 2008 Sep 15;198(6):827-35
18657037
BMC Med Res Methodol. 2009 Feb 10;9:9
19208234
Lancet. 2007 Feb 24;369(9562):657-66
17321311
Stat Med. 1994 Nov 15;13(21):2233-47
7846422
Stat Med. 2000 Dec 30;19(24):3309-24
11122497
J Infect Dis. 2010 May 15;201(10):1463-9
20370481
Stat Med. 2008 Jun 30;27(14):2665-83
17910008