Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multivariate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

A comparison of computational models for predicting yield sooting index

Sooting propensity, a measurement of how much particulate matter is produced when a fuel is burned, is a property of significant interest among researchers who are striving to discover the next generation of cleaner, more efficient fuels and fuel additives. Many compounds are not viable as fuels and/or fuel additives, and as a result, designing cleaner-burning biofuels using only experimental techniques is inefficient. Predictive models have been instrumental in reducing this inherent difficulty, providing researchers with a tool to preemptively screen compounds before production and testing. The present work compares the accuracies and interpretabilities of existing models used to predict a particular measure of sooting propensity, Yield Sooting Index (YSI). These models include artificial neural networks, graph neural networks, and multivariate equations. A novel equation for predicting YSI based on atom path count and bond order is proposed, which can highlight key structural components that contribute to YSI. It was found that artificial neural networks slightly outperform graph neural networks and greatly outperform multivariate equations in blind (test set) prediction accuracy; however, graph neural networks and multivariate equations provide significantly more interpretability as to how compound structure relates to YSI. Predictions of YSI are compared to experimental measurements for previously un-tested compounds with cetane numbers comparable to diesel fuel (50-60) (butyl decanoate, ethyl decanoate, 1,4-bis(ethenoxymethyl)cyclohexane, and 5-heptyloxolan-2-one), and it was found that these compounds produce significantly less soot compared to diesel fuel.

09 BIOMASS FUELS↗

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nonlinear Decoupling Control with ANFIS-based Unmodeled Dynamics Compensation for A Class of Complex Industrial Processes

Complex industrial processes are multivariable and generally exhibit strong coupling among their control loops with heavy nonlinear nature that make it very difficult to obtain an accurate model. As a result, the conventional and data driven control methods are difficulty to apply. Using a twin-tank level control system as an example, a novel multivariable decupling control algorithm with adaptive neural fuzzy inference system (ANFIS) based unmodeled dynamics (UD) compensation is proposed in this paper for a class of complex industrial processes. At first, a nonlinear multivariable decupling controller with UD compensation is introduced. Different from the general method, the decomposition estimation algorithm by ANFIS is used to estimate the UD, and the desired estimating and decoupling control effects are achieved. Secondly, the proposed method does not require the complex switching mechanism which has been commonly used in the literature. This simplifies the obtained decoupling algorithm and its realization. Thirdly, based on some new lemmas and theorems, the conditions on the stability and convergence of the closed-loop system are analyzed to show the uniform boundedness of all the variables. Finally, experimental tests on a heavily coupled nonlinear twin-tank system are carried out to demonstrate the effectiveness and the practicability of the proposed method.

Complex industrial process, ANFIS, multivariable a↗

SU-F-R-07: Radiomics of CT Features and Associations and Correlation with Outcomes Following Lung SBRT

Objective: To investigate the prognostic significance of image gradients and in predicting clinical outcomes in a patients with non-small cell lung cancer treated with stereotactic body radiotherapy (SBRT) on 71 patients with 83 treated lesions. Methods: The records of patients treated with lung SBRT were retrospectively reviewed. When applicable, SBRT target volumes were modified to exclude any overlap with pleura, chestwall, or mediastinum. The ITK software package was utilized to generate quantitative measures of image intensity, inhomogeneity, shape morphology and first and second-order CT textures. Multivariate and univariate models containing CT features were generated to assess associations with clinicopathologic factors. Results: On univariate analysis, tumor size (HR 0.54, p=0.045) sumHU (HR 0.31, p=0.044) and short run grey level emphasis STD (HR 0.22, p=0.019) were associated with regional failure-free survival; meanHU (HR 0.30, p=0.035), long run emphasis (HR 0.21, p=0.011) and long run low grey level emphasis (HR 0.14, p=0.005) was associated with distant failure-free survival (DFFS). No features were significant on multivariate modeling however long run low grey level emphasis had a hazard ratio of 0.12 (p=0.061) for DFFS. Adenocarcinoma and squamous cell carcinoma differed with respect to long run emphasis STD (p=0.024), short run low grey level emphasis STD (p<0.001), and long run low grey level emphasis STD (p=0.024). Multivariate modeling of texture features associated with tumor histology was used to estimate histologies of 18 lesions treated without histologic confirmation. Of these, MVA suggested the same histology as a prior metachronous lung malignancy in 3/7 patients. Conclusion: Extracting radiomics features on clinical datasets was feasible with the ITK package with minimal effort to identify pre-treatment quantitative CT features with prognostic factors for distant control after lung SBRT.

