Skip to main content

Risk prediction models with incomplete data with application to prediction of estrogen receptor-positive breast cancer: prospective data from the Nurses' Health Study



A number of breast cancer risk prediction models have been developed to provide insight into a woman's individual breast cancer risk. Although circulating levels of estradiol in postmenopausal women predict subsequent breast cancer risk, whether the addition of estradiol levels adds significantly to a model's predictive power has not previously been evaluated.


Using linear regression, the authors developed an imputed estradiol score using measured estradiol levels (the outcome) and both case status and risk factor data (for example, body mass index) from a nested case-control study conducted within a large prospective cohort study and used multiple imputation methods to develop an overall risk model including both risk factor data from the main cohort and estradiol levels from the nested case-control study.


The authors evaluated the addition of imputed estradiol level to the previously published Rosner and Colditz log-incidence model for breast cancer risk prediction within the larger Nurses' Health Study cohort. The follow-up was from 1980 to 2000; during this time, 1,559 invasive estrogen receptor-positive breast cancer cases were confirmed. The addition of imputed estradiol levels significantly improved risk prediction; the age-specific concordance statistic increased from 0.635 ± 0.007 to 0.645 ± 0.007 (P < 0.001) after the addition of imputed estradiol.


Circulating estradiol levels in postmenopausal women appear to add to other lifestyle factors in predicting a woman's individual risk of breast cancer.


Breast cancer risk prediction models have been developed for use as an entry criterion into breast cancer chemoprevention trials (for example, National Surgical Adjuvant Breast and Bowel Project tamoxifen trial and the Study of Tamoxifen and Raloxifene), in counseling women on the potential use of chemopreventives, and to provide insight into a woman's individual breast cancer risk [14]. The initial Gail model incorporated a subset of breast cancer risk factors, namely age, age at menarche, age at first birth, family history of breast cancer or of atypical hyperplasia, and history of breast biopsies [5, 6]. Subsequently, several groups have developed more extensive statistical models that incorporate a greater number of breast cancer risk factors [1, 4].

In postmenopausal women, circulating levels of estradiol predict subsequent breast cancer risk [710], particularly for estrogen receptor (ER)-positive disease [9]. However, plasma estradiol levels have not previously been evaluated within risk prediction models, and whether their addition would add to the model's predictive power is unknown. Plasma estradiol is available only from the Nurses Health Study nested case-control data set. We initially attempted to evaluate the addition of estradiol concentrations to the Rosner and Colditz risk prediction model using data from the Nurses' Health Study nested case-control data set. However, with the relatively modest size of the nested case-control data set and the large number of parameters to be estimated, a number of the risk factor parameters did not adequately reflect those from the parent cohort. Thus, development of an accurate risk prediction tool requires a large sample size as in the main Nurses' Health Study cohort. The purpose of this paper is to describe a methodology for developing a risk prediction rule when one or more predictors are incompletely observed and to apply it to assess the predictive power of plasma estradiol after adjusting for standard breast cancer risk factors as well as the effect of other risk factors adjusted for plasma estradiol.

Materials and methods


The Nurses' Health Study cohort was established in 1976, when 121,700 female US nurses (30 to 55 years old) responded to a mailed questionnaire that inquired about reproductive history and a range of lifestyle factors in addition to disease diagnoses [11, 12]. Follow-up questionnaires have been mailed biennially to update exposure information and any major medical events. Deaths are reported by family members or the postal service or are identified by a search of the National Death Index. We estimate that mortality ascertainment is 98% complete [13, 14]. This investigation was approved by the Brigham and Women's Hospital institutional review board.

Identification of breast cancer cases

On each questionnaire, we inquired whether breast cancer had been diagnosed and, if so, the date of diagnosis. All women who reported breast cancer (or the next of kin for decedents) were contacted for permission to review medical records to confirm the diagnosis. We include only invasive cases of breast cancer confirmed by the pathology report. ER and progesterone receptor status of the tumor was determined from the medical record. In this report, we evaluate ER+ cases only, as we previously observed that plasma estradiol most strongly predicted this tumor subtype [9].

Population for analysis

The population of women used in this analysis has been described in detail in several previous publications [1, 15]. Briefly, we excluded women with unknown, inconsistent, or out-of-range reports for height, weight in 1976 or at age 18, age at menarche or menopause or each pregnancy, parity, and duration or type of postmenopausal hormone (PMH) use (n = 42,886). Additionally, women with a simple hysterectomy (and hence unknown age at menopause) (n = 10,301) were excluded. Participants who were ineligible for the study (for example, prevalent cancer in 1976) or no follow-up after 1978 (n = 2,360) were excluded. In the current analysis, women who were premenopausal throughout follow-up were excluded (n = 6,342), but once they became postmenopausal, they could contribute person-time. Overall, 59,812 participants remained for this analysis. These women contributed 750,086 person-years from 1980 to 2000, during which 1,559 incident invasive ER+ breast cancer cases occurred.

Blood subcohort and nested case-control study

From 1989 to 1990, 32,826 cohort members provided blood samples. Informed consent was obtained from each participant; details about the blood collection methods have been published previously [16, 17]. Briefly, women arranged to have their blood drawn and shipped with an icepack via overnight courier to our laboratory, where it was processed and archived in liquid nitrogen freezers. Estradiol is stable in cooled whole blood for 24 to 48 hours [18]. At blood collection, women completed a short questionnaire that included questions on recent use of PMH (within the last 3 months). Follow-up of the blood study cohort was 99% in 2000.

In the current analyses, we used a previously described nested case-control study of sex steroids and breast cancer risk with cases diagnosed after blood collection through 31 May 1998 [9, 16]. In addition, cases diagnosed up through 31 May 2000 and their matched controls (that is, a 2-year extension of the published report [9]) are included. At blood collection, cases and controls were postmenopausal, were not recent users of PMH, and had no prior diagnosed cancer (except nonmelanoma skin cancer). Control subjects were matched by age, month/year and time of day of blood collection, and fasting status and had not been diagnosed with breast cancer before the diagnosis date of their matched case. To mimic the larger population used in the risk prediction modeling, only cases and controls meeting the inclusion criteria described above were included (for example, no prior simple hysterectomy). Women were considered postmenopausal if they reported having a natural menopause (for example, no menstrual cycles during the previous 12 months) or had a bilateral oophorectomy. In all, 164ER+ cases and 346 controls were included.

Laboratory assays

