Local Recurrence After Excision of Merkel Cell Carcinoma: What Do We Actually Know, and What Should We Do Next?

2026
Merkel Cell Carcinoma
Radiation Therapy
Local Recurrence
This Perspectives on the Science article summarizes a multidisciplinary Society of Cutaneous Oncology Journal Club discussion of Kavanagh et al., Local Recurrence Following Complete Surgical Excision of Primary Merkel Cell Carcinoma. The discussion examines uncertainty in local recurrence risk after margin-negative excision, variation in postoperative radiation therapy practice, and future evidence-generation approaches that integrate curated multi-institutional data, causal inference, and probabilistic modeling to support decision-making under uncertainty.
Authors
Affiliations

David Michael Miller, MD, PhD

Massachusetts General Hospital

Harvard Medical School

Kevin Scott Emerick, MD

Mass Eye and Ear

Harvard Medical School

Howard Lane Kaufman, MD

Mass Eye and Ear

Harvard Medical School

Aubriana M. McEvoy, MD, MS

Washington University School of Medicine in St. Louis

Siteman Cancer Center

Lisa Zaba, MD, PhD

Stanford University School of Medicine

Stanford Cancer Center

Ajay Sharma, MD

Massachusetts General Hospital

Harvard Medical School

Adewunmi Adelaja, MD, PhD

Harvard Combined Dermatology Residency Program

Harvard Medical School

Ann W. Silk, MD

Dana-Farber Cancer Institute

Harvard Medical School

Nikhil I. Khushalani, MD

H Lee Moffitt Cancer Center and Research Institute

Isaac Brownell, MD, PhD

National Institutes of Health

Daniel G. Coit, MD

Memorial Sloan Kettering Cancer Center

Paul Nghiem, MD, PhD

University of Washington

Fred Hutch Cancer Center

Vernon K Sondak, MD

H Lee Moffitt Cancer Center and Research Institute

Published

July 16, 2026

Doi
Keywords

Merkel Cell Carcinoma, Local Recurrence, Adjuvant Radiation, Postoperative Radiation, Causal Inference

Equal senior authorship: Daniel G. Coit, Paul Nghiem, and Vernon K. Sondak contributed equally as senior authors. Within this designation, names are listed alphabetically.

Featured Article

Kavanagh, F. G. et al. Local Recurrence Following Complete Surgical Excision of Primary Merkel Cell Carcinoma. Annals of Surgical Oncology (2026).1

Introduction

The natural history of local recurrence after excision of localized Merkel cell carcinoma (MCC) remains incompletely defined. This uncertainty is not merely descriptive: it directly informs one of the most debated decisions in localized MCC management, namely whether adjuvant radiation therapy (RT) should be recommended after complete excision. Although postoperative RT is frequently considered for patients with adverse clinicopathologic features, substantial uncertainty persists regarding the baseline risk of local recurrence, the clinicopathologic factors that meaningfully modify that risk, the patients most likely to derive clinically meaningful benefit from treatment, and the outcomes of patients who develop local recurrence and undergo salvage therapy24. More fundamentally, this debate raises a broader scientific question: how should clinicians make evidence-based decisions when definitive evidence is unlikely to emerge in the foreseeable future?

Reported rates of local recurrence following excision of MCC vary substantially across the literature. Historical studies reinforced the perception that MCC carries a high risk of local recurrence at the primary site; in their 1991 review of 234 patients, Shaw and Rumball reported local recurrence after excision alone in 39% of patients5. By contrast, several more contemporary institutional cohorts have reported substantially lower local recurrence rates, often in the range of approximately 1–8%1,68. However, these estimates are not directly interchangeable. Published series differ with respect to patient selection, margin status, pathologic nodal evaluation, use or omission of adjuvant RT, definitions of local versus satellite, in-transit, or regional recurrence, and duration of follow-up. Consequently, clinicians are often asked to make treatment recommendations using recurrence estimates that differ by more than an order of magnitude, while recognizing that those estimates arise from different clinical eras, staging practices, treatment strategies, and endpoint definitions. The challenge, therefore, is not simply to identify a single “correct” recurrence rate, but to determine how best to reason and make decisions in the face of persistent uncertainty.

Against this background, Kavanagh and colleagues recently reported a large single-institution experience focused specifically on local recurrence following complete surgical excision of primary MCC1. The study is best understood as a natural-history analysis rather than a comparative trial of postoperative RT. Its clinical relevance, however, is direct: because primary-site RT is intended to prevent local recurrence, a focused estimate of that event inevitably informs radiation decision-making.

The clinical implications of this uncertainty are considerable. Current NCCN guidelines recommend consideration of postoperative RT for patients with adverse risk factors including lymphovascular invasion (LVI), immunosuppression, primary tumor size greater than 1 cm, close or positive margins, positive sentinel lymph node biopsy (SLNB), and head and neck primary site2. Yet the evidentiary basis supporting these recommendations remains heterogeneous (Table 1). Many of the risk factors incorporated into modern treatment algorithms were identified through retrospective studies designed to evaluate recurrence broadly, nodal positivity, disease-specific survival, or overall survival, rather than local recurrence specifically1,926. As a result, relatively few studies have been designed or adequately powered to determine which clinicopathologic features independently predict local recurrence after complete excision.

Table 1. Selected Evidence Informing NCCN Adverse Risk Factors in Merkel Cell Carcinoma
Studies informing adverse clinicopathologic features commonly used in postoperative radiation therapy decision-making.
Key Study Design / N Principal Finding Relevant to Risk Factor
Lymphovascular Invasion
Fields et al. (MSKCC) (2011) Retrospective, n = 153 Patients with LVI were more likely to have a positive SLNB. No patient without LVI experienced recurrence or MCC-specific death in this cohort.
Conic et al. (NCDB) (2019) Population-based, n = 3,048 LVI was associated with greater odds of SLN positivity (OR 3.45, p < 0.001) and was the strongest predictor evaluated.
Harounian et al. (2021) Retrospective, n = 122 Compared with patients without LVI, those with LVI had worse RFS (HR 2.3, p = 0.04), worse DSS (HR 5.2, p = 0.007), and more frequent SLN positivity (62.0% vs. 25.5%, p = 0.003).
Bartlett et al. (MSKCC) (2026) Retrospective, n = 342 Tumors with LVI and tumors >2 cm were more likely to have a positive SLN on multivariable analysis.
Immunosuppression
Liang et al. (U. Wisconsin) (2015) Retrospective, n = 87 Immunosuppressed patients had worse survival than immunocompetent patients (HR 2.01, 95% CI 1.1–3.7).
Yusuf et al. (NCDB) (2021) Population-based, n = 3,882 Immunosuppressed patients had worse 3-year OS than immunocompetent patients (44.6% vs. 68.7%, p < 0.0001; adjusted mortality HR 2.36, 95% CI 2.03–2.75). Among immunosuppressed patients, solid-organ transplant recipients had the lowest 3-year OS (32.7%).
McEvoy et al. (Seattle) (2022) Prospective cohort, n = 618 Immunosuppressed patients had a higher risk of recurrence than immunocompetent patients (HR 2.4, 95% CI 1.7–3.3, p < 0.001).
Primary Tumor >1 cm
Stokes et al. (VA) (2009) VA database, n = 346 Regional lymph-node metastasis was less frequent for tumors ≤1 cm than for tumors >1 cm (4% vs. 24%, p < 0.0001).
Smith et al. (Moffitt) (2015) Retrospective, n = 375 Increasing tumor diameter was associated with greater SLN positivity (p = 0.007), worse OS (p = 0.003), and worse DSS (p = 0.008).
Bartlett et al. (MSKCC) (2026) Retrospective, n = 342 Even among patients with tumors ≤1 cm and no LVI, 17% had a positive SLN, indicating that small size alone did not exclude occult nodal disease.
Close/Positive Margins
Andruska et al. (Wash U) (2021) Retrospective, n = 79 stage I–II Patients with close or positive margins ≤3 mm had worse distant RFS (HR 6.83, 95% CI 1.80–25.9). Each additional 1 cm of WLE margin was associated with improved regional RFS, distant RFS, and DSS.
Tarabadkar et al. (Seattle) (2021) Retrospective, n = 188 Among patients receiving adjuvant RT, local recurrence was approximately 1% regardless of margin width. Without RT, local recurrence was 20% after margins ≤1 cm versus 0% after margins >1 cm (p = 0.049).
Lodde et al. (DeCOG) (2025) Multicenter, n = 1,049 Margins of 1–2 cm were associated with the most favorable outcomes; margins >2 cm provided no additional observed benefit.
Kavanagh et al. (MSKCC) (2026) Retrospective, n = 447 Among patients treated with negative-margin surgery alone, the estimated 1-year local recurrence rate was 1.8%.
Positive SLNB
Servy et al. (2016) Retrospective, n = 87 successful SLNB Patients with a negative SLNB had longer DFS and OS than patients with a positive SLNB on multivariable analysis.
Harounian et al. (2021) Retrospective, n = 122 Patients with a positive SLN had worse RFS than patients with a negative SLN (HR 3.6, p = 0.007).
Ahn et al. (Westmead) (2026) Retrospective, n = 75 SLNB was positive in 37% of patients. Among patients with nodal disease, the overall false-negative rate was approximately 20%.
Bartlett et al. (MSKCC) (2026) Retrospective, n = 342 Patients with a positive SLN had a higher 5-year recurrence rate than patients with a negative SLN (39% vs. 12%, p < 0.001); SLN status was the only independent recurrence predictor identified.
Head/Neck Primary Site
Cullison et al. (SEER) (2021) Population-based, n = 9,407 Among patients with localized MCC, scalp or neck primary site was associated with the highest 5-year cumulative incidence of MCC-specific mortality, approximately 26%.
McEvoy et al. (Seattle) (2022) Prospective cohort, n = 618 Anatomic site was not significantly associated with recurrence in this cohort (p = 0.44), in which adjuvant RT was used frequently.
Nayak et al. (Pan-Canadian) (2023) Multicenter, n = 400 head/neck MCC In this head and neck MCC cohort, 5-year OS was 42%, and approximately one-half of patients experienced recurrence within 5 years.
Brenner et al. (Israel) (2025) Multicenter, n = 191 Patients with head and neck primaries had worse 5-year OS than patients with primaries at other anatomic sites (51.6% vs. 65.2%, p = 0.025; adjusted mortality HR 1.77, 95% CI 1.10–2.84, p = 0.018).
Wang et al. (2026) Single-center cohort, n = 86 For head and neck MCC, SLNB sensitivity was 52.6%. Among patients with nodal disease who underwent SLNB, the false-negative rate was 26.7%; SLNB mapping or identification failed in 29.0%.
The table is organized by the adverse clinicopathologic feature each study helps inform. Selected studies evaluated recurrence, sentinel lymph-node positivity, disease-specific survival, overall survival, surgical margins, or anatomic site; not all evaluated primary-site local recurrence as a discrete endpoint.
Abbreviations: CI, confidence interval; DeCOG, Dermatologic Cooperative Oncology Group; DFS, disease-free survival; DSS, disease-specific survival; HR, hazard ratio; LN, lymph node; LR, local recurrence; LVI, lymphovascular invasion; MCC, Merkel cell carcinoma; MSKCC, Memorial Sloan Kettering Cancer Center; NCDB, National Cancer Database; OR, odds ratio; OS, overall survival; RFS, recurrence-free survival; RT, radiation therapy; SLN, sentinel lymph node; SLNB, sentinel lymph node biopsy; VA, Veterans Affairs; WLE, wide local excision.

Thus, while the NCCN-listed adverse features used to guide postoperative RT recommendations are clinically reasonable and widely used, the evidence linking them specifically to local recurrence remains largely indirect. The challenge, therefore, is not simply that more evidence is needed. Rather, clinicians must make consequential treatment decisions using evidence that is incomplete, variably applicable, and unlikely to become definitive in the near future.

To explore how these evidence gaps are interpreted in practice, the Society of Cutaneous Oncology (SoCO) convened a multidisciplinary Journal Club bringing together dermatologists, medical oncologists, surgical oncologists, radiation oncologists, and trainees from across the MCC community. Here, we review the findings of Kavanagh et al., summarize survey results and expert discussion generated during that session, and consider how future evidence generation in MCC can support more rigorous decision-making under persistent uncertainty.

Evolution of the Evidence Base

Early Observations and the Emergence of Adjuvant Radiation