60 APPLIED LIFE SCIENCES↗

TU-FG-201-05: Varian MPC as a Statistical Process Control Tool

Purpose: Quality assurance in radiotherapy requires the measurement of various machine parameters to ensure they remain within permitted values over time. In Truebeam release 2.0 the Machine Performance Check (MPC) was released allowing beam output and machine axis movements to be assessed in a single test. We aim to evaluate the Varian Machine Performance Check (MPC) as a tool for Statistical Process Control (SPC). Methods: Varian’s MPC tool was used on three Truebeam and one EDGE linac for a period of approximately one year. MPC was commissioned against independent systems. After this period the data were reviewed to determine whether or not the MPC was useful as a process control tool. Analyses on individual tests were analysed using Shewhart control plots, using Matlab for analysis. Principal component analysis was used to determine if a multivariate model was of any benefit in analysing the data. Results: Control charts were found to be useful to detect beam output changes, worn T-nuts and jaw calibration issues. Upper and lower control limits were defined at the 95% level. Multivariate SPC was performed using Principal Component Analysis. We found little evidence of clustering beyond that which might be naively expected such as beam uniformity and beam output. Whilst this makes multivariate analysis of little use it suggests that each test is giving independent information. Conclusion: The variety of independent parameters tested in MPC makes it a sensitive tool for routine machine QA. We have determined that using control charts in our QA programme would rapidly detect changes in machine performance. The use of control charts allows large quantities of tests to be performed on all linacs without visual inspection of all results. The use of control limits alerts users when data are inconsistent with previous measurements before they become out of specification. A. Carver has received a speaker’s honorarium from Varian.

07 ISOTOPE AND RADIATION SOURCES↗

Ensembles of radial basis function networks for spectroscopic detection of cervical precancer

The mortality related to cervical cancer can be substantially reduced through early detection and treatment. However, current detection techniques, such as Pap smear and colposcopy, fail to achieve a concurrently high sensitivity and specificity. In vivo fluorescence spectroscopy is a technique which quickly, noninvasively and quantitatively probes the biochemical and morphological changes that occur in precancerous tissue. A multivariate statistical algorithm was used to extract clinically useful information from tissue spectra acquired from 361 cervical sites from 95 patients at 337-, 380-, and 460-nm excitation wavelengths. The multivariate statistical analysis was also employed to reduce the number of fluorescence excitation-emission wavelength pairs required to discriminate healthy tissue samples from precancerous tissue samples. The use of connectionist methods such as multilayered perceptrons, radial basis function (RBF) networks, and ensembles of such networks was investigated. RBF ensemble algorithms based on fluorescence spectra potentially provide automated and near real-time implementation of precancer detection in the hands of nonexperts. The results are more reliable, direct, and accurate than those achieved by either human experts or multivariate statistical algorithms.

Cervix Neoplasms/pathology↗

External Influences on Modeled and Observed Cloud Trends

