<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="research-article"><?properties manuscript?><front><journal-meta><journal-id journal-id-type="nlm-journal-id">9009644</journal-id><journal-id journal-id-type="pubmed-jr-id">1090</journal-id><journal-id journal-id-type="nlm-ta">Epidemiology</journal-id><journal-id journal-id-type="iso-abbrev">Epidemiology</journal-id><journal-title-group><journal-title>Epidemiology (Cambridge, Mass.)</journal-title></journal-title-group><issn pub-type="ppub">1044-3983</issn><issn pub-type="epub">1531-5487</issn></journal-meta><article-meta><article-id pub-id-type="pmid">29746368</article-id><article-id pub-id-type="pmc">6126532</article-id><article-id pub-id-type="doi">10.1097/EDE.0000000000000848</article-id><article-id pub-id-type="manuscript">NIHMS964991</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Contrasting causal effects of workplace interventions</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Izano</surname><given-names>Monika A.</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Brown</surname><given-names>Daniel M.</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Neophytou</surname><given-names>Andreas M.</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Garcia</surname><given-names>Erika</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib><contrib contrib-type="author"><name><surname>Eisen</surname><given-names>Ellen A.</given-names></name><xref ref-type="aff" rid="A1">1</xref></contrib></contrib-group><aff id="A1">
<label>1</label>School of Public Health, University of California, Berkeley, CA</aff><author-notes><corresp id="FN1">Corresponding Author: Ellen A. Eisen, Sc.D.; Professor and Head; Division of Environmental Health Sciences; School of Public Health; University of California, Berkeley; 50 University Hall #7360; Berkeley, CA. 94720; Phone: (510) 643-5310; <email>eeisen@berkeley.edu</email></corresp></author-notes><pub-date pub-type="nihms-submitted"><day>18</day><month>5</month><year>2018</year></pub-date><pub-date pub-type="ppub"><month>7</month><year>2018</year></pub-date><pub-date pub-type="pmc-release"><day>01</day><month>7</month><year>2019</year></pub-date><volume>29</volume><issue>4</issue><fpage>542</fpage><lpage>546</lpage><!--elocation-id from pubmed: 10.1097/EDE.0000000000000848--><abstract><p id="P1">Occupational exposure guidelines are ideally based on estimated effects of static interventions that assign constant exposure over a working lifetime. Static effects are difficult to estimate when follow-up extends beyond employment because their identifiability requires additional assumptions. Effects of dynamic interventions that assign exposure while at work, allowing subjects to leave and become unexposed thereafter, are more easily identifiable but result in different estimates. Given the practical implications of exposure limits, we explored the drivers of the differences between static and dynamic interventions in a simulation study where workers could terminate employment due to an intermediate adverse health event that functions as a time-varying confounder. The two effect estimates became more similar with increasing strength of the health event and outcome relationship, and with increasing time between health event and employment termination. Estimates were most dissimilar when the intermediate health event occurred early in employment, providing an effective screening mechanism.</p></abstract><kwd-group><kwd>Healthy-worker survivor effect</kwd><kwd>workplace interventions</kwd><kwd>dynamic interventions</kwd><kwd>static interventions</kwd></kwd-group></article-meta></front><body><sec sec-type="intro" id="S1"><title>INTRODUCTION</title><p id="P2">Assessments of workplace risk are generally based on observational studies of occupational cohorts. Estimates from these studies are often subject to bias due to the healthy worker survivor effect (HWSE), a ubiquitous process that results in the healthiest workers accruing the most exposure over their lifetimes.<sup><xref rid="R1" ref-type="bibr">1</xref>&#x02013;<xref rid="R4" ref-type="bibr">4</xref></sup> The potential outcomes framework defines causal effects as contrasts between the distributions of counterfactual outcomes under hypothetical interventions.<sup><xref rid="R5" ref-type="bibr">5</xref></sup> Counterfactuals may be defined under static regimens that assign exposure independently of a worker&#x02019;s characteristics, or under dynamic regimens that assign exposure according to a worker&#x02019;s observed past.</p><p id="P3">A static workplace intervention corresponds to a target trial<sup><xref rid="R6" ref-type="bibr">6</xref>,<xref rid="R7" ref-type="bibr">7</xref></sup> in which workers hired into an industry at the start of their working lives are assigned a fixed exposure level until retirement age. The effects of static interventions reflect the biologic effect of exposure that would have been observed if workers did not select out of the workforce. Since workers that experience exposure-related adverse health outcomes tend to leave work, such strategies do not correspond to anything commonly observed in the real world. Given the paucity of observed data to estimate such effects, most estimates necessarily rely on the assumption of no unmeasured confounding between leaving work and the outcome.<sup><xref rid="R8" ref-type="bibr">8</xref></sup></p><p id="P4">A dynamic intervention also corresponds to a target trial in which workers hired at the start of their working lives are assigned a fixed exposure level while at work; however, workers may leave and become unexposed for the remainder of follow-up. These strategies more closely resemble the real-world observed data, and their effects are more straightforward to estimate.</p><p id="P5">The US Occupational Safety and Health Administration (OSHA) is mandated to establish standards for occupational hazards that assure, &#x0201c;on the basis of the best available evidence that no employee will suffer material impairment of health or functional capacity, even if such employee has a regular exposure to the hazard dealt with by such standard for the period of his working life.&#x0201d;<sup><xref rid="R9" ref-type="bibr">9</xref></sup> This directive is typically interpreted as requiring that OSHA risk assessments be based on static effects of 45 years of exposure.<sup><xref rid="R10" ref-type="bibr">10</xref></sup> Since static effects are difficult to estimate <sup><xref rid="R11" ref-type="bibr">11</xref></sup>, we conducted a simulation study to investigate the factors that drive differences between the static and dynamic measures in workplace studies.</p></sec><sec sec-type="methods" id="S2"><title>METHODS</title><p id="P6">Motivated by typical aspects of occupational studies, we simulated 500 cohorts of 50,000 workers, under four scenarios, for a maximum follow-up of 20 years. For each worker and year <italic>t</italic> of follow-up, we simulated the indicators: <italic>W</italic>(<italic>t</italic>), for active employment; <italic>E</italic>(<italic>t</italic>) for workplace exposure; <italic>Y</italic>(<italic>t</italic>), for the outcome; and <italic>H</italic>(<italic>t</italic>), for the intermediate adverse health event, one of the mechanisms through which exposure affects the outcome. Exposure had a direct effect on the outcome in addition to that mediated by <italic>H</italic>(<italic>t</italic>). We additionally simulated <italic>S</italic>, an indicator of a worker&#x02019;s susceptibility to <italic>H</italic>(<italic>t</italic>) and <italic>Y</italic>(<italic>t</italic>). Only susceptible workers were at risk of experiencing <italic>H</italic>(<italic>t</italic>) and <italic>Y</italic>(<italic>t</italic>) and, once they experienced <italic>H</italic>(<italic>t</italic>), workers maintained that status. Workers could only terminate employment if in poor health.</p><sec id="S3"><title>Defining the Interventions</title><p id="P7">We generated outcomes by setting the nodes <italic>E</italic>(<italic>t</italic>) and <italic>W</italic>(<italic>t</italic>) in the data generating system sequentially for each worker and year of follow-up. Our static interventions {<italic>d<sub>s</sub></italic><sub>,1</sub>, <italic>d<sub>s</sub></italic><sub>,0</sub>} set exposure <italic>E</italic>(<italic>t</italic>) to 1 or 0 respectively, and <italic>W</italic>(<italic>t</italic>) to 1 for all years. The counterfactual outcomes <italic>Y</italic><sub><italic>i,d</italic><sub><italic>s</italic>,1</sub></sub>(<italic>t</italic>) and <italic>Y</italic><sub><italic>i,d</italic><sub><italic>s</italic>,0</sub></sub>(<italic>t</italic>) denote the outcomes that worker <italic>i</italic> would experience in year <italic>t</italic> if they were always at work, and respectively always exposed or unexposed.</p><p id="P8">Our dynamic interventions assigned exposure according to employment status. Intervention <italic>d<sub>d</sub></italic><sub>,1</sub> set exposure to 1 while a worker was actively employed and to 0 once the worker terminated employment. Intervention <italic>d<sub>d</sub></italic><sub>,0</sub> assigned workers to no exposure before and after termination. The counterfactual outcomes <italic>Y</italic><sub><italic>i,d</italic><sub><italic>d</italic>,1</sub></sub>(<italic>t</italic>) and <italic>Y</italic><sub><italic>i,d</italic><sub><italic>d</italic>,0</sub></sub>(<italic>t</italic>) correspond to worker <italic>t</italic>&#x02019;s outcome in year <italic>t</italic> if they were always exposed and unexposed while at work, respectively.</p><p id="P9">The causal relationships between the variables under the two interventions are presented in the directed acyclic graph (DAG) in <xref rid="F1" ref-type="fig">Figure 1</xref>.</p></sec><sec id="S4"><title>Simulation Scenarios</title><p id="P10">The <xref rid="SD1" ref-type="supplementary-material">eAppendix</xref> includes the data-generating equations, the <italic>R</italic> code used to create the datasets and results, and the distribution of covariates under intervention <italic>d<sub>d</sub></italic><sub>,1</sub> (<xref rid="SD1" ref-type="supplementary-material">eTable 1</xref>), for each scenario. In Scenario 1, <italic>H</italic>(<italic>t</italic>) was predicted by cumulative exposure, and workers who experienced <italic>H</italic>(<italic>t</italic>) terminated employment in the following year. <italic>H</italic>(<italic>t</italic>) had a moderate effect on the probability of developing the outcome, as might be the case with asthma [<italic>H</italic>(<italic>t</italic>)] and chronic obstructive pulmonary disease [COPD, <italic>Y</italic>(<italic>t</italic>)]. Scenarios 2, 3, and 4 were compared to Scenario 1.</p><p id="P11">In Scenario 2, <italic>H</italic>(<italic>t</italic>) was a stronger predictor of the outcome than in Scenario 1, as might be the case with abnormal lung function and COPD.</p><p id="P12">In Scenario 3 we altered Scenario 1 by increasing the time between <italic>H</italic>(<italic>t</italic>) and leaving work. For example, years may pass before an asthmatic worker terminates employment due to worsening of symptoms. To reflect this temporal relationship, employment status was predicted by adverse health events occurring 10 years prior.</p><p id="P13">Finally, in Scenario 4 we altered the temporal relationship between the exposure and <italic>H</italic>(<italic>t</italic>). Here <italic>H</italic>(<italic>t</italic>) is an acute response that occurs in the first year of exposure, as may be the case if workers exposed to an irritant immediately experience severe asthma.</p></sec><sec id="S5"><title>Comparing static and dynamic interventions</title><p id="P14">For each scenario, dataset, and worker, we generated four counterfactual outcomes: exposed and unexposed for both the static and dynamic interventions. The reported scenario-specific survival (<xref rid="F2" ref-type="fig">Figures 2</xref>&#x02013;<xref rid="F4" ref-type="fig">4</xref>) and cumulative incidence (<xref rid="T1" ref-type="table">Table</xref>) estimates under each exposure level and intervention were averages across the 500 datasets. Static and dynamic exposure effects were measured by the risk differences (<italic>RD</italic>) contrasting the cumulative incidence of disease between exposed and unexposed workers. The ratio (<italic>R</italic>) of static and dynamic RDs represents the factor by which the risk among the exposed was reduced by early employment termination.</p><p id="P15">All datasets were simulated using the <italic>simcausal</italic> R package.<sup><xref rid="R12" ref-type="bibr">12</xref></sup> All analyses were performed in the R programming language, version 3.3.2.<sup><xref rid="R13" ref-type="bibr">13</xref></sup></p></sec></sec><sec sec-type="results" id="S6"><title>RESULTS</title><p id="P16">In the <xref rid="T1" ref-type="table">Table</xref> we report static and dynamic exposure effects for each scenario, and their ratios over time. In Scenario 1, a static effect of 0.42 indicates that the 20-year risk of disease if always exposed would be 42% greater than if unexposed. A dynamic effect of 0.07 indicates that, when workers terminated employment for health-related reasons, the 20-year risk was only 7% greater if exposed while at work versus unexposed. The corresponding static-to-dynamic effect ratio of 6.00 in year 20 [<italic>R</italic><sub>1</sub>(20) = 6.00] indicates that the effect of exposure under a static intervention was six times stronger than under a dynamic intervention.</p><p id="P17"><xref rid="F2" ref-type="fig">Figures 2</xref> &#x02013; <xref rid="F4" ref-type="fig">4</xref> present comparisons of Scenarios 2, 3, and 4 with Scenario 1. For each scenario, we present survival curves for exposed and unexposed workers according to static and dynamic interventions. Survival probabilities were identical for the &#x0201c;never expose&#x0201d; regimen under both static and dynamic interventions; survival probabilities were smaller under the static &#x0201c;always expose&#x0201d; versus the dynamic &#x0201c;expose while at work&#x0201d; interventions (<xref rid="F2" ref-type="fig">Figures 2</xref>&#x02013;<xref rid="F4" ref-type="fig">4</xref>).</p><p id="P18">The first comparison assesses how the static-to-dynamic effect ratio changes as a function of the strength of the intermediate health events and outcome relationship (<xref rid="F2" ref-type="fig">Figure 2</xref>). Static and dynamic effects were more alike when intermediate health had a stronger effect on the outcome (Scenario 2), reflected by smaller ratios in Scenario 2 than in 1 (<xref rid="T1" ref-type="table">Table</xref>, <italic>R</italic><sub>1</sub>(20) = 6.00, <italic>R</italic><sub>2</sub>(20) = 2.00).</p><p id="P19">In <xref rid="F3" ref-type="fig">Figure 3</xref> we contrast Scenarios 1 and 3. Leaving work was determined by health status in the previous year in Scenario 1, but 10 years prior in Scenario 3. Static effects were the same in both scenarios (0.42). However, because exposed workers remained at work longer in Scenario 3, the dynamic effect was larger (0.07 vs. 0.32) and the ratio was smaller (<italic>R</italic><sub>1</sub>(20) = 6.00, <italic>R</italic><sub>3</sub>(20) = 1.31).</p><p id="P20">In <xref rid="F4" ref-type="fig">Figure 4</xref> we contrast Scenarios 3 and 4, where workers only experienced the adverse health event during their very first year of exposure. While static effects were similar in both scenarios (0.42 vs. 0.43), dynamic effects were smaller in Scenario 4 (0.07 vs. 0.05) and the ratio increased [<italic>R</italic><sub>1</sub>(20) = 6.00, <italic>R</italic><sub>4</sub>(20) = 8.60].</p></sec><sec sec-type="discussion" id="S7"><title>DISCUSSION</title><p id="P21">Static and dynamic interventions are expected to result in different exposure effect estimates.<sup><xref rid="R11" ref-type="bibr">11</xref></sup> We aimed to examine key features of the data generating distribution likely to drive the differences between static and dynamic intervention effects, since both can answer causal questions. This matter has practical implications in the context of the OSHA mandate: It is important to note that not all causal effects of the same exposure level, on the same disease, in the same population, over the same follow-up period, are necessarily equivalent.</p><p id="P22">Dynamic parameters better approximate static estimates if the consequences of exposure that cause leaving work are strongly related to the outcome, or if symptomatic workers remain at work despite declining health. In both cases, the amount of exposure accrued while actively employed already places workers at high risk for the outcome. Dynamic parameters generally provide conservative estimates of the disease risk associated with long-term exposure when the consequences of exposure that lead to leaving work occur early in employment. In this case, leaving work provides a screening mechanism for susceptible workers, while a larger proportion of unsusceptible survivors remain at work (<xref rid="SD1" ref-type="supplementary-material">eTable 1</xref>).</p></sec><sec sec-type="supplementary-material" id="S8"><title>Supplementary Material</title><supplementary-material content-type="local-data" id="SD1"><label>Supplemental Digital Content_ eAppendix</label><media xlink:href="NIHMS964991-supplement-Supplemental_Digital_Content__eAppendix.docx" orientation="portrait" xlink:type="simple" id="d36e582" position="anchor"/></supplementary-material></sec></body><back><ack id="S9"><p><bold>Funding Sources:</bold> Preparation of this manuscript was supported by grants # 1R01OH010228 (P.I. Ellen A. Eisen), and #T42OH008429 (M.A.I.), by the National Institute for Occupational Safety and Health (NIOSH), Centers for Disease Control and Prevention, Department of Health and Human Services.</p></ack><fn-group><fn fn-type="COI-statement" id="FN2"><p><bold>Conflict of Interests:</bold> The authors have no conflict of interest to disclose.</p></fn><fn id="FN3"><p><bold>Data and Code Availability:</bold> The code for simulating datasets and conducting analyses is provided in the <xref rid="SD1" ref-type="supplementary-material">eAppendix</xref>.</p></fn></fn-group><glossary id="GL"><title>Abbreviations</title><def-list><def-item><term id="G1">HWSE</term><def><p>Healthy-worker survivor effect</p></def></def-item><def-item><term id="G2">RD</term><def><p>risk difference</p></def></def-item></def-list></glossary><ref-list><ref id="R1"><label>1</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Fox</surname><given-names>AJ</given-names></name><name><surname>Collier</surname><given-names>PF</given-names></name></person-group><article-title>Low mortality rates in industrial cohort studies due to selection for work and survival in the industry</article-title><source>Br J Prev Soc Med</source><volume>30</volume><fpage>225</fpage><lpage>30</lpage><year>1976</year><pub-id pub-id-type="pmid">1009272</pub-id></element-citation></ref><ref id="R2"><label>2</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gilbert</surname><given-names>ES</given-names></name></person-group><article-title>Some confounding factors in the study of mortality and occupational exposures</article-title><source>Am J Epidemiol</source><volume>116</volume><fpage>177</fpage><lpage>88</lpage><year>1982</year><pub-id pub-id-type="pmid">7102651</pub-id></element-citation></ref><ref id="R3"><label>3</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Monson</surname><given-names>RR</given-names></name></person-group><article-title>Observations on the healthy worker effect</article-title><source>J Occup Med</source><volume>28</volume><fpage>425</fpage><lpage>433</lpage><year>1986</year><pub-id pub-id-type="pmid">3723215</pub-id></element-citation></ref><ref id="R4"><label>4</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Arrighi</surname><given-names>HM</given-names></name><name><surname>Hertz-Picciotto</surname><given-names>I</given-names></name></person-group><article-title>The Evolving Concept of the Healthy Worker Survivor Effect</article-title><source>Epidemiology</source><volume>5</volume><fpage>189</fpage><lpage>196</lpage><year>1994</year><pub-id pub-id-type="pmid">8172994</pub-id></element-citation></ref><ref id="R5"><label>5</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rubin</surname><given-names>DB</given-names></name></person-group><article-title>Estimating causal effects of treatments in randomized and nonrandomized studies</article-title><source>J Educ Psychol</source><volume>66</volume><fpage>688</fpage><year>1974</year></element-citation></ref><ref id="R6"><label>6</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hern&#x000e1;n</surname><given-names>MA</given-names></name><etal/></person-group><article-title>Observational studies analyzed like randomized experiments: an application to postmenopausal hormone therapy and coronary heart disease</article-title><source>Epidemiology</source><volume>19</volume><fpage>766</fpage><year>2008</year><pub-id pub-id-type="pmid">18854702</pub-id></element-citation></ref><ref id="R7"><label>7</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hern&#x000e1;n</surname><given-names>MA</given-names></name><name><surname>Robins</surname><given-names>JM</given-names></name></person-group><article-title>Using big data to emulate a target trial when a randomized trial is not available</article-title><source>Am J Epidemiol</source><volume>183</volume><fpage>758</fpage><lpage>764</lpage><year>2016</year><pub-id pub-id-type="pmid">26994063</pub-id></element-citation></ref><ref id="R8"><label>8</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>DM</given-names></name><etal/></person-group><article-title>The healthy worker survivor effect: target parameters and target populations</article-title><source>Curr Environ Heal Reports</source><fpage>1</fpage><lpage>9</lpage><year>2017</year></element-citation></ref><ref id="R9"><label>9</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Act</surname><given-names>WS</given-names></name></person-group><article-title>Occupational Safety and Health Act of 1970</article-title><source>Public Law</source><volume>91</volume><fpage>596</fpage><year>1970</year></element-citation></ref><ref id="R10"><label>10</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rice</surname><given-names>FL</given-names></name><etal/></person-group><article-title>Crystalline silica exposure and lung cancer mortality in diatomaceous earth industry workers: a quantitative risk assessment</article-title><source>Occup Environ Med</source><volume>58</volume><fpage>38</fpage><lpage>45</lpage><year>2001</year><pub-id pub-id-type="pmid">11119633</pub-id></element-citation></ref><ref id="R11"><label>11</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Petersen</surname><given-names>M</given-names></name><etal/></person-group><article-title>Targeted maximum likelihood estimation for dynamic and static longitudinal marginal structural working models</article-title><source>J Causal Inference</source><volume>2</volume><fpage>147</fpage><lpage>185</lpage><year>2014</year><pub-id pub-id-type="pmid">25909047</pub-id></element-citation></ref><ref id="R12"><label>12</label><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Sofrygin</surname><given-names>O</given-names></name><name><surname>van der Laan</surname><given-names>MJ</given-names></name><name><surname>Neugebauer</surname><given-names>R</given-names></name></person-group><source>{simcausal}: Simulating Longitudinal Data with Causal Inference Applications</source><year>2015</year></element-citation></ref><ref id="R13"><label>13</label><element-citation publication-type="book"><collab>Team, R. C</collab><source>R: A language and environment for statistical computing</source><publisher-name>R Foundation for Statistical Computing</publisher-name><publisher-loc>Vienna, Austria</publisher-loc><year>2013</year><comment>2014</comment></element-citation></ref></ref-list></back><floats-group><fig id="F1" orientation="portrait" position="float"><label>Figure 1</label><caption><title>Directed acyclic graphs (DAGs) representing the relationships between variables in scenarios 1 under the static (a) and dynamic (b) regimens of interest</title><p><italic>E</italic>(<italic>t</italic>), <italic>H</italic>(<italic>t</italic>), <italic>W</italic>(<italic>t</italic>), <italic>Y</italic>(<italic>t</italic>) denote workplace exposure, intermediate adverse health event, active employment status, and the outcome in year <italic>t. S</italic> is a baseline indicator of a worker&#x02019;s susceptibility to <italic>H</italic>(<italic>t</italic>) and <italic>Y</italic>(<italic>t</italic>)that modifies the effect of the exposure. The active employment (<italic>W</italic>(2)) and exposure (<italic>E</italic>(1), <italic>E</italic>(2)) nodes are set rather than predicted by the past under the static regimen. Under the dynamic intervention, <italic>H</italic>(1) predicts subsequent employment status (<italic>W</italic>(2)), which determines future exposure (<italic>E</italic>(2)).</p></caption><graphic xlink:href="nihms964991f1"/></fig><fig id="F2" orientation="portrait" position="float"><label>Figure 2</label><caption><title>Evaluating the role of the intermediate health event - outcome relationship</title><p>Counterfactual survival curves among cohorts of exposed (dashed line) and unexposed (dotted line) workers following static (grey line) and dynamic (black line) interventions under Scenarios 1 and 2.</p></caption><graphic xlink:href="nihms964991f2"/></fig><fig id="F3" orientation="portrait" position="float"><label>Figure 3</label><caption><title>Evaluating the role of the temporal relationship between the intermediate health event and leaving work</title><p>Counterfactual survival curves among cohorts of exposed (dashed line) and unexposed (dotted line) workers following static (grey line) and dynamic (black line) interventions under Scenarios 1 and 3.</p></caption><graphic xlink:href="nihms964991f3"/></fig><fig id="F4" orientation="portrait" position="float"><label>Figure 4</label><caption><title>Evaluating the role of the temporal relationship between exposure and the intermediate health event</title><p>Counterfactual survival curves among cohorts of exposed (dashed line) and unexposed (dotted line) following static (grey line) and dynamic (black line) interventions under Scenarios 1 and 4.</p></caption><graphic xlink:href="nihms964991f4"/></fig><table-wrap id="T1" position="float" orientation="landscape"><label>Table</label><caption><p>Risk difference (RD) estimates for static and dynamic interventions, as well as the ratios (<italic>R</italic>) of static and dynamic risk differences, over time (<italic>t</italic>).</p></caption><table frame="hsides" rules="groups"><thead><tr><th valign="bottom" rowspan="3" align="left" colspan="1"/><th colspan="10" valign="bottom" align="center" rowspan="1">Year of follow-up</th></tr><tr><th colspan="10" valign="bottom" align="center" rowspan="1">
<hr/></th></tr><tr><th valign="bottom" align="center" rowspan="1" colspan="1">t =2</th><th valign="bottom" align="center" rowspan="1" colspan="1">t =4</th><th valign="bottom" align="center" rowspan="1" colspan="1">t =6</th><th valign="bottom" align="center" rowspan="1" colspan="1">t =8</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 10</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 12</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 14</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 16</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 18</th><th valign="bottom" align="center" rowspan="1" colspan="1">t = 20</th></tr></thead><tbody><tr><td colspan="11" align="left" valign="top" rowspan="1"><bold>Scenario 1</bold><xref rid="TFN1" ref-type="table-fn">a</xref></td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Static RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td><td align="center" valign="top" rowspan="1" colspan="1">0.12</td><td align="center" valign="top" rowspan="1" colspan="1">0.19</td><td align="center" valign="top" rowspan="1" colspan="1">0.28</td><td align="center" valign="top" rowspan="1" colspan="1">0.36</td><td align="center" valign="top" rowspan="1" colspan="1">0.42</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Dynamic RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.02</td><td align="center" valign="top" rowspan="1" colspan="1">0.03</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.05</td><td align="center" valign="top" rowspan="1" colspan="1">0.06</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1"><italic>R</italic><sub>1</sub>(<italic>t</italic>)</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">2.00</td><td align="center" valign="top" rowspan="1" colspan="1">2.33</td><td align="center" valign="top" rowspan="1" colspan="1">3.00</td><td align="center" valign="top" rowspan="1" colspan="1">3.80</td><td align="center" valign="top" rowspan="1" colspan="1">4.67</td><td align="center" valign="top" rowspan="1" colspan="1">5.14</td><td align="center" valign="top" rowspan="1" colspan="1">6.00</td></tr><tr><td colspan="11" align="left" valign="top" rowspan="1"><bold>Scenario 2</bold><xref rid="TFN2" ref-type="table-fn">b</xref></td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Static RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.02</td><td align="center" valign="top" rowspan="1" colspan="1">0.06</td><td align="center" valign="top" rowspan="1" colspan="1">0.13</td><td align="center" valign="top" rowspan="1" colspan="1">0.22</td><td align="center" valign="top" rowspan="1" colspan="1">0.31</td><td align="center" valign="top" rowspan="1" colspan="1">0.38</td><td align="center" valign="top" rowspan="1" colspan="1">0.41</td><td align="center" valign="top" rowspan="1" colspan="1">0.41</td><td align="center" valign="top" rowspan="1" colspan="1">0.40</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Dynamic RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td><td align="center" valign="top" rowspan="1" colspan="1">0.11</td><td align="center" valign="top" rowspan="1" colspan="1">0.14</td><td align="center" valign="top" rowspan="1" colspan="1">0.16</td><td align="center" valign="top" rowspan="1" colspan="1">0.18</td><td align="center" valign="top" rowspan="1" colspan="1">0.19</td><td align="center" valign="top" rowspan="1" colspan="1">0.20</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1"><italic>R</italic><sub>2</sub>(<italic>t</italic>)</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">2.00</td><td align="center" valign="top" rowspan="1" colspan="1">1.50</td><td align="center" valign="top" rowspan="1" colspan="1">1.86</td><td align="center" valign="top" rowspan="1" colspan="1">2.00</td><td align="center" valign="top" rowspan="1" colspan="1">2.21</td><td align="center" valign="top" rowspan="1" colspan="1">2.38</td><td align="center" valign="top" rowspan="1" colspan="1">2.28</td><td align="center" valign="top" rowspan="1" colspan="1">2.16</td><td align="center" valign="top" rowspan="1" colspan="1">2.00</td></tr><tr><td colspan="11" align="left" valign="top" rowspan="1"><bold>Scenario 3</bold><xref rid="TFN3" ref-type="table-fn">c</xref></td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Static RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td><td align="center" valign="top" rowspan="1" colspan="1">0.12</td><td align="center" valign="top" rowspan="1" colspan="1">0.19</td><td align="center" valign="top" rowspan="1" colspan="1">0.28</td><td align="center" valign="top" rowspan="1" colspan="1">0.36</td><td align="center" valign="top" rowspan="1" colspan="1">0.42</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Dynamic RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.07</td><td align="center" valign="top" rowspan="1" colspan="1">0.12</td><td align="center" valign="top" rowspan="1" colspan="1">0.18</td><td align="center" valign="top" rowspan="1" colspan="1">0.24</td><td align="center" valign="top" rowspan="1" colspan="1">0.29</td><td align="center" valign="top" rowspan="1" colspan="1">0.32</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1"><italic>R</italic><sub>3</sub>(<italic>t</italic>)</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">1.06</td><td align="center" valign="top" rowspan="1" colspan="1">1.17</td><td align="center" valign="top" rowspan="1" colspan="1">1.24</td><td align="center" valign="top" rowspan="1" colspan="1">1.31</td></tr><tr><td colspan="11" align="left" valign="top" rowspan="1"><bold>Scenario 4</bold><xref rid="TFN4" ref-type="table-fn">d</xref></td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Static RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.03</td><td align="center" valign="top" rowspan="1" colspan="1">0.05</td><td align="center" valign="top" rowspan="1" colspan="1">0.08</td><td align="center" valign="top" rowspan="1" colspan="1">0.13</td><td align="center" valign="top" rowspan="1" colspan="1">0.20</td><td align="center" valign="top" rowspan="1" colspan="1">0.29</td><td align="center" valign="top" rowspan="1" colspan="1">0.37</td><td align="center" valign="top" rowspan="1" colspan="1">0.43</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1">Dynamic RD</td><td align="center" valign="top" rowspan="1" colspan="1">0.00</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.01</td><td align="center" valign="top" rowspan="1" colspan="1">0.02</td><td align="center" valign="top" rowspan="1" colspan="1">0.02</td><td align="center" valign="top" rowspan="1" colspan="1">0.03</td><td align="center" valign="top" rowspan="1" colspan="1">0.03</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.04</td><td align="center" valign="top" rowspan="1" colspan="1">0.05</td></tr><tr><td align="left" valign="top" rowspan="1" colspan="1"><italic>R</italic><sub>4</sub>(<italic>t</italic>)</td><td align="center" valign="top" rowspan="1" colspan="1">--</td><td align="center" valign="top" rowspan="1" colspan="1">1.00</td><td align="center" valign="top" rowspan="1" colspan="1">3.00</td><td align="center" valign="top" rowspan="1" colspan="1">2.50</td><td align="center" valign="top" rowspan="1" colspan="1">4.00</td><td align="center" valign="top" rowspan="1" colspan="1">4.33</td><td align="center" valign="top" rowspan="1" colspan="1">6.67</td><td align="center" valign="top" rowspan="1" colspan="1">7.25</td><td align="center" valign="top" rowspan="1" colspan="1">9.25</td><td align="center" valign="top" rowspan="1" colspan="1">8.60</td></tr></tbody></table><table-wrap-foot><fn id="TFN1"><label>a</label><p>Scenario 1: The intermediate health had a moderate effect on the outcome, it could occur in any year <italic>t</italic>, and the intermediate health event in year <italic>t</italic> &#x02212; 1 predicted leaving work in year <italic>t</italic>.</p></fn><fn id="TFN2"><label>b</label><p>Scenario 2: The intermediate health event had a stronger effect on the outcome than in Scenario 1. All other relationships between variables were as in Scenario 1.</p></fn><fn id="TFN3"><label>c</label><p>Scenario 3: The intermediate health event in year <italic>t</italic> &#x02212; 10 predicted leaving work in year <italic>t</italic>. All other relationships between variables were as in Scenario 1.</p></fn><fn id="TFN4"><label>d</label><p>Scenario 4: The intermediate health event could only occur during the first year, <italic>t</italic> = 1. All other relationships between variables were as in Scenario 1.</p></fn></table-wrap-foot></table-wrap></floats-group></article>