Estradiol was measured by radioimmunoassay following extraction and celite column chromatography, as previously described [9]. The coefficient of variation was less than or equal to 11%.

Description of the risk prediction model

We fit the log-incidence model of breast cancer to incident ER+cases, as previously described [1, 15]. We assume that incidence at time, t(I t ), is proportional to the number of cell divisions, C t , accumulated throughout life up to age t; that is,

I t = kC t

The cumulative number of breast cell divisions is factored as follows:

C t = C o x i = 0 t 1 ( C i + 1 / C i ) = C o x i = 0 t 1 λ i MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaee4qam0aaSbaaSqaaiabbsha0bqabaGccqGH9aqpcqqGdbWqdaWgaaWcbaGaee4Ba8gabeaakiabbIha4naarahabaGaeiikaGIaee4qam0aaSbaaSqaaiabbMgaPjabgUcaRiabigdaXaqabaGccqGGVaWlcqqGdbWqdaWgaaWcbaGaeeyAaKgabeaakiabcMcaPaWcbaGaeeyAaKMaeyypa0JaeGimaadabaGaeeiDaqNaeyOeI0IaeGymaedaniabg+GivdGccqGH9aqpcqqGdbWqdaWgaaWcbaGaee4Ba8gabeaakiabbIha4naarahabaGaeq4UdW2aaSbaaSqaaiabbMgaPbqabaaabaGaeeyAaKMaeyypa0JaeGimaadabaGaeeiDaqNaeyOeI0IaeGymaedaniabg+Givdaaaa@5983@

Thus, λ i = C i+1/C i represents the rate of increase of breast cell divisions from age i to age i+1. Log (λ i ) is assumed to be a linear function of risk factors that are relevant at age i. The set of risk factors and their magnitude may vary according to the stage of reproductive life. Details of the representation of the C i are given in [1, 15]. The overall model is given by:

log I = α + β 0 ( t m t 0 ) + β 1 b + β 2 ( t 1 t 0 ) b 1 , t 1 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGceaqabeaacyGGSbaBcqGGVbWBcqGGNbWzcqWGjbqscqGH9aqpcqaHXoqycqGHRaWkcqaHYoGydaWgaaWcbaGaeGimaadabeaakiabcIcaOiabdsha0naaBaaaleaacqWGTbqBaeqaaOGaeyOeI0IaemiDaq3aaSbaaSqaaiabicdaWaqabaGccqGGPaqkaeaacqGHRaWkcqaHYoGydaWgaaWcbaGaeGymaedabeaakiabdkgaIjabgUcaRiabek7aInaaBaaaleaacqaIYaGmaeqaaOGaeiikaGIaemiDaq3aaSbaaSqaaiabigdaXaqabaGccqGHsislcqWG0baDdaWgaaWcbaGaeGimaadabeaakiabcMcaPiabdkgaInaaBaaaleaacqaIXaqmcqGGSaalcqWG0baDcqGHsislcqaIXaqmaeqaaaaaaa@58E6@
+ γ 1 ( t t m ) m A + γ 2 ( t t ) m m B + δ 1 p m h A + δ 2 p m h B + δ 3 p m h C + δ 4 p m h c u r , t + ( δ 4 + δ 5 ) p m h p a s t , t + β 3 B M I 1 + β 3 * B M I 2 + β 4 h 1 + β 4 * h 2 + α 1 b b d + α 2 b b d t o + α 3 b b d ( t m t o ) + α 4 b b d ( t t m ) m t + φ f h x + β 5 a l c 1 + β 5 a l c 2 + β 5 a l c 3 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGceaqabeaacqGHRaWkcqaHZoWzdaWgaaWcbaGaeGymaedabeaakiabcIcaOiabdsha0jabgkHiTiabdsha0naaBaaaleaacqWGTbqBaeqaaOGaeiykaKIaemyBa02aaSbaaSqaaiabdgeabbqabaGccqGHRaWkcqaHZoWzdaWgaaWcbaGaeGOmaidabeaakiabcIcaOiabdsha0jabgkHiTiabdsha0naaBeaaleaacqWGTbqBaeqaaOGaeiykaKIaemyBa02aaSbaaSqaaiabdkeacbqabaaakeaacqGHRaWkcqaH0oazdaWgaaWcbaGaeGymaedabeaakiabbccaGiabdchaWjabd2gaTjabdIgaOnaaBaaaleaacqWGbbqqaeqaaOGaey4kaSIaeqiTdq2aaSbaaSqaaiabikdaYaqabaGccqqGGaaicqWGWbaCcqWGTbqBcqWGObaAdaWgaaWcbaGaemOqaieabeaakiabgUcaRiabes7aKnaaBaaaleaacqaIZaWmaeqaaOGaeeiiaaIaemiCaaNaemyBa0MaemiAaG2aaSbaaSqaaiabdoeadbqabaaakeaacqGHRaWkcqaH0oazdaWgaaWcbaGaeGinaqdabeaakiabbccaGiabdchaWjabd2gaTjabdIgaOnaaBaaaleaacqWGJbWycqWG1bqDcqWGYbGCcqGGSaalcqWG0baDaeqaaOGaey4kaSIaeiikaGIaeqiTdq2aaSbaaSqaaiabisda0aqabaGccqGHRaWkcqaH0oazdaWgaaWcbaGaeGynaudabeaakiabcMcaPiabbccaGiabdchaWjabd2gaTjabdIgaOnaaBaaaleaacqWGWbaCcqWGHbqycqWGZbWCcqWG0baDcqGGSaalcqWG0baDaeqaaOGaey4kaSIaeqOSdi2aaSbaaSqaaiabiodaZaqabaGccqWGcbGqcqWGnbqtcqWGjbqsdaWgaaWcbaGaeGymaedabeaakiabgUcaRiabek7aInaaDaaaleaacqaIZaWmaeaacqGGQaGkaaGccqWGcbGqcqWGnbqtcqWGjbqsdaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabek7aInaaBaaaleaacqaI0aanaeqaaOGaemiAaG2aaSbaaSqaaiabigdaXaqabaGccqGHRaWkcqaHYoGydaqhaaWcbaGaeGinaqdabaGaeiOkaOcaaOGaemiAaG2aaSbaaSqaaiabikdaYaqabaaakeaacqGHRaWkcqaHXoqydaWgaaWcbaGaeGymaedabeaakiabbccaGiabdkgaIjabdkgaIjabdsgaKjabgUcaRiabeg7aHnaaBaaaleaacqaIYaGmaeqaaOGaeeiiaaIaemOyaiMaemOyaiMaemizaqMaeeiiaaIaemiDaq3aaSbaaSqaaiabd+gaVbqabaGccqGHRaWkcqaHXoqydaWgaaWcbaGaeG4mamdabeaakiabbccaGiabdkgaIjabdkgaIjabdsgaKjabbccaGiabcIcaOiabdsha0naaBaaaleaacqWGTbqBaeqaaOGaeyOeI0IaemiDaq3aaSbaaSqaaiabd+gaVbqabaGccqGGPaqkcqGHRaWkcqaHXoqydaWgaaWcbaGaeGinaqdabeaakiabbccaGiabdkgaIjabdkgaIjabdsgaKjabbccaGiabcIcaOiabdsha0jabgkHiTiabdsha0naaBaaaleaacqWGTbqBaeqaaOGaeiykaKIaemyBa02aaSbaaSqaaiabdsha0bqabaaakeaacqGHRaWkcqaHgpGAcqqGGaaicqWGMbGzcqWGObaAcqWG4baEcqGHRaWkcqaHYoGydaWgaaWcbaGaeGynaudabeaakiabbccaGiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqaIXaqmaeqaaOGaey4kaSIaeqOSdi2aa0baaSqaaiabiwda1aqaaiabgEHiQaaakiabbccaGiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqaIYaGmaeqaaOGaey4kaSIaeqOSdi2aa0baaSqaaiabiwda1aqaaiabgEHiQiabgEHiQaaakiabbccaGiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqaIZaWmaeqaaaaaaa@0B39@

