Article Text

Original research
Change in prevalence rates of physical and sexual intimate partner violence against women: data from two cross-sectional studies in New Zealand, 2003 and 2019
  1. Janet Fanslow1,
  2. Ladan Hashemi1,
  3. Zarintaj Malihi1,
  4. Pauline Gulliver1,
  5. Tracey McIntosh2
  1. 1 School of Population Health, The University of Auckland Faculty of Medical and Health Sciences, Auckland, New Zealand
  2. 2 School of Māori Studies and Pacific Studies, Faculty of Arts, The University of Auckland Faculty of Arts, Auckland, New Zealand
  1. Correspondence to Dr Janet Fanslow; J.Fanslow{at}auckland.ac.nz

Abstract

Objectives To explore changes in reported prevalence of physical and sexual intimate partner violence (IPV) between 2003 and 2019. The impact of sociodemographic differences between the two samples and between group differences were also examined. Changes in attitudes supportive of violence and in help-seeking behaviour following disclosure were also explored.

Design Two cross-sectional studies.

Setting and participants Cross-sectional studies on family violence conducted in New Zealand in 2003 and 2019. Ever-partnered female respondents aged 18–64 years old were included (2003 n=2674, 2019 n=944).

Main outcome measures Prevalence rates of lifetime and past 12-month physical and sexual IPV, attitudes towards gender roles and acceptability of a man hitting his wife, help sought and received following disclosure were compared between the study years.

Results Lifetime prevalence of physical IPV was unchanged between 2003 and 2019 (AOR=0.89; 95% CI 0.73 to 1.08). There was a significant decrease in the proportion of women who reported experiencing 12-month physical IPV (AOR=0.53; 95% CI 0.29 to 0.97). Small reductions in rates for lifetime sexual IPV were also observed (AOR=0.74; 95% CI 0.59 to 0.95). In 2019, fewer women agreed with one or more statements supportive of traditional gender roles (48.1% (95% CI 45.7% to 50.5%) in 2003; 38.4% (95% CI 33.8% to 43.2%) in 2019). A significant decrease was noted in the proportion of women who sought help from informal sources (from 71.3% (95% CI 68.1% to 74.2%) in 2003 to 64.6% (95% CI 58.7% to 70.1%) in 2019). No significant changes in seeking help from formal sources, or perceived helpfulness from any source were noted.

Conclusion While the reductions in 12-month physical and lifetime sexual IPV are positive, prevention efforts need to be established, maintained and strengthened to address the substantial lifetime prevalence of IPV. Efforts to strengthen responses from formal and informal sources continue to be needed.

  • epidemiology
  • public health
  • social medicine
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.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, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

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

  • The current investigation used large, representative samples of women from population-based surveys in 2003 and 2019.

  • Regular and comparable surveys of violence exposure, agreement with attitudes supportive of violence and help-seeking behaviours provide an understanding of the effectiveness of population-based policies and programmes.

  • True prevalence estimates may be higher as it is expected that women in severely abusive relationships would be unable or unwilling to participate in both surveys.

  • Observed changes may reflect societal changes or environmental factors not considered in this investigation.

  • Regular and comparable surveys of violence exposure are required to determine if the observed changes are sustained and represent a trend.

Introduction

Intimate partner violence (IPV) has been reported by the UN Secretary-General (2006) as ‘the most common form of violence experienced by women globally’.1 IPV includes physical and sexual violence, as well as psychological abuse, controlling behaviour and economic abuse.

Efforts to respond to IPV in high-income countries include the introduction of legislation or national action plans, and strengthening the non-for-profit sector to respond to the violence experienced.2 However, the effectiveness of these strategies is not clear, as there is a lack of consistent and reliable data available to monitor changes in the prevalence of IPV over time.

The limited research available tends to rely on analysis of IPV homicide data, or other forms of administrative data from agencies such as health providers, police or courts.2 While providing useful insights, these data do not reflect the magnitude of the problem at the population level, as many who experience IPV frequently do not present to services, or the underlying cause of their presentation may not be identified or recorded.2 3