Understanding the cloud response to external forcing is a major challenge for climate science. This crucial goal is complicated by intermodel differences in simulating present and future cloud cover and by observational uncertainty. This is the first formal detection and attribution study of cloud changes over the satellite era. Presented herein are CMIP5 (Coupled Model Intercomparison Project - Phase 5) model-derived fingerprints of externally forced changes to three cloud properties: the latitudes at which the zonally averaged total cloud fraction (CLT) is maximized or minimized, the zonal average CLT at these latitudes, and the height of high clouds at these latitudes. By considering simultaneous changes in all three properties, the authors define a coherent multivariate fingerprint of cloud response to external forcing and use models from phase 5 of CMIP (CMIP5) to calculate the average time to detect these changes. It is found that given perfect satellite cloud observations beginning in 1983, the models indicate that a detectable multivariate signal should have already emerged. A search is then made for signals of external forcing in two observational datasets: ISCCP (International Satellite Cloud Climatology Project) and PATMOS-x (Advanced Very High Resolution Radiometer (AVHRR) Pathfinder Atmospheres - Extended). The datasets are both found to show a poleward migration of the zonal CLT pattern that is incompatible with forced CMIP5 models. Nevertheless, a detectable multivariate signal is predicted by models over the PATMOS-x time period and is indeed present in the dataset. Despite persistent observational uncertainties, these results present a strong case for continued efforts to improve these existing satellite observations, in addition to planning for new missions.

clouds (meteorology)↗

A Novel Approach to Detect Accelerated Aged and Surface-Mediated Degradation in Explosives by UPLC-ESI-MS.

A new approach was created for studying energetic material degradation. This approach involved detecting and tentatively identifying non-volatile chemical species by liquid chromatography-mass spectrometry (LC-MS) with multivariate statistical data analysis that form as the CL-20 energetic material thermally degraded. Multivariate data analysis showed clear separation and clustering of samples based on sample group: either pristine or aged material. Further analysis showed counter-clockwise trends in the principal components analysis (PCA), a type of multivariate data analysis, Scores plots. These trends may indicate that there was a discrete shift in the chemical markers as the went from pristine to aged material, and then again when the aged CL-20 mixed with a potentially incompatible material was thermally aged for 4, 6, or 9 months. This new approach to studying energetic material degradation should provide greater knowledge of potential degradation markers in these materials.

36 MATERIALS SCIENCE↗

Flexible classification with spatial quantile comparison and novel statistical features applied to spent nuclear fuel analysis

Multivariate classification algorithms are a common tool in contemporary data science that have begun to be applied to nuclear forensic analyses. Using a multivariate generalization of quantile–quantile plots for comparing unknown statistical distributions, a new model-free multivariate classifier called the Quantile–Quantile Comparator has been developed and tested on the analysis of simulated irradiated nuclear fuel. A method to estimate the sample-specific probability of having made an incorrect classification decision by analyzing the distribution of quantile comparisons has been developed, and methods have been developed to analyze the partition of the known data library used to analyze unknown data. By analyzing this partition, the analyst may identify classes that are likely to be confused, identify test data that are inconsistent with the data library and return a “None-of-the-Above” classification decision, and optimize the class boundaries. In this paper, the development of these capabilities is described, and a case study using real experimental assay data collected from nuclear fuel extracted from the newly-released In conclusion, Spent Fuel Isotopic Composition Database to demonstrate the efficacy of these capabilities is presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

pN0(i+) Breast Cancer: Treatment Patterns, Locoregional Recurrence, and Survival Outcomes

