A cohort of 636 patients treated between June 2018 and April 2021 using the helical TomoTherapy® platform (Radixact) was collected, 512 planned using the Precision TPS while the remaining 124 using the RayStation TPS.
Three groups of complexity metrics were extracted using the in-house developed Matlab® (The MathWorks Inc, Natick, MA, USA) library TCoMX [1]:
Group A: 12 typical delivery parameters (e.g. Gantry Period, Treatment Time, etc)
Group B: 14 metrics proposed by Santos and colleagues [2]
Group C: 42 newly developed metrics describing the field geometry and beam modulation [3].
Additionally, 174 radiomics features (Group D) were extracted from the 2D sinograms using the IBSI compliant software S-IBEX [4].
Separate linear models were realized for the two TPSs to predict the 3%,2mm gamma index passing rate computed with global normalization and a 10% threshold [5]. Four different models were created considering four different sets of variables obtained by progressively adding one group at a time to Group A.
A variables pre selection based on the Pearson’s Correlation Coefficient (PCC) was applied to remove strongly correlated variables (|r|>0.8). Each model was trained 1000 times by randomly splitting the dataset into 75% training and 25% validation at each try. In each training phase, variables could be added (removed) progressively to (from) the model through a stepwise approach based on the p-value of an F-statistics. The starting point was always a constant model.
The mean absolute error (MAE) and PCC between the true and predicted gamma indexes were computed to measure the performance.