Much of the early literature describing recurrence patterns in MCC was generated across clinical eras that differed substantially in staging, pathologic evaluation, surgical management, use of adjuvant RT, and follow-up. Many historical series also predated routine SLNB-based nodal staging, which was incorporated into NCCN recommendations by 2006 and became increasingly integrated into practice after the first MCC-specific AJCC staging system was introduced in 20102,27,28. Recurrence estimates therefore varied widely across studies and often reflected mixed endpoints that combined local, regional, and distant failures. These reports, including the 1991 review by Shaw and Rumball, helped establish the perception that MCC was associated with a substantial risk of local recurrence after excision alone5.

In the absence of prospective randomized data, several institutional series and population-based analyses subsequently suggested that postoperative RT could improve local control and, in some reports, overall survival6,2931. These studies contributed to the increasing adoption of adjuvant RT and established local control as a major therapeutic objective in localized MCC. However, interpretation of this literature remains constrained by retrospective design, inconsistent definitions of treatment and recurrence, variable follow-up, incomplete details regarding RT field, dose, and timing, and limited physician-curated review of recurrence events. As noted above, population-based registries provide complementary scale, but generally lack information on the date and site of first recurrence, salvage treatment, disease-specific survival, and competing causes of death. As a result, associations between adjuvant RT and improved overall survival in registry-based analyses should be interpreted cautiously and should not be assumed to establish a causal survival benefit.

Risk Stratification and Selective Radiation

As clinical practice became increasingly risk-adapted, the question evolved. Rather than asking only whether RT could improve local control, clinicians increasingly asked which patients were most likely to benefit after complete excision. During this period, adverse clinicopathologic features including lymphovascular invasion, immunosuppression, larger tumor size, positive SLNB, close or positive margins, and head and neck location emerged as important prognostic factors and were progressively incorporated into guideline recommendations3,7,13. Importantly, however, many of these features were identified as predictors of recurrence broadly, nodal involvement, systemic progression, disease-specific survival, or overall survival, rather than local recurrence specifically. Their use as indications for adjuvant RT therefore rests, in part, on the assumption that features associated with worse overall prognosis also identify patients at sufficiently increased risk of local recurrence to justify treatment.

At the same time, several institutional cohorts began reporting low rates of local recurrence among selected patients treated with complete excision and selective use or omission of adjuvant RT. Fields and colleagues reported a local recurrence rate of approximately 3% following complete resection, although 23% of patients received adjuvant local RT, making this a selective-therapy rather than surgery-alone recurrence estimate7. Frohm and colleagues observed only two true local or satellite recurrences among 104 patients managed without local RT, including one true local recurrence and one satellite recurrence, although the cohort was restricted to tumors smaller than 2 cm and reflected care at a high-volume center8. More recently, recurrence analyses from large prospective cohorts have reinforced that MCC recurrence risk varies substantially by stage and time from diagnosis, underscoring the need to interpret recurrence estimates in relation to patient and disease context13. Similarly, the observation arm of ADMEC-O provides a useful reminder that outcomes in completely resected MCC cohorts may differ from historical expectations, although its endpoint was disease-free survival rather than isolated local recurrence, and patients could receive other standard local therapies32. In this prospective randomized trial, local-regional relapse, not further subdivided, occurred in 10 of 61 patients on the observation arm after a minimum follow-up of 2 years for surviving patients, even though 74% of patients had received postoperative RT. Nevertheless, the investigators concluded that “the spontaneous course of MCC was better than historical data suggested”32.

These observations must be interpreted alongside a substantial body of retrospective evidence associating adjuvant RT with improved outcomes. In a systematic review and meta-analysis of 29 studies including 17,179 patients, Petrelli and colleagues found that adjuvant RT was associated with improved local disease-free survival, locoregional disease-free survival, and overall survival4. Similarly, Bhatia and colleagues, using the National Cancer Data Base, reported an association between adjuvant RT and improved overall survival in stage I and II MCC, although not in stage III disease3. However, registry-based overall survival analyses should be interpreted cautiously. Population-based registries generally lack information on the date and site of first recurrence, salvage treatment, and disease-specific survival, and overall survival is vulnerable to competing mortality and treatment-selection bias in early-stage MCC. Thus, reported associations between adjuvant RT and improved overall survival should not be assumed to establish a causal survival benefit.

More recent registry-based studies have also suggested that the timing of postoperative or adjuvant RT may be clinically important. In a prospective registry analysis of 124 patients with stage I/II MCC treated with postoperative RT, Alexander and colleagues found that initiation of RT more than 8 weeks after surgery was associated with increased locoregional recurrence and MCC-specific death33. Similarly, in the DeCOG multicenter registry of 1,049 patients, Lodde, Becker, and colleagues reported that delayed adjuvant RT was associated with poorer progression-free probability and disease-specific survival, while expansion of RT to the lymph node bed did not provide additional benefit21. These observational data do not eliminate concerns regarding treatment selection and confounding, but they suggest that timing, field selection, and patient risk may all matter.

Consequently, the available literature presents a persistent tension. Meta-analytic and registry-based studies report associations between adjuvant RT and improved outcomes, yet selected institutional cohorts report very low rates of local recurrence after complete excision, including among patients with adverse clinicopathologic features. Reconciling these observations requires greater clarity about the baseline risk of local recurrence, the clinical features that meaningfully modify that risk, the clinical consequences of local recurrence, and the absolute benefit that RT provides in different patient subgroups. This tension provides the immediate context for the Kavanagh study.

The Kavanagh Study

Study Overview

Kavanagh and colleagues addressed this question through a retrospective review of patients treated at Memorial Sloan Kettering Cancer Center over a four-decade period1. The investigators identified 447 patients with clinically localized MCC (no clinical evidence of nodal or distant metastasis at presentation) who underwent complete surgical excision with pathologically negative margins. Of these, 393 patients were managed with surgery alone, while 54 received adjuvant RT to the primary site. The authors focused specifically on local recurrence outcomes among patients managed without postoperative primary-site RT.

The study addressed a notable limitation of much of the existing literature. While numerous reports have evaluated recurrence broadly or examined the association between RT and clinical outcomes, relatively few have focused specifically on local recurrence after complete excision. By restricting the analysis to patients with pathologically negative margins and clearly defined local recurrence outcomes, Kavanagh and colleagues provided a more direct estimate of the event most relevant to decisions regarding postoperative primary-site RT.

Principal Findings

Among the 447 patients included in the study, only eight local recurrences were identified. All occurred among patients managed with surgery alone. Seven of the eight local recurrences occurred within the first year following treatment, yielding an estimated one-year local recurrence rate of 1.8% among the 393 patients managed without postoperative RT1.

Several observations deserve emphasis. First, local recurrence was uncommon. Second, local recurrences tended to occur early, consistent with broader MCC recurrence patterns in which recurrence risk is highest soon after diagnosis. Third, the clinical significance of local recurrence appeared to depend strongly on whether it occurred as an isolated event or in association with regional or distant disease. Among the six patients who developed isolated local recurrence, all achieved durable local control after additional local therapy. By contrast, the two patients who developed local recurrence with synchronous regional and distant metastatic disease died of disease, suggesting a biologically aggressive phenotype rather than a failure of local salvage alone. These findings raise important questions regarding the balance between recurrence prevention, salvageability, treatment burden, and RT-associated morbidity.

The study also highlights the challenges inherent in MCC outcomes research. Despite representing one of the largest institutional series focused specifically on local recurrence, only eight local recurrence events were observed. Given this low event rate, the study had limited ability to identify clinicopathologic predictors of local recurrence or to precisely estimate recurrence risk within individual subgroups.

Persistent Uncertainty in Current Practice

If the available evidence supported a single, compelling interpretation, one might expect practice patterns to converge. Instead, uncertainty surrounding local recurrence risk has translated into substantial variation in how clinicians estimate recurrence risk, recommend postoperative RT, and individualize treatment decisions. To better understand how clinicians interpret the current evidence, participants in the June 2026 Society of Cutaneous Oncology Journal Club completed a structured survey examining practice patterns, perceptions of recurrence risk, and attitudes toward postoperative RT.

Figure 1: Primary role of survey respondents. Respondents represented a multidisciplinary group, with medical oncology and medical dermatology comprising the largest groups, followed by surgical oncology, Mohs surgery, radiation oncology, trainees, and research.
Figure 2: Monthly MCC patient volume among survey respondents. Most respondents reported personally evaluating or managing at least one patient with MCC in a typical month; 28% reported managing more than 10 MCC patients per month.

Variation in Baseline Risk Estimates and Radiation Practice

Despite broad agreement that local recurrence is an important outcome in MCC, survey responses revealed substantial variation in both perceived recurrence risk and approaches to postoperative RT (Figures 3–8). Baseline estimates of local recurrence after margin-negative excision ranged from less than 2% to greater than 10%, with most respondents estimating the risk to be between 6% and 10% (Figure 3). Notably, nearly one-half of respondents reported being minimally familiar or not familiar at all with the Kavanagh et al. study before the session (Figure 4), highlighting the limited dissemination of newer data specifically focused on local recurrence. Together, these findings suggest that variation in clinical practice reflects not simply differences in physician preference, but differences in how clinicians interpret an incomplete and evolving evidence base.

This uncertainty was also reflected in treatment recommendations. Equal proportions of respondents (20%) reported routinely recommending postoperative RT for nearly all patients or rarely recommending it at all, while the largest group (44%) indicated that their recommendations were individualized according to patient and tumor-specific risk factors (Figure 5). The factors most frequently influencing RT recommendations included tumor size greater than 1 cm, close or positive margins, lymphovascular invasion, immunosuppression, and positive sentinel lymph node biopsy findings (Figure 6). These factors closely mirror those incorporated into current NCCN recommendations, suggesting that variation in practice reflects not disregard for existing guidelines, but differences in how clinicians weigh the strength of the evidence supporting those risk factors and their implications for individual patients.

Figure 3: Estimated local recurrence risk after margin-negative excision alone for clinically localized MCC prior to the Journal Club discussion. Respondents were asked: “Before this Journal Club, what local recurrence risk would you generally estimate after margin-negative excision alone for a clinically localized MCC primary?” Most respondents estimated local recurrence risk at 6–10%, although estimates ranged from <2% to >10%.
Figure 4: Familiarity with the Kavanagh et al. study before participation in the Journal Club. Respondents were asked: “Before this Journal Club, how familiar were you with the Kavanagh et al. study on local recurrence after complete surgical excision of primary Merkel cell carcinoma?” Familiarity varied substantially, with nearly one-half of respondents reporting that they were minimally familiar or not familiar at all.
Figure 5: Frequency with which respondents recommend or support adjuvant RT to the primary tumor site after complete excision of clinically localized MCC with negative margins and no clinically evident nodal or distant disease. Respondents were asked how often they recommend or support adjuvant primary-site RT in this setting. Most respondents reported recommending RT often depending on risk factors, while smaller groups reported almost always or rarely/never recommending RT.
Figure 6: Factors reported to most strongly influence decisions regarding adjuvant primary-site RT. Respondents could select multiple factors. Tumor size >1 cm, close or positive margins, lymphovascular invasion, immunosuppression, and positive sentinel lymph node biopsy were the most commonly selected factors.
Figure 7: General approach to adjuvant primary-site RT when guideline-defined adverse features are present after margin-negative excision. Respondents were asked how they generally approach adjuvant primary-site RT when guideline-defined adverse features are present. Responses were split between generally following guideline-based recommendations and individualizing decisions based on overall recurrence risk and treatment burden.
Figure 8: Evidence most likely to change current practice regarding use or omission of adjuvant primary-site RT after complete excision of clinically localized MCC. “Prospective registry or clinical trial data” was the option most frequently selected.

Evidence Thresholds and Areas of Equipoise

Survey responses also demonstrated that uncertainty in MCC management extends beyond the binary decision of whether postoperative RT should be recommended. Respondents reported a broad range of preferred dose and fractionation approaches, including conventional fractionation, moderate hypofractionation, and single-fraction 8 Gy, with substantial variation across clinical scenarios (Figures 9–10). Similarly, respondents differed in the degree to which they individualize dose and fractionation based on perceived recurrence risk, with equal proportions reporting minimal, somewhat, or substantial variation in approach.

Practice variation was therefore not limited to whether RT was used. Rather, uncertainty extended across multiple dimensions of care: baseline recurrence risk estimates, treatment thresholds, risk factors used in decision-making, dose/fractionation strategies, and assumptions about the detectability and salvageability of local recurrence. This pattern is consistent with meaningful clinical equipoise and supports the need for prospective evidence generation capable of informing not only whether RT should be used, but for whom, in what form, under what clinical circumstances, and with what expected benefit compared with observation and salvage therapy if local recurrence occurs.