Purpose: To examine treatment patterns, recurrence, and survival outcomes in patients with pN0(i+) breast cancer. Methods and Materials: Subjects were 5999 women with AJCC (6th edition) pT1-3, pN0-N1a, M0 breast cancer diagnosed between 2003 and 2006. Of these, 4342 (72%) had pN0, 96 (2%) had pN0(i+), 349 (6%) had pNmic (micrometastases >0.2 mm to ≤2 mm), and 1212 (20%) had pN1a (1-3 positive macroscopic nodes) disease. Treatment characteristics and 5-year Kaplan-Meier local recurrence, regional recurrence (RR), locoregional recurrence (LRR), and overall survival were compared between nodal subgroups. Multivariable analysis was performed using Cox regression modeling. A 1:3 case-match analysis examined outcomes in pN0(i+) cases compared with pN0 controls matched for similar tumor and treatment characteristics. Results: Median follow-up was 4.8 years. Adjuvant systemic therapy use increased with nodal stage: 81%, 92%, 95%, and 94% in pN0, pN0(i+), pNmic, and pN1a disease, respectively (P<.001). Nodal radiation therapy (RT) use also increased with nodal stage: 1.7% in pN0, 27% in pN0(i+), 33% in pNmic, and 63% in pN1a cohorts (P<.001). Five-year Kaplan-Meier outcomes in pN0 versus pN0(i+) cases were as follows: local recurrence 1.7% versus 3.7% (P=.20), RR 0.5% versus 2.2% (P=.02), and LRR 2.1% versus 5.8% (P=.02). There were no RR events in 26 patients with pN0(i+) disease who received nodal RT and 2 RR events in 70 patients who did not receive nodal RT. On multivariable analysis, pN0(i+) was not associated with worse locoregional control or survival. On case-match analysis, LRR and overall survival were similar between pN0(i+) and matched pN0 counterparts. Conclusions: Nodal involvement with isolated tumor cells is not a significant prognostic factor for LRR or survival in this study's multivariable and case-match analyses. These data do not support the routine use of nodal RT in the setting of pN0(i+) disease. Prospective studies are needed to define optimal locoregional management for women with pN0(i+) breast cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Preliminary Toxicity Analysis of 3-Dimensional Conformal Radiation Therapy Versus Intensity Modulated Radiation Therapy on the High-Dose Arm of the Radiation Therapy Oncology Group 0126 Prostate Cancer Trial

Purpose: To give a preliminary report of clinical and treatment factors associated with toxicity in men receiving high-dose radiation therapy (RT) on a phase 3 dose-escalation trial. Methods and Materials: The trial was initiated with 3-dimensional conformal RT (3D-CRT) and amended after 1 year to allow intensity modulated RT (IMRT). Patients treated with 3D-CRT received 55.8 Gy to a planning target volume that included the prostate and seminal vesicles, then 23.4 Gy to prostate only. The IMRT patients were treated to the prostate and proximal seminal vesicles to 79.2 Gy. Common Toxicity Criteria, version 2.0, and Radiation Therapy Oncology Group/European Organization for Research and Treatment of Cancer late morbidity scores were used for acute and late effects. Results: Of 763 patients randomized to the 79.2-Gy arm of Radiation Therapy Oncology Group 0126 protocol, 748 were eligible and evaluable: 491 and 257 were treated with 3D-CRT and IMRT, respectively. For both bladder and rectum, the volumes receiving 65, 70, and 75 Gy were significantly lower with IMRT (all P<.0001). For grade (G) 2+ acute gastrointestinal/genitourinary (GI/GU) toxicity, both univariate and multivariate analyses showed a statistically significant decrease in G2+ acute collective GI/GU toxicity for IMRT. There were no significant differences with 3D-CRT or IMRT for acute or late G2+ or 3+ GU toxicities. Univariate analysis showed a statistically significant decrease in late G2+ GI toxicity for IMRT (P=.039). On multivariate analysis, IMRT showed a 26% reduction in G2+ late GI toxicity (P=.099). Acute G2+ toxicity was associated with late G3+ toxicity (P=.005). With dose–volume histogram data in the multivariate analysis, RT modality was not significant, whereas white race (P=.001) and rectal V70 ≥15% were associated with G2+ rectal toxicity (P=.034). Conclusions: Intensity modulated RT is associated with a significant reduction in acute G2+ GI/GU toxicity. There is a trend for a clinically meaningful reduction in late G2+ GI toxicity with IMRT. The occurrence of acute GI toxicity and large (>15%) volumes of rectum >70 Gy are associated with late rectal toxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Effectiveness of Chemoradiation for Head and Neck Cancer in an Older Patient Population

