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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A NOVEL, FULLY AUTOMATED PIPELINE FOR PERIOD ESTIMATION IN THE EROS 2 DATA SET

We present a new method to discriminate periodic from nonperiodic irregularly sampled light curves. We introduce a periodic kernel and maximize a similarity measure derived from information theory to estimate the periods and a discriminator factor. We tested the method on a data set containing 100,000 synthetic periodic and nonperiodic light curves with various periods, amplitudes, and shapes generated using a multivariate generative model. We correctly identified periodic and nonperiodic light curves with a completeness of ∼90% and a precision of ∼95%, for light curves with a signal-to-noise ratio (S/N) larger than 0.5. We characterize the efficiency and reliability of the model using these synthetic light curves and apply the method on the EROS-2 data set. A crucial consideration is the speed at which the method can be executed. Using a hierarchical search and some simplification on the parameter search, we were able to analyze 32.8 million light curves in ∼18 hr on a cluster of GPGPUs. Using the sensitivity analysis on the synthetic data set, we infer that 0.42% of the sources in the LMC and 0.61% of the sources in the SMC show periodic behavior. The training set, catalogs, and source code are all available at http://timemachine.iic.harvard.edu.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Imaging nanoscale lattice variations by machine learning of x-ray diffraction microscopy data

In this paper, we present a novel methodology based on machine learning to extract lattice variations in crystalline materials, at the nanoscale, from an x-ray Bragg diffraction-based imaging technique. By employing a full-field microscopy setup, we capture real space images of materials, with imaging contrast determined solely by the x-ray diffracted signal. The data sets that emanate from this imaging technique are a hybrid of real space information (image spatial support) and reciprocal lattice space information (image contrast), and are intrinsically multidimensional (5D). By a judicious application of established unsupervised machine learning techniques and multivariate analysis to this multidimensional data cube, we show how to extract features that can be ascribed physical interpretations in terms of common structural distortions, such as lattice tilts and dislocation arrays. Finally, we demonstrate this 'big data' approach to x-ray diffraction microscopy by identifying structural defects present in an epitaxial ferroelectric thin-film of lead zirconate titanate.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A new MR-SAD algorithm for the automatic building of protein models from low-resolution X-ray data and a poor starting model

Determining macromolecular structures from X-ray data with resolution worse than 3Å remains a challenge. Even if a related starting model is available, its incompleteness or its bias together with a low observation-to-parameter ratio can render the process unsuccessful or very time-consuming. Yet, many biologically important macromolecules, especially large macromolecular assemblies, membrane proteins and receptors, tend to provide crystals that diffract to low resolution. A new algorithm to tackle this problem is presented that uses a multivariate function to simultaneously exploit information from both an initial partial model and low-resolution single-wavelength anomalous diffraction data. The new approach has been used for six challenging structure determinations, including the crystal structures of membrane proteins and macromolecular complexes that have evaded experts using other methods, and large structures from a 3.0 Å resolution F1-ATPase data set an a 4.5Å resolution SecYEG–SecA complex data set. All of the models were automatically built by the method to Rfree values of between 28.9 and 39.9% and were free from the initial model bias.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Towards the extraction of the crystal cell parameters from pair distribution function profiles

The approach based on atomic pair distribution function (PDF) has revolutionized structural investigations by X-ray/electron diffraction of nano or quasi-amorphous materials, opening up the possibility of exploring short-range order. However, the ab initio crystal structural solution by the PDF is far from being achieved due to the difficulty in determining the crystallographic properties of the unit cell. A method for estimating the crystal cell parameters directly from a PDF profile is presented, which is composed of two steps: first, the type of crystal cell is inferred using machine-learning approaches applied to the PDF profile; second, the crystal cell parameters are extracted by means of multivariate analysis combined with vector superposition techniques. The procedure has been validated on a large number of PDF profiles calculated from known crystal structures and on a small number of measured PDF profiles. The lattice determination step has been benchmarked by a comprehensive exploration of different classifiers and different input data. The highest performance is obtained using the k -nearest neighbours classifier applied to whole PDF profiles. Descriptors calculated from the PDF profiles by recurrence quantitative analysis produce results that can be interpreted in terms of PDF properties, and the significance of each descriptor in determining the prediction is evaluated. The cell parameter extraction step depends on the cell metric rather than its type. Monometric, dimetric and trimetric cells have top-1 estimates that are correct 40, 20 and 5% of the time, respectively. Promising results were obtained when analysing real nanocrystals, where unit cells close to the true ones are found within the top-1 ranked solution in the case of monometric cells and within the top-6 ranked solutions in the case of dimetric cells, even in the presence of a crystalline impurity with a weight fraction up to 40%.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Forced Oscillation Source Location in Power Systems Using MVMD-assisted DEF in TF Plane