Figure 9: Preferred dose and fractionation when respondents recommend or support adjuvant RT to the primary tumor site after complete excision of clinically node-negative MCC. Respondents were asked what dose/fractionation they most commonly use when they recommend or support adjuvant radiation. Responses included conventional fractionation, moderate hypofractionation, single-fraction 8 Gy, and scenario-dependent approaches.
Figure 10: Degree to which preferred dose and fractionation change based on perceived risk of local recurrence. Respondents varied in whether they generally use the same regimen, individualize within a limited range, or use substantially different approaches for lower- versus higher-risk cases.

Clinical Scenarios Reveal Areas of Consensus and Disagreement

When respondents were presented with specific clinical scenarios, a consistent pattern emerged. Broad consensus was observed at the extremes of risk, whereas greater variation appeared in intermediate-risk settings where the balance between expected benefit and treatment burden was less certain (Figures 11–18). Observation was generally favored for a low-risk extremity primary lacking adverse features (Figure 11), whereas most respondents recommended RT for a patient with multiple guideline-defined risk factors, including head and neck location, lymphovascular invasion, close margins, and tumor size greater than 1 cm (Figure 13). These responses suggest that clinicians largely agree on which scenarios appear higher or lower risk, but differ in how they translate uncertain estimates of absolute benefit into treatment recommendations.

The scenarios also demonstrated that patient factors may influence the preferred approach to radiation delivery more than the binary decision to treat. In a scenario incorporating advanced age, frailty, travel distance, and treatment burden, most respondents still recommended RT (Figure 15), but single-fraction RT emerged as the dominant preferred regimen among those who would treat (Figure 16). Similarly, when a patient expressed a strong preference to minimize treatment burden, recommendations were distributed across observation and abbreviated or conventional RT approaches (Figure 17). Finally, in the setting of microscopic nodal involvement, respondents frequently selected nodal-basin RT and, in many cases, also selected RT to the local excision bed (Figure 18). This finding suggests that a positive SLNB often increases perceived treatment intensity and broadens consideration of RT fields, although whether microscopic nodal disease should influence treatment of the local excision bed remains a separate and unresolved question.

Taken together, these scenario-based responses suggest that practice variation in localized MCC is not random. Rather, clinicians appear to integrate tumor risk, anatomic site, nodal status, patient frailty, treatment burden, and feasibility of abbreviated radiation approaches. Future evidence generation should therefore capture not only recurrence outcomes, but also salvage treatment, disease-specific outcomes, and the patient- and clinician-level factors that shape treatment decisions. Such data are essential if future research is to improve decision-making rather than simply produce more precise recurrence estimates.

Figure 11: Responses to a low-risk localized MCC scenario. Respondents were asked about a 74-year-old immunocompetent patient with a 0.8 cm forearm MCC, negative margins, 0.7 cm closest pathologic margin, absent lymphovascular invasion, and negative SLNB. Most respondents favored observation without adjuvant RT or discussion of RT without generally recommending treatment.
Figure 12: Preferred dose and fractionation among respondents who would recommend RT in the low-risk scenario. Among the small number of respondents recommending RT, single-fraction 8 Gy was the most commonly selected approach.
Figure 13: Responses to a higher-risk head and neck MCC scenario. Respondents were asked about a 76-year-old immunocompetent patient with a 1.6 cm cheek MCC, negative margins with a closest pathologic margin of 0.2 cm, lymphovascular invasion present, and negative SLNB. Most respondents recommended adjuvant RT.
Figure 14: Preferred dose and fractionation among respondents recommending RT in the higher-risk head and neck scenario shown above. Conventional fractionation, moderate hypofractionation, and single-fraction RT were all represented, suggesting lack of a dominant dose/fractionation strategy.
Figure 15: Responses to a higher-risk MCC scenario incorporating advanced age, frailty, comorbidity, and substantial travel burden. Respondents were asked about an 86-year-old patient who lives far from the treatment center and has a 1.8 cm lower leg MCC, negative margins, 0.4 cm closest pathologic margin, lymphovascular invasion present, and negative SLNB. Most respondents continued to recommend adjuvant RT despite competing clinical considerations.
Figure 16: Preferred dose and fractionation among respondents recommending RT in the frailty and treatment-burden scenario above. Single-fraction 8 Gy was the most commonly selected approach, suggesting that patient factors may influence RT delivery even when the decision to treat is retained.
Figure 17: Responses to a preference-sensitive scenario involving a patient interested in minimizing treatment burden. Respondents were asked about a 72-year-old immunocompetent patient with a 1.2 cm upper arm MCC, negative margins, 0.6 cm closest pathologic margin, absent lymphovascular invasion, and negative SLNB, who was open to radiation only if it meaningfully reduced risk. Responses were distributed across observation, single-fraction RT, hypofractionation, and conventional fractionation.
Figure 18: Management preferences for a patient with microscopic nodal involvement identified on SLNB. Respondents were asked about a 70-year-old immunocompetent patient with a 1.4 cm forearm MCC, negative margins, 0.8 cm closest pathologic margin, absent lymphovascular invasion, and a positive SLNB without clinical or radiographic nodal or distant disease. Nodal-basin RT and primary-site RT were both commonly selected.

Defining the Clinical Question

Although the Kavanagh study focused on local recurrence after complete surgical excision of primary MCC, the broader debate surrounding adjuvant RT may be driven as much by the questions being asked as by the available evidence. Clinicians may appear to disagree about the role of RT while, in reality, attempting to answer fundamentally different clinical questions. Disagreement may arise not only from how the available data are interpreted, but also from assumptions regarding baseline recurrence risk, the mechanisms by which RT might provide benefit, the clinical significance of the events being prevented, the effectiveness of salvage therapy, and the harms associated with treatment. Clarifying these questions is essential before considering what type of future evidence is needed—and whether additional retrospective studies are likely to resolve the controversy.

What Is the True Baseline Risk?

A natural instinct is to ask, “What is the true baseline risk of local recurrence?” Yet the available evidence suggests that this question may not have a single definitive answer. Among patients treated with margin-negative excision alone, Kavanagh et al. observed a 1-year local recurrence rate of just 1.8%, substantially lower than estimates reported in many historical series and lower than recurrence risks commonly cited during clinical decision-making1.

Interpretation of this estimate remains challenging. Approximately 12% of patients in the Kavanagh cohort received adjuvant RT, and these patients differed from those managed with surgery alone: they were younger, had larger tumors, more frequent LVI or tumor emboli, and a higher prevalence of NCCN-listed high-risk features, although they also had wider clinical margins1. These differences suggest nonrandom treatment selection, but they do not establish that the features prompting RT were independent predictors of local recurrence. Indeed, 75% of patients managed with surgery alone had at least one NCCN-listed high-risk feature, and the low number of local recurrence events precluded reliable identification of clinicopathologic predictors. As a result, the observed local recurrence rate among surgery-alone patients may not represent the risk that would have been observed in an entirely unselected population, but the Kavanagh data also challenge the assumption that currently used adverse features reliably discriminate patients at high risk for local recurrence.

A key area of uncertainty is whether the real-world risk of local recurrence is closer to 2%, as suggested by Kavanagh et al., or meaningfully higher in broader clinical practice. Several prior institutional and multicenter series report higher rates of recurrence or progression, but many were smaller, included mixed-treatment cohorts, used heterogeneous recurrence endpoints, or did not isolate local recurrence after margin-negative excision with the same methodological specificity. At the same time, framing the debate as a choice between 2% and 10% may oversimplify a more nuanced reality. The risk of local recurrence is unlikely to be a single fixed quantity applicable to all patients; rather, it almost certainly varies according to patient characteristics, tumor biology, treatment selection, and other measured and unmeasured factors. The more meaningful question may therefore be not whether one study has identified the “true” recurrence rate, but how uncertainty about baseline risk should influence clinical decision-making. Ultimately, the implications of this uncertainty depend not only on the absolute recurrence risk, but also on the expected benefit of RT, the effectiveness of salvage therapy, and the clinical consequences of the events being prevented.

What Does Primary-Site RT Prevent?

A central unresolved question is whether primary-site RT should be understood narrowly as an intervention to prevent local recurrence, or more broadly as a treatment that may sterilize residual microscopic dermal or lymphatic disease before it declares itself as in-transit, nodal, or distant progression.

The Kavanagh data most directly address the first question: whether complete excision without primary-site RT is followed by clinically consequential local failure. After negative-margin surgery, local recurrence was uncommon and occurred early. Among the eight patients who developed local recurrence, six had local recurrence alone, while two had local recurrence with simultaneous regional and distant metastatic disease. All six patients with isolated local recurrence received additional local therapy, achieved durable local control, and none subsequently died of MCC. By contrast, the two patients with local recurrence accompanied by synchronous regional and distant metastases progressed rapidly despite systemic therapy, a pattern the authors interpreted as more consistent with aggressive disease biology than failure of initial local treatment. Thus, within the limits of an analysis focused on local recurrence, these data provide little evidence that omission of primary-site RT commonly leads to unsalvageable local failure or to a substantial burden of metastatic relapse clearly arising from uncontrolled local disease.

The small number of recurrence events limits definitive conclusions regarding whether RT might prevent some events beyond local recurrence alone. Nevertheless, the available data more strongly support a model in which most isolated local recurrences behave as local events that can be effectively salvaged, rather than as precursors of later regional or distant progression. This distinction matters. If primary-site RT primarily prevents rare and salvageable local recurrences, then many patients may receive treatment without meaningful clinical benefit. If, instead, RT interrupts a broader pathway of locoregional or systemic progression, then local recurrence alone may underestimate its clinical value.

At present, the Kavanagh data provide stronger support for the former hypothesis than the latter, while also underscoring the need for future studies to capture not only whether recurrence occurs, but what type of recurrence occurs, how it is treated, and whether it affects disease-specific outcomes.

How Salvageable Is the Event Being Prevented?

The clinical value of adjuvant RT depends not only on how often recurrence occurs, but also on what happens when recurrence develops. A local recurrence that can be detected early and successfully salvaged through surgery, RT, or systemic therapy carries different implications than a recurrence that serves as a precursor to uncontrolled regional or distant dissemination.

An implicit assumption in many discussions of local recurrence is that all recurrences carry equivalent clinical importance. In practice, however, the consequences of recurrence may vary substantially according to timing, location, patient fitness, and the availability of effective salvage options. This consideration is particularly relevant in MCC, where advances in immune checkpoint blockade have altered the therapeutic landscape and improved outcomes for many patients with recurrent disease.

The Kavanagh experience is notable in this regard. Among patients with isolated local recurrence, durable disease control was achieved with additional local therapy, and none subsequently died of MCC. In contrast, the two patients who developed local recurrence together with regional and distant metastatic disease experienced poor outcomes despite systemic therapy. This pattern suggests that the clinical significance of recurrence may depend less on the presence of a local recurrence itself and more on the biological context in which it occurs.

If most isolated local recurrences remain highly salvageable, then the threshold for recommending routine adjuvant RT may be higher. Conversely, if local recurrence occurs in association with regional or distant progression, it is important to distinguish whether the local event is a cause of subsequent spread, a marker of aggressive tumor biology, or one manifestation of disease that was already systemic. These scenarios have different implications for adjuvant RT. Preventing a local recurrence would have broader clinical value only if that event would otherwise contribute meaningfully to subsequent progression. If, instead, local recurrence occurs simultaneously with regional or distant disease, then treatment of the primary site addresses only one component of a more aggressive biological process. Thus, the value of adjuvant RT cannot be judged solely by counting recurrences. It must also be judged by the clinical significance of the events prevented, the effectiveness of salvage therapy, and the extent to which prevention of local recurrence changes disease-specific outcomes.

What Are the Harms of Treatment?

Discussions regarding adjuvant RT often focus on recurrence prevention while devoting comparatively less attention to treatment-related harms. Yet the risk-benefit calculation depends critically on both sides of the equation.

Conventional postoperative RT may be associated with acute toxicity, fibrosis, lymphedema, alopecia, wound-healing complications, treatment burden, financial cost, and, in rare circumstances, serious long-term sequelae. These considerations become increasingly important when the anticipated absolute reduction in recurrence risk is small. Even if RT is highly effective at preventing local recurrence, the clinical value of treatment depends on whether the expected benefit outweighs these competing harms.

Importantly, RT should not be viewed as a single intervention. Contemporary practice encompasses a spectrum of approaches ranging from conventional fractionation to abbreviated hypofractionated regimens and single-fraction RT. Survey responses demonstrated substantial variation in preferred dose and fractionation schedules, particularly in scenarios involving advanced age, frailty, or logistical barriers to treatment. These findings suggest that the debate may not simply be whether RT should be delivered, but how it should be delivered and at what cost in terms of patient burden.