t = age

t o = age at menarche

t m = age at menopause

s t = parity at age t

t i = age at ith birth, i = 1, ..., s t

b = birth index = i = 1 s ( t m t i ) b i t MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaWaaabCaeaacqGGOaakcqWG0baDdaWgaaWcbaGaemyBa0gabeaakiabgkHiTiabdsha0naaBaaaleaacqWGPbqAaeqaaaqaaiabdMgaPjabg2da9iabigdaXaqaaiabdohaZXGaemiDaqhaniabggHiLdGccqGGPaqkcqWGIbGydaWgaaWcbaGaemyAaKMaemiDaqhabeaaaaa@4334@ for parous women, = 0 for nulliparous women

b it = 1 if parity is greater than or equal to i at age t (otherwise, b it = 0)

m A = 1 if natural menopause (otherwise, m A = 0)

m B = 1 if bilateral oophorectomy (otherwise,m B = 0)

bbd = 1 if breast disease is benign (otherwise, bbd = 0)

fhx = 1 if there is a family history of breast cancer in mother or sister (otherwise, fhx = 0)

pmh A = number of years on oral estrogen

pmh B = number of years on oral estrogen and progestin

pmh C = number of years on other types of PMHs

pmh cur,t = 1 if current user of postmenopausal hormones at age t (otherwise, pmh cur,t = 0)

pmh past,t = 1 if past user of postmenopausal hormones at age t (otherwise, pmh past,t = 0)

BMI j = body mass index (BMI) at age j (kg/m2)

alc j = alcohol consumption (grams) at age j

h = height (inches)