Fast and accurate forced oscillation source location (FOSL) is critical for mitigating forced oscillations in power systems. Here, to overcome the shortcomings of traditional dissipating energy flow (DEF) in decomposing the forced oscillation components from the measurement responses, a multivariate variational mode decomposition (MVMD)-assisted FOSL is proposed in this paper to locate the forced oscillation sources from measurements in time-frequency (TF) plane. The multi-channel measurement matrix of each generator is formed, then the multi-channel modes are decomposed from the formed multi-channel measurement matrix using the MVMD. Further, the intrinsic mode functions (IMFs) associated with the forced oscillation mode are distinguished from the decomposed modes according to the relative energy weights. With the distinguished IMFs, the MVMD-assisted DEF calculation method is developed to locate the forced oscillation source. Finally, the performance of the proposed method is evaluated using the simulation data of the WECC 179-bus test system and IEEE-NASPI Oscillation Source Location Contest as well as the actual PMU data of ISO New England. The results validated the effectiveness and accuracy of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Time Series Classification for Locating Forced Oscillation Sources

Here, this article presents a machine learning based time-series classification method for using synchrophasor measurements to locate the source of forced oscillation (FO) for fast disturbance removal. First, multivariate time series (MTS) matrices are constructed by the most informative measurements selected by sequential feature selection from each power plant. Then, the Mahalanobis matrix is trained such that the Mahalanobis distance between the MTSs from the same class (i.e., with the same FO source location) are minimized and from different classes (i.e., with different FO source locations) are maximized. This allows MTSs to be classified by classifiers with class membership corresponding to the location of each FO source. To meet the runtime requirements of online matching, class templates are constructed to reduce data size and improve matching efficiency. To account for uncertainty in identifying the exact beginning of an FO event, dynamic time warping is used to align the out-of-sync MTSs. IEEE 39bus and WECC 179bus systems are used for algorithm development and validation. Simulation results demonstrate that the algorithm meets online operation runtime requirement with high accuracy using misaligned data sets.

42 ENGINEERING↗
DOE OSTI2013

A Model-Based Framework for Air Quality Indices and Population Risk Evaluation, with an Application to the Analysis of Scottish Air Quality Data

Summary The paper is devoted to the development of a statistical framework for air quality assessment at the country level and for the evaluation of the ambient population exposure and risk with respect to airborne pollutants. The framework is based on a multivariate space–time model and on aggregated indices defined at different levels of aggregation in space and time. The indices are evaluated, uncertainty included, by considering both the model outputs and the information on the population spatial distribution. The framework is applied to the analysis of air quality data for Scotland for 2009 referring to European and Scottish air quality legislation.

Unbalanced networks↗
DOE OSTI2014

Predicting error in detecting mammographic masses among radiology trainees using statistical models based on BI-RADS features

Purpose: The purpose of this study is to explore Breast Imaging-Reporting and Data System (BI-RADS) features as predictors of individual errors made by trainees when detecting masses in mammograms. Methods: Ten radiology trainees and three expert breast imagers reviewed 100 mammograms comprised of bilateral medial lateral oblique and craniocaudal views on a research workstation. The cases consisted of normal and biopsy proven benign and malignant masses. For cases with actionable abnormalities, the experts recorded breast (density and axillary lymph nodes) and mass (shape, margin, and density) features according to the BI-RADS lexicon, as well as the abnormality location (depth and clock face). For each trainee, a user-specific multivariate model was constructed to predict the trainee's likelihood of error based on BI-RADS features. The performance of the models was assessed using area under the receive operating characteristic curves (AUC). Results: Despite the variability in errors between different trainees, the individual models were able to predict the likelihood of error for the trainees with a mean AUC of 0.611 (range: 0.502–0.739, 95% Confidence Interval: 0.543–0.680,p < 0.002). Conclusions: Patterns in detection errors for mammographic masses made by radiology trainees can be modeled using BI-RADS features. These findings may have potential implications for the development of future educational materials that are personalized to individual trainees.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