Taken together, these considerations suggest that disagreements regarding adjuvant RT in localized MCC extend well beyond estimates of local recurrence risk alone. Clinicians may arrive at different treatment recommendations because they hold different assumptions regarding baseline risk, treatment effect, salvageability, and toxicity. These questions are closely intertwined, and each may require different forms of evidence to answer. The challenge is therefore not simply to generate more data, but to generate evidence that meaningfully informs clinical decision-making.

Are We Asking More of Retrospective Studies Than They Can Answer?

Taken together, these considerations suggest that many of the questions surrounding adjuvant RT in localized MCC may be difficult to answer using conventional retrospective analyses alone. Decisions regarding RT are influenced by tumor characteristics, surgical technique, patient factors, physician judgment, and institutional practice patterns, creating substantial potential for confounding by indicationa. At the same time, local recurrence remains a relatively uncommon event, limiting statistical power even in large institutional series.

Confounding by indication occurs when the reason a treatment is given is also related to the outcome being studied. In the context of adjuvant RT for MCC, patients selected for radiation may have larger tumors, closer margins, lymphovascular invasion, head and neck primary sites, immunosuppression, or other features that make recurrence more likely.

As a result, comparing patients who received RT with those who did not can be misleading. If the RT group has similar outcomes despite having higher-risk disease, RT may be beneficial. If the RT group has worse outcomes, this may reflect baseline risk rather than harm from radiation. Statistical adjustment can reduce measured confounding, but uncertainty remains whenever treatment selection is influenced by clinical judgment or other factors that are incompletely measured.

Retrospective studies remain indispensable for understanding the natural history of disease, describing patterns of care, estimating recurrence risk, and generating clinically important hypotheses. Their principal limitation is not that they are retrospective, but that they are often asked to answer causal questions for which they were not designed. This distinction is particularly important when the objective shifts from describing recurrence to estimating which patients derive meaningful net benefit from treatment. Such questions require consideration of treatment selection, competing risks, and patient-centered outcomes—including toxicity, functional outcomes, cosmetic outcomes, travel burden, financial toxicity, and quality of life—which are often incompletely captured in retrospective datasets.

The rarity of local recurrence further compounds these challenges. If the true recurrence risk after margin-negative excision is approximately 2–8%, even large institutional cohorts may contain relatively few events. Small differences in event counts can therefore produce meaningfully different recurrence estimates, while subgroup analyses may remain unstable despite seemingly large sample sizes. Consequently, simply accumulating larger retrospective cohorts may not fully resolve the underlying uncertainty.

Viewed through this lens, the central challenge is no longer simply obtaining a more precise estimate of recurrence risk. Rather, it is developing an evidence-generation framework capable of estimating individualized treatment benefit while explicitly accounting for treatment selection, competing risks, patient-centered outcomes, and uncertainty. Addressing that challenge will likely require complementary approaches, including prospective data collection, harmonized analyses of existing multicenter cohorts, and modern causal and probabilistic methods capable of estimating treatment benefit under uncertainty. The goal is not to replace retrospective research, but to integrate multiple sources of evidence in a framework that supports better clinical decision-making.

What Should Future Evidence Generation Look Like?

The limitations of the existing evidence base do not diminish the importance of the clinical question. Rather, they clarify how future evidence generation should be designed. The central limitation is not simply that many studies are retrospective. It is that much of the available evidence comes from administrative datasets or retrospective institutional series in which treatment rationale, recurrence definitions, radiation fields and fractionation, salvage treatment, toxicity, treatment burden, patient preferences, and disease-specific outcomes are incompletely or inconsistently captured.

If the central goal is to identify which patients derive meaningful benefit from adjuvant primary-site RT, then future studies must capture more than whether recurrence occurs. They must also define baseline risk, treatment strategy, recurrence pattern, salvageability, toxicity, treatment burden, patient preferences, and the clinical rationale underlying treatment selection.

Survey responses supported the need for more structured multi-institutional evidence generation. Almost all respondents rated a prospective multicenter study of radiation after surgical excision of MCC as either very important or moderately important (Figure 19). This finding suggests that, despite differences in current practice, there is broad agreement that the existing evidence base remains insufficient to guide patient selection and treatment approach. Importantly, such a study should not simply catalogue whether adjuvant RT was used. Its highest value would be to more clearly define the natural history of local recurrence after complete excision, using standardized recurrence definitions, physician-curated event review, detailed treatment information, salvage treatment data, and disease-specific outcomes. In this respect, the Kavanagh study provides a useful methodological template: future datasets should seek to replicate, refute, or refine its findings across institutions and practice settings.

Figure 19: Priority assigned to prospective multicenter evidence generation after surgical excision of MCC. Respondents were asked: “How important would a prospective multicenter study be to better inform the role of radiation following surgical excision of MCC?” Most respondents rated such a study as very important or moderately important, supporting prospective evidence generation to address unresolved questions regarding patient selection, recurrence risk, treatment approach, salvageability, and disease-specific outcomes.

Prospective evidence generation should not be equated with the feasibility of a conventional randomized trial. A randomized trial would provide the cleanest estimate of treatment effect, but feasibility may be limited by the rarity of MCC, low local recurrence event rates, cost, heterogeneity in radiation approaches, and the need for genuine clinician and patient equipoise. The Journal Club discussion highlighted this challenge clearly: some clinicians favor adjuvant RT for nearly all patients with localized MCC, whereas others use it rarely or selectively. Conversely, an unstructured registry may be easier to implement but risks reproducing many of the same interpretive limitations that affect existing retrospective studies. The goal, therefore, is not to identify a single ideal study design, but to build an evidence-generation strategy capable of answering clinically meaningful questions despite the practical constraints of studying a rare cancer.

The most useful approach may lie between these extremes: a highly curated, multi-institutional observational effort, ideally prospective in design but also capable of incorporating protocolized retrospective case capture when key data elements can be reliably abstracted. Such efforts are feasible when organized around prespecified data elements, standardized definitions, transparent governance, equitable academic credit, and a shared commitment to data quality. A study designed in this manner would standardize recurrence definitions, capture the clinical rationale for treatment decisions, characterize radiation fields and fractionation in detail, record toxicity and patient-centered outcomes, and include salvage treatment and disease-specific outcomes as core endpoints. It should also be designed from the outset to support causal inference and hierarchical probabilistic modeling. Rather than asking whether retrospective or prospective evidence is inherently superior, this framework recognizes that rigorous science comes from combining complementary sources of evidence within a principled analytic approach.

Why Causal Inference Matters

A central challenge in studying adjuvant RT is that treatment is not assigned randomly. Clinicians selectively recommend RT to patients perceived to be at higher risk based on tumor size, lymphovascular invasion, immunosuppression, margin status, surgical technique, anatomic location, nodal status, frailty, and provider or institutional practice patterns. Consequently, treatment assignment itself contains information about recurrence risk.

Traditional analyses often compare outcomes among patients who received RT with those who did not. The difficulty is that these groups are often fundamentally different from the outset. Patients selected for RT are frequently those clinicians believe are most likely to recur. As a result, observed outcome differences may reflect underlying baseline risk rather than treatment effect. This is the central problem of confounding by indication.

A causal inference framework begins by explicitly modeling this decision process. Rather than treating treatment selection as an inconvenient source of bias, causal inference asks why treatment was selected, which factors influenced that decision, and how those factors also relate to recurrence. Directed acyclic graphs (DAGs) provide one useful approach by requiring investigators to specify the assumed relationships among tumor characteristics, patient characteristics, treatment selection, recurrence outcomes, salvage therapy, and survival before data analysis begins3436.

In the hypothetical DAG shown in Figure 20, several tumor- and patient-level variables influence both the likelihood of receiving RT and the risk of recurrence. Tumor size, lymphovascular invasion, margin status, primary site, SLNB status, immunosuppression, frailty, and surgical technique may all function as confounders if they affect both treatment selection and recurrence risk. Provider or institutional practice patterns may also influence treatment assignment, radiation fields, fractionation, surgical technique, and surveillance intensity. The red arrow represents the causal effect of interest: the effect of adjuvant RT on recurrence.

Figure 20: Hypothetical directed acyclic graph illustrating potential relationships among tumor-level risk factors, patient-level factors, surgical technique, provider or institutional practice patterns, treatment selection, radiation delivery, recurrence, salvage therapy, and disease-specific survival in localized MCC. Tumor size, lymphovascular invasion, margin status, primary site, SLNB status, immunosuppression, frailty, and surgical technique may influence both the likelihood of receiving adjuvant RT and the risk of recurrence, creating confounding by indication. Provider or institutional practice patterns may also influence surgical technique, RT fields, RT fractionation, and surveillance intensity. The red arrow represents the causal effect of interest: the effect of adjuvant RT on recurrence.

DAGs are useful because they clarify the estimand and help identify adjustment sets: the variables that should be conditioned on to estimate a causal effect under the assumptions encoded in the graph37,38. For example, if tumor size, lymphovascular invasion, primary site, margin status, and immunosuppression influence both RT selection and recurrence, then a naïve comparison of recurrence rates between patients who did and did not receive RT will be biased. Adjustment for these common causes is necessary to estimate the total effect of RT on recurrence. By contrast, variables that occur downstream of RT, such as salvage therapy or post-recurrence survival, should not be adjusted for when estimating the total effect of RT on recurrence, because doing so may block part of the causal pathway or introduce bias.

This distinction is particularly important for multi-institutional MCC evidence generation. A DAG does not prove that the assumed causal structure is correct. Rather, it makes assumptions explicit, transparent, and open to critique. By specifying the causal model before data analysis begins, investigators can determine which variables must be collected, which variables belong in the primary adjustment set, and which variables should be treated as mediators, colliders, post-treatment variables, or downstream outcomes rather than confounders.

Importantly, this approach shifts the focus from asking whether RT and recurrence are associated to estimating the causal effect of RT after accounting for the factors that influenced treatment selection in the first place. For a disease in which treatment decisions are highly individualized, this distinction is critical. Causal inference does not eliminate uncertainty, but it provides a disciplined framework for making that uncertainty explicit and for asking whether the available data can support the clinical question being posed.

Why Probabilistic Models May Be Useful in a Rare Disease

Even with prospective data collection, MCC presents a fundamental inferential challenge: local recurrence is uncommon. Sparse events limit precision, and multicenter studies introduce additional heterogeneity because referral patterns, surgical practice, RT thresholds, fractionation schedules, surveillance intensity, and salvage approaches differ across institutions. Consequently, even thoughtfully designed studies may not produce definitive answers to every clinically important question.

The challenge is therefore not simply to estimate whether RT works on average. Rather, it is to estimate how likely RT is to provide a clinically meaningful benefit for an individual patient. Probabilistic modelsb align naturally with this decision problem because they estimate probabilities that map directly onto clinical decisions. Rather than asking only whether a treatment effect crosses an arbitrary threshold for statistical significance, they estimate the probability that benefit exceeds thresholds clinicians and patients consider meaningful.

In practical terms, this shifts the question from “Was RT statistically significant?” to “How likely is RT to help enough to justify treatment?” A model could estimate, for example, the probability that RT reduces recurrence risk by more than 1%, 3%, or 5%. Those probabilities can then be considered alongside treatment burden, toxicity, salvageability, patient preferences, and competing risks to support individualized treatment decisions.

In a rare disease, a study may not produce enough events to answer every question with a definitive yes or no.

A Bayesian analysis can still be useful because it estimates probabilities that map directly onto clinical decisions. Instead of asking only whether RT is “statistically significant,” it asks: What is the probability that RT lowers recurrence risk by enough to matter?

This does not eliminate uncertainty. It quantifies uncertainty in a way that is transparent, clinically interpretable, and directly relevant to patient care.

A multicenter MCC study naturally includes several levels of variation: patients differ from one another, and institutions differ in case mix, referral patterns, and treatment practice. A hierarchical Bayesian model provides a way to represent this structure in the analysis. Rather than assuming that all institutions are identical, the model can allow baseline recurrence risk, treatment patterns, and, where data permit, treatment-associated effects to vary across centers.

This approach is particularly useful in rare diseases, where many centers or clinically important subgroups may contribute relatively few patients or events. Through partial pooling, estimates from smaller or sparser groups are informed by the broader multicenter experience, while still allowing larger centers or more data-rich groups to contribute their own information. This helps avoid overfitting to sparse subsets and can produce estimates that are more stable than analyzing each center separately, while avoiding the oversimplification of collapsing all patients into a single undifferentiated dataset.

