Article Text

Are familial factors underlying the association between socioeconomic position and prescription medicine? A register-based study on Danish twins
  1. Mia Madsen1,2,
  2. Per Kragh Andersen3,
  3. Mette Gerster4,
  4. Anne-Marie Nybo Andersen5,
  5. Merete Osler2,6,
  6. Kaare Christensen1
  1. 1Unit of Epidemiology, University of Southern Denmark, Odense, Denmark
  2. 2Danish Aging Research Center, Odense, Denmark
  3. 3Section of Biostatistics, University of Copenhagen, Copenhagen, Denmark
  4. 4Unit of Biostatistics, University of Southern Denmark, Odense, Denmark
  5. 5Section of Social Medicine, University of Copenhagen, Copenhagen, Denmark
  6. 6Glostrup University Hospital, Research Center for Prevention and Health, Glostrup, Denmark
  1. Correspondence to Mia Madsen; mimadsen{at}


Objectives Although well established, the association between socioeconomic position and health and health behaviour is not clearly understood, and it has been speculated that familial factors, for example, dispositional factors or exposures in the rearing environment, may be underlying the association. The objective was to compare prescription fillings within twin pairs who are partly or fully genetically identical and share childhood exposures.

Design Twin cohort study.

Setting Denmark.

Participants Data from the Danish Twin Registry were linked to registers in Statistics Denmark and the Danish Registry of Medicinal Product statistics. A total of 8582 monozygotic (MZ) and 15 788 dizygotic same sex (DZSS) twins were included.

Outcome measures Number of prescription fillings during follow-up (1995–2005) was analysed according to education and income. Results of unpaired and intrapair analyses were compared.

Results An inverse social gradient in filling of prescriptions for all-purpose and system-specific drugs was observed in the unpaired analyses. In the intrapair analyses, associations were attenuated some in DZSS and more in MZ twins. Filling of drugs targeting the nervous system was still strongly associated with income in the intrapair analyses.

Conclusions Familial factors seem to account for part of the observed social inequality in filling of prescription medicine.

  • Epidemiology
  • Genetics
  • Social Medicine

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See:

Statistics from

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • This study is based on nation-wide register data covering more than a decade, and we present results from a largely unselected population of more than 8000 monozygotic and 15 000 dizygotic same sex twins.

  • This is the first twin study on health inequalities to report such clear-cut findings with sufficient power for meaningful interpretation of intrapair analyses and for comparison of zygosity-specific results.

  • Prescription fillings are not a perfect measure of health status, but we still argue that it is an inventive and fairly valid indicator in a Danish setting where prescription medicine has to be prescribed by a medical doctor and accounts for 90% of the total drug sale.

  • The twin study is clearly a powerful approach, but when interpreting findings, alternative explanations for the pattern of attenuation should be considered, including measurement error and unshared confounding. However, in this case we find it unlikely that these potential sources of bias could fully explain our findings.


Social inequalities in health and health behaviour have been widely demonstrated in many societies, including the Scandinavian countries.1 Often social disparities in health are not confined to the most marginalised groups of society but are expressed as a gradient over the entire spectrum of social stratification.2 The inverse relationship between socioeconomic position (SEP) and health has consistently been observed across different social indicators and for numerous types of health outcomes.3 ,4 Use of prescription medicine is one health-related outcome that has been found to be socially patterned, where Danish data have shown a greater use among people with low SEP compared with people with high SEP.5 ,6 Denmark has a tax-financed decentralised healthcare system with a partly need-dependent reimbursement system for outpatient prescription drugs. All prescription drugs have to be prescribed by a medical doctor and, compared with a number of other countries, prescription medicine constitute the majority of the total drug sales in Denmark (approximately 90%).7 Many studies have shown substantially reduced social differences in drug use when health status have been taken into account,6 ,8 ,9 which may indicate that filling of prescription medicine to a large extent reflects health status. This is supported by the fact that comorbidity indices, based on prescription drug dispensings, have been established as strong predictors of mortality,10 and that drug use has proven to be a valid indicator of self-rated health.11 ,12 However, other factors such as healthcare seeking behaviour and access to healthcare13 ,14 may also influence drug use, but the social gradient in these factors would, theoretically, be expected to be opposite that of health status, that is, greater drug use among people with high SEP compared with low SEP. Such a positive association is generally not supported by the existing Danish literature.5 ,6 ,15