WE-E-17A-05: Complementary Prognostic Value of CT and 18F-FDG PET Non-Small Cell Lung Cancer Tumor Heterogeneity Features Quantified Through Texture Analysis

Purpose: Previous studies have shown that CT or 18F-FDG PET intratumor heterogeneity features computed using texture analysis may have prognostic value in Non-Small Cell Lung Cancer (NSCLC), but have been mostly investigated separately. The purpose of this study was to evaluate the potential added value with respect to prognosis regarding the combination of non-enhanced CT and 18F-FDG PET heterogeneity textural features on primary NSCLC tumors. Methods: One hundred patients with non-metastatic NSCLC (stage I–III), treated with surgery and/or (chemo)radiotherapy, that underwent staging 18F-FDG PET/CT images, were retrospectively included. Morphological tumor volumes were semi-automatically delineated on non-enhanced CT using 3D SlicerTM. Metabolically active tumor volumes (MATV) were automatically delineated on PET using the Fuzzy Locally Adaptive Bayesian (FLAB) method. Intratumoral tissue density and FDG uptake heterogeneities were quantified using texture parameters calculated from co-occurrence, difference, and run-length matrices. In addition to these textural features, first order histogram-derived metrics were computed on the whole morphological CT tumor volume, as well as on sub-volumes corresponding to fine, medium or coarse textures determined through various levels of LoG-filtering. Association with survival regarding all extracted features was assessed using Cox regression for both univariate and multivariate analysis. Results: Several PET and CT heterogeneity features were prognostic factors of overall survival in the univariate analysis. CT histogram-derived kurtosis and uniformity, as well as Low Grey-level High Run Emphasis (LGHRE), and PET local entropy were independent prognostic factors. Combined with stage and MATV, they led to a powerful prognostic model (p<0.0001), with median survival of 49 vs. 12.6 months and a hazard ratio of 3.5. Conclusion: Intratumoral heterogeneity quantified through textural features extracted from both CT and FDG PET images have complementary and independent prognostic value in NSCLC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-P-18: Intensity-Modulated Radiation Therapy for Cervical Esophageal Squamous Cell Carcinoma

Purpose: To retrospectively analyze the outcomes and prognostic factors of cervical esophageal squamous cell carcinoma (SCC) treated with intensity modulated radiation therapy (IMRT). Methods: Thirty-seven patients with cervical esophageal SCC treated with IMRT were analyzed retrospectively. They received 54–66 Gy in 27–32 fractions. Nineteen patients received concurrent (n=12) or sequential (n=7) platinum-based two drugs chemoradiotherapy. Overall survival (OS), local control rates (LCR) and prognostic factors were evaluated. Acute toxicities and patterns of first failures were observed. Results: The median follow-up was 46 months for alive patients. The l-, 3-, 4- and 5-year OS of the all patients were 83.8%, 59.1%, 47.5% and 32.6% respectively. The median survival time was 46 months. The l-, 3-,4- and 5-year LCR were 82.9%, 63.0%, 54.5% and 54.5%, respectively. Univariate and Multivariate analysis all showed that size of GTV was an independent prognostic factor (p=0.033, p=0.039). There were no patients with Grade 3 acute radiation esophagitis and Grade 2–4 acute pneumonitis. The local failure accounted for 70.0% of all treatment-related failures. Conclusion: IMRT is safe and effective in the treatment of cervical esophageal squamous cell carcinoma. Size of GTV is an independent prognostic factor. Local failure still remains the main reason of treatment failures. The authors declare no conflicts of interest in preparing this article.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-253: The Radiomics Toolbox in the Computational Environment for Radiological Research (CERR)