Hierarchical Bayesian models can also incorporate prior information, including from existing observational studies, and update estimates as new data accrue. This feature is particularly attractive in MCC, where decades of institutional experience already exist but uncertainty remains substantial. At the same time, prior information must be incorporated transparently. Posterior estimates may be sensitive to prior assumptions, particularly when historical evidence arises from studies affected by treatment-selection bias. Accordingly, Bayesian analyses should prespecify prior distributions, justify their selection, and evaluate sensitivity to alternative assumptions.

Importantly, Bayesian modeling does not eliminate confounding, selection bias, or sparse data. Rather, its value is that it expresses the remaining uncertainty in a form that directly informs clinical decision-making. Instead of asking whether uncertainty has been eliminated, clinicians can ask whether the probability of clinically meaningful benefit is sufficient to justify treatment for a particular patient. In this way, probabilistic models help translate imperfect evidence into transparent, individualized decisions.

Visualizing Probabilistic Benefit

This framework addresses a practical question that arose repeatedly throughout our review and discussion: how can rigorous clinical inference continue when important outcomes are uncommon and definitive studies may remain difficult to conduct? One answer is to report decision-relevant probabilities rather than relying solely on a single treatment-effect estimate or a statistical-significance threshold.

The example in Figure 21 is hypothetical rather than derived from observed data. It demonstrates how a prospective hierarchical Bayesian model could estimate the posterior distribution of absolute recurrence reduction with adjuvant RT and summarize the probability that benefit exceeds clinically meaningful thresholds such as 1%, 3%, or 5%. These probabilities can then be interpreted alongside treatment burden, toxicity, salvageability, patient preferences, and other patient-specific considerations. The goal is not to eliminate uncertainty, but to quantify it honestly so that clinicians and patients can make better-informed decisions.

Figure 21: Illustrative posterior distributions and threshold probabilities, derived from hypothetical rather than actual data, showing how a Bayesian model could estimate clinically meaningful benefit from adjuvant RT. The upper panel shows posterior distributions for the absolute reduction in recurrence risk across lower- and higher-risk patient profiles, with dashed vertical lines marking prespecified benefit thresholds of 1%, 3%, and 5%. The lower panel shows the posterior probability that benefit exceeds each of these thresholds. This framework shifts the question from whether RT is statistically significant on average to the probability that RT provides a benefit large enough to matter clinically for a given patient risk profile. ARR, absolute risk reduction.

Bayesian Assurance and Study Feasibility

The same probabilistic framework can also be applied before a study begins to evaluate whether a proposed design is likely to produce clinically meaningful evidence. For the illustrative simulations in Table 2, we asked how large a prospective cohort would need to be for a Bayesian analysis to have at least 80% assurance of reaching an actionable conclusion, defined as a greater than 80% posterior probability that the absolute recurrence reduction (ARR) exceeded 3%.

This framing differs from conventional power calculations. Statistical power asks how often a study will detect an effect of a specified size under repeated sampling. Bayesian assurance instead asks how often a proposed study is expected to produce a clinically actionable posterior conclusion, explicitly incorporating uncertainty in the true treatment effect before the study begins. As a result, assurance shifts the emphasis from achieving statistical significance to generating evidence that is sufficiently informative to support clinical decision-making.

Table 2. Bayesian Assurance for Clinically Specified Local-Recurrence Scenarios
Approximate total cohort size needed for each separate comparison design to achieve 80% Bayesian assurance that ARR exceeds 3%.
No RT
1-year LR
RT
1-year LR
ARR Relative risk
reduction
Approximate total prospective cohort size, patients
(separate comparison designs; do not sum across columns)
Any RT (65%)
vs no RT (35%)
Conventional or
hypofractionated RT
(each 25%) vs no RT (35%)
Single-fraction RT (15%)
vs no RT (35%)
5% 1% 4% 80% 5,000 8,000 6,000
10% 5% 5% 50% 3,000 3,000 5,000
10% 2% 8% 80% 400 400 600
10% 1% 9% 90% 200 300 400
20% 10% 10% 50% 400 500 800
20% 5% 15% 75% 200 200 200
Values in the three rightmost columns are alternative design scenarios for separate pairwise comparisons and should not be summed across columns. Each value represents the approximate total prospective cohort size in patients needed to achieve 80% Bayesian assurance for that comparison design. Candidate total cohort sizes were evaluated from 100 to 10,000 patients.
Allocation assumptions are shown in the column headers and are expressed as percentages of the full prospective capture cohort. The conventional/hypofractionated column models conventional RT and hypofractionated RT as separate arms with identical 25% allocation and operating characteristics; it is not a pooled 50% multi-fraction comparison. For arm-specific comparisons, patients receiving other RT regimens are captured in the total cohort but do not contribute to that pairwise contrast; therefore, the reported values should be read as total cohort sizes, not per-arm sample sizes.
Scenarios with ARR ≤3% are not shown because a rule requiring Pr(ARR >3%) >80% is not expected to achieve 80% assurance when the true effect is at or below 3%. High-benefit scenarios, including 20% → 5%, are illustrative data-generating assumptions rather than predictions.

Table 2. Scenario-based Bayesian assurance calculations using local recurrence as the endpoint. Values represent the approximate total prospective cohort size required to achieve at least 80% Bayesian assurance under each clinically specified local-recurrence scenario. Bayesian assurance was defined as the proportion of simulated future studies that met a prespecified posterior decision rule. For these simulations, the clinically meaningful effect threshold was an absolute risk reduction (ARR) greater than 3%, and the posterior decision rule required greater than 80% posterior probability that ARR exceeded this threshold.

For each scenario, we specified assumed 1-year local recurrence risks with and without primary-site RT, which together define the corresponding ARR and relative risk reduction. Simulated studies were then generated across a prespecified grid of total cohort sizes ranging from 100 to 10,000 patients. Within each simulated study, patients were allocated according to prespecified comparison-specific proportions: 35% no RT, 65% any RT pooled, 25% conventional RT or hypofractionated RT, each modeled separately with identical allocation and operating characteristics, and 15% single-fraction RT. These allocation assumptions are part of the design being evaluated; different treatment distributions would produce different assurance estimates and sample-size requirements.

For each simulated study, beta-binomial posterior updating was used to estimate the probability that ARR exceeded the clinically meaningful threshold of 3%. A simulated study was considered clinically actionable if this posterior probability exceeded 80%. The table reports the smallest total cohort size at which at least 80% of simulated studies satisfied this decision rule.

These values should be interpreted as approximate operating characteristics of a proposed study design rather than exact sample-size requirements. Assurance summarizes how often a proposed design is expected to yield a clinically actionable posterior conclusion across many possible future studies, not whether any single realized study will necessarily do so. Individual studies conducted under the same design may still produce stronger or weaker posterior evidence because of sparse events, random variation, or differences in the observed distribution of patients across treatment groups.

Scenarios in which the true ARR is at or below the clinically meaningful threshold of 3% are not shown because increasing sample size would be expected to concentrate the posterior at or below that threshold, rather than increase the probability that ARR exceeds it. This observation illustrates an important principle of study design. When the untreated baseline risk of local recurrence is very low—as suggested by the Kavanagh cohort—a large absolute benefit from RT may simply not be biologically possible. In that setting, the scientific objective shifts. Rather than attempting to demonstrate a clinically meaningful absolute benefit, the more informative question becomes whether the available evidence can confidently exclude such a benefit for selected patients.

Framed this way, Bayesian assurance helps determine not only when a study is likely to demonstrate a clinically meaningful reduction in local recurrence, but also when it is likely to provide sufficiently informative evidence to support treatment de-escalation. However, because these calculations use local recurrence as the endpoint, they do not by themselves establish whether primary-site RT improves disease-specific survival or overall survival. This distinction is clinically important: if local recurrences are uncommon and effectively salvageable, then a reduction in recurrence may not translate into improved disease-specific outcomes. In that setting, the clinical value of routine adjuvant treatment should be judged not only by the magnitude of recurrence reduction, but also by salvageability, toxicity, treatment burden, competing risks, and patient preferences.

Although the scenario-based assurance calculations focus on local recurrence, the broader framework extends to any clinically relevant outcome. A curated multi-institutional evidence-generation effort should evaluate not only local recurrence, but also regional and locoregional recurrence, recurrence-free or progression-free survival, overall survival, salvageability, toxicity, treatment burden, and patient preferences. Local recurrence serves as a useful modeled endpoint because it is directly relevant to the Kavanagh study and to the intended effect of primary-site RT. More broadly, however, study feasibility for any endpoint depends on the underlying event rate, the expected absolute treatment effect, treatment allocation, and the degree of uncertainty investigators and clinicians are willing to accept.

This perspective differs from conventional sample-size or power calculations. Traditional power calculations ask how many patients are required to reject a null hypothesis under prespecified assumptions about the true treatment effect and type I error. Bayesian assurance instead asks how often a proposed study design is expected to produce evidence that is sufficiently informative to support a prespecified clinical decision. In this framework, the objective shifts from achieving statistical significance to generating clinically actionable evidence.

This distinction is particularly relevant in rare diseases. If the true absolute benefit of RT is large—for example, a reduction in 1-year local recurrence from 10% to 1%—a prospective study may produce informative evidence with only a few hundred patients. If the true ARR is closer to 5%, several thousand patients may be required, particularly for less commonly used treatment approaches such as single-fraction RT. Conversely, if the true ARR is closer to 1%, even a large study may be unlikely to demonstrate that benefit exceeds commonly accepted clinical thresholds, despite supporting the existence of a small treatment effect.

This illustrates why accurately characterizing the natural history of local recurrence is so important. The maximum achievable absolute benefit of RT is constrained by the baseline recurrence risk. If the untreated 1-year local recurrence risk is approximately 2%, as reported by Kavanagh et al., then an ARR greater than 3% is mathematically impossible for local recurrence, even if RT completely eliminated all local failures. In that setting, the scientific objective shifts from demonstrating a large absolute benefit to more precisely estimating baseline risk, identifying patients most likely to benefit, and determining how confidently clinically meaningful benefit can be excluded for selected low-risk patients. Even when a prespecified benefit threshold is unlikely to be exceeded, posterior probabilities remain informative by clarifying what is likely, unlikely, and still uncertain.

From Recurrence Rates to Decision Quality

Ultimately, the goal of future evidence generation is not simply to estimate local recurrence more precisely, but to improve clinical decision-making. Patients and clinicians need to know which individuals can safely omit RT, which patients are likely to benefit from treatment, whether abbreviated RT provides an acceptable balance of benefit and burden, and how recurrence risk should be weighed alongside salvageability, toxicity, competing risks, and patient preferences.

No single study is likely to eliminate uncertainty, particularly in a rare disease characterized by low event rates, evolving treatment practices, and heterogeneous clinical decision-making. That should not be interpreted as failure. The value of a study lies not only in producing definitive answers, but also in reducing uncertainty, refining clinically important estimates, and improving the decisions that clinicians and patients make together. By making assumptions explicit and quantifying uncertainty transparently, causal inference and probabilistic modeling provide a framework for learning even when definitive answers remain out of reach.

Conclusion

The study by Kavanagh et al. has prompted an important reexamination of adjuvant primary-site RT following complete excision of localized MCC. The low local recurrence rates reported in their cohort challenge prevailing assumptions, but the broader controversy extends beyond any single recurrence estimate. It reflects unresolved questions regarding baseline risk, treatment selection, clinically meaningful benefit, recurrence salvageability, treatment burden, and patient preferences.

The SoCO Journal Club demonstrated that these uncertainties are reflected in contemporary practice. Although clinicians broadly agree on many of the factors that influence risk, substantial variation remains in how recurrence risk is estimated, how treatment benefit is judged, and what level of uncertainty is acceptable when recommending RT.

Future evidence generation should therefore focus not simply on collecting more data, but on generating more informative evidence. Curated multi-institutional data collection, standardized recurrence definitions, explicit causal frameworks, and probabilistic models offer complementary approaches for estimating individualized treatment benefit while acknowledging the realities of studying a rare disease.

Perhaps the most important contribution of the Kavanagh study is not that it resolves the controversy surrounding adjuvant RT, but that it clarifies the uncertainty that remains. For rare diseases such as MCC, the perfect study may never arrive. Scientific rigor, then, is not the pursuit of certainty; it is the disciplined quantification of uncertainty. Progress comes from measuring uncertainty honestly, communicating it transparently, and incorporating it into better clinical decisions.

Materials and Methods

Platform and Writing Tools

