Background {#Sec1}
==========

Identifying patients at increased risk for perioperative complications imperative for determining candidates for surgical intervention and proper preoperative patient counselling. This is particularly relevant in the head and neck oncology discipline where complications such as fistula development and free flap necrosis can result in significantly extended length of stay and decreased quality of life. Surgery-specific risk calculators are available for various surgical practices \[[@CR1]--[@CR4]\], and risk factors for cardiovascular and head and neck-specific complications for head and neck procedures including duration of anaesthesia and ASA are known \[[@CR5], [@CR6]\]; however, there are no head and neck-specific risk calculators available.

The universal surgical risk calculator created by the American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP) was developed as both a surgical aid and informed consent tool to improve the overall decision-making process. The risk calculator is an open-access online tool that uses an algorithm and validated data from over 500 hospitals and 2.7 million operations performed in the Unites States to predict the likelihood of 12 postoperative outcomes \[[@CR7]--[@CR11]\]. This online program accepts the input of 20 comorbidity and demographic-related, patient-specific variables, (such as age group, gender, smoking status, functional status, body-mass index, etc.) in conjunction with a surgery-specific Current Procedure Terminology (CPT) code, to predict patients' risk of 12 postoperative outcomes within 30 days after surgery \[[@CR12], [@CR13]\]. Patient-specific variables chosen for the risk calculator are broadly applicable for patients undergoing one of the 1887 procedures covered \[[@CR11], [@CR13], [@CR14]\].

Recently, this system has been evaluated in several surgical subspecialties including head and neck surgery \[[@CR8], [@CR15]--[@CR17]\], urology \[[@CR9]\], neurology \[[@CR10]\], pulmonology \[[@CR18]\], gynecology \[[@CR19]\], gastroenterology \[[@CR12], [@CR20]--[@CR23]\], orthopedics \[[@CR24]\] and general surgery with variable results. Despite these numerous publications, very few studies have examined the risk calculator's efficacy in a Canadian setting \[[@CR13], [@CR14]\] and none, so far, have evaluated head and neck surgery in a Canadian setting. We aim to provide the first comprehensive evaluation of the risk calculator for head and neck surgery in a Canadian setting by assessing its accuracy in predicting complications following six major head and neck procedures.

Our study evaluates the ACS NSQIP surgical risk calculator for predictive accuracy by comparing forecasted 30-day postoperative outcomes to observed incidences for patients who underwent one of six head and neck surgeries: total thyroidectomy, total laryngectomy, hemiglossectomy, partial glossectomy, laryngopharyngectomy and composite resection. These representative procedures were chosen for review because they are frequently performed by head and neck surgeons (thyroidectomy) or are associated with increased complexity and resultant higher complication rates and average length of stay (composite resection, laryngopharyngectomy). Separate statistical analyses were completed for the overall study population, and for each individual procedure. Furthermore, stratification and statistical analysis by free flap utilization, where patients were either separated into the Free Flap reconstructed (FF) group or Non-Free Flap reconstructed (NFF) group, was completed. This sub-stratification was warranted due to the risk calculator lacking separate CPT codes for head and neck procedures that included free flap reconstruction. This is an important distinction because it is well-established that free flap reconstruction increases patients' time under anesthesia, and it increases risk of surgical site infection and donor site complications \[[@CR5], [@CR25], [@CR26]\].