Purpose: To present an open source and free platform to facilitate radiomics research — The “Radiomics toolbox” in CERR. Method: There is scarcity of open source tools that support end-to-end modeling of image features to predict patient outcomes. The “Radiomics toolbox” strives to fill the need for such a software platform. The platform supports (1) import of various kinds of image modalities like CT, PET, MR, SPECT, US. (2) Contouring tools to delineate structures of interest. (3) Extraction and storage of image based features like 1st order statistics, gray-scale co-occurrence and zonesize matrix based texture features and shape features and (4) Statistical Analysis. Statistical analysis of the extracted features is supported with basic functionality that includes univariate correlations, Kaplan-Meir curves and advanced functionality that includes feature reduction and multivariate modeling. The graphical user interface and the data management are performed with Matlab for the ease of development and readability of code and features for wide audience. Open-source software developed with other programming languages is integrated to enhance various components of this toolbox. For example: Java-based DCM4CHE for import of DICOM, R for statistical analysis. Results: The Radiomics toolbox will be distributed as an open source, GNU copyrighted software. The toolbox was prototyped for modeling Oropharyngeal PET dataset at MSKCC. The analysis will be presented in a separate paper. Conclusion: The Radiomics Toolbox provides an extensible platform for extracting and modeling image features. To emphasize new uses of CERR for radiomics and image-based research, we have changed the name from the “Computational Environment for Radiotherapy Research” to the “Computational Environment for Radiological Research”.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-119: Analysis the Efficacy of Different Radiotherapy Methods and Failure Mode in No-Metastasis Esophageal Squamous Cell Carcinoma

Purpose: To evaluate the curative effect of radio (chemo) therapy and mode of treatment failure in no-metastasis and lesion length ≤ 5.0cm esophageal squamous cell carcinoma (ESCC). Methods: There were 158 eligible patients were retrospectively analyzed, to analysis the curative effect of radio (chemo) therapy, prognosis factors, toxicity and prognostic index model. Results: To all patients the 1, 3, 5 overall survival rate were 83.54%, 52.53%, 32.58%, the local recurrence rate were 15.08%, 33.60% and 38.14%; distant metastasis rate were 10.64%, 25.21% and 36.06%; tumor specific survival rate were 76.64%, 54.07% and 44.51%. Multivariate analysis showed that patients with ECOG grade (χ2=13.945, P=0.000), short-term effect (χ2=19.360, P=0.000) and different radiotherapy methods (χ2=9.866, P=0.002) as the independent prognostic factors. Prognostic index model showed that the survival rate was significantly higher in the lower value of PI group than in the larger value of PI group (χ2=49.19, P=0.0000). In our whole group, there were simple locoregional recurrence (LR) 40 cases (25.3%), simple Distant metastasis (DM) 31 cases (19.6%), LR and DM in 14 cases (8.9%) after treatment. The chi-square test showed that there were no significant difference in the incidence of Elective Nodal Irradiation (ENI )and Involved Field Irradiation (IFI) patients with LR and DM ( χ2=2.363, 2.950, P=0.124, 0.085). Conclusion: Radio (chemo) therapy has a good curative effect in no-metastasis and lesion length ≤ 5.0cm ESCC patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-CD-BRB-04: Automated Radiomic Features Complement the Prognostic Value of VASARI in the TCGA-GBM Dataset

Purpose: To compare the complementary prognostic value of automated Radiomic features to that of radiologist-annotated VASARI features in TCGA-GBM MRI dataset. Methods: For 96 GBM patients, pre-operative MRI images were obtained from The Cancer Imaging Archive. The abnormal tumor bulks were manually defined on post-contrast T1w images. The contrast-enhancing and necrotic regions were segmented using FAST. From these sub-volumes and the total abnormal tumor bulk, a set of Radiomic features quantifying phenotypic differences based on the tumor intensity, shape and texture, were extracted from the post-contrast T1w images. Minimum-redundancy-maximum-relevance (MRMR) was used to identify the most informative Radiomic, VASARI and combined Radiomic-VASARI features in 70% of the dataset (training-set). Multivariate Cox-proportional hazards models were evaluated in 30% of the dataset (validation-set) using the C-index for OS. A bootstrap procedure was used to assess significance while comparing the C-Indices of the different models. Results: Overall, the Radiomic features showed a moderate correlation with the radiologist-annotated VASARI features (r = −0.37 – 0.49); however that correlation was stronger for the Tumor Diameter and Proportion of Necrosis VASARI features (r = −0.71 – 0.69). After MRMR feature selection, the best-performing Radiomic, VASARI, and Radiomic-VASARI Cox-PH models showed a validation C-index of 0.56 (p = NS), 0.58 (p = NS) and 0.65 (p = 0.01), respectively. The combined Radiomic-VASARI model C-index was significantly higher than that obtained from either the Radiomic or VASARI model alone (p = <0.001). Conclusion: Quantitative volumetric and textural Radiomic features complement the qualitative and semi-quantitative annotated VASARI feature set. The prognostic value of informative qualitative VASARI features such as Eloquent Brain and Multifocality is increased with the addition of quantitative volumetric and textural features from the contrast-enhancing and necrotic tumor regions. These results should be further evaluated in larger validation cohorts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-CD-BRB-05: Radiation Damage Signature of White Matter Fiber Bundles Using Diffusion Tensor Imaging (DTI)