Other attempts to measure changes in IPV occurrence over time have relied on data from general crime victimisation surveys,4 but the overall framing of these questionnaires (ie, surveys about ‘crime’) tends to lower the reporting of the violent behaviours within intimate relationships.2 3 Surveys conducted for other purposes (eg, health surveys) which include a dedicated module on family violence provide some information, but can also be problematic, as space limitations for specific modules means that they might not be able to include questions that canvas the full range of violent experiences.5

The emerging consensus is that ‘population-based stand-alone surveys are the instruments of choice’ for collecting statistics on violence against women.6 To date, specific violence against women surveys have been carried out in several high-income countries (eg, USA,7 Canada,8 Australia,9 European Union,10 Finland,11 12 Spain,13 New Zealand14). However, with the exception of Australia and Finland, the surveys have generally been one-off efforts and thus do not allow for time-related comparisons. Without regular, comparable surveys, it is not possible to determine if there are overall changes in the occurrence of IPV, or if there are differential patterns of change for specific subgroups within the population.

According to the WHO, violence results from the complex interplay between individual, relationship, social, cultural and environmental factors.15 The ecological model has been important in helping determine risk and protective factors associated with violence occurrence, but also holds promise for prevention, as it carries the assumption that changes in contributing factors can potentially lead to changes in prevalence.16 To date, the limited research that has explored differences in the prevalence of IPV over time has suggested that population-level changes in demographic factors, such as shifts in age, education, relationship status, and socioeconomic factors may contribute to the observed prevalence changes.4 6 17 18 However, changes in environmental and social norms that may condone or help perpetuate violence, and associated effects on violence occurrence have received scant attention in the research.

Community-level norms, such as acceptance of ‘traditional’ gender roles and beliefs in the justification of ‘circumstances in which it is acceptable for a man to hit his wife’ are associated with perpetration of IPV.19 In some countries, women’s acceptance of these attitudes has been found to be associated with increased risk of IPV victimisation.20 For these reasons, attitudes have been a key target of community education campaigns aimed at preventing violence against women.21 However, to date, there has been little examination of the effectiveness of these initiatives at changing attitudes, or on any associated changes in violence rates.21–23

New Zealand is one of few high-income countries where more than one comprehensive population-based survey of violence against women has been conducted: the first survey was conducted in 2003, and the second survey in 2019. Between the two surveys, a series of actions were taken to address family violence including; legislation (eg, amendments to family violence law and protection for victims act), and prevention campaigns (eg, the Family Violence: It’s not ok national campaign, and the Accident Compsensation Corporation (ACC)-funded Mates and Dates high schools programme on healthy relationships). Many of these initiatives have focused on addressing physical and sexual violence and have included strong messaging about the importance of help-seeking by those experiencing violence. Comparable surveys on attitudes supporting violence over time may provide evidence about the impact of such campaigns at the population level.

In the current study, using data from two New Zealand cross-sectional population-based surveys we aimed to: (1) describe changes in the reported prevalence rates of physical and sexual IPV between 2003 and 2019, (2) examine whether changes in women’s sociodemographic characteristics were associated with changes in IPV prevalence rates, and (3) determine whether changes in the reported prevalence rates were consistent across population subgroups. We also sought to determine if there were (4) changes in attitudes supportive of violence and (5) changes in help-seeking by those who reported experiencing IPV.

Method

Procedure and participants

Data were drawn from two cross-sectional studies on family violence conducted in New Zealand in 2003 and 2019. A comprehensive description of the methods used in the 2003 and 2019 surveys has been previously presented.14 24 A brief description of the two surveys is presented here.

The 2003 study was conducted in Auckland and Waikato regions. For the 2019 study, Northland was also included in the sampling.

Sampling strategies were similar in both surveys. A population-based cluster sampling scheme with a fixed number of dwellings per cluster was used for both studies. Primary sampling units (PSUs) were based on meshblock boundaries which contain between 50 and 100 dwellings. The starting point consisted of a randomly selected street and street number within each PSU. Interviewers made up to seven visits to each selected household to identify and recruit study participants. Non-residential, aged-care and short-term residential properties were excluded from both surveys. Interviewer training and support procedures were comparable across survey waves.

Eligibility

To be eligible to participate in the survey, household members needed to be able to speak conversational English, have lived in the household for at least 1 month and slept in the house for four or more nights a week.