Purpose: The purpose of this study was to compare chemoradiation therapy (CRT) with radiation therapy (RT) only in an older patient population with head and neck squamous cell carcinoma (HNSCC). Methods and Materials: Using the Surveillance, Epidemiology, and End Results (SEER)-Medicare linked database (1992-2007), we identified a retrospective cohort of nonmetastatic HNSCC patients and divided them into treatment groups. Comparisons were made between CRT and RT cohorts. Propensity scores for CRT were estimated from covariates associated with receipt of treatment using multivariable logistic regression. Standardized mortality ratio weights (SMRW) were created from the propensity scores and used to balance groups on measured confounders. Multivariable and SMR-weighted Cox proportional hazard models were used to estimate the hazard ratio (HR) of death for receipt of CRT versus RT among the whole group and for separate patient and tumor categories. Results: The final cohort of 10,599 patients was 68% male and 89% white. Median age was 74 years. Seventy-four percent were treated with RT, 26% were treated with CRT. Median follow-up points for CRT and RT survivors were 4.6 and 6.3 years, respectively. On multivariable analysis, HR for death with CRT was 1.13 (95% confidence interval [CI]: 1.07-1.20; P<.01). Using the SMRW model, the HR for death with CRT was 1.08 (95% CI: 1.02-1.15; P=.01). Conclusions: Although the addition of chemotherapy to radiation has proven efficacious in many randomized controlled trials, it may be less effective in an older patient population treated outside of a controlled trial setting.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Fractionation of Palliative Radiation Therapy for Bone Metastases in Ontario: Do Practice Guidelines Guide Practice?

Purpose: To evaluate the effect of a provincial practice guideline on the fractionation of palliative radiation therapy for bone metastases (PRT.B) in Ontario. Methods and Materials: The present retrospective study used electronic treatment records linked to Ontario's population-based cancer registry. Hierarchical multivariable regression analysis was used to evaluate temporal trends in the use of single fractions (SFs), controlling for patient-related factors associated with the use of SFs. Results: From 1984 to 2012, 43.9% of 161,835 courses of PRT.B were administered as SFs. The percentage of SF courses was greater for older patients (age <50 years, 39.8% vs age >80 years, 52.5%), those with a shorter life expectancy (survival >12 months, 36.9% vs < 1 month, 53.6%), and those who lived farther from a radiation therapy center (<10 km, 42.1% vs > 50 km, 47.3%). The percentage of SFs to spinal fields was lower than that to other skeletal sites (31.5% vs 57.1%). The percentage of SFs varied among the cancer centers (range, 26.0%-67.8%). These differences were all highly significant in the multivariable analysis (P<.0001). In 2004, Cancer Care Ontario released a practice guideline endorsing the use of SFs for uncomplicated bone metastases. The rate of use of SFs increased from 42.3% in the pre-guideline period (1999-2003) to 52.6% in the immediate post-guideline period (2004-2007). However, it subsequently decreased again to 44.0% (2009-2012). These temporal trends were significant after controlling for patient-related factors in the multivariable analysis (P<.0001). Large intercenter variations in the use of SFs persisted after publication of the guideline. Conclusions: The publication of an Ontario practice guideline endorsing the use of SF PRT.B was associated with only a transient increase in the use of SFs in Ontario and did little to reduce intercenter variations in fractionation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Radiation Therapy Noncompliance and Clinical Outcomes in an Urban Academic Cancer Center

Purpose: To examine associations between radiation therapy (RT) noncompliance and clinical outcomes. Methods and Materials: We reviewed all patients who completed courses of external beam RT with curative intent in our department from the years 2007 to 2012 for cancers of the head and neck, breast, lung, cervix, uterus, or rectum. Patients who missed 2 or more scheduled RT appointments (excluding planned treatment breaks) were deemed noncompliant. Univariate, multivariable, and propensity-matched analyses were performed to examine associations between RT noncompliance and clinical outcomes. Results: Of 1227 patients, 266 (21.7%) were noncompliant. With median follow-up of 50.9 months, 108 recurrences (8.8%) and 228 deaths (18.6%) occurred. In univariate analyses, RT noncompliance was associated with increased recurrence risk (5-year cumulative incidence 16% vs 7%, P<.001), inferior recurrence-free survival (5-year actuarial rate 63% vs 79%, P<.001), and inferior overall survival (5-year actuarial rate 72% vs 83%, P<.001). In multivariable analyses that were adjusted for disease site and stage, comorbidity score, gender, ethnicity, race, and socioeconomic status (SES), RT noncompliance was associated with inferior recurrence, recurrence-free survival, and overall survival rates. Propensity score–matched models yielded results nearly identical to those seen in univariate analyses. Low SES was associated with RT noncompliance and was associated with inferior clinical outcomes in univariate analyses, but SES was not associated with inferior outcomes in multivariable models. Conclusion: For cancer patients being treated with curative intent, RT noncompliance is associated with inferior clinical outcomes. The magnitudes of these effects demonstrate that RT noncompliance can serve as a behavioral biomarker to identify high-risk patients who require additional interventions. Treatment compliance may mediate the associations that have been observed linking SES and clinical outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Adjuvant Radiation Therapy for Margin-Positive Vulvar Squamous Cell Carcinoma: Defining the Ideal Dose-Response Using the National Cancer Data Base