Methods {#Sec2}
=======

Study design {#Sec3}
------------

A retrospective review of all patients that underwent head and neck surgery at a single Canadian tertiary care referral center, between July 2015 and March 2017 was completed following research ethics board approval. Procedures included for review were total thyroidectomy, total laryngectomy, hemiglossectomy, partial glossectomy, laryngopharyngectomy and composite resection. Each patient's demographic-, procedure-, and complication- related data were collected and maintained in a secure database.

Demographic information was manually entered into the ACS NSQIP surgical risk calculator \[[@CR25]\] to fill 20 patient-specific risk factors including age group, gender, functional status, BMI, ASA class, hypertension, smoking status and COPD. The most relevant Current Procedural Terminology (CPT) codes were then selected based on the type, extent, and attributes of the procedure; the following guidelines were used, hemiglossectomy = 41135, partial glossectomy = 41120, laryngectomy = 31365, laryngopharyngectomy with free flap = 31395, composite resection with free flap = 41153 or 41155, thyroidectomy with central neck dissection = 60252, thyroidectomy without neck dissection = 60240. To maintain consistency, "Surgeon Adjustment of Risk" was not altered.

Patient-specific estimates of postoperative risk, determined by the risk calculator, including serious complications, any complications, pneumonia, cardiac complications, surgical site infection (SSI), urinary tract infection (UTI), venous thromboembolism (VTE), renal failure, readmission, discharge to nursing or rehabilitation facility, return to operating room (ROR), and death, were recorded. Complications were simply tallied, with no additional weight given to multiple complications in the same patient. Observed incidences of complication as determined by database review were compared  with the risk calculator's predictions.

Statistical analysis {#Sec4}
--------------------

Statistical analyses were completed for the overall study population, as well as for subset analysis by procedure type, and by presence or absence of free flap reconstruction (FF group and NFF group, respectively). Predicted and observed rates of incidences were then compared within each stratification. Brier score, ROC curves, likelihood ratios and Wilcoxon signed--ranked tests for nonparametric data were generated with SPSS, version 24.0.

Discrimination - ROC curves {#Sec5}
---------------------------

Evaluation of the ACS NSQIP surgical risk calculator for discrimination ability was completed using area under (AUC) the receiver operating characteristics (ROC) curve, sometimes called c-statistic. Receiver operating characteristics curve graphs model sensitivity (true positive predictions) against 1 -- specificity (false positive predictions) to evaluate how well a model distinguishes between higher and lower risk units within a population. The area under this graph is considered an accurate representation of a model's discrimination and is scored between 1.0 to 0.5, where the former is considered perfectly predictive and the latter is equivalent to chance. C-statistic values \> 0.7 are considered adequately discriminative whereas values of \> 0.8 are considered strongly discriminative \[[@CR9], [@CR10], [@CR26], [@CR27]\].

Calibration - brier score {#Sec6}
-------------------------

Calibration is a measure of how well a model's predictions 'fit' an observed incidence rate over a collection of predictions \[[@CR10]\]. The Brier score, which is the sum of the mean squared differences between predicted values and binary outcomes, was used to evaluate the ACS NSQIP surgical risk calculator for calibration. Incidences were given a value of 1 and non-incidences were given a value of 0. Brier score values range from 0.0 to 1.0, where smaller values represent increased calibration and accuracy. Several thresholds, varying from 0.01 to 0.16, have been suggested to quantify sufficient accuracy for the Brier score \[[@CR16], [@CR17], [@CR19], [@CR21]\]; 0.09 was used as it is considered accurate for outcomes with low incidences \[[@CR28], [@CR29]\].

Results {#Sec7}
=======

Patient characteristics {#Sec8}
-----------------------

Patients that underwent operations at a single, tertiary, Canadian teaching hospital between July 2015 and March 2017 were retrospectively reviewed for the purposes of this study. A total of 131 patients underwent operations, and 24 patients lacking adequate records or without 30-day follow-up information were excluded from the study, leaving a final population size of 107 individuals. Twenty-seven postoperative complications were observed within 30 days of surgery for the patients in this study (25%). The patient group had a mean age of 61 years, was composed approximately evenly of males and females (54, 53 years, respectively) and was, on average, overweight (Body Mass Index = 27.6 mg/kg^2^). Most patients were non-diabetic (87%) and exhibited a high incidence of hypertension (40%). There were no emergency cases (0%), and one patient was on dialysis before their operation (1%). The most commonly observed postoperative complications were ROR (12 incidences, 11%), SSI (11 incidences, 10%) and readmission (5 incidences, 5%). A full outline of demographic information can be found in Table [1](#Tab1){ref-type="table"}.Table 1Demographic and Medical InformationCharacteristicValue (*n* = 107)Age, mean (SD), y61.1 (13.4)Male gender, \# (%)54 (50)Height, mean (SD), m1.65 (0.14)Weight, mean (SD), kg74.1 (17.2)BMI, mean (SD)27.6 (6.65) Underweight, \# (%)5 (5) Normal weight, \# (%)29 (27) Overweight, \# (%)29 (27) Obese, \# (%)23 (21) Unknown, \# (%)21 (20)Functional Status, \# (%) Independent106 (99) Partially Dependent1 (1) Dependent0 (0)ASA Class, \# (%) I26 (24) II54 (50) III27 (25) IV0 (0)Diabetes, \# (%) No93 (87) Oral Medication4 (4) Insulin10 (9)Hypertension requiring medication, \# (%) Yes43 (40) No64 (60)Dyspnea, \# (%) No92 (86) With moderate exertion12 (11) At rest3 (3)Current Smoker (within 1 year), \# (%) Yes23 (21) No84 (79)History of severe COPD, \# (%) Yes13 (12) No94 (88)Dialysis, \# (%) Yes1 (1) No106 (99)Steroid use for chronic condition, \# (%) Yes8 (7) No99 (93)Emergency case: Yes, \# (%)0 (0)Ascites within 30 days prior to surgery: Yes, \# (%)0 (0)Sepsis within 48 h: Yes, \# (%)0 (0)Acute renal failure: Yes, \# (%)0 (0)Use of mechanical ventilation: Yes, \# (%)0 (0)Congestive Heart Failure 30 days prior to surgery: Yes, \# (%)0 (0)

Predicted vs. observed outcomes {#Sec9}
-------------------------------

Summary analyses demonstrating total observed and predicted incidences of complications without stratification were calculated (Fig. [1](#Fig1){ref-type="fig"}). The risk calculator exhibited 100% accuracy in predicting renal failure and adequate performance in predicting readmission (predicted = 6, observed = 5) and ROR (predicted = 9, observed = 12). However, it underestimated serious complications (35%), any complication (30%), cardiac complications (100%), SSI (36%), UTI (50%), and ROR (25%). Overestimations were equally common with pneumonia (50%), VTE, readmission (20%), death, and discharge to nursing and rehabilitation (133%), all scoring below predicted.Fig. 1Predicted vs. observed incidences for total study population

Stratification of the data based on use of free flap reconstruction revealed that the NSQIP- risk calculator was accurate in patients that did not receive microvascular reconstruction. Figure [2](#Fig2){ref-type="fig"} depicts the stratification of predicted complications as calculated by the ACS NSQIP surgical risk calculator compared to observed incidences for those patients who underwent free flap reconstruction (*n* = 58) and those that did not (*n* = 49). The risk calculator correctly predicted the number of complications in the NFF patient group for 11 of 12 outcomes, but exhibited a 100% overestimation for the remaining outcome---readmission (predicted = 2, observed = 1). For the patient cohort that underwent free flap reconstruction, the risk calculator demonstrated more limited accuracy, predicting only renal failure (predicted = 0, observed = 0), readmission (predicted = 4, observed = 4) and death (predicted = 0, observed = 0), within acceptable error. Six of the remaining nine outcomes were underestimated, while the rest were overestimated.Fig. 2Predicted vs. observed incidences by free flap reconstruction

Free flap vs. no free flap {#Sec10}
--------------------------

In our population, 58 (54%) patients underwent free flap reconstruction as part of their procedure. C-statistic and Brier score analyses were calculated for each of these stratifications. Overall, the risk calculator demonstrated improved calibration and comparable discrimination in predicting complications for the non-free flap reconstructed group when compared to the free flap reconstructed group (Table [2](#Tab2){ref-type="table"}). Values shown in bold are considered to meet or exceed the threshold values for their respective test statistics. Twelve of twelve Brier score values in the NFF group were below threshold (0.09) and 3 of 12 c-statistic values reached adequacy (0.70). Conversely, in the FF group, only 8 of 12 Brier score values showed acceptable calibration and 4 of 12 c-statistic values scored above 0.70. Specifically, SSI (Brier = 0.019, ROC = 0.885), readmission (Brier = 0.018, ROC = 1.00), and discharge to nursing or rehabilitation (Brier = 0.023, ROC = 0.854) exhibited adequate discrimination and calibration in the NFF group, and pneumonia (Brier = 0.017, ROC = 0.912), cardiac complication (Brier = 0.051, ROC = 0.767), readmission (Brier = 0.063, ROC = 0.713), and discharge to nursing or rehabilitation (Brier = 0.028, ROC = 0.987) scored adequately in the FF group. C-statistic values shown as 'N/A' were incalculable due to an absence of observed complications for those outcomes.Table 2Brier score and ROC AUC by Risk Calculator outcomeWithout Free Flap^\*,†,a^With Free Flap^\*,†,a^Combined^\*,†,a^OutcomeBrier ScoreROCBrier ScoreROCBrier ScoreROCSerious Complications**0.0838**0.52200.25080.60600.1743**0.7070**Any Complications**0.0861**0.55000.24740.62800.1735**0.7190**Pneumonia**0.0208**0.3230**0.01710.91200.0188**0.5430Cardiac Complications**0.0000**N/A**0.05090.76700.02760.8380**SSI**0.01930.8850**0.14530.6830**0.08760.7900**UTI**0.0001**N/A**0.0339**0.4420**0.0184**0.4980VTE**0.0001**N/A**0.0002**N/A**0.0001**N/ARenal Failure**0.0000**N/A**0.0000**N/A**0.0000**N/AReadmission**0.01831.00000.06270.71300.0445**0.6040ROR**0.0234**0.21900.14060.6880**0.08690.7500**Death**0.0000**N/A**0.0001**N/A**0.0001**N/AD/C SNF/Rehab**0.02320.85400.02810.98700.02590.8510**\*Brier scores of \< 0.09 and ROC \> 0.7 (bold values) were considered accurate^a^N/A indicates no ROC score could be calculated due to the absence of complications^†^Bold values meet or exceed their respective accuracy thresholds

When considering the entire study population without stratification by free flap, the risk calculator benefited from the effects of increased heterogeneity on discrimination, with 6 of 12 outcomes scoring above 0.70. Outcomes that scored satisfactorily for both discrimination and calibration were: cardiac complication (Brier = 0.028, ROC = 0.838), SSI (Brier = 0.088, ROC = 0.790), ROR (Brier = 0.087, ROC = 0.750), and discharge to nursing or rehabilitation (Brier = 0.026, ROC = 0.851). Two of the remaining outcomes, serious complications and any complications, scored acceptably for discrimination (ROC = 0.707 and 0.719, respectively), but inadequately for calibration (Brier = 0.174 and 0.174, respectively). Conversely, pneumonia, UTI, and readmission met calibration requirements (Brier = 0.019, 0.018, and 0.045, respectively), but did not score above the threshold for discrimination (ROC = 0.543, 0.498, and 0.604, respectively). Incidences of death, venous thromboembolism or renal failure were not observed; thus, c-statistics for these outcomes could not be calculated.

Stratification by procedure {#Sec11}
---------------------------

Evaluation of the risk calculator was also stratified by procedure (Table [3](#Tab3){ref-type="table"}) where patients were separated into one of the following groups: hemiglossectomy with free flap, partial glossectomy, laryngectomy, laryngopharyngectomy with free flap, composite resection with free flap, and thyroidectomy. Thyroidectomy was further subdivided into procedures with central neck dissection and those without. Table [3](#Tab3){ref-type="table"} indicates the number of patients in each stratum, the percentage of patients with stage IV cancer for cases with available staging information (*n* = 89), and the instances of preoperative radiation. One patient was excluded from analysis as a result of procedural characteristics not adequately fitting into the above groups. The occurrences of preoperative radiation (Likelihood Ratio = 21.689, *p* = 0.003) and the incidences of stage IV cancer (Likelihood Ratio = 49.314, *p* = 0.001) were significantly different between groups. This is to be expected considering the different standards of care associated with different disease locations and extensions. The risk calculator demonstrated adequate calibration and discrimination for laryngopharyngectomy with free flap (Brier = 0.073, ROC = 0.842), composite resection with free flap (Brier = 0.081, ROC = 0.842), thyroidectomy overall (Brier = 0.024, ROC = 0.817), and thyroidectomy with central neck dissection (Brier = 0.023, ROC = 0.816) or without central neck dissection (Brier = 0.027, ROC = 0.845). The risk calculator's ability to predict complications secondary to laryngectomy showed excellent discrimination (ROC = 0.956) but fell just short of acceptable calibration (Brier = 0.095). The remaining procedures, hemiglossectomy with free flap and partial glossectomy, exhibited acceptable calibration (Brier = 0.003 and 0.001, respectively) but, due to lack of observed complications, had incalculable c-statistics. Overall, the risk calculator performed exceptionally well in predicting complications by procedure; however, this is biased by the few patients in most of the groups.Table 3Brier Score, ROC AUC, and demographic information by procedure typeProcedureNBrier Score^\*,†^ROC^a,†^% T4aRadiated PatientsCPT Code(s)Hemiglossectomy with Free Flap4**0.0032**N/A25%041,135Partial Glossectomy12**0.0015**N/A0%041,120Laryngectomy50.0954**0.956**100%331,365Laryngopharyngectomy with Free Flap16**0.07320.842**81%231,395Composite Resection with Free Flap39**0.08070.842**30%441,153 & 41,155Thyroidectomy Overall30**0.02410.817**7%060,252 & 60,240Thyroidectomy with neck dissection21**0.02290.816**7%060,252Thyroidectomy without neck dissection9**0.02730.845**0%060,240\*Brier scores of \< 0.09 and ROC \> 0.7 (bold values) were considered accurate^a^N/A indicates no ROC score could be calculated due to the absence of complications^†^Bold values meet or exceed their respective accuracy thresholds

Overall {#Sec12}
-------

The ACS NSQIP surgical risk calculator's overall discrimination and calibration scores were calculated by evaluating all predictions of complications (12 predictions for each patient) together; this was completed for both FF and NFF groups, as well as, for the entire study population (Table [4](#Tab4){ref-type="table"}). In all three instances the risk calculator exhibited accurate discrimination and calibration: FF (Brier = 0.081, ROC = 0.85), NFF (Brier = 0.023, ROC = 0.78), and combined (Brier = 0.055, ROC = 0.86).Table 4Brier score and ROC AUC for all complications calculated togetherWithout Free Flap^a,b^With Free Flap^a,b^Combined^a,b^OutcomeBrier ScoreROCBrier ScoreROCBrier ScoreROCOverall**0.02300.77900.08140.85200.05480.8590**^a^Brier scores of \< 0.09 and ROC \> 0.7 (bold values) were considered accurate^b^Bold values meet or exceed their respective accuracy thresholds

Length of stay {#Sec13}
--------------

Length of stay analyses were completed for the overall study population and for both NFF and FF stratifications. Analyses were completed using the Wilcoxon signed--ranked test for nonparametric data \[[@CR12], [@CR22]\] and one patient was excluded from analysis due to ongoing hospital stay at time of manuscript development. Figure [3](#Fig3){ref-type="fig"} compares the box plots of observed and predicted lengths of stay for each sub-stratification and the total study population. Predicted lengths of stay for the total study population and for the FF stratified group were found to be significantly different than observed lengths of stay (*p* = 0.001). Conversely, predicted length of stay for the NFF stratified group was not found to be significantly different from observed length of stay (*p* = 0.764).Fig. 3Comparison between predicted and observed length of hospital stay for study population and free flap stratifications. Footnote: \*Observed and predicted lengths of stay were found to significantly differ

Discussion {#Sec14}
==========

Improving the informed consent process, providing surgeons with essential preoperative information, and acting as a quality control measure, are all potential benefits obtainable from an accurate surgical risk calculator. For these reasons, and others, various subspecialties have evaluated the efficacy of the ACS NSQIP surgical risk calculator in practice. Results from these studies are markedly inconsistent, with conclusions for the risk calculator ranging from entirely non-predictive for all outcomes \[[@CR9]\] to useful within judicious parameters \[[@CR22]\].

One potential cause for these varying results is lack of agreement on a threshold value for the Brier Score. Several studies consider scores below 0.01 as predictive \[[@CR9], [@CR16]\], referencing Bilimoria et al. and Cohen et al. as justification for this decision. Bilimoria et al., however, does not use a Brier score threshold of 0.01 in their publication and interestingly, considered a value of 0.069 as sufficiently accurate \[[@CR30]\]. Additionally, Cohen et al. does not at all use, nor refer to, the Brier score test in their analysis \[[@CR31]\]. Other published threshold values are much higher, including 0.16 \[[@CR21]\] and 0.14 \[[@CR22]\]. Due to these discrepancies, we decided to use 0.09, a relatively intermediate value that has been shown to indicate accuracy for models with low incidence rates \[[@CR28], [@CR29]\].

Some studies have found that the surgical risk calculator systemically underestimates the true risk of postoperative complications \[[@CR20], [@CR32]\]. Lack of surgery-specific metrics and predicted outcomes inherent to the risk calculator are often implicated as a cause of this inaccuracy \[[@CR32]\]. For this reason, many studies have tested the validity of supplementary specialty-specific preoperative metrics to support the 21 already captured by the ACS NSQIP surgical risk calculator program \[[@CR12], [@CR23], [@CR33]\]. As opposed to evaluating additional head and neck --specific preoperative metrics, we decided stratification by surgical complexity, free flap reconstruction versus no free flap reconstruction, and surgical type, as was previously unstudied, could delineate the risk calculator\'s accuracy as a function of surgical complexity. Initial evaluation of the risk calculator for observed versus predicted events, without stratification, exhibited results similar to other studies with 5 of 12 and 6 of 12 outcomes overestimated \[[@CR9], [@CR23]\] and underestimated, respectively. However, upon stratification by free flap utilization, the primary driver of these inaccuracies was elucidated. In the FF group, the risk calculator over- or under-predicted 9 of 12 outcomes. In the NFF group, however, outcomes were exceptionally well predicted, with 11 of 12 complications accurately projected by the risk calculator. We believe increased technical difficulty, potential sacrifice of vital structures (i.e. internal jugular vein, spinal accessory nerve, phrenic nerve, sympathetic trunk, and the thoracic duct), longer duration of anaesthesia, and extended post-operative stay, are potential contributors to the disparity between these stratifications.

Statistical results similar to that of predicted and observed outcomes was demonstrated with LOS; length of stay predictions for the overall population and FF sub-group were found to significantly differ from observed postoperative stays. However, when LOS analysis was completed for the NFF stratified group, predicted and observed stays were not found to significantly differ. Furthermore, Fig. [3](#Fig3){ref-type="fig"} clearly depicts the dissimilarity between predicted and observed LOS for the FF group and overall population, when compared to that of the NFF group. Our overall population and the FF sub-group results are consistent with many other publications, which found little to no association between observed and predicted length of stay \[[@CR12], [@CR15]--[@CR17], [@CR21], [@CR22]\]. Massoumi et al. implicated the level of care offered by a hospital (primary vs. tertiary) and its practices as potential generators of these inconsistencies \[[@CR20]\].

Our results follow a general trend present in the literature: increasingly complex procedures, where postoperative outcomes are closely influenced by the case-by-case features of the procedure, the skill level of the surgeon and postoperative care practices, are more likely miss-projected by the ACS NSQIP risk calculator, which cannot consider surgical intricacies, surgeon experience or hospital-specific acute care standards \[[@CR8], [@CR20]\]. An essential first step in mitigating this inaccuracy, and a function currently missing from the ACS NSQIP surgical risk calculator, is the combination of CPT codes to build a case more representative of the actual procedure \[[@CR32]\].

Our analysis of the risk calculator for calibration and discrimination followed a similar trend as LOS with all 12 outcomes properly calibrated in the NFF group and 8 of 12 calibrated in the FF group. As mentioned by Winoker et al., it is important to emphasize that few incidences of observed complications can contribute to low Brier scores \[[@CR9]\]. This can perhaps be seen in outcomes such as renal failure, VTE, and UTI in our cohort. Interestingly, the risk calculator exhibited slightly better discrimination for outcomes in the FF group (4 of 12) than in the NFF group (3 of 12). One possible explanation for this phenomenon is the improved performance of models that exhibited higher levels of heterogeneity. As shown by Cohen et al., ACS NSQIP surgical risk calculator estimations, projected for heterogeneous populations of patients (those that undergo many different procedures and have a diversity of preoperative morbidities), score better for discrimination than risk projections for homogenous patient populations \[[@CR9]\]. The risk calculator scored well for both calibrated and discrimination for 6 of 12 outcomes in the total study population. Our results demonstrating the risk calculator's accuracy for a portion of complications predicted are consistent with that of previous publications \[[@CR21], [@CR22], [@CR33]\]; however, they are inconsistent with other evaluations specific to the head and neck discipline. Two potential reasons for this discrepancy are differences in Brier score threshold values used and patient population heterogeneity. The former, explained in the discussion above, has resulted in previous studies considering acceptable scores as non-predictive \[[@CR16], [@CR17]\]. The latter refers to previous studies evaluating only one head and neck surgical procedure and thus exhibiting lower levels of heterogeneity \[[@CR15], [@CR16]\].

Analysis by procedure type further substantiated the risk calculator's accuracy in our patient cohort, with 7 of 8 evaluated procedures showing good calibration and all calculable procedures meeting adequate discrimination levels. Potential reasons why the risk calculator could have shown poor calibration for laryngectomy, which scored just above the Brier threshold, includes relatively few patients (*n* = 5) having undergone this procedure and thus a potential for sampling error, and the relatively high incidences of T4 cancer (100%) and prior radiation (60%). Interestingly, of the procedures analyzed that met the Brier score threshold, composite resection with free flap and laryngopharyngectomy with free flap had the highest values. This relative decrease in calibration can potentially be attributed to the increased complexity associated with these surgeries. This is certainly supported considering the results of the free flap stratified analyses.

For the overall population analyses, in which all surgical procedures and complications were pooled together, discrimination and calibration computations met their respective thresholds this provides the strongest evidence in our analysis of the risk calculator's effectiveness in predicting complications secondary to head and neck surgery.

Primary strengths {#Sec15}
-----------------

We present the first comprehensive evaluation of the ACS NSQIP surgical risk calculator in the head and neck discipline in a Canadian setting. Furthermore, to the authors' knowledge we are the first publication on the risk calculator to stratify by microvascular reconstruction, a surgical procedure which is associated with longer OR times and increased complication rates. Finally, our findings demonstrate good performance of the ACS NSQIP surgical risk calculator in predicting postoperative complications in patients undergoing head and neck surgery without free flap reconstruction. Given these results, we regularly utilize the risk calculator to determine individual patient's risk of complications as a decision-support tool to evaluate the safety of the procedure and to provide more comprehensive preoperative counselling to patients and their families.

Limitations {#Sec16}
-----------

We are primarily limited by the size of our study cohort (*n* = 107), having few patients with complications (27), and being a single institution retrospective review. The latter two limitations have been implicated by Cohen et al. as potential sources of inaccuracy when evaluating the risk calculator \[[@CR11]\]. Specifically, they concluded that three criteria must be met to fairly evaluate the risk calculator: a patient cohort acquired from more than one care center, a minimum of 100 patients with an incidence, and sufficient cohort heterogeneity. Based on the diverse list of procedures reviewed, we believe the level of heterogeneity required has been met in our cohort.

Next steps {#Sec17}
----------

Specialty-specific surgical risk calculators have shown improved predictive accuracy over the ACS NSQIP risk calculator \[[@CR12], [@CR23], [@CR32]\], suggesting development of a head and neck-specific risk calculator is warranted. One potential caveat of these additional metrics, as was pointed out by Cohen et al., is the potential for overfitting \[[@CR11]\]. Application of surgery-specific risk calculators should be used only on the procedures and settings which they have been examined. To this end, a multicenter collaboration as part of the Canadian Association of Head and Neck Surgeons is planned.

The purpose of a risk calculator is to determine the potential complications our patients may endure based on patient comorbidities and the procedures planned. The ACS NSQIP surgical risk calculator predicts the major complications of any surgery; however, it does not predict complications specific to our patient population that impair quality of life and increase LOS including fistula development and free flap failure requiring a second free flap or locoregional flap. A specialty-specific risk calculator would address this shortcoming.

Conclusion {#Sec18}
==========

The ACS NSQIP surgical risk calculator has the potential to act as a quality improvement metric and aid in the informed consent process through preoperative planning and postoperative prevention of potential mortality and morbidity. Though data collected for the risk calculator was amassed entirely from American hospitals with the intended use of improving the standard of care within the American health care system, our results suggest the potential utility of the risk calculator for Canadian Head and Neck oncology. Specifically, judicious application of the risk calculator for head and neck surgical procedures that do not involve microvascular reconstruction is supported.

ACS NSQIP

:   American College of Surgeons National Surgical Quality Improvement Program

AUC

:   Area under the curve

COPD

:   Chronic obstructive pulmonary disease

CPT

:   Current Procedural Terminology

FF

:   Free flap reconstructed patient group

LOS

:   Length of stay

NFF

:   Non-free flap reconstructed patient group

REB

:   Research ethics board

ROC

:   Receiver operating characteristics

ROR

:   Return to operating room

SSI

:   Surgical site infection

UTI

:   Urinary tract infection

VTE

:   Venous thromboembolism

The authors would like to acknowledge Michael Hong, BSc, for his contributions to the manuscript.

Funding {#FPar1}
=======

No forms of funding were acquired for this research.

Availability of data and materials {#FPar2}
==================================

Supporting data can be requested through contact with the corresponding author.

All authors were involved with the conception and design of the study, analysis and interpretation of data, revising the manuscript, and have given final approval. Further to this PV was involved with the design, direction and writing of the manuscript, including figures and statistical analyses. MO was involved in developing and editing the manuscript, completing the statistical analysis and figures, and editing the manuscript. KH and DE were the attending surgeons of the patients analyzed in this study.

Ethics approval and consent to participate {#FPar3}
==========================================

Research ethics approval was obtained from the Sunnybrook Health Sciences Center Research Ethics Board (REB). The following REB number applies 263--2017.

Consent for publication {#FPar4}
=======================

Not Applicable.

Competing interests {#FPar5}
===================

The authors declare that they have no competing interests.

Publisher's Note {#FPar6}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