This Perspectives on the Science article was produced using Quarto®, an open-source scientific and technical publishing system39. OpenAI GPT-5.5 was used to assist with manuscript drafting, narrative structure, editing, reference reconciliation, and code troubleshooting40. OpenEvidence was used to support literature review and verification of selected clinical and scientific concepts41. The authors reviewed, edited, and verified all AI-assisted content, analyses, citations, and interpretations, and take full responsibility for the accuracy, integrity, and final content of the manuscript.

Software Environment

Analyses and visualizations were conducted using R and RStudio42,43. Data manipulation and plotting used the tidyverse suite of packages44, including dplyr, tidyr, stringr, purrr, tibble, forcats, readr, and ggplot245. Tables were generated using gt46. Figures were exported as high-resolution PNG files and incorporated into the Quarto manuscript using knitr47.

Survey Analysis

Survey data were collected using REDCap electronic data capture tools hosted at Mass General Brigham48,49. Survey responses were summarized descriptively using counts and proportions. Survey visualizations were generated using custom R functions built around ggplot2 to maintain consistent formatting across figures. No formal hypothesis testing was performed for the survey analyses; results were interpreted descriptively as measures of practice variation, perceived uncertainty, and areas of clinical equipoise.

Literature Synthesis and Evidence Table

The manuscript summarizes evidence relevant to local recurrence after excision of localized Merkel cell carcinoma, with particular focus on the study by Kavanagh et al. and on the adverse clinicopathologic features used in postoperative radiation decision-making. The evidence table was developed to summarize selected studies informing NCCN adverse risk factors, including lymphovascular invasion, immunosuppression, primary tumor size, surgical margins, positive sentinel lymph node biopsy, and head and neck primary site. The table was intended as a narrative evidence synthesis rather than a formal systematic review or meta-analysis. Studies were included when they informed recurrence, sentinel lymph node positivity, disease-specific survival, overall survival, surgical margins, or anatomic risk, recognizing that not all studies evaluated local recurrence as a discrete endpoint.

Bayesian Assurance Calculations

Bayesian assurance calculations were performed to illustrate how prospective or highly curated multi-institutional evidence generation could be evaluated before study launch. These simulations were not intended to provide definitive sample-size requirements. Rather, they were designed to show how assumptions about baseline recurrence risk, true absolute risk reduction, treatment allocation, and posterior decision thresholds influence the feasibility of obtaining clinically actionable evidence.

The simulations considered baseline 1-year local recurrence risks without radiation therapy of 2%, 5%, and 10%, and assumed true ARR scenarios of 1%, 5%, and 9% where mathematically possible. For each scenario, the implied 1-year local recurrence risk with RT was calculated as the baseline risk minus the assumed ARR. Scenarios in which the assumed ARR exceeded the baseline recurrence risk were labeled “Not possible,” because they would imply a negative recurrence risk with RT.

For each candidate total cohort size, simulated studies were generated using prespecified allocation assumptions. No RT allocation was fixed at 35%. RT allocation varied by comparison: 65% for any RT pooled, 25% for conventional or hypofractionated RT, and 15% for single-fraction RT. These allocation assumptions are part of the design being evaluated; required total cohort sizes would differ under alternative allocation ratios. For example, if only 15% of patients were managed without RT, larger total cohort sizes would generally be required to estimate baseline recurrence risk and absolute benefit with comparable assurance.

For each simulated study, recurrence counts were generated under the assumed data-generating scenario, and beta-binomial posterior updating was used to estimate the posterior probability that ARR exceeded a clinically meaningful threshold. The primary illustrative decision rule required greater than 80% posterior probability that ARR exceeded 3%. Bayesian assurance was defined as the proportion of simulated future studies that met this posterior decision rule. Candidate total cohort sizes were evaluated up to 10,000.

Values in the assurance table should be interpreted as approximate design operating characteristics rather than exact sample-size requirements. A scenario that does not reach 80% assurance at a given total cohort size may still generate some individual simulated studies in which the posterior probability of clinically meaningful benefit is high. Conversely, other simulated studies under the same assumptions may generate weaker or inconclusive posterior evidence because of sparse events, random variation, or limited representation of one comparison group. Thus, assurance summarizes how often a proposed design would be expected to yield an actionable posterior conclusion across many possible future studies, not whether such a conclusion is possible in any single realized study.

# bayesian_assurance_mcc.R
#
# Bayesian assurance simulation for MCC local recurrence
# ------------------------------------------------------
# Revised Table 2 version.
#
# Main revision:
#   Table 2 is now organized around clinically specified recurrence-risk
#   scenarios rather than abstract baseline-risk / true-ARR combinations.
#   Rows show actual 1-year local recurrence risks with no RT and with RT,
#   along with ARR and relative risk reduction. Columns show approximate
#   total prospective cohort size needed for 80% Bayesian assurance under
#   prespecified treatment-allocation assumptions.
#
# Base decision rule for Table 2:
#   Pr(ARR > 3%) > 80%
#
# Assurance definition:
#   Across simulated future studies generated under each assumed true scenario,
#   assurance is the proportion of studies that meet the posterior decision rule.


# ============================================================
# Step 1: setup
# ============================================================

library(tidyverse)
library(glue)
library(scales)
library(gt)

theme_set(
  theme_minimal(base_size = 12) +
    theme(
      plot.title = element_text(face = "bold", hjust = 0.5),
      plot.subtitle = element_text(hjust = 0.5)
    )
)

set.seed(20260617)

# Helper for percent labels
pct_label <- function(x, accuracy = 1) {
  scales::percent(x, accuracy = accuracy)
}

# Candidate total prospective cohort sizes.
# Values are total patients in the prospective capture cohort, not per-arm N.
# If 80% assurance is not reached by the largest value in this grid, the table
# reports the result as >max(n_total_values).
n_total_values <- c(
  100, 200, 300, 400, 500, 600, 800, 1000, 1200,
  1600, 2000, 3000, 4000, 5000, 6000, 8000, 10000
)

# Primary posterior decision rule for Table 2.
arr_threshold_base <- 0.03
posterior_threshold_base <- 0.80
target_assurance_base <- 0.80

# Simulation settings.
# Increase these for final manuscript-quality runs if desired.
n_sim_table <- 200
n_draws_table <- 3000


# ============================================================
# Step 2: posterior and assurance functions
# ============================================================

# ------------------------------------------------------------
# Simulate one future study with unequal group sizes
# ------------------------------------------------------------

simulate_one_study_unequal <- function(
    baseline_lr,
    rt_lr,
    n_no_rt,
    n_rt
) {
  
  events_no_rt <- rbinom(
    n = 1,
    size = n_no_rt,
    prob = baseline_lr
  )
  
  events_rt <- rbinom(
    n = 1,
    size = n_rt,
    prob = rt_lr
  )
  
  tibble(
    n_no_rt = n_no_rt,
    n_rt = n_rt,
    events_no_rt = events_no_rt,
    events_rt = events_rt,
    observed_lr_no_rt = events_no_rt / n_no_rt,
    observed_lr_rt = events_rt / n_rt,
    observed_arr = observed_lr_no_rt - observed_lr_rt
  )
}


# ------------------------------------------------------------
# Bayesian posterior probability that ARR exceeds threshold
# ------------------------------------------------------------

posterior_prob_arr_gt <- function(
    events_no_rt,
    n_no_rt,
    events_rt,
    n_rt,
    arr_threshold = arr_threshold_base,
    n_draws = n_draws_table
) {
  
  # Weakly informative Jeffreys priors:
  # p_no_rt ~ Beta(0.5, 0.5)
  # p_rt    ~ Beta(0.5, 0.5)
  
  p_no_rt_draws <- rbeta(
    n_draws,
    shape1 = events_no_rt + 0.5,
    shape2 = n_no_rt - events_no_rt + 0.5
  )
  
  p_rt_draws <- rbeta(
    n_draws,
    shape1 = events_rt + 0.5,
    shape2 = n_rt - events_rt + 0.5
  )
  
  arr_draws <- p_no_rt_draws - p_rt_draws
  
  mean(arr_draws > arr_threshold)
}


# ------------------------------------------------------------
# Analyze one simulated study
# ------------------------------------------------------------

analyze_one_study <- function(
    study_data,
    arr_threshold = arr_threshold_base,
    posterior_threshold = posterior_threshold_base,
    n_draws = n_draws_table
) {
  
  post_prob <- posterior_prob_arr_gt(
    events_no_rt = study_data$events_no_rt,
    n_no_rt = study_data$n_no_rt,
    events_rt = study_data$events_rt,
    n_rt = study_data$n_rt,
    arr_threshold = arr_threshold,
    n_draws = n_draws
  )
  
  study_data |>
    mutate(
      arr_threshold = arr_threshold,
      posterior_prob_arr_gt_threshold = post_prob,
      posterior_threshold = posterior_threshold,
      actionable = posterior_prob_arr_gt_threshold > posterior_threshold
    )
}


# ------------------------------------------------------------
# Estimate assurance for unequal no-RT vs RT comparison
# ------------------------------------------------------------

estimate_assurance_unequal <- function(
    baseline_lr,
    rt_lr,
    n_no_rt,
    n_rt,
    arr_threshold = arr_threshold_base,
    posterior_threshold = posterior_threshold_base,
    n_sim = n_sim_table,
    n_draws = n_draws_table
) {
  
  sim_data <- purrr::map_dfr(
    seq_len(n_sim),
    \(i) {
      simulate_one_study_unequal(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        n_no_rt = n_no_rt,
        n_rt = n_rt
      ) |>
        analyze_one_study(
          arr_threshold = arr_threshold,
          posterior_threshold = posterior_threshold,
          n_draws = n_draws
        ) |>
        mutate(sim_id = i)
    }
  )
  
  sim_data |>
    summarize(
      assurance = mean(actionable),
      mean_observed_arr = mean(observed_arr),
      mean_posterior_prob = mean(posterior_prob_arr_gt_threshold),
      .groups = "drop"
    ) |>
    mutate(
      baseline_lr = baseline_lr,
      rt_lr = rt_lr,
      true_arr = baseline_lr - rt_lr,
      n_no_rt = n_no_rt,
      n_rt = n_rt,
      arr_threshold = arr_threshold,
      posterior_threshold = posterior_threshold,
      .before = 1
    )
}


# ------------------------------------------------------------
# Find minimum total cohort size for target assurance
# ------------------------------------------------------------

find_min_n_for_assurance <- function(
    baseline_lr,
    rt_lr,
    arr_threshold = arr_threshold_base,
    allocation_no_rt = 0.35,
    allocation_rt = 0.65,
    candidate_n_total = n_total_values,
    posterior_threshold = posterior_threshold_base,
    target_assurance = target_assurance_base,
    n_sim = n_sim_table,
    n_draws = n_draws_table
) {
  
  true_arr <- baseline_lr - rt_lr
  
  if (rt_lr < 0 || rt_lr > baseline_lr) {
    return(
      tibble(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        true_arr = true_arr,
        arr_threshold = arr_threshold,
        min_n_total = NA_real_,
        n_no_rt_at_min = NA_real_,
        n_rt_at_min = NA_real_,
        max_n_total = max(candidate_n_total),
        assurance_at_max_n = NA_real_,
        result_type = "Invalid scenario"
      )
    )
  }
  
  # If the assumed true ARR is at or below the decision threshold, increasing
  # sample size should concentrate the posterior at or below the threshold.
  # Such scenarios are not included in the main Table 2, but this guardrail
  # prevents misleading finite-sample artifacts from being reported as if they
  # were meaningful sample-size results.
  if (true_arr <= arr_threshold) {
    return(
      tibble(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        true_arr = true_arr,
        arr_threshold = arr_threshold,
        min_n_total = NA_real_,
        n_no_rt_at_min = NA_real_,
        n_rt_at_min = NA_real_,
        max_n_total = max(candidate_n_total),
        assurance_at_max_n = NA_real_,
        result_type = "ARR at or below threshold"
      )
    )
  }
  
  assurance_by_n <- purrr::map_dfr(
    candidate_n_total,
    \(n_total) {
      
      n_no_rt <- round(n_total * allocation_no_rt)
      n_rt <- round(n_total * allocation_rt)
      
      estimate_assurance_unequal(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        n_no_rt = n_no_rt,
        n_rt = n_rt,
        arr_threshold = arr_threshold,
        posterior_threshold = posterior_threshold,
        n_sim = n_sim,
        n_draws = n_draws
      ) |>
        mutate(n_total = n_total)
    }
  )
  
  reached <- assurance_by_n |>
    filter(assurance >= target_assurance) |>
    arrange(n_total) |>
    slice_head(n = 1)
  
  if (nrow(reached) == 0) {
    assurance_at_max_n <- assurance_by_n |>
      filter(n_total == max(candidate_n_total)) |>
      pull(assurance)
    
    return(
      tibble(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        true_arr = true_arr,
        arr_threshold = arr_threshold,
        min_n_total = NA_real_,
        n_no_rt_at_min = NA_real_,
        n_rt_at_min = NA_real_,
        max_n_total = max(candidate_n_total),
        assurance_at_max_n = assurance_at_max_n,
        result_type = "Not reached within range"
      )
    )
  }
  
  tibble(
    baseline_lr = baseline_lr,
    rt_lr = rt_lr,
    true_arr = true_arr,
    arr_threshold = arr_threshold,
    min_n_total = reached$n_total,
    n_no_rt_at_min = reached$n_no_rt,
    n_rt_at_min = reached$n_rt,
    max_n_total = max(candidate_n_total),
    assurance_at_max_n = assurance_by_n |>
      filter(n_total == max(candidate_n_total)) |>
      pull(assurance),
    result_type = "Reached"
  )
}