Purpose: Positive surgical margins after radical vulvectomy for vulvar cancer portend a high risk for local relapse, which may be challenging to salvage. We assessed the impact of adjuvant radiation therapy (aRT) on overall survival (OS) and the dose-response relationship using the National Cancer Data Base. Methods and Materials: Patients with vulvar squamous cell carcinoma who underwent initial extirpative surgery with positive margins from 1998 to 2012 were included. Factors associated with aRT and specific dose levels were analyzed using logistic regression. Log-rank and multivariable Cox proportional hazards modeling were used for OS analysis. Results: We identified 3075 patients with a median age of 66 years (range, 22-90 years); the median follow-up time was 36.4 months (interquartile range [IQR] 15.4-71.0 months). Stage IA/B disease represented 41.2% of the cohort. Sixty-three percent underwent lymph node assessment, with a 45% positivity rate. In total, 1035 patients (35.3%) received aRT, with a median dose of 54.0 Gy (IQR 48.6-60.0 Gy). The 3-year OS improved from 58.5% to 67.4% with aRT (P<.001). On multivariable analysis, age, Charlson-Deyo score ≥1, stage ≥II, tumors ≥4 cm, no aRT, and adverse nodal characteristics led to inferior survival. Dose of aRT was positively associated with OS as a continuous variable on univariate analysis (P<.001). The unadjusted 3-year OS for dose subsets 30.0 to 45.0 Gy, 45.1 to 53.9 Gy, 54.0 to 59.9 Gy, and ≥60 Gy was 54.3%, 55.7%, 70.1%, and 65.3%, respectively (P<.001). Multivariable analysis using a 4-month conditional landmark revealed that the greatest mortality reduction occurred in cumulative doses ≥54 Gy: 45.1 to 53.9 Gy (hazard ratio [HR] 0.94, P=.373), 54.0 to 59.9 Gy (HR 0.75, P=.024), ≥60 Gy (HR 0.71, P=.015). No survival benefit was seen with ≥60 Gy compared with 54.0 to 59.9 Gy (HR 0.95, P=.779). Conclusions: Patients with vulvar squamous cell carcinoma and positive surgical margins derive an OS benefit from aRT with a seemingly optimal dose in the range of 54.0 to 59.9 Gy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Stereotactic Body Radiation Therapy and the Influence of Chemotherapy on Overall Survival for Large (≥5 Centimeter) Non-Small Cell Lung Cancer