Of the households invited, 88.3% in 2003 and 78% in 2019 agreed to participate. Of the eligible women, 75.8% in 2003 and 63.7% in 2019 participated, yielding an overall response rate of 66.9% in 2003 and 63.7% in 2019. Figure 1 demonstrates the number of people invited and those who were interviewed and included in the analyses for each survey year.

Figure 1

Flow diagram of female participants in the 2003 and 2019 population-based studies on family violence in New Zealand.

Participants of the 2003 study were 2855 women aged 18–64 years. In 2019, the eligible population was expanded to include women and men aged 16 years and older resulting in 2888 completed interviews (n=1464 women, n=1423 men, n=1 other). For the purpose of this paper, only ever-partnered women aged 18–64 years from each sample were included, equivalent to almost 94% of all women aged 18–64 years surveyed in both waves (2003, n=2674; 2019, n=944).

Representativeness

In both surveys, the ethnicity, marital status, and area-level deprivation distribution of the samples were closely comparable to the general population; however, both samples were under-represented for younger women (ages 20–29 in 2003, 16–29 in 2019).14 24 Demographic characteristics of ever-partnered women aged 18–64 years in the 2003 and 2019 surveys are presented in table 1.

Table 1

Demographic characteristics of ever-partnered women aged 18–64 years in 2003 and 2019 surveys

Safety and ethics considerations

Ethics and safety recommendations for research on violence against women were followed throughout the research.25 One individual was randomly selected from each household for the interview. In households with more than one eligible resident, the participant was randomly selected. Interviews were conducted in privacy with no one over the age of 2 years present. At the completion of the interview, interviewers provided all respondents with a list of approved support agencies regardless of disclosure status. Written informed consent was obtained from all participants.

Patient and public involvement

No patients or members of the public were involved in the design, conduct or reporting or dissemination plans of our research.

Study instrument and measures

To collect data, the WHO Multi-Country Study on Women’s Health and Domestic Violence Against Women26 was used in both surveys.

‘Intimate partners’ included male current or ex-partners that the women were married to or had lived with, or current regular male sexual partners. Definitions are presented in online supplemental table 1 for: physical and sexual IPV; sociodemographic characteristics; attitudes towards gender roles, acceptance of attitudes justifying a man hitting his wife, and sources of help sought (who was told about the IPV) and help received (sources who provided help). All questions used for analyses were identical in the two surveys.

Supplemental material

Analytic procedure

To explore whether there were any underlying differences in demographic characteristics of the respondents at the two time periods, the 2003 and 2019 samples were compared in terms of age, relationship status, education attainment, access to an independent source of income, and area-level deprivation using χ2 tests.

Then, the prevalence rates of physical and sexual IPV were compared between two samples with results presented as percentages with 95% CIs. As the results for ‘moderate’ and ‘severe’ physical IPV showed similar patterns to any physical IPV, in the following analyses, only the results for any physical IPV are presented. Any act of sexual IPV was considered as severe. To identify evidence of differences in the estimated prevalence over time, ORs and 95% CIs for reported experience of physical and sexual IPV were calculated using univariate logistic regression models, with the study year as the predictor. The same procedure was followed for assessing differences in women’s attitudes towards gender roles, attitudes towards acceptability of a man hitting his wife, help sought, and help received between the study years. For help-seeking variables, the analyses were restricted to women who reported lifetime experience of physical or sexual IPV only.

Then, to determine if the noted differences in the prevalence rates of IPV between the two study years found in the univariate analyses remained significant after controlling for sociodemographic characteristics, the following steps were taken:

  • First, the association between each sociodemographic characteristic and each type of IPV (lifetime or 12 month physical and sexual IPV) was explored using univariate logistic regression models with results presented as unadjusted ORs with 95% CIs.

  • Second, multivariate analyses were conducted, with the study year and sociodemographic characteristics included, and results were presented as adjusted ORs (AOR) with 95% CIs.

Finally, to determine whether the noted changes in the reported prevalence rates were consistent across population subgroups, multivariate logistic regression models with interaction terms (between each sociodemographic characteristic and the study year) were tested. Potential confounders (eg, age, education, relationship status, independent income, and area-level deprivation) and the study year were included in these analyses.