Although well established, the inverse association between SEP and health and health behaviour is not clearly understood, but many factors during the life course may play a role in the production of these social health differences. Suggested pathways include: material deprivation, behavioural factors and psychological factors.16 It has also been suggested that it is not SEP in adulthood per se that influences health, but underlying familial factors already at play in early life,17 ,18 for example, rearing environment or genetic make-up are associated both with the obtainment of adult SEP and health status later in life.

In a recent commentary, Gilman and Loucks19 discuss why childhood environment may be such an important underlying factor in the production of social health differences. They argue that because childhood encompasses developmentally sensitive periods for the acquisition of cognitive, psychological and social skills that have a profound impact on an individual's capacity for long-term health and provide opportunities (and constraints) for adult socioeconomic conditions, it can potentially influence the association between any adult risk factor and later life health.19 The potential role of genetic setup in explaining health inequalities is reviewed by Mackenbach,20 who suggests that the most plausible set of hypotheses relates to genetic influences of personal attributes such as cognitive ability, personality and bodily and mental fitness that may influence both SEP and adult health status.

With the discordant twin pair design, it is possible to investigate the effect of adult SEP on prescription fillings eliminating confounding from factors shared by a pair of twins, exploiting that twins are partly or fully genetically identical and have experienced a similar childhood environment due to their common upbringing. This approach has been applied in a number of previous twin studies investigating different health outcomes, including mortality,21 hospitalisations,22 breast cancer23 and also softer endpoints such as depression scores and self-rated health.24 ,25 However, the majority of these studies suffer from power limitations either due to small sample sizes and/or rare outcomes, making interpretations difficult.

In this study we aimed to investigate if there was an effect of SEP in adulthood on the number of prescription fillings above and beyond the effect of shared familial factors. Filling of prescription medicine is a quantitative health measure rendering sufficient power for the intrapair comparisons. Differences in genetic relatedness between monozygotic (MZ; genetically identical) and dizygotic (DZ) twins (share on average 50% of segregating genes) may further give an indication of the possible types of confounding (ie, genetic or environmental factors). In the theoretical situation of no measurement error and no confounding from unshared factors, genetic confounding would be indicated if an observed social gradient in prescription medicine was partly attenuated in DZ twins and fully attenuated in MZ twins, when shared familial factors had been taken into account. Likewise, an attenuation of similar size in DZ and MZ twins would be compatible with shared environmental confounding. Finally, to support a SEP effect on prescription fillings not due to confounding from shared familial factors, the association would have to persist in both DZ and MZ twins.

To support the utility of prescription medicine as a health indicator we also analysed its association with mortality.


Study population

Data from the Danish Twin Registry were linked with information in the Danish Registry of Medicinal Products and administrative registers in Statistics Denmark. The study population consisted of 24 370 twins, including 8582 MZ and 15 788 dizygotic same sex (DZSS) twins, who were born during 1921–1965 and alive and resident in Denmark at the beginning of follow-up (1 January 1995), after exclusions due to triplets (n=384) unknown zygosity (UZ; n=3731), and dizygotic opposite sex pairs (n=16 710), missing data on education or income (n=845) and non-intact twin pairs in which one twin was dead, emigrated or had missing information (n=2847). To evaluate the representativeness of the results from the unpaired twin analyses, we also analysed a 5% population-based sample of the Danish population (n=137 300).

Social indicators