Purpose: Stereotactic body radiation therapy (SBRT) for ≥5 cm lesions is poorly defined, largely owing to the low sample sizes in existing studies. The present analysis examined the SBRT outcomes and assessed the effect of chemotherapy in this population. Methods and Materials: The National Cancer Data Base was queried for primary non-small cell lung cancer ≥5 cm treated with SBRT (≤10 fractions). Patient, tumor, and treatment parameters were extracted. The primary outcome was overall survival (OS). Statistical methods involved Kaplan-Meier analysis and multivariable Cox proportional hazards modeling. Results: From 2004 to 2012, data from 201 patients were analyzed. The median follow-up was 41.1 months. The median tumor size was 5.5 cm (interquartile range 5.0-6.0), with cT2a, cT2b, and cT3 disease in 24.9%, 53.2%, and 21.9%, respectively. The median total SBRT dose and fractionation was 50 Gy in 4 fractions, and 92.5% of the patients underwent SBRT with ≤5 fractions. The median OS was 25.1 months. Of the 201 patients, 15% received chemotherapy. The receipt of chemotherapy was associated with longer OS (median 30.6 vs 23.4 months; P=.027). On multivariable analysis, worse OS was seen with increasing age (hazard ratio [HR] 1.03; P=.012), poorly differentiated tumors (HR 2.06; P=.049), and T3 classification (HR 2.13; P=.005). On multivariable analysis, chemotherapy remained independently associated with improved OS (HR 0.57; P=.039). Conclusions: SBRT has utility in the setting of tumors ≥5 cm, with chemotherapy associated with improved OS in this subset. These hypothesis-generating data now raise the necessity of performing prospective analyses to determine whether chemotherapy confers outcome benefits after SBRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Understanding Competing Risks

Highlights: • Clinicians and researchers need to be cognizant of competing risks (ie, the occurrence of an event that precludes future development of the outcome of interest). • Censoring is appropriate when a patient being censored at a particular time (eg, owing to loss to follow-up) is deemed to have approximately the same future risk of developing the outcome of interest as patients still being followed, and when the reason for censoring is not related to events of interest. • Censoring is not appropriate in many scenarios; for example, a dead patient can never develop a future cancer recurrence or late effect. • Kaplan-Meier (KM) analysis is commonly used for a time-to-event analysis when the primary outcome (eg, overall survival) has no meaningful competing risk and censoring is appropriate. Cox proportional hazards modeling can be used for multivariable analysis to examine the association between covariates and the primary outcome. • In the presence of competing events, incorrect use of the KM analysis will likely lead to an overestimation of the risk of the outcome of interest and biased estimates from both univariable and multivariable analyses. Therefore, KM usually should not be used to determine local recurrence risk because humanity is always at risk for dying, which is a competing risk. • A competing risk analysis is therefore most appropriate when there are competing events to the primary outcome of interest. The most commonly used multivariable model in these situations is the Fine-Gray model. • Many fundamental questions in radiation oncology—local and regional recurrence, treatment-related toxicity—are best addressed using competing risk methodologies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Practical Guide to Chemometric Analysis of Optical Spectroscopic Data

The methodology and mathematical treatment of several classic multivariate methods for the analysis of spectroscopic data is demonstrated in a straightforward way that can be used as a basis for teaching an undergraduate introductory course on chemometric analysis. The multivariate techniques of classical least squares (CLS), principal component regression (PCR), and partial least squares (PLS), as well as the univariate Beer’s law method have been described and compared, building students’ understanding by starting with the univariate method and progressing step by step into the multivariate methods. Equations for the production of regression vectors from training set spectral data is described and their use demonstrated for the prediction of constituent concentrations on a separate validation set of spectra. Extreme care is taken to ensure consistency in variable formatting of data matrices. This provides a key foundation to understanding how spectral data are manipulated using these different mathematical approaches for building quantitative regression models. Each method is applied to a real-world data set, and the results are discussed to show students the types of information that can be gleaned from each method. A training set comprised of 20 infrared absorbance spectra containing 3 constituents (benzene, polystyrene, and gasoline) of known composition are used to demonstrate the matrix operations for each regression method. A separate set of 12 real-world napalm samples (containing benzene, polystyrene and gasoline) are used as a validation set to demonstrate the ability to utilize the regression models on an unknown dataset. A toolbox (PNNL Chemometric Toolbox) written in MATLAB language is supplied in the Supplemental Information file and can be used as a companion for understanding the development and deployment of the chemometric algorithms described in this paper. The datasets of the infrared spectra are also supplied, allowing users to build and inspect the chemometric models on their own. Finally, the Toolbox includes scripts to assist users in loading their own datasets into MATLAB and performing CLS, PCR, and PLS on their data.

Upper-Division Undergraduate, Analytical Chemistry↗