B M I 1 = j = t o t m 1 ( B M I j 21.8 ) + j = t m t 1 ( B M I j 24.4 ) p m h c u r , j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOqaiKaemyta0KaemysaK0aaSbaaSqaaiabigdaXaqabaGccqGH9aqpdaaeWbqaaiabcIcaOiabdkeacjabd2eanjabdMeajnaaBaaaleaacqWGQbGAaeqaaOGaeyOeI0IaeGOmaiJaeGymaeJaeiOla4IaeGioaGJaeiykaKcaleaacqWGQbGAcqGH9aqpcqWG0baDdaWgaaadbaGaem4Ba8gabeaaaSqaaiabdsha0naaBaaameaacqWGTbqBaeqaaSGaeyOeI0IaeGymaedaniabggHiLdGccqGHRaWkdaaeWbqaaiabcIcaOiabdkeacjabd2eanjabdMeajnaaBaaaleaacqWGQbGAaeqaaOGaeyOeI0IaeGOmaiJaeGinaqJaeiOla4IaeGinaqJaeiykaKIaemiCaaNaemyBa0MaemiAaG2aaSbaaSqaaiabdogaJjabdwha1jabdkhaYjabcYcaSiabdQgaQbqabaaabaGaemOAaOMaeyypa0JaemiDaq3aaSbaaWqaaiabd2gaTbqabaaaleaacqWG0baDcqGHsislcqaIXaqma0GaeyyeIuoaaaa@6D4B@
B M I 2 = j = t m t 1 ( B M I j 24.4 ) ( 1 p m h c u r . j ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOqaiKaemyta0KaemysaK0aaSbaaSqaaiabikdaYaqabaGccqGH9aqpdaaeWbqaaiabcIcaOiabdkeacjabd2eanjabdMeajnaaBaaaleaacqWGQbGAaeqaaOGaeyOeI0IaeGOmaiJaeGinaqJaeiOla4IaeGinaqJaeiykaKIaeeiiaaIaeiikaGIaeGymaeJaeyOeI0IaemiCaaNaemyBa0MaemiAaG2aaSbaaSqaaiabdogaJjabdwha1jabdkhaYjabc6caUiabdQgaQbqabaGccqGGPaqkaSqaaiabdQgaQjabg2da9iabdsha0naaBaaameaacqWGTbqBaeqaaaWcbaGaemiDaqNaeyOeI0IaeGymaedaniabggHiLdaaaa@58F1@
h 1 = ( h 64.5 ) ( t m t o ) + ( h 64.4 ) j = t m t 1 p m h c u r , j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiAaG2aaSbaaSqaaiabigdaXaqabaGccqGH9aqpcqGGOaakcqWGObaAcqGHsislcqaI2aGncqaI0aancqGGUaGlcqaI1aqncqGGPaqkcqGGOaakcqWG0baDdaWgaaWcbaGaemyBa0gabeaakiabgkHiTiabdsha0naaBaaaleaacqWGVbWBaeqaaOGaeiykaKIaey4kaSIaeiikaGIaemiAaGMaeyOeI0IaeGOnayJaeGinaqJaeiOla4IaeGinaqJaeiykaKYaaabCaeaacqWGWbaCcqWGTbqBcqWGObaAdaWgaaWcbaGaem4yamMaemyDauNaemOCaiNaeiilaWIaemOAaOgabeaaaeaacqWGQbGAcqGH9aqpcqWG0baDdaWgaaadbaGaemyBa0gabeaaaSqaaiabdsha0jabgkHiTiabigdaXaqdcqGHris5aaaa@606E@
h 2 = ( h 64.4 ) j = t m t 1 ( 1 p m h c u r , j ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiAaG2aaSbaaSqaaiabikdaYaqabaGccqGH9aqpcqGGOaakcqWGObaAcqGHsislcqaI2aGncqaI0aancqGGUaGlcqaI0aancqGGPaqkdaaeWbqaaiabcIcaOiabigdaXiabgkHiTiabdchaWjabd2gaTjabdIgaOnaaBaaaleaacqWGJbWycqWG1bqDcqWGYbGCcqGGSaalcqWGQbGAaeqaaOGaeiykaKcaleaacqWGQbGAcqGH9aqpcqWG0baDdaWgaaadbaGaemyBa0gabeaaaSqaaiabdsha0jabgkHiTiabigdaXaqdcqGHris5aaaa@52B7@
a l c 1 = j = 18 t m 1 a l c j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaemiBaWMaem4yam2aaSbaaSqaaiabigdaXaqabaGccqGH9aqpdaaeWbqaaiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqWGQbGAaeqaaaqaaiabdQgaQjabg2da9iabigdaXiabiIda4aqaaiabdsha0naaBaaameaacqWGTbqBaeqaaSGaeyOeI0IaeGymaedaniabggHiLdaaaa@4535@
a l c 2 = j = t m t 1 a l c j p m h c u r , j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaemiBaWMaem4yam2aaSbaaSqaaiabikdaYaqabaGccqGH9aqpdaaeWbqaaiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqWGQbGAaeqaaOGaemiCaaNaemyBa0MaemiAaG2aaSbaaSqaaiabdogaJjabdwha1jabdkhaYjabcYcaSiabdQgaQbqabaaabaGaemOAaOMaeyypa0JaemiDaq3aaSbaaWqaaiabd2gaTbqabaaaleaacqWG0baDcqGHsislcqaIXaqma0GaeyyeIuoaaaa@4F81@
a l c 3 = j = t m t 1 a l c j ( 1 p m h c u r , j ) . MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyyaeMaemiBaWMaem4yam2aaSbaaSqaaiabiodaZaqabaGccqGH9aqpdaaeWbqaaiabdggaHjabdYgaSjabdogaJnaaBaaaleaacqWGQbGAaeqaaOGaeiikaGIaeGymaeJaeyOeI0IaemiCaaNaemyBa0MaemiAaG2aaSbaaSqaaiabdogaJjabdwha1jabdkhaYjabcYcaSiabdQgaQbqabaGccqGGPaqkaSqaaiabdQgaQjabg2da9iabdsha0naaBaaameaacqWGTbqBaeqaaaWcbaGaemiDaqNaeyOeI0IaeGymaedaniabggHiLdGccqGGUaGlaaa@5415@

The general rationale for a log-incidence model is that the number of precancerous cells increases multiplicatively with time but that historical exposures differentially affect the rate of increase. Specifically, for breast cancer, the number of precancerous cells is assumed to increase annually at the rate of exp(β 0) prior to menopause for nulliparous women, at the rate of exp(β 0 + β 1s) prior to menopause for parous women with parity = s, and so forth. Finally, the number of precancerous cells increases immediately after the first birth by exp [β 2(t1 - t0)]. The incidence rate of breast cancer is assumed to be approximately proportional to the number of precancerous cells.

The log-incidence model was fit using iteratively reweighted least squares with PROC NLIN in SAS (SAS version 6.12; SAS Institute Inc., Cary, NC, USA) (1996). The parameters of the model are readily interpretable in a relative risk (RR) context. For example, exp(-β 0) = RR for a 1-year increase in age at menarche among nulliparous women, exp [-(β 0 + β 2)] = RR for a 1-year increase in age at menarche among parous women, and so forth. In this analysis, women were followed until they had an event (ER+ breast cancer) or were censored if they developed (a) ER- breast cancer, (b) breast cancer in which ER status is unknown, or (c) other types of cancer except nonmelanoma skin cancer or (d) if they died.

Imputation and inclusion of estradiol in the risk prediction model

Ideally, we would have estradiol levels measured on each main study participant at several points in time. However, since this was not possible, we used an indirect approach to impute estradiol. Let x = estradiol and z = other covariates in the risk prediction model.

From the main study, we can obtain Pr(D|z) given under the rare disease assumption by:

Pr(D| z ¯ ) = exp( α + β ¯ z ¯ )/[1 + exp( α + β ¯ z ¯ )] exp ( α + β ˜ z ˜ ) . MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeeiuaaLaeeOCaiNaeeikaGIaeeiraqKaeeiFaW3aaWaaaeaacqqG6bGEaaGaeeykaKcccaGae8xpa0JaeeyzauMaeeiEaGNaeeiCaaNaeeikaGIaeqySdeMae83kaSYaaWaaaeaacqaHYoGyaaWaaWaaaeaacqqG6bGEaaGaeeykaKIaee4la8Iaee4waSLaeeymaeJae83kaSIaeeyzauMaeeiEaGNaeeiCaaNaeeikaGIaeqySdeMaey4kaSYaaWaaaeaacqaHYoGyaaWaaWaaaeaacqqG6bGEaaGaeeykaKIaeeyxa0LaeyyrIaKagiyzauMaeiiEaGNaeiiCaaNaeiikaGIaeqySdeMaey4kaSIafqOSdiMba4bacuWG6bGEgaGhaiabcMcaPiabc6caUaaa@61BB@

We want to estimate Pr(D|x,z), where under the rare disease assumption

Pr(D|x,z) exp(α* + β*z + δ*x).

From the blood study, we can estimate δ* based on conditional logistic regression. Indeed, in principle, we could also estimate β* from the blood study, but the estimates will be very imprecise due to the small sample size. Therefore, we used the main study population to estimate the parameters in Equation 5 by estimating x for all subjects in the main study based on a linear regression derived from the blood study:

x = α o + γ o y + γ z imp + e,

where x = ℓn (estradiol) as a continuous variable, y = 1 if case and 0 if control, and Z imp = a subset of the other covariates Z in the risk prediction model. Z imp was ascertained by first forcing in y and then using stepwise-up regression to determine the subset of components of Z in the main study which were significantly associated with x at the 5% level.

In the blood study, estradiol levels on average were higher for cases than controls. The rationale for including y as a covariate in Equation 6 is to account for this relationship in the main study as well. In addition, because there is substantial overlap between the estradiol distribution of cases and controls, we used an imputation strategy to estimate x by adding error to the prediction such that for each main study participant we obtain

x ^ i = α o + γ o y i + γ ¯ z ¯ i m p , i + e i , MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiEaGNbaKaalmaaBaaameaacqWGPbqAaeqaaOGaeyypa0JaeqySde2aaSbaaSqaaiabd+gaVbqabaGccqGHRaWkcqaHZoWzdaWgaaWcbaGaem4Ba8gabeaakiabdMha5TWaaSbaaWqaaiabdMgaPbqabaGccqGHRaWkdaadaaqaaiabeo7aNbaadaqfqaqabSqaaiabdMgaPjabd2gaTjabdchaWjabcYcaSiabdMgaPbqab4qaamaamaaabaGaemOEaOhaaaaakiabgUcaRiabdwgaLTWaaSbaaWqaaiabdMgaPbqabaWccqGGSaalaaa@4C4B@

where (a) e i = an error term that is normally distributed with mean 0 and variance σ2, (b) y i = 1 if a breast cancer case and = 0 otherwise, (c) σ2 is estimated from Equation 6, and e is obtained by the RANNOR function of SAS so as to add error to the estimate of x for individual women. We then fit the model in Equation 5 using x ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiEaGNbaKaaaaa@2F59@ instead of x, thus obtaining the model

Pr ( D | x , z ) exp ( α * * + β * * z + δ * * x ^ ) . MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xI8qiVKIOFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGagiiuaaLaeiOCaiNaeiikaGIaemiraqKaeiiFaWNaemiEaGNaeiilaWIaemOEaONaeiykaKIaeyyrIaKagiyzauMaeiiEaGNaeiiCaaNaeiikaGIaeqySdeMaeiOkaOIaeiOkaOIaey4kaSIaeqOSdiMaeiOkaOIaeiOkaOIaemOEaONaey4kaSIaeqiTdqMaeiOkaOIaeiOkaOIafmiEaGNbaKaacqGGPaqkcqGGUaGlaaa@4F9C@

Since the parameter estimates in Equation 8 may be influenced by the random error introduced in Equation 7, we repeated this imputation approach four additional times and used multiple imputation [19] to combine estimates from the separate imputations to obtain an overall estimate.

To assess the additional predictive power of serum estradiol, we computed age-specific (5-year age groups) deciles of the risk function without estradiol (model A) as well as including imputed estradiol (model B). From the cross-classification of risk decile model A × risk decile model B, we then compared the observed number of cases in specific risk deciles of model B with the expected number of cases within strata defined by model A risk decile. Specifically, let Xij = the number of breast cancer cases, Nij = the number of person-years, and pij = Xij/Nij, which is the estimated incidence rate within the ith age-specific risk decile for model A and the jth age-specific risk decile for model B, and let ln(pij) = αi + β(j - 1). 100% × [exp( β MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuqqRPxAKvMB6bYrY9gDLn3AGiuraeXatLxBI9gBaebbnrfifHhDYfgasaacPi6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafqOSdiMbambaaaa@2F8B@ )-1] is an estimate of the percentage increase in breast cancer incidence for an increase of one model B risk decile, holding the model A risk decile constant [20]. We wish to test the hypothesis H0: β = 0 versus H1: β ≠ 0. This approach of cross-classifying individuals by two different risk prediction rules is similar to the reclassification table approach used to compare risk prediction rules in the Framingham Heart Study [21]. In addition, to assess the predictive ability of our risk prediction models, we used the area under the receiver operating characteristic curve (that is, the concordance or C statistic). This statistic ranges from 0.5 to 1.0 and represents the probability that, for a randomly selected pair of women, one with ER+ breast cancer and one without breast cancer, the woman with ER+ breast cancer has the higher estimated disease probability. Also, we compared the C statistic for different risk prediction rules [22]. In our primary analysis, we evaluated the addition of imputed estradiol levels to risk prediction models in the entire cohort. As a secondary approach, we calculated Rosner and Colditz model risk scores in the entire cohort and then, in the nested case-control data set, assessed the impact of adding this score to the plasma estradiol and breast cancer model.


Within the nested case-control data set, we observed a significant association between plasma estradiol and risk of breast cancer (P trend < 0.001), with an RR for the top (Q4) versus bottom (Q1) quartile category of 3.3 (95% confidence interval [CI] = 1.8 to 6.0) for ER+ breast cancer (Table 1). Each of the variables in the Rosner and Colditz risk prediction model was considered as a potential predictor of plasma estradiol. BMI was most strongly related to estradiol level; in addition, the birth index, case status, and duration of postmenopause each contributed modestly but significantly (Table 2). Other variables, including family history of breast cancer, alcohol intake, and history of BBD, did not contribute significantly to the model and thus were dropped from further consideration. The r 2, for the regression model, was 0.219.

Table 1 Relative risk of estrogen receptor-positive breast cancer by quartile of postmenopausal plasma estradiol concentration (164 cases and 346 controls)a
Table 2 Imputation of log estradiol using data from the Nurses' Health Study breast cancer nested case-control studya

Most variables in the Rosner and Colditz model incorporate a time component (for example, postmenopausal BMI = average BMI postmenopause × duration postmenopause). Because exposure status at the time of blood draw might be most strongly correlated with estradiol, we also evaluated each of the variables at the time of blood draw or, for variables ascertained only on the main study questionnaire (for example, alcohol intake), within 2 years of blood draw. All results were similar.

In Table 3, we present the standard risk prediction model without estradiol and an enhanced model with imputed estradiol based on an average of five imputations. The regression coefficient for loge estradiol varied over the five imputations (0.400 to 0.582) with an average of 0.477. For a one-unit increase in loge estradiol, the RR was 1.61 (95% CI = 1.35 to 1.93). Its inclusion also caused changes in some of the other model parameters. As expected, the regression coefficient for postmenopausal BMI decreased from 0.0038 (P < 0.001) to 0.0023 (P = 0.002). For example, if we compare two 70-year-old postmenopausal women with age at menopause of 50 years, no PMH use, and constant BMI of 20 and 30, respectively, from age 50 to 70, then the RRs of breast cancer at age 70 for the woman with BMI = 30 versus the woman with BMI = 20 would be exp(0.0038 × 10 × 20) = 2.1 (95% CI = 1.7 to 2.7) without controlling for estradiol and exp(0.0023 × 10 × 20) = 1.6 (95% CI = 1.2 to 2.1) after controlling for estradiol. Also, the regression coefficient for the birth index decreased from -0.0030 ± 0.0007 (P < 0.001) to -0.0020 ± 0.0008 (P = 0.02). To interpret this difference, we compare two postmenopausal women with age at menarche of 13 and age at menopause of 50, one of whom had four births at ages 20, 23, 26, and 29 whereas the other was nulliparous. The birth indices for these two women are 102 [(50 - 20) + (50 - 23) + (50 - 26) + (50 - 29)] and 0, respectively. The RRs for the parous versus nulliparous woman are exp [-0.0030(102)] = 0.74 (95% CI = 0.64 to 0.85) without and exp [-0.0020(102)] = 0.82 (95% CI = 0.69 to 0.96) with controlling for loge estradiol. Finally, the regression coefficient for both the duration after natural menopause and the duration after bilateral oophorectomy increased after adjusting for loge estradiol. To interpret this finding, we compare two postmenopausal women with age at natural menopause of 45 and 55 years, respectively. The RRs of breast cancer for the second compared with the first woman are exp [10(0.102-0.049)] = 1.7 (95% CI = 1.4 to 2.0) without and exp [10(0.101-0.056)] = 1.6 (95% CI = 1.3 to 1.9) with adjusting for loge estradiol. It appears that part of the effect of BMI, parity, and late menopause on the incidence of breast cancer is mediated in part by changes in loge estradiol caused by obesity (increase), multiparity with first birth at an early age (decrease), and delayed menopause (increase), respectively.

Table 3 Risk prediction model for breast cancer with and without estradiol (loge)

We now present the cross-classification of model A × model B risk decile in Table 4. It is clear from Table 4 that, within most model A risk deciles, there are important differences in estimated incidence according to model B risk decile (often twofold). Overall, for a given model A risk decile, the observed number of cases was higher than expected when the model B decile was high and lower than expected when the model B decile was low. The overall slope was β = 0.511 ± 0.034 (P < 0.001), indicating that there is a significant estimated 67% increase in breast cancer incidence for an increase of one model B age-specific risk decile, holding the age-specific model A risk decile constant. This indicated that there is substantial increased predictive power upon adding loge estradiol to the risk prediction model.

Table 4 Cross-classification of model A risk decilea × model B risk decileb

In addition, we compared the age-specific C statistics between model A versus model B. We found C statistics of 0.635 ± 0.007 for model A and 0.645 ± 0.007 for model B (C statistic model A versus C statistic model B; P < 0.001). Using our secondary approach (applying the population-based risk scores to the subset of women in the nested case-control study), the RRs of breast cancer by plasma hormone level were similar, though slightly attenuated, compared with those in Table 1. For example, the RRs of ER+ breast cancer with increasing quartile of estradiol were 1.0, 1.5, 1.4, and 2.5 (95% CI = 1.5 to 4.2). Finally, in Table 5, we present the 5-year incidence of breast cancer by age and model B risk decile after adjusting for competing mortality risks [23]. The RR of breast cancer comparing women at the highest versus the lowest age-specific decile ranges from 5.0 to 8.5. For example, for 60- to 64-year-old women, the absolute 5-year risk of breast cancer is 436/105(0.4%) for women in the first decile and 2,982/105 (3.0%) for women in the 10th decile (RR = 6.8), indicating substantial differences in absolute risk according to the model B risk equation.

Table 5 Five-year risk of breast cancer by age and model B risk decile


In the Nurses' Health Study, we found that estradiol levels, as imputed from a nested case-control study within the same cohort, added significantly to the Rosner and Colditz risk prediction model, which already includes most confirmed breast cancer risk factors. There was an increase of 67% in incidence per increase of one model B risk decile, holding model A risk decile constant. The increase in the C statistic was also statistically significant.

Strengths of this study include the large size of the cohort and the large number of available questionnaire-based breast cancer risk factors. Additionally, prospectively assessed estradiol levels were available in a subset of the same women. Through the use of both risk factors and case status in the linear regression, our imputed values as applied to the larger cohort accounted for both the association between hormone level and breast cancer and the correlation between hormone levels and other risk factors already in the risk prediction model.

One limitation of the study was that we did not have measured estradiol levels on all cohort members; however, this is a limitation of all large prospective studies because of the high cost of the assays. In addition, due to our desire to have consistent eligibility criteria throughout, only 164 cases and 346 controls in the nested case-control study met all criteria for the model (for example, known age at menopause). Thus, it was not possible within this small data set to provide a sufficiently precise evaluation of the Rosner and Colditz model (which contains 22 beta coefficients); in our initial attempt to evaluate the model, all beta coefficients had wide CIs. In our secondary analyses, in which we used the risk score within the nested case-control data set, plasma estradiol again contributed significantly to the model. However, the RR for estradiol in the secondary analysis for Q4 versus Q1 was 2.5 (95% CI = 1.5 to 4.2) when controlling for other risk factors using the Rosner-Colditz risk scores versus 3.3 (95% CI = 1.8 to 6.0) when individual risk factors were used within the case-control study. One would expect that the effects of other risk factors are more accurately measured by a single risk score derived from a large cohort study than individual risk factors derived from a relatively small nested case-control data set. More generally, this may indicate better control for confounding in small case-control studies based on risk scores derived from large cohorts versus internal control for confounding based on individual risk factors whose regression coefficients are poorly estimated in small case-control studies. With further follow-up and within a collaboration across several cohorts at other institutions, we plan to re-evaluate the case-control approach.

In imputing estradiol levels, only BMI, the birth index, and duration of postmenopause (in addition to case status) were significant predictors of loge (estradiol). The correlation between BMI and estradiol was expected given that aromatization of androgens to estrogens in postmenopausal women occurs in adipose tissue [24]. The association with the birth index, a summary variable representing the number and spacing of pregnancies, has been evaluated less frequently and data are not as consistent [2528]. The association with duration of postmenopause may be due to declines in estrogen after menopause. Alcohol intake, which previously has been found to correlate with estrogen levels in several studies [29], did not contribute significantly here, which is consistent with our previous report from a subset of the current population showing no correlation with estradiol [17]. To our knowledge, no other lifestyle factors have consistently been shown to predict postmenopausal estradiol levels. The correlation between measured and imputed estradiol was 0.47.

With the inclusion of imputed estradiol, the C statistic increased from 0.635 to 0.645, which is a modest improvement but suggests reasonable discriminatory ability overall. However, the reclassification table approach (Table 4) indicated that a substantial difference in incidence is explained by including imputed estradiol, suggesting that the C statistic may be relatively insensitive to additions of single predictors to risk prediction models [21, 30]. However, the relationship of one or a combination of risk factors with disease must be very strong – RRs on the order of 100 to 200 between exposed and unexposed – to serve as a screening tool at the individual level [3133]. Continued expansion of current models with other risk factors (for example, genetic factors, mammographic density, or cytology from nipple aspirate fluid [34]) may further improve the C statistic. In addition, if chemopreventive agents were developed with few risks (and an acceptable cost-benefit ratio), the need to minimize the false-positive rate would decrease, similar to the use of cholesterol-lowering agents for the prevention of heart disease.


In summary, our data indicate that circulating estradiol levels in postmenopausal women may contribute significantly to current risk prediction models. Further assessment of estradiol in other studies and of other biomarkers that predict risk is needed to continue to improve our ability to predict breast cancer risk and inform prevention strategies. Similar approaches can be used to incorporate other breast cancer biomarkers in overall risk prediction models.



benign breast disease


body mass index

C statistic:

concordance statistic


confidence interval


estrogen receptor


postmenopausal hormone


Q2, Q3, and Q4 = first, second, third, and fourth quartile


relative risk.


  1. Colditz GA, Rosner B: Cumulative risk of breast cancer to age 70 years according to risk factor status: data from the Nurses' Health Study. Am J Epidemiol. 2000, 152: 950-964. 10.1093/aje/152.10.950.

    CAS  Article  PubMed  Google Scholar 

  2. Fisher B, Costantino JP, Wickerham DL, Redmond CK, Kavanah M, Cronin WM, Vogel V, Robidoux A, Dimitrov N, Atkins J, Daly M, Wieand S, Tan-Chiu E, Ford L, Wolmark N: Tamoxifen for prevention of breast cancer: report of the National Surgical Adjuvant Breast and Bowel Project P-1 Study. J Natl Cancer Inst. 1998, 90: 1371-1388. 10.1093/jnci/90.18.1371.

    CAS  Article  PubMed  Google Scholar 

  3. Freedman AN, Seminara D, Gail MH, Hartge P, Colditz GA, Ballard-Barbash R, Pfeiffer RM: Cancer risk prediction models: a workshop on development, evaluation, and application. J Natl Cancer Inst. 2005, 97: 715-723.

    Article  PubMed  Google Scholar 

  4. Tyrer J, Duffy SW, Cuzick J: A breast cancer prediction model incorporating familial and personal risk factors. Stat Med. 2004, 23: 1111-1130. 10.1002/sim.1668.

    Article  PubMed  Google Scholar 

  5. Costantino JP, Gail MH, Pee D, Anderson S, Redmond CK, Benichou J, Wieand HS: Validation studies for models projecting the risk of invasive and total breast cancer incidence. J Natl Cancer Inst. 1999, 91: 1541-1548. 10.1093/jnci/91.18.1541.

    CAS  Article  PubMed  Google Scholar 

  6. Gail MH, Costantino JP, Bryant J, Croyle R, Freedman L, Helzlsouer K, Vogel V: Weighing the risks and benefits of tamoxifen treatment for preventing breast cancer. J Natl Cancer Inst. 1999, 91: 1829-1846. 10.1093/jnci/91.21.1829.

    CAS  Article  PubMed  Google Scholar 

  7. Kaaks R, Rinaldi S, Key TJ, Berrino F, Peeters PH, Biessy C, Dossus L, Lukanova A, Bingham S, Khaw KT, Allen NE, Bueno-de-Mesquita HB, van Gils CH, Grobbee D, Boeing H, Lahmann PH, Nagel G, Chang-Claude J, Clavel-Chapelon F, Fournier A, Thiébaut A, González CA, Quirós JR, Tormo MJ, Ardanaz E, Amiano P, Krogh V, Palli D, Panico S, Tumino R, et al: Postmenopausal serum androgens, oestrogens and breast cancer risk: the European prospective investigation into cancer and nutrition. Endocr Relat Cancer. 2005, 12: 1071-1082. 10.1677/erc.1.01038.

    CAS  Article  PubMed  Google Scholar 

  8. Key T, Appleby P, Barnes I, Reeves G: Endogenous sex hormones and breast cancer in postmenopausal women: reanalysis of nine prospective studies. J Natl Cancer Inst. 2002, 94: 606-616.

    CAS  Article  PubMed  Google Scholar 

  9. Missmer SA, Eliassen AH, Barbieri RL, Hankinson SE: Endogenous estrogen, androgen, and progesterone concentrations and breast cancer risk among postmenopausal women. J Natl Cancer Inst. 2004, 96: 1856-1865.

    CAS  Article  PubMed  Google Scholar 

  10. Zeleniuch-Jacquotte A, Shore RE, Koenig KL, Akhmedkhanov A, Afanasyeva Y, Kato I, Kim MY, Rinaldi S, Kaaks R, Toniolo P: Postmenopausal levels of oestrogen, androgen, and SHBG and breast cancer: long-term results of a prospective study. Br J Cancer. 2004, 90: 153-159. 10.1038/sj.bjc.6601517.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. Colditz GA: The nurses' health study: a cohort of US women followed since 1976. J Am Med Womens Assoc. 1995, 50: 40-44.

    CAS  PubMed  Google Scholar 

  12. Colditz GA, Hankinson SE: The Nurses' Health Study: lifestyle and health among women. Nat Rev Cancer. 2005, 5: 388-396. 10.1038/nrc1608.

    CAS  Article  PubMed  Google Scholar 

  13. Rich-Edwards JW, Corsano KA, Stampfer MJ: Test of the National Death Index and Equifax Nationwide Death Search. Am J Epidemiol. 1994, 140: 1016-1019.

    CAS  PubMed  Google Scholar 

  14. Stampfer MJ, Willett WC, Speizer FE, Dysert DC, Lipnick R, Rosner B, Hennekens CH: Test of the National Death Index. Am J Epidemiol. 1984, 119: 837-839.

    CAS  PubMed  Google Scholar 

  15. Colditz GA, Rosner BA, Chen WY, Holmes MD, Hankinson SE: Risk factors for breast cancer according to estrogen and progesterone receptor status. J Natl Cancer Inst. 2004, 96: 218-228.

    CAS  Article  PubMed  Google Scholar 

  16. Hankinson SE, Willett WC, Manson JE, Colditz GA, Hunter DJ, Spiegelman D, Barbieri RL, Speizer FE: Plasma sex steroid hormone levels and risk of breast cancer in postmenopausal women. J Natl Cancer Inst. 1998, 90: 1292-1299. 10.1093/jnci/90.17.1292.

    CAS  Article  PubMed  Google Scholar 

  17. Hankinson SE, Willett WC, Manson JE, Hunter DJ, Colditz GA, Stampfer MJ, Longcope C, Speizer FE: Alcohol, height, and adiposity in relation to estrogen and prolactin levels in postmenopausal women. J Natl Cancer Inst. 1995, 87: 1297-1302. 10.1093/jnci/87.17.1297.

    CAS  Article  PubMed  Google Scholar 

  18. Hankinson SE, London SJ, Chute CG, Barbieri RL, Jones L, Kaplan LA, Sacks FM, Stampfer MJ: Effect of transport conditions on the stability of biochemical markers in blood. Clin Chem. 1989, 35: 2313-2316.

    CAS  PubMed  Google Scholar 

  19. Rubin D: Multiple Imputation for Nonresponse in Surveys. 1987, New York: Wiley

    Book  Google Scholar 

  20. Rosner B: Fundamentals of Biostatistics. 2007, Belmont, CA: Wadsworth, 6

    Google Scholar 

  21. Pencina MJ, D' Agostino RBS, D' Agostino RBJ, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27: 157-172. 10.1002/sim.2929.

    Article  PubMed  Google Scholar 

  22. Rosner B, Glynn RJ: Power and sample size estimation for the Wilcoxon rank sum test with application to comparisons of C statistics from alternative prediction models. Biometrics. 2008, May 28.

  23. Arias E: United States life tables, 2004. Natl Vital Stat Rep. 2007, 56: 1-39.

    PubMed  Google Scholar 

  24. Key TJ, Appleby PN, Reeves GK, Roddam A, Dorgan JF, Longcope C, Stanczyk FZ, Stephenson HE, Falk RT, Miller R, Schatzkin A, Allen DS, Fentiman IS, Key TJ, Wang DY, Dowsett M, Thomas HV, Hankinson SE, Toniolo P, Akhmedkhanov A, Koenig K, Shore RE, Zeleniuch-Jacquotte A, Berrino F, Muti P, Micheli A, Krogh V, Sieri S, Pala V, Venturelli E, et al: Body mass index, serum sex hormones, and breast cancer risk in postmenopausal women. J Natl Cancer Inst. 2003, 95: 1218-1226.

    CAS  Article  PubMed  Google Scholar 

  25. Bernstein L, Ross RK: Endogenous hormones and breast cancer risk. Epidemiol Rev. 1993, 15: 48-65.

    CAS  PubMed  Google Scholar 

  26. Lamar CA, Dorgan JF, Longcope C, Stanczyk FZ, Falk RT, Stephenson HE: Serum sex hormones and breast cancer risk factors in postmenopausal women. Cancer Epidemiol Biomarkers Prev. 2003, 12: 380-383.

    CAS  PubMed  Google Scholar 

  27. Madigan MP, Troisi R, Potischman N, Dorgan JF, Brinton LA, Hoover RN: Serum hormone levels in relation to reproductive and lifestyle factors in postmenopausal women (United States). Cancer Causes Control. 1998, 9: 199-207. 10.1023/A:1008838412423.

    CAS  Article  PubMed  Google Scholar 

  28. Verkasalo PK, Thomas HV, Appleby PN, Davey GK, Key TJ: Circulating levels of sex hormones and their relation to risk factors for breast cancer: a cross-sectional study in 1092 pre- and postmenopausal women (United Kingdom). Cancer Causes Control. 2001, 12: 47-59. 10.1023/A:1008929714862.

    CAS  Article  PubMed  Google Scholar 

  29. Singletary KW, Gapstur SM: Alcohol and breast cancer: review of epidemiologic and experimental evidence and potential mechanisms. JAMA. 2001, 286: 2143-2151. 10.1001/jama.286.17.2143.

    CAS  Article  PubMed  Google Scholar 

  30. Cook NR: Use and misuse of the receiver operating characteristic curve in risk prediction. Circulation. 2007, 115: 928-935. 10.1161/CIRCULATIONAHA.106.672402.

    Article  PubMed  Google Scholar 

  31. Rockhill B, Kawachi I, Colditz GA: Individual risk prediction and population-wide disease prevention. Epidemiol Rev. 2000, 22: 176-180.

    CAS  Article  PubMed  Google Scholar 

  32. Rockhill B, Spiegelman D, Byrne C, Hunter DJ, Colditz GA: Validation of the Gail et al. model of breast cancer risk prediction and implications for chemoprevention. J Natl Cancer Inst. 2001, 93: 358-366. 10.1093/jnci/93.5.358.

    CAS  Article  PubMed  Google Scholar 

  33. Wald NJ, Hackshaw AK, Frost CD: When can a risk factor be used as a worthwhile screening test?. BMJ. 1999, 319: 1562-1565.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Tice JA, Miike R, Adduci K, Petrakis NL, King E, Wrensch MR: Nipple aspirate fluid cytology and the Gail model for breast cancer risk assessment in a screening population. Cancer Epidemiol Biomarkers Prev. 2005, 14: 324-328. 10.1158/1055-9965.EPI-04-0289.

    Article  PubMed  Google Scholar 

Download references


Support for this project was provided by National Institutes of Health grants P01 CA87969, CA49449, and CA089393 (Specialized Programs of Research Excellence in Breast Cancer). GAC was supported in part by a clinical research professorship from the American Cancer Society. We acknowledge Marion McPhee for programming assistance and Jessica Bugg for preparation assistance.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bernard Rosner.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BR, GAC, and SEH contributed to the design of the study and to the analysis and interpretation of the data. JDI contributed to the revision of the manuscript and added important clinical insight. All authors read and approved the final manuscript.

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Rosner, B., Colditz, G.A., Iglehart, J.D. et al. Risk prediction models with incomplete data with application to prediction of estrogen receptor-positive breast cancer: prospective data from the Nurses' Health Study. Breast Cancer Res 10, R55 (2008).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI:


  • Breast Cancer
  • Estradiol
  • Breast Cancer Risk
  • Estradiol Level
  • Breast Cancer Risk Factor