We used information on the highest attained education in 1995 at age of 30+ of the study population, categorised according to the International Standard Classification of Education of 199726 into: primary/lower secondary education (basic compulsory education), upper secondary education (secondary) and tertiary education (Bachelor's degree and above). Income was measured as an average of the equivalised gross household income (Danish currency) in 1994–1995, grouped into quartiles within each birth cohort. ‘Equivalised’ income is adjusted to take into account that households can share resources. Inspired by the standard Organisation for Economic Co-operation and Development equivalence scale,27 the household-denominator was constructed according to an equivalence scale that assigned a value of 1 to the first household member, of 0.7 to the second household member and of 0.5 to each additional household member. Calculating quartiles by birth cohort is aiming at taking into account that income is not fully comparable at different ages during the life course. An implication of this is that the income measure reflects a relative rather than absolute income level.

Prescription medicine

The Danish Registry of Medicinal Products was established in 1995 and contains information on the total sale of prescription medicine in Denmark registered by the pharmacies via a computerised accounting system. The drug information is recorded according to the Anatomical Therapeutic Chemical classification system (ATC)28 where the active substances in drugs are divided into different groups, primarily according to the organ system on which they act. Other information includes date of dispensing and daily defined doses (DDD), defined as the assumed average maintenance dose per day for a drug used for its main indication in adults.28 Outcomes of interest were defined as the number of prescription fillings during follow-up for all-purpose medication (ACT-groups: A, B, C, D, G, H, J, L, M, N, P, Q, R, S, V, excluding Q which is for veterinary use, and G02B and G03A, which are for contraceptive use) and four groups of system-specific medication: ATC-C (heart/circulatory system), ATC-N (nervous system), ATC-R (respiratory system) and ATC-A (alimentary tract/metabolism). If no ATC code was assigned a prescription, it was excluded (1%). DDD was used as a supplementary outcome measure.

Other covariates

Sex and age were included in the analyses. Analyses were performed separately for MZ and DZSS twins. Zygosity had previously been determined by questionnaire,29 a method that has proven valid with an overall misclassification frequency of less than 5% validated against a classification based on genetic markers.30

Data analysis

Analyses were carried out using Cox regression for recurrent events (proportional means regression), thereby accommodating multiple events (prescription fillings) per individual.31 Age was used as underlying time-variable and participants were followed from age at baseline to the age at the end of follow-up (31 December 2005) or death or emigration, and the timing of each prescription filling during this period was recorded. Both unpaired and intrapair analyses were performed. In the unpaired analyses robust SEs were used to account for the interdependence of observations within pairs and within each individual. Results of the unpaired analyses were compared with the intrapair analyses, which were performed using a stratified Cox model: ENij(t)=μ0j(t)eβ'z_ij. Here ENij(t) denotes the mean number of prescription fillings at time t for an individual i belonging to twin pair j, μ0j(t) the pair-specific baseline mean, and eβ the effect of the covariates adjusted for shared familial factors. Thus, in this analysis the mean number of prescription fillings was compared within twin pairs, thereby controlling for familial factors shared by a pair of twins. In the stratified model only SEP-discordant twin pairs contribute information to the β-parameter for SEP, although concordant pairs may contribute information to other β-parameters.

As counting process models do not allow simultaneous events, only one prescription of interest per day was included in the analysis.

Analyses were performed separately for zygosity and sex, and potential interactions between the social indicators and zygosity, sex and age, respectively, were evaluated by including interaction terms in the models and using Wald test statistics. To analyse the association between prescriptions and mortality, we used a Cox regression model, continuously updating the number of prescription fillings over time. Descriptive intrapair correlations were calculated as polychoric correlation coefficients.


Since DDD may reflect the actual drug use more precisely than the number of prescription fillings, mean DDD were also analysed by means of a fixed-effects linear regression model taking censoring into account.32 This was treated as a subanalysis due to the highly right-skewed distribution of DDD, making linear regression a less optimal choice of model.

Permission for linkage and use of data for this study was obtained from the Danish Data Protection Board (2000-54-0047).


According to table 1, showing the distribution of education and income in the population, men were better educated than women. Generally, the education of the twin population was comparable to that of the population-based sample. Only UZ twins had a markedly lower educational status and a higher proportion of missing information. The table also shows the mean income within each quartile. Twins have a somewhat higher mean income as compared to the population-based sample.

Table 1

Descriptives on educational status and income in a 5% population-based sample (n=144 749) and a population of Danish twins (n=48 887), including intact and non-intact twin pairs, stratified on sex and zygosity

The proportion of education-discordant twin pairs was 30% for MZ twins and 42% for DZSS twins, while 61% MZ twins and 68% DZSS twins were discordant on income (table 2).

Table 2

Descriptives on intrapair discordances on educational status and income, stratified on zygosity in a population of Danish twins (n=24 370)

In the twin population, 1.9% did not redeem any drug prescriptions during the total follow-up period and 25% redeemed less than one prescription/year, on average (results not shown). Overall, women filled more prescriptions than men for all age groups and the number of prescription fillings increased with age in both sexes (table 3).

Table 3

Mean number of prescription fillings after 1 and 10 years of follow-up in a 5% population-based sample and a population of Danish twins, stratified on age and sex

ATC-N drug prescriptions were the most frequently redeemed in the youngest age group. The mean number of prescription fillings in the twin population was similar to that in the population-based sample, although slightly higher in the population-based sample.

In order to get a better idea of the clinical outcomes underlying the observed prescription patterns, we have included a more detailed description of the most incident types of drugs at the 3-digit and 4-digit ATC level (see online supplementary table S1). From the table it appears that the ATC-C category primarily comprises drugs indicated for hypertension and ischaemic heart disease, while the ATC-N category includes drugs indicated for a wide range of different diagnoses, including chronic pain disorders, insomnia, anxiety, psychoses, depression, etc.28 In the ATC-R category the main types of drugs are drugs for obstructive airway disease such as asthma and chronic obstructive lung disorder (COPD).28 Finally, drugs for acid-related disorders and diabetes are the main constituents of the ATC-A category. The results of an analysis of the association between SEP and prescription fillings for all-purpose and system-specific drugs are presented in table 4.

Table 4

Ratio (MR) of the mean number of prescription fillings (1995–2005) for all-purpose and system-specific drugs according to educational status and income in a 5% population-based sample (137 300) and a population of Danish twins (n=24 370), showing results from unpaired and intrapair analyses stratified on zygosity

For all-purpose drugs, MZ twins in the highest income quartile had a 30% lower mean number of prescription fillings than MZ twins in the lowest income quartile (MR=0.70 (95% CI 0.65 to 0.77)). Similarly, MZ twins with a tertiary education had 17% lower mean number of prescription fillings than twins with a primary education (MR=0.82 (95% CI 0.75 to 0.89)). The association was slightly less strong for ATC-C drugs while stronger associations could be observed for the other system-specific drugs, particularly for ATC-N. In the unpaired analyses, no large differences could be observed between MZ and DZSS twins and there was no strong evidence for an interaction between SEP and zygosity in the unpaired analysis (education p=0.81, income p=0.95). In contrast, the effect of SEP was moderated by zygosity in the intrapair analyses (education p=0.02, income p=0.05), suggesting a more attenuated effect of SEP in MZ twins than in DZSS twins. In MZ twins, generally no statistically significant effect of education could be observed, while income was still somewhat associated with all-purpose medication and quite strongly with ACT-N drugs. In DZSS twins, a minor or moderate attenuation could be observed, and all effect parameters were statistically significant, except for ATC-C drugs. In the population-based sample, the associations were similar to that in the overall twin population.

When using DDD (for all-purpose drugs) as outcome measure (table 5), results similar to the results based on the number of prescription fillings were obtained, although the attenuation patterns in MZ twins were somewhat more clear in this analysis, showing no statistically significant associations with DDD of neither education nor income. In the analysis of number of prescriptions as a predictor of mortality, we found a 60% increase in mortality rate during follow-up for a person who, on average, filled one prescription per month during the follow-up period compared with a person without any prescription fillings (HR=1.60 (95% CI 1.25 to 2.05)).

Table 5

Results from a linear regression model showing number of daily defined doses (DDD) for all-purpose drugs according to educational status and income in a population of Danish twins (n=24 370) stratified on zygosity


Summary of findings

In this large Danish twin population we demonstrated an inverse social gradient in fillings of prescription medicine for all-purpose and system-specific drugs. That is, a greater use of prescription medicine among those with a low SEP thus suggesting a worse health status of this group. Furthermore, the number of prescription fillings was shown to be associated with a greater risk of mortality. In the intrapair analyses, where shared familial factors were controlled per design, the association between SEP and prescription fillings attenuated to a wide extent within MZ twins and statistically significantly less in DZSS twins. For education, there was generally no association with prescription fillings in the intrapair analyses of MZ twins, except for drugs targeting the nervous system, while there was still some association of income with all-purpose medication and most system-specific drugs. However, no effect of income was seen in the intrapair analyses when DDD were used as outcome measure. In spite of some variations between system-specific drugs, the overall pattern was the same as for all-purpose drugs, except for ATC-N drugs and to a lesser degree ATC-A drugs, which still showed some effect of SEP in the intrapair analyses.


The fact that the association was attenuated in the intrapair analyses suggests that underlying factors shared by a pair of twins account for some of the observed association between SEP and prescription fillings. However, there may be a number of alternative or complementary explanations for this finding: first, part of it is likely caused by measurement error. Measurement error is random and therefore not shared by twins. In the intrapair analysis, variance due to shared factors is removed, and a relatively larger proportion of exposure variance is thus due to measurement error. This will result in increasing attenuation of parameter estimates in the intrapair analyses.33 However, since data come from official registers it can be expected that education and income have limited measurement error. Second, the behaviour of one twin could be affected by the behaviour of the other twin, that is, ‘social interaction’ either by ‘imitation’ or ‘differentiation’. The implications of a potential imitation could be smaller within-pair differences, which could be misinterpreted as confounding from underlying familial factors. This may be of some concern in this study since fillings of prescription medicine is likely to have a behavioural component. Another implication of the discordant twin pair approach is that only twin pairs who are discordant on exposure contribute with statistical information to the estimation of SEP effects in the intrapair analysis. This could potentially induce bias in the comparison of the unpaired and intrapair analysis, since the two analyses are based on only partly overlapping populations. Since MZ twins are less likely to be discordant, the selection bias may be most severe in MZ twins.34 We examined this concern by restricting the unpaired analysis only to include groups with exposure-discordant twins. Results were similar to those of the main analysis in both MZ and DZ twins. A further issue relates to the fact that twins who are discordant on exposure are special in the sense that they, in spite of overlapping genes and rearing environment, differ in their SEP. Thus, unshared factors may explain these differences. This implies that the intrapair analysis may still be confounded by unshared confounding factors, and that this confounding may be enforced in the intrapair analysis.33 Generally, confounding from unshared factors is thought to be more severe when confounders are less correlated than exposure.33 This may or may not be true in our case, since we can only hypothesise on the correlation of unmeasured confounders. Nevertheless, the correlation of education and income (MZ=0.79, DZ=0.58 and MZ=0.49, DZ=0.33) are more or less of the same magnitude as phenotypes like cognitive ability35 and body mass index36 that have previously been investigated using similar designs. In any case, as shown by Frisell et al35 in order for increased confounding to explain the attenuation of the association observed in this study, the potential non-shared confounders would have to create a positive association between education and prescription fillings. In most cases, however, we would expect a negative association, although filling of prescription medicine is indeed a complex health outcome and it is theoretically possible for a confounder to have a dual effect on prescription fillings by (1) increasing the capability of getting the relevant prescriptions, for example, due to more regular heath checks or better understanding of symptoms (2) at the same time reducing the need for prescription fillings due to better health status. Nevertheless, in a Danish context health status does seem to be the main determinant of prescription medicine10 ,37 ,38 and was clearly associated with mortality in this data. Thus, in the presence of non-shared confounding we would expect the net effect of any confounding to produce a negative association between education and prescription fillings resulting in a bias away from the null,33 and therefore it could not explain the observed attenuation of effect in the intrapair analyses observed in this study.

The finding of greatest attenuation in MZ twins is compatible with the scenario of genetic confounding described in the introduction. Since the genetic overlap is greater in MZ twins than DZSS twins, one would expect a larger attenuation in MZ twins, if genetic factors gave rise to the confounding. However, a number of other factors could produce a greater attenuation in MZ twins: first of all, the aforementioned measurement error may be most pronounced in MZ twins, due to the fact that increasing correlation in exposure will result in a greater attenuation of intrapair parameter estimates.33 In addition, the comparison of strength of association between MZ and DZ twins to indicate the source of familial confounding relies on the equal environment assumption: that is, MZ and DZ twins share their environment equally. However, this assumption may be challenged, and if it does not hold, that is, if MZ twins not only share their genes but also their environment more closely than DZ twins,39 shared environmental factors could also produce a greater attenuation of effect in MZ twins. However, the interplay between environment and genes is complex, and it is conceivable that the environment to a certain extent is shaped in response to the genetic dispositions of an individual. If, for instance, a genetically determined similar appearance of MZ twins also elicits a more similar treatment from their peers, then these similar environmental stimuli can also be ascribed to genetic factors.40

Thus, although several factors may contribute to the observed attenuation of association between SEP and filling of prescription drugs in the intrapair analyses and to the differences in strength of association between MZ and DZ twins, the findings suggest that a part of the observed social inequality in prescription fillings is explained by either genetic or environmental factors shared by twins.

The finding of a persistent association of particularly income with fillings of prescriptions for ATC-N drugs clearly stood out from the general results. We therefore carried out a post hoc subanalysis of specific categories of ATC-N drugs to test whether it were certain types of drugs driving the association. The analysis revealed that it was mainly for analgesics and psycholeptics that the association persisted, while the association was somewhat attenuated for psychoanaleptics (primarily antidepressants). This finding may support the notion of reverse causation, that is, that health status influences SEP rather than the other way around, since it is likely that diagnoses such as chronic pain disorders, psychoses, anxiety and bipolar disorders associated with these types of drugs may affect the work ability and thus the income of a person. Income is much more sensitive to reverse causation than education, which is usually obtained early in life and remains constant once it is obtained.

Existing literature

In this study we argue that fillings of all-purpose and system-specific drug prescriptions to some extent serve as an indicator of general health status and of conditions such as ischaemic heart disease, hypertension, COPD and diabetes. The finding of weaker effects of SEP, when underlying familial factors were taken into account, accords with a large part of the life-course literature that supports the importance of early life factors in the aetiology of chronic disease in adulthood.41 ,42

In addition, a number of previous twin and sibling studies have investigated the association between various SEP indicators and health outcomes. While some studies have shown strong or some evidence of a causal effect of adult SEP4347 other suggest no or little effect of SEP, when familial factors have been adjusted for.22 ,4851 Other studies show mixed results.21 ,23–25 ,47 ,52–54 However, the studies vary with respect to a number of important parameters hampering overall comparisons. Most importantly, investigated health outcomes differ and may have different aetiologies and show different social gradients and patterns of familial confounding, but even for specific-health outcomes there is no clear picture of evidence. For example, for mortality, two sibling45 ,46 and two twin studies21 ,44 have investigated the association with education. Although all studies support some effect of education after controlling for shared family factors, important differences exist: the twin study by Madsen et al21 demonstrated an effect of education within twins only in older cohorts of men, but not in younger men and not in women. The influence of shared family factors was thus suggested in the latter, but not in the first subgroup of this study. The two sibling studies also suggested some effect of shared family factors, while the twin study by Lundborg43 did not support any such effect. For self-rated health, the majority of the existing twin studies24 ,25 ,53 have shown mixed results, although Osler et al47 and Lundborg43 did find a statistically significant association between social class respective education and self-rated health in the intrapair analyses. One important reason for the inconsistencies may be that many studies have suffered from power limitations in the intrapair analyses due to a small number of events.21 ,23 ,25 ,47 ,53 This is not a concern in this study, where almost all of the exposure discordant twins contributed with information to the intrapair analyses, thus providing some of the most powerful twin results on social inequality in health to date.

Strengths and limitations

The twin approach is clearly a strength of this study, since it makes it possible to adjust for unmeasured background factors shared by twins that may underlie the association between SEP and health. Furthermore, the large sample size and the long follow-up allow a great number of events to be analysed. In addition, the use of population-based registers has the advantage of relatively little selection and measurement error in the study population. Finally, this study not only focuses on prescription medicine in general, but also on drugs targeting specific organ systems which may give us a better idea of potential mechanisms.

Among study limitations is the fact that filling of prescription medicine is not a perfect measure of health status. We do, however, argue that the observed association with mortality and the inverse relationship with education and income likely reflect a greater need for medicine among people of lower SEP. In case that access to healthcare also plays a role, presumably favouring people of high SEP, the social gradient in health status in this study could be underestimated. Although data did not contain information on over-the-counter medication, it only constitutes a minor part of the total drug consumption in Denmark,7 and since we were interested in drugs as an indicator of health status, it could be argued that the most valid approach is to focus on drugs that require a formal medical indication (ie, prescription).

As models for counting processes only allow one event at a time, the number of prescription fillings was defined as the main outcome of interest. Therefore, it was reassuring that the subanalysis on DDD led to similar results.

Although the distribution of SEP indicators and prescription fillings were not entirely identical in twins and the population-based sample, the associations between SEP and prescription fillings were indeed, which supports generalisability to the general population in line with previous studies that show that health and mortality in twins are similar to the background population.55 ,56 Furthermore, the fact that the analysis based only on discordant twin pairs show marginal results similar to those based on the total twin sample was reassuring.

Conclusion and perspectives

In this study, based on a comprehensive twin cohort, we observed a social inequality in fillings of prescription medicine, showing fewer prescription fillings among people with high SEP. This association attenuated within twin pairs for most types of prescription fillings, except those targeting the nervous system. The attenuation was most pronounced in MZ twins. Measurement error can contribute to the observed attenuation of association, but it is likely to account only for a minor fraction, considering the exposure and outcome measures in this study. As for potential non-shared confounding we would, if anything, expect it to bias the estimates away from the null.

We therefore interpret the findings as an indication that at least a part of the observed social inequality in prescription fillings is explained by shared familial factors. The implications of these findings are by no means a disqualification of SEP as a determinant of health, but rather an indication that family matters in the production of social health inequalities; most likely through early environment, and possibly also genetic factors. Recognising biological components and early life exposures as important explanatory factors may improve our understanding of the multifaceted aetiology of health inequalities over the life course and increase the potential for effective and targeted interventions.


Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Files in this Data Supplement:


  • Contributors KC and MO formulated the original study idea and provided funding. MM, KC, MO and A-MNA further developed the study idea and design of the study. MM carried out the data analyses in collaboration with PKA and MG. All authors participated in interpretation and discussion of the findings. MM drafted the first version of the manuscript, but all authors have contributed to the revisions and have approved the final version of the paper.

  • Funding This work is supported by a grant from the VELUX Foundation (95-103-11419), a grant from the Danish Cancer Society (SU08001) and from the Danish Council for Independent Research (social sciences, 09-06008). The funding bodies had no influence on the contents of this paper.

  • Competing interests None.

  • Ethics approval The Danish Data Protection Board (2000-54-0047).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.