All analyses were performed on a pooled dataset of the two samples. Missing data including: do not know, do not remember, and no responses were excluded from all analyses. Less than 4% of any variable had missing data in both surveys. All analyses were conducted using Stata/SE V.15.127 survey commands to allow for stratification by sample location (region), clustering by PSUs, and weighting of data to account for the number of eligible participants in each household.

Results

Differences between two study samples in terms of sociodemographic characteristics are presented in table 1. In general, there were more women over 45 years in 2019 (51.4%) compared with 2003 (39.3%). Additionally, a higher proportion of the sample had attained tertiary education in 2019 (65.1%) compared with 44.8% in 2003. A smaller proportion of women in 2019 reporting having an independent source of income (72.5%) compared with 79.5% in 2003.

Physical IPV

Changes in physical IPV prevalence rates

Lifetime physical IPV prevalence: the lifetime prevalence of physical IPV remained relatively unchanged between 2003 and 2019, with almost 30% of ever-partnered women aged 18–64 reporting having experienced at least one episode of physical violence (table 2). After controlling for sociodemographic factors, adjusted AORs showed no significant difference in the reported prevalence rates of lifetime physical IPV between the two study years (AOR=0.89; 95% CI 0.73 to 1.08).

Table 2

Characteristics of women reporting lifetime and past-12 month Physical IPV in the pooled database from two cross-sectional studies in New Zealand

12-month physical IPV prevalence: the 12-month prevalence of physical IPV decreased from 5% in 2003 to 2·4% in 2019 (OR=0.46; 95% CI 0.27 to 0.79). The AOR showed that, after controlling for sociodemographic factors, the decrease in 12-month physical IPV was attenuated but still remained significant (AOR=0.53; 95% CI 0.29 to 0.97).

Characteristics of women reporting lifetime and past-12 month physical IPV

Lifetime physical IPV: All sociodemographic factors were significantly associated with reporting lifetime physical IPV in the multivariate model, with the exception of ‘access to independent income’ and ‘educational attainment’. Women aged 25 years and above were more likely to report having experienced at least one act of lifetime physical IPV. Compared with married women, a higher proportion of women who were cohabiting, divorced or widowed reported experiencing lifetime physical IPV. Similarly, those who were living in moderately or most deprived areas were more likely to report the experience of lifetime physical IPV compared with those living in the least deprived areas (table 2).

Past 12-month physical IPV: at the multivariate level, age and relationship status were significantly associated with reports of experiencing past 12-month physical IPV. A lower proportion of women aged 45 years and older reported experiencing past 12-month physical IPV compared with those younger than 45 years. A higher proportion of those who were cohabiting or divorced compared with those who were married reported this experience (table 2).

Sexual IPV

Changes in sexual IPV prevalence rates

Lifetime prevalence: a significant decrease in the reported lifetime prevalence of sexual IPV was found in the univariate analysis, from 16.9% in 2003 to 13.1% in 2019 (OR=0.74; 95% CI 0.59 to 0.92). After controlling for sociodemographic variables, the significant decrease in the reported experience of lifetime sexual IPV remained unchanged (AOR=0.74; 95% CI 0.59 to 0.95).

12-month prevalence: no significant differences in the 12-month prevalence rates of sexual IPV between the two study years was found in the univariate analysis (approximately 1% in both study years) (OR=0.50, 95% CI 0.23 to 1.10). After controlling for sociodemographic factors, the non-significant difference in 12-month sexual IPV between two study years remained unchanged (AOR=0.50; 95% CI to 0.19 to 1.35).

Characteristics of women reporting lifetime and past 12 month sexual IPV

Lifetime sexual IPV: at the multivariate level, age, relationship status, education attainment and area-deprivation level were significantly associated with lifetime sexual IPV. Women were more likely to report having experienced lifetime sexual IPV if they were: aged 25 and over; cohabiting, divorced or separated, or widowed; or living in the most deprived areas. Those who had some tertiary education were less likely to report lifetime experience of sexual IPV compared with those with primary or secondary education (table 3).

Table 3

Characteristics of women with lifetime and past 12 month Sexual IPV in the pooled database from two cross-sectional studies in New Zealand

Past 12-month sexual IPV: those who were divorced/separated were more likely to report having experienced 12-month sexual IPV compared with married women. Those living in the most deprived areas were also more likely to report 12-month sexual IPV. Women aged 55 years and above were less likely to report having experienced sexual IPV in the past 12 months compared with younger women (table 3).

No significant interaction was found between study year and sociodemographic factors (data not shown).

Changes in women’s attitudes

In 2003, 48.1% agreed with at least one of the statements indicating agreement with traditional gender roles, compared with 38.4% in 2019. While not common in 2003, it was even less common in 2019 for women to agree with the justifications for a man to hit his wife if he finds out she has been unfaithful (3.8% agreement in 2003, 1.8% agreement in 2019; table 4).

Table 4

Prevalence rates and changes in women’s attitudes towards traditional gender roles in relationships and attitudes towards acceptability of a man hitting his wife

Changes in help-seeking behaviours

There was an overall reduction in the proportion of women who had sought help from formal or informal sources, with three-quarters (77%) of women who had experienced IPV reporting that they had told someone about the violence in 2003 compared with 70% in 2019. This reduction appears to be driven by the significant reduction in the proportion of women who sought help from informal sources (from 71.3% in 2003 to 64.6% in 2019). There was no change in the proportion of women who sought help from ‘formal’ sources between the two study years. Similarly, there was no significant change in the proportion of women who reported that they received help from formal sources (table 5).

Table 5

Prevalence rates and changes in help sought and help received between 2003 and 2019 by those who reported at least one type of sexual or physical IPV

Discussion

Changes in prevalence of physical and sexual IPV between 2003 and 2019 were explored using two population-based surveys. Our findings indicated that the lifetime prevalence of physical IPV remained relatively unchanged between 2003 and 2019, with almost one-third (30%) of women in both surveys reporting having experienced at least one act of physical IPV in their lifetime. This is similar to reported prevalence rates from the EU 28-countries study (33%),28 and the USA (30.6%),29 and is comparable to the global average.30 While lifetime prevalence of physical IPV was unchanged, there was a significant decrease in the proportion of women who reported experiencing 12-month physical IPV. Small reductions in rates for lifetime sexual IPV were also observed. Population changes in sociodemographic characteristics did not fully explain the decreases in IPV prevalence over time, and the noted changes were consistent across sub-groups of the population.

In 2003, 48.1% of women agreed with one or more of the statements supportive of traditional gender roles, compared with 38.4% in 2019. These were low percentages of agreement compared with women in low-income and middle-income countries.31–33 Agreement with attitudes supportive of justifications for a man hitting his wife was low in both the 2003 (0.2%–3.8%) and 2019 surveys (0.2%–2.3%), and extremely low compared with results reported from low-and middle-income countries.34 35 but comparable with high income countries.36 Even with this low rate of agreement, change was still observed, with a significant reduction in agreement with the statement that ‘it is acceptable for a man to hit his wife if he found out she was unfaithful’, from 3.8% in 2003 to 1.8% in 2019.

Overall, among women who experienced IPV, the rates of disclosure (telling someone about the violence) were high (77% in 2003, 70% in 2019), compared with findings from low-income and middle-income countries,37 38 and comparable with high-income countries.39 It should be noted, however, that most disclosures were made to informal sources, such as family or friends. There was no change in ‘help received’ from formal sources (21.1% in 2003, 19.4% in 2019). This warrants further attention, to determine if this is due to limited service capacity, or limits in the quality of help currently available.

Possible explanations for the study findings include: actual changes in perpetrator behaviour over time; or changes due to differences in methods, measurement or samples.

There is some evidence that changes in perpetrator behaviour may have occurred, as the reduction in the 12-month prevalence of physical and lifetime sexual IPV between 2003 and 2019 is consistent with a reduction in 12-month prevalence of psychological IPV noted in the same sample.40 Changes in perpetrator behaviour are possible, as there have been a series of strategies and campaigns implemented between the two study years. These included: changes in legislation (eg, amendments to family violence law), and the introduction of prevention campaigns and programmes (eg, the Family Violence: It’s not ok national campaign,41 and the Accident Compensation Corporation-funded Mates and Dates high schools programmes on healthy relationships42). These actions may have contributed to changes in societal awareness and understandings of attitudes supportive of violence against women as there is some evidence that these initiatives had wide population reach.40 This interpretation is supported by our findings on the reduction in women’s agreement with attitudes towards traditional gender roles and reduction in women’s agreement with the acceptability of a man hitting his wife if she was unfaithful. Other studies have also noted the relationship between attitudes to violence and victimisation.43 44

An additional feature of these societal actions was the call for those experiencing violence to reach out for help.41 Our findings suggest that there has been no change in women contacting formal sources of help, and a small but significant reduction in talking with informal sources. As help seeking can be related to the severity of violence experienced, it is possible that the lack of change in accessing formal help among women is related to the reduction of current physical, and lifetime sexual IPV between the studied years and a possible decrease of high severity cases. However, it is also possible that activities designed to encourage community engagement in violence prevention may need additional resourcing to ensure a sustained response and appropriate access to necessary services. Further research with larger sample sizes will be important to verify this finding.

The alternate explanation of the observed changes being due to differences in study methods or sample difference seem less likely. Specifically, the comparability of methods across the two surveys, including use of identical questions in the two survey waves, lends strength to the interpretation that the prevalence changes noted are real. Additionally, while there were some differences in the characteristics of the two samples, the AOR showed that after controlling for all sociodemographic factors, the observed differences in prevalence still remained significant.

The observed reduction in 12-month prevalence of physical IPV is positive, and parallels overall reductions in crime rates reported by crime and victimisation surveys,45 and is similar to reductions in prevalence of IPV documented in Australia between 1996 and 2005.46 It may be the result of more women recognising abusive behaviour and taking their own actions to leave abusive relationships. However, further efforts and investment are needed to ensure that those who ask for help actually receive help. Importantly, the stability of the lifetime prevalence of physical IPV should heighten efforts to develop and implement comprehensive and sustained prevention work with those who use violence in relationships.

Strengths

Strengths include: the representativeness of the samples obtained, and the use of comparable methods and comparable questions across the two survey waves. Additionally, the 15-year time gap between the two survey waves is sufficient to determine if real change occurred.12

Limitations and recommendations for future studies

Changes between two time points are not sufficient to determine if the change represents a trend, so caution is needed when interpreting the changes observed. Overall, the prevalence estimates obtained may under-report what is happening in the population as a whole, either because of stigma,47 or because of the overall response rate for the study. While we successfully surveyed over 63% of eligible women, those with greater levels of exposure to violence may be less likely to have participated. Future studies would benefit from larger sample sizes, which would improve the chance of detecting real changes in low base rate phenomena, such as 12-month prevalence of sexual IPV.

Conclusion

The observed reduction in 12-month physical and lifetime sexual IPV prevalence rates, changes in attitudes about the acceptability of violence, and the increases in help seeking are positive. However, work is still needed to address the substantial problem of IPV, as the lifetime prevalence rate of one in three women experiencing IPV remained stable over the 15-year time interval. This means that prevention efforts must be increased and sustained, and that adequate structures and resources must be available to respond to those seeking help.

Acknowledgments

The authors gratefully acknowledge participants, the interviewers and the study project team, led by Patricia Meagher-Lundberg. They also acknowledge the representatives from the Ministry of Justice, the Accident Compensation Corporation, the New Zealand Police, and the Ministry of Education, who were part of the Governance Group for Family and Sexual Violence at the inception of the study. The two cross-sectional studies from which this study used data are based on the WHO Violence Against Women Instrument as developed for use in the WHO Multi-Country Study on Women’s Health and Domestic Violence and has been adapted from the version used in Asia and the Pacific by kNOwVAWdata (V.12.03).

References

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.

Footnotes

  • Contributors JF and PG contributed to the conception and design of the study. TM contributed to the application for funding of 2019 study. LH managed the data cleaning, and conducted the analyses, with contributions from ZM. LH, JF, ZM and PG interpreted the data, drafted the article and revised it. All authors contributed to the manuscript and approved the final version.

  • Funding This work was supported by the Health Research Council of New Zealand (Grant 02/207) for the 2003 study and the New Zealand Ministry of Business, Innovation and Employment, Contract number CONT-42799-HASTR-UOA for the 2019 study.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval Ethics approval was granted by the University of Auckland Human Participants Ethics Committee in 2003 (Ref number: 2002/199) and 2019 (Reference number 2015/ 018244).

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

  • Data availability statement No data are available. Data are unavailable due to the confidentiality and sensitivity of the data and Māori data sovereignty.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.

Linked Articles