Purpose: To develop an automated and scalable approach and identify temporal, spatial and dosimetric patterns of radiation damage of white matter (WM) fibers following partial brain irradiation. Methods: An automated and scalable approach was developed to extract DTI features of 22 major WM fibers from 33 patients with low-grade/benign tumors treated by radiation therapy (RT). DTI scans of the patients were performed pre-RT, 3- and 6-week during RT, and 1, 6 and 18 months after RT. The automated tractography analysis was applied to 198 datasets as: (1) intra-subject registration of longitudinal DTI, (2) spatial normalization of individual-patient DTI to the Johns Hopkins WM Atlas, (3) automatic fiber tracking regulated by the WM Atlas, and (4) segmentation of WM into 22 major tract profiles. Longitudinal percentage changes in fractional anisotropy (FA), and mean, axial and radial diffusivity (MD/AD/RD) of each tract from pre-RT were quantified and correlated to 95%, 90% and 80% percentiles of doses and mean doses received by the tract. Heatmaps were used to identify clusters of significant correlation and reveal temporal, spatial and dosimetric signatures of WM damage. A multivariate linear regression was further carried out to determine influence of clinical factors. Results: Of 22 tracts, AD/MD changes in 12 tracts had significant correlation with doses, especially at 6 and 18 months post-RT, indicating progressive radiation damage after RT. Most interestingly, the DTI-index changes in the elongated tracts were associated with received maximum doses, suggesting a serial-structure behavior; while short association fibers were affected by mean doses, indicating a parallel-structure response. Conclusion: Using an automated DTI-tractography analysis of whole brain WM fibers, we reveal complex radiation damage patterns of WM fibers. Damage in WM fibers that play an important role in the neural network could be associated with late neurocognitive function declines after brain irradiation. NIH NS064973.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Computer-aided diagnosis for classifying benign versus malignant thyroid nodules based on ultrasound images: A comparison with radiologist-based assessments

Purpose: To develop a semiautomated computer-aided diagnosis (CAD) system for thyroid cancer using two-dimensional ultrasound images that can be used to yield a second opinion in the clinic to differentiate malignant and benign lesions. Methods: A total of 118 ultrasound images that included axial and longitudinal images from patients with biopsy-confirmed malignant (n = 30) and benign (n = 29) nodules were collected. Thyroid CAD software was developed to extract quantitative features from these images based on thyroid nodule segmentation in which adaptive diffusion flow for active contours was used. Various features, including histogram, intensity differences, elliptical fit, gray-level co-occurrence matrixes, and gray-level run-length matrixes, were evaluated for each region imaged. Based on these imaging features, a support vector machine (SVM) classifier was used to differentiate benign and malignant nodules. Leave-one-out cross-validation with sequential forward feature selection was performed to evaluate the overall accuracy of this method. Additionally, analyses with contingency tables and receiver operating characteristic (ROC) curves were performed to compare the performance of CAD with visual inspection by expert radiologists based on established gold standards. Results: Most univariate features for this proposed CAD system attained accuracies that ranged from 78.0% to 83.1%. When optimal SVM parameters that were established using a grid search method with features that radiologists use for visual inspection were employed, the authors could attain rates of accuracy that ranged from 72.9% to 84.7%. Using leave-one-out cross-validation results in a multivariate analysis of various features, the highest accuracy achieved using the proposed CAD system was 98.3%, whereas visual inspection by radiologists reached 94.9% accuracy. To obtain the highest accuracies, “axial ratio” and “max probability” in axial images were most frequently included in the optimal feature sets for the authors’ proposed CAD system, while “shape” and “calcification” in longitudinal images were most frequently included in the optimal feature sets for visual inspection by radiologists. The computed areas under curves in the ROC analysis were 0.986 and 0.979 for the proposed CAD system and visual inspection by radiologists, respectively; no significant difference was detected between these groups. Conclusions: The use of thyroid CAD to differentiate malignant from benign lesions shows accuracy similar to that obtained via visual inspection by radiologists. Thyroid CAD might be considered a viable way to generate a second opinion for radiologists in clinical practice.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Statistical-learning strategies generate only modestly performing predictive models for urinary symptoms following external beam radiotherapy of the prostate: A comparison of conventional and machine-learning methods

Purpose: Given the paucity of available data concerning radiotherapy-induced urinary toxicity, it is important to ensure derivation of the most robust models with superior predictive performance. This work explores multiple statistical-learning strategies for prediction of urinary symptoms following external beam radiotherapy of the prostate. Methods: The performance of logistic regression, elastic-net, support-vector machine, random forest, neural network, and multivariate adaptive regression splines (MARS) to predict urinary symptoms was analyzed using data from 754 participants accrued by TROG03.04-RADAR. Predictive features included dose-surface data, comorbidities, and medication-intake. Four symptoms were analyzed: dysuria, haematuria, incontinence, and frequency, each with three definitions (grade ≥ 1, grade ≥ 2 and longitudinal) with event rate between 2.3% and 76.1%. Repeated cross-validations producing matched models were implemented. A synthetic minority oversampling technique was utilized in endpoints with rare events. Parameter optimization was performed on the training data. Area under the receiver operating characteristic curve (AUROC) was used to compare performance using sample size to detect differences of ≥0.05 at the 95% confidence level. Results: Logistic regression, elastic-net, random forest, MARS, and support-vector machine were the highest-performing statistical-learning strategies in 3, 3, 3, 2, and 1 endpoints, respectively. Logistic regression, MARS, elastic-net, random forest, neural network, and support-vector machine were the best, or were not significantly worse than the best, in 7, 7, 5, 5, 3, and 1 endpoints. The best-performing statistical model was for dysuria grade ≥ 1 with AUROC ± standard deviation of 0.649 ± 0.074 using MARS. For longitudinal frequency and dysuria grade ≥ 1, all strategies produced AUROC>0.6 while all haematuria endpoints and longitudinal incontinence models produced AUROC<0.6. Conclusions: Logistic regression and MARS were most likely to be the best-performing strategy for the prediction of urinary symptoms with elastic-net and random forest producing competitive results. The predictive power of the models was modest and endpoint-dependent. New features, including spatial dose maps, may be necessary to achieve better models.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-103: Analysis of Hepatitis B Virus Reactivation After Conformal Radiotherapy in Patients with Hepatocellular Carcinoma Using the Lyman NTCP Model

Purpose: The aim of this research was to investigate the feasibility of Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) model in analyzing hepatitis B virus (HBV) reactivation in patients receiving conformal radiotherapy for patients with hepatocellular carcinoma (HCC). Methods: Between June 2009 and June 2012, 108 HBV-related HCC patients (90 were specifically selected and 18 patients were excluded) treated with conformal RT at three centers were enrolled in this retrospective study. They were all diagnosed as HCC by pathology or cytology. All 90 patients were followed up to September 2013 with a median follow-up time of 25.2 months. The parameters (TD50 (1), n, and m) of the modified LKB NTCP model were derived using maximum likelihood estimation. Bootstrap and leave-one-out were employed to test the generalizability of the results for use in a general population. Results: The incidences of complications in the study population were as follows: radiation-induced liver diseases (RILD) were 17.6%, HBV reactivation was 24.8%, and HBV reactivation-induced hepatitis was 22.7%, respectively. In multivariate analysis, the NTCP (p<0.001), and V20 were associated with HBV reactivation. TD50 (1), m and n were 42.9Gy (95% CI) (38.2–46.8), 0.14 (0.12–0.15) and 0.30 (0.2–0.33), respectively, for HBV reactivation. Bootstrap and leave-one-out results showed that the HBV parameter fits were extremely robust. Conclusion: A modified LKB NTCP model has been established to predict HBV reactivation for patients with HCC receiving conformal RT. The finding derives parameters set to predict potential endpoints of HBV reactivation.

60 APPLIED LIFE SCIENCES↗