# ============================================================
# Step 3: clinically specified recurrence-risk scenarios
# ============================================================

# These scenarios intentionally use the recurrence risks clinicians discuss:
#   1-year LR with no RT versus 1-year LR with RT.
#
# The main table includes only scenarios in which the assumed true ARR exceeds
# the Table 2 decision threshold of ARR >3%. Lower-ARR scenarios remain in this
# tibble for transparency and optional sensitivity analyses, but they are not
# expected to support a decision rule requiring Pr(ARR >3%) >80%.

clinical_scenarios <- tibble::tribble(
  ~scenario_id, ~baseline_lr, ~rt_lr, ~include_in_main_table, ~clinical_note,
  "2_to_1",    0.02, 0.01, FALSE, "Low baseline-risk scenario; ARR is below the 3% decision threshold.",
  "5_to_2",    0.05, 0.02, FALSE, "Borderline scenario; ARR equals the 3% decision threshold.",
  "5_to_1",    0.05, 0.01, TRUE,  "Lower event-rate scenario with a 4% ARR.",
  "10_to_5",   0.10, 0.05, TRUE,  "Vern lower-bound scenario: 50% relative reduction and 5% ARR.",
  "10_to_2",   0.10, 0.02, TRUE,  "Higher-benefit scenario with 8% ARR.",
  "10_to_1",   0.10, 0.01, TRUE,  "Optimistic scenario with 9% ARR.",
  "20_to_10",  0.20, 0.10, TRUE,  "High-risk scenario with 50% relative reduction and 10% ARR.",
  "20_to_5",   0.20, 0.05, TRUE,  "Very high-benefit selected-risk scenario with 15% ARR."
) |>
  mutate(
    true_arr = baseline_lr - rt_lr,
    relative_reduction = true_arr / baseline_lr,
    scenario_label = glue(
      "No RT {pct_label(baseline_lr)} → RT {pct_label(rt_lr)}"
    )
  )

main_assurance_scenarios <- clinical_scenarios |>
  filter(include_in_main_table)


# ============================================================
# Step 4: treatment-allocation assumptions
# ============================================================

# Values are proportions of the total prospective capture cohort.
# For arm-specific comparisons, patients in other RT groups remain part of the
# total captured cohort but do not contribute to that specific pairwise contrast.

comparison_allocation <- tibble(
  comparison = c(
    "Any RT (65%) vs no RT (35%)",
    "Conventional or hypofractionated RT (each 25%) vs no RT (35%)",
    "Single-fraction RT (15%) vs no RT (35%)"
  ),
  allocation_no_rt = c(0.35, 0.35, 0.35),
  allocation_rt = c(0.65, 0.25, 0.15)
)

comparison_levels <- comparison_allocation$comparison


# ============================================================
# Step 5: calculate Table 2 assurance results
# ============================================================

assurance_input <- tidyr::crossing(
  scenario_id = main_assurance_scenarios$scenario_id,
  comparison = comparison_allocation$comparison
) |>
  left_join(main_assurance_scenarios, by = "scenario_id") |>
  left_join(comparison_allocation, by = "comparison") |>
  select(
    scenario_id,
    baseline_lr,
    rt_lr,
    clinical_note,
    true_arr,
    relative_reduction,
    scenario_label,
    comparison,
    allocation_no_rt,
    allocation_rt
  )

min_n_results_clinical <- assurance_input |>
  pmap_dfr(
    \(scenario_id,
      baseline_lr,
      rt_lr,
      clinical_note,
      true_arr,
      relative_reduction,
      scenario_label,
      comparison,
      allocation_no_rt,
      allocation_rt) {
      
      find_min_n_for_assurance(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        arr_threshold = arr_threshold_base,
        allocation_no_rt = allocation_no_rt,
        allocation_rt = allocation_rt,
        candidate_n_total = n_total_values,
        posterior_threshold = posterior_threshold_base,
        target_assurance = target_assurance_base,
        n_sim = n_sim_table,
        n_draws = n_draws_table
      ) |>
        mutate(
          scenario_id = scenario_id,
          scenario_label = scenario_label,
          clinical_note = clinical_note,
          relative_reduction = relative_reduction,
          comparison = comparison,
          allocation_no_rt = allocation_no_rt,
          allocation_rt = allocation_rt,
          .before = 1
        )
    }
  )


format_min_n_cell <- function(
    result_type,
    min_n_total,
    n_no_rt_at_min,
    n_rt_at_min,
    max_n_total
) {
  # Reached cells show the minimum tested total cohort size.
  # Not-reached cells show >max_n_total because 80% assurance was not
  # achieved within the tested grid.
  case_when(
    result_type == "Invalid scenario" ~ "Invalid scenario",
    result_type == "ARR at or below threshold" ~ "ARR ≤ threshold",
    result_type == "Not reached within range" ~ glue("<strong>&gt;{comma(max_n_total)}</strong>"),
    TRUE ~ glue("<strong>{comma(min_n_total)}</strong>")
  )
}


table2_clinical_assurance_table <- min_n_results_clinical |>
  mutate(
    no_rt_lr_label = pct_label(baseline_lr),
    rt_lr_label = pct_label(rt_lr),
    true_arr_label = pct_label(true_arr),
    relative_reduction_label = pct_label(relative_reduction),
    table_cell = format_min_n_cell(
      result_type = result_type,
      min_n_total = min_n_total,
      n_no_rt_at_min = n_no_rt_at_min,
      n_rt_at_min = n_rt_at_min,
      max_n_total = max_n_total
    )
  ) |>
  select(
    scenario_id,
    no_rt_lr_label,
    rt_lr_label,
    true_arr_label,
    relative_reduction_label,
    comparison,
    table_cell
  ) |>
  pivot_wider(
    names_from = comparison,
    values_from = table_cell
  ) |>
  arrange(
    factor(scenario_id, levels = main_assurance_scenarios$scenario_id)
  ) |>
  select(-scenario_id)


# ============================================================
# Step 6: build revised Table 2 gt object
# ============================================================

arm_specific_gt <- table2_clinical_assurance_table |>
  gt() |>
  tab_header(
    title = md("**Table 2. Bayesian Assurance for Clinically Specified Local-Recurrence Scenarios**"),
    subtitle = md(
      "Approximate total prospective cohort size needed for 80% Bayesian assurance that ARR exceeds 3%."
    )
  ) |>
  cols_label(
    no_rt_lr_label = md("**No RT<br>1-year LR**"),
    rt_lr_label = md("**RT<br>1-year LR**"),
    true_arr_label = md("**ARR**"),
    relative_reduction_label = md("**Relative risk<br>reduction**"),
    `Any RT (65%) vs no RT (35%)` = md("**Any RT (65%)<br>vs no RT (35%)**"),
    `Conventional or hypofractionated RT (each 25%) vs no RT (35%)` = md("**Conventional or<br>hypofractionated RT<br>(each 25%) vs no RT (35%)**"),
    `Single-fraction RT (15%) vs no RT (35%)` = md("**Single-fraction RT (15%)<br>vs no RT (35%)**")
  ) |>
  tab_spanner(
    label = md("**Approximate total prospective cohort size, patients**"),
    columns = all_of(comparison_levels)
  ) |>
  fmt_markdown(columns = everything()) |>
  cols_align(
    align = "center",
    columns = everything()
  ) |>
  cols_width(
    no_rt_lr_label ~ gt::pct(9),
    rt_lr_label ~ gt::pct(9),
    true_arr_label ~ gt::pct(7),
    relative_reduction_label ~ gt::pct(10),
    all_of(comparison_levels) ~ gt::pct(25)
  ) |>
  tab_source_note(
    source_note = md(
      glue(
        "Values show the approximate **total prospective cohort size in patients** needed to achieve 80% Bayesian assurance. Candidate total cohort sizes were evaluated from {comma(min(n_total_values))} to {comma(max(n_total_values))} patients."
      )
    )
  ) |>
  tab_source_note(
    source_note = md(
      "Allocation assumptions are shown in the column headers and are expressed as percentages of the full prospective capture cohort. The conventional/hypofractionated column models conventional RT and hypofractionated RT as separate arms with identical 25% allocation and operating characteristics; it is not a pooled 50% multi-fraction comparison. For arm-specific comparisons, patients receiving other RT regimens are captured in the total cohort but do not contribute to that pairwise contrast; therefore, the reported values should be read as total cohort sizes, not per-arm sample sizes."
    )
  ) |>
  tab_source_note(
    source_note = md(
      "Scenarios with ARR ≤3% are not shown because a rule requiring Pr(ARR >3%) >80% is not expected to achieve 80% assurance when the true effect is at or below 3%. High-benefit scenarios, including 20% → 5%, are illustrative data-generating assumptions rather than predictions."
    )
  ) |>
  tab_options(
    table.font.size = px(10.5),
    data_row.padding = px(3),
    column_labels.font.weight = "bold",
    table.width = gt::pct(100)
  )

arm_specific_gt


# ============================================================
# Step 7: optional heatmap for any RT pooled vs no RT
# ============================================================

# This heatmap uses the same clinically specified scenarios as Table 2.
# It shows assurance across candidate total cohort sizes for any RT pooled
# versus no RT under the decision rule Pr(ARR >3%) >80%.

n_sim_heatmap <- 300
n_draws_heatmap <- 1000

heatmap_grid <- tidyr::crossing(
  scenario_id = main_assurance_scenarios$scenario_id,
  n_total = n_total_values
) |>
  left_join(main_assurance_scenarios, by = "scenario_id") |>
  select(
    scenario_id,
    baseline_lr,
    rt_lr,
    clinical_note,
    true_arr,
    relative_reduction,
    scenario_label,
    n_total
  )

four_group_heatmap_results <- heatmap_grid |>
  pmap_dfr(
    \(scenario_id,
      baseline_lr,
      rt_lr,
      clinical_note,
      true_arr,
      relative_reduction,
      scenario_label,
      n_total) {
      
      n_no_rt <- round(n_total * 0.35)
      n_rt <- round(n_total * 0.65)
      
      estimate_assurance_unequal(
        baseline_lr = baseline_lr,
        rt_lr = rt_lr,
        n_no_rt = n_no_rt,
        n_rt = n_rt,
        arr_threshold = arr_threshold_base,
        posterior_threshold = posterior_threshold_base,
        n_sim = n_sim_heatmap,
        n_draws = n_draws_heatmap
      ) |>
        mutate(
          scenario_id = scenario_id,
          scenario_label = glue(
            "No RT {pct_label(baseline_lr)} → RT {pct_label(rt_lr)} ",
            "(ARR {pct_label(true_arr)})"
          ),
          n_total = n_total
        )
    }
  )

scenario_order <- main_assurance_scenarios |>
  mutate(
    scenario_label = glue(
      "No RT {pct_label(baseline_lr)} → RT {pct_label(rt_lr)} ",
      "(ARR {pct_label(true_arr)})"
    )
  ) |>
  pull(scenario_label)

four_group_assurance_heatmap <- four_group_heatmap_results |>
  mutate(
    scenario_label = factor(
      scenario_label,
      levels = scenario_order
    ),
    n_total_label = factor(
      n_total,
      levels = n_total_values
    )
  ) |>
  ggplot(aes(x = n_total_label, y = scenario_label, fill = assurance)) +
  geom_tile(color = "white") +
  geom_text(aes(label = percent(assurance, accuracy = 1)), size = 3.1) +
  scale_fill_gradient(
    labels = percent_format(accuracy = 1),
    limits = c(0, 1)
  ) +
  labs(
    x = "Total prospective cohort size",
    y = "Assumed true scenario",
    fill = "Assurance",
    title = "Bayesian assurance for any RT versus no RT",
    subtitle = "Allocation: 35% no RT, 65% any RT; decision rule Pr(ARR >3%) >80%"
  )

four_group_assurance_heatmap


# ============================================================
# Step 8: optional scenario-context table
# ============================================================

# This table is not intended as main Table 2. It documents the clinical
# scenarios that were considered but excluded from the main table because the
# assumed ARR is at or below the ARR >3% decision threshold.

below_threshold_scenarios <- clinical_scenarios |>
  filter(!include_in_main_table) |>
  transmute(
    `No RT 1-year LR` = pct_label(baseline_lr),
    `RT 1-year LR` = pct_label(rt_lr),
    ARR = pct_label(true_arr),
    `Relative reduction` = pct_label(relative_reduction),
    Interpretation = clinical_note
  )

below_threshold_scenarios_gt <- below_threshold_scenarios |>
  gt() |>
  tab_header(
    title = md("**Context Scenarios Not Included in Main Table 2**"),
    subtitle = md("These scenarios are clinically relevant but do not exceed the ARR >3% decision threshold used in Table 2.")
  ) |>
  cols_align(
    align = "center",
    columns = c(`No RT 1-year LR`, `RT 1-year LR`, ARR, `Relative reduction`)
  ) |>
  cols_align(
    align = "left",
    columns = Interpretation
  ) |>
  tab_options(
    table.font.size = px(11),
    data_row.padding = px(3),
    column_labels.font.weight = "bold",
    table.width = gt::pct(100)
  )

below_threshold_scenarios_gt

Bibliography

1.
Kavanagh, F. G. et al. Local recurrence following complete surgical excision of primary merkel cell carcinoma. Annals of Surgical Oncology (2025) doi:10.1245/s10434-025-18670-2.
2.
3.
4.
Petrelli, F. et al. Adjuvant radiotherapy for merkel cell carcinoma: A systematic review and meta-analysis. Radiotherapy and Oncology 134, 211–219 (2019).
5.
Shaw, J. H. F. & Rumball, E. Merkel cell tumour: Clinical behaviour and treatment. British Journal of Surgery 78, 138–142 (1991).
6.
Allen, P. J. et al. Merkel cell carcinoma: Prognosis and treatment of patients from a single institution. Journal of Clinical Oncology 23, 2300–2309 (2005).
7.
8.
9.
10.
11.
Conic, R. R. Z. et al. Sentinel lymph node biopsy in merkel cell carcinoma: Predictors of sentinel lymph node positivity and association with overall survival. Journal of the American Academy of Dermatology 81, 364–372 (2019).
12.
Bartlett, E. K. et al. Outcomes of patients undergoing sentinel lymph node biopsy for merkel cell carcinoma. Annals of Surgical Oncology (2025) doi:10.1245/s10434-025-18792-7.
13.
14.
Yusuf, M. B. et al. Immune status in merkel cell carcinoma: Relationships with clinical factors and independent prognostic value. Annals of Surgical Oncology 28, 6154–6165 (2021).
15.
Liang, E. et al. Merkel cell carcinoma analysis of outcomes: A 30-year experience. PLOS ONE 10, e0129476 (2015).
16.
17.
18.
Cullison, C. R., Zheng, D. X., Levoska, M. A., Scott, J. F. & Bordeaux, J. S. Tumor primary site as a prognostic factor for merkel cell carcinoma disease-specific death. Journal of the American Academy of Dermatology 85, 1259–1266 (2021).
19.
20.
21.
Lodde, G. C. et al. Clinical course of merkel cell carcinoma: A DeCOG multicenter study of 1049 patients. European Journal of Cancer 221, 115406 (2025).
22.
23.
Ahn, J., Howle, J. R. & Veness, M. J. Sentinel lymph node biopsy in clinically node-negative merkel cell carcinoma: An update of the westmead hospital experience. ANZ Journal of Surgery (2026) doi:10.1111/ans.70692.
24.
25.
26.
Wang, S. Y. et al. Accuracy of sentinel lymph node biopsy to predict nodal metastasis in head and neck merkel cell carcinoma. JAMA Otolaryngology–Head & Neck Surgery 152, 313–322 (2026).
27.
Messina, J. L. et al. Selective lymphadenectomy in patients with Merkel cell (cutaneous neuroendocrine) carcinoma. Annals of Surgical Oncology 4, 389–395 (1997).
28.
Gunaratne, D. A., Howle, J. R. & Veness, M. J. Sentinel lymph node biopsy in Merkel cell carcinoma: A 15-year institutional experience and statistical analysis of 721 reported cases. British Journal of Dermatology 174, 273–281 (2016).
29.
Strom, T. et al. Radiation therapy is associated with improved outcomes in merkel cell carcinoma. Annals of Surgical Oncology 23, 3572–3578 (2016).
30.
Mojica, P., Smith, D. & Ellenhorn, J. D. I. Adjuvant radiation therapy is associated with improved survival in merkel cell carcinoma of the skin. Journal of Clinical Oncology 25, 1043–1047 (2007).
31.
Lewis, K. G., Weinstock, M. A., Weaver, A. L. & Otley, C. C. Adjuvant local irradiation for merkel cell carcinoma. Archives of Dermatology 142, 693–700 (2006).
32.
33.
Alexander, N. A. et al. Increased risk of recurrence and disease-specific death following delayed postoperative radiation for merkel cell carcinoma. Journal of the American Academy of Dermatology 90, 261–268 (2024).
34.
Greenland, S., Pearl, J. & Robins, J. M. Causal diagrams for epidemiologic research. Epidemiology 10, 37–48 (1999).
35.
Pearl, J. Causality: Models, reasoning, and inference. (Cambridge University Press, 2009). doi:10.1017/CBO9780511803161.
36.
Hernan, M. A. & Robins, J. M. Causal inference: What if. (Chapman & Hall/CRC, 2020).
37.
Textor, J., Hardt, J. & Knuppel, S. DAGitty: A graphical tool for analyzing causal diagrams. Epidemiology 22, 745 (2011).
38.
Textor, J., Zander, B. van der, Gilthorpe, M. S., Liskiewicz, M. & Ellison, G. T. H. Robust causal inference using directed acyclic graphs: The R package dagitty. International Journal of Epidemiology 45, 1887–1894 (2016).
39.
Posit Software, PBC. Quarto. (2026).
40.
OpenAI. ChatGPT: OpenAI GPT-5.5. (2026).
41.
OpenEvidence. OpenEvidence. (2026).
42.
R Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, 2026).
43.
44.
Wickham, H. et al. Welcome to the tidyverse. Journal of Open Source Software 4, 1686 (2019).
45.
Wickham, H. ggplot2: Elegant graphics for data analysis. (Springer-Verlag New York, 2016). doi:10.1007/978-3-319-24277-4.
46.
47.
Xie, Y. Dynamic documents with r and knitr. (Chapman; Hall/CRC, 2015).
48.
49.
Harris, P. A. et al. The REDCap consortium: Building an international community of software platform partners. Journal of Biomedical Informatics 95, 103208 (2019).

NCCN Disclaimer

NCCN makes no warranties of any kind whatsoever regarding their content, use, or application and disclaims any responsibility for their application or use in any way.

Appendix

Abbreviations

AJCC, American Joint Committee on Cancer; ADMEC-O, Adjuvant Avelumab in Merkel Cell Carcinoma trial; APP, advanced practice provider; ARR, absolute risk reduction; CI, confidence interval; DAG, directed acyclic graph; DeCOG, Dermatologic Cooperative Oncology Group; DFS, disease-free survival; DSS, disease-specific survival; ENT, ear, nose, and throat; HR, hazard ratio; LN, lymph node; LR, local recurrence; LVI, lymphovascular invasion; MCC, Merkel cell carcinoma; MSKCC, Memorial Sloan Kettering Cancer Center; NCDB, National Cancer Database; NCCN, National Comprehensive Cancer Network; OR, odds ratio; OS, overall survival; Pr, probability; RFS, recurrence-free survival; RT, radiation therapy; SEER, Surveillance, Epidemiology, and End Results; SLN, sentinel lymph node; SLNB, sentinel lymph node biopsy; SoCO, Society of Cutaneous Oncology; VA, Veterans Affairs; WLE, wide local excision.

Acknowledgments

The authors thank the members of the Society of Cutaneous Oncology community who participated in the June 12, 2026 Journal Club. Their thoughtful discussion, clinical perspectives, and survey responses directly informed the questions, analyses, and recommendations developed in this Perspectives article.

Citation

For attribution, please cite this work as:

Miller DM, Emerick KS, Kaufman HL, McEvoy AM, Zaba L, Sharma A, Adelaja A, Silk AW, Khushalani NI, Brownell I, Coit DG, Nghiem P, Sondak VK. Local Recurrence After Excision of Merkel Cell Carcinoma: What Do We Actually Know, and What Should We Do Next? Journal of Cutaneous Oncology. 2026;4(2). https://doi.org/10.59449/joco.2026.07.16 Copied!

BibTeX citation:

@article{miller2026localrecurrence,
  author = {Miller, David Michael and Emerick, Kevin Scott and Kaufman, Howard Lane and McEvoy, Aubriana M. and Zaba, Lisa and Sharma, Ajay and Adelaja, Ade and Silk, Ann W. and Khushalani, Nikhil I. and Brownell, Isaac and Coit, Daniel G. and Nghiem, Paul and Sondak, Vernon K.},
  title = {Local Recurrence After Excision of Merkel Cell Carcinoma: What Do We Actually Know, and What Should We Do Next?},
  journal = {Journal of Cutaneous Oncology},
  year = {2026},
  volume = {4},
  number = {2},
  url = {https://journalofcutaneousoncology.io/perspectives/Vol_4_Issue_2/Local_Recurrence_After_Excision_MCC/},
  doi = {10.59449/joco.2026.07.16},
  issn = {2837-1933},
  publisher = {Society of Cutaneous Oncology},
  langid = {en}
}

Disclaimer

This site represents our opinions only. See our full Disclaimer

Disclosures

Conflict of Interests
Dr. Miller has received honoraria for serving as a consultant or participating on advisory boards for Almirall, Bristol Myers Squibb, Merck, EMD Serono, Regeneron, Sanofi Genzyme, Pfizer, Castle Biosciences, Checkpoint Therapeutics, and Sun Pharma. He serves on the Unloxcyt steering committee for Sun Pharma. He has stock options from Avstera Therapeutics. He has received institutional research funding from Regeneron, Kartos Therapeutics, NeoImmune Tech, Inc., Project Data Sphere, ECOG-ACRIN, and the American Skin Association.

Dr. Khushalani owns stock in Asensus Surgical, and Bellicum Pharmaceuticals. He participates in data and safety monitoring for AstraZeneca and Incyte Corporation. He has served as a consultant for Bristol Myers Squibb, Castle Biosciences, IDEAYA Biosciences, Immunocore, Instil Bio, IO Biotech, Iovance Biotherapeutics, Merck, Mural Oncology, MyCareGorithm, Nektar, Novartis, Regeneron Pharmaceuticals, Replimune, T-Knife Therapeutics, Delcath, and Sun Pharmaceuticals.

Dr. Nghiem reports compensation/support from UpToDate (honoraria), Almirall (advisory role), Incyte (institutional research funding), and has a patent pending for high-affinity T-cell receptors that target the Merkel polyomavirus, Patent filed: “Merkel cell polyomavirus T antigen-specific TCRs and uses thereof” (institution).

Dr. Kaufman is an employee of Ankyra Therapeutics. He has served in advisory roles for Castle Biosciences, Marengo Therapeutics, Tatum Biosciences, and Virogin Biotech; serves on the Board of Crichton Biosciences; holds stock in Replimune, Inc.; and has received honoraria from the Society for Immunotherapy of Cancer.

Dr. Sondak is a compensated consultant for Bristol Myers Squibb, Merck, and Mural Oncology, and receives research funding from Neogene Therapeutics and Skyline.

Dr. Emerick reports advisory board honoraria from Regeneron, royalties from UpToDate, and research funding from Haystack.

Dr. Silk reports receiving grants/research support to her institution from Regeneron and Merck, advisory board fees from Merck and Regeneron, and royalties from UpToDate, Inc.

Drs. Coit, McEvoy, Zaba, Sharma, Adelaja, and Brownell report no relevant disclosures.

License

This work is licensed under a creative commons BY-NC-ND license

Creative Commons License

Publication Stage

  • Published