Robust, fast and accurate mapping of diffusional mean kurtosis.

Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.

Generalized interacting multiple model Kalman filtering algorithm for maneuvering target tracking under non-Gaussian noises.

The traditional interacting multiple model Kalman filtering algorithm (IMM-KF) can deal with the maneuvering target problem under Gaussian noise by soft switching among possible motion models. In practice, its performance is likely to degrade when handling non-Gaussian noise. We introduce the Gaussian mixture model (GMM) into the IMM-KF, and the GMM is utilized to model the non-Gaussian measurement noise as a mixture of multiple Gaussian probability densities with a certain probability. Then, a GIMM framework is proposed that enables accurate switching and fusion among multiple possible motion and noise models. And combined with Kalman filtering (KF), a GIMM-KF algorithm is proposed that enables accurate state estimation of maneuvering targets under non-Gaussian noise conditions. Subsequently, the provided simulations and experiments validate that the GIMM-KF algorithm outperforms existing methods in terms of accuracy, stability and robustness.

A novel class of non-Gaussian system performance assessment and controller parameter tuning methods.

Traditional variance-based control performance assessment (CPA) and controller parameter tuning (CPT) methods tend to ignore non-Gaussian external disturbances. To address this limitation, this study proposes a novel class of CPA and CPT methods for non-Gaussian single-input single-output systems, denoted as data Gaussianization (inverse) transformation methods. The idea of quantile transformation is used to transform the non-Gaussian data with the goal of maximizing mutual information into virtual Gaussian data. In addition, optimal system data for the virtual loop are mapped back to the actual non-Gaussian system using quantile inverse transformation. Furthermore, a CARMA model-based recursive extended least square algorithm and a CARMA model-based least absolute deviation iterative algorithm are used to identify virtual Gaussian and non-Gaussian system process models, respectively, while implementing the CPT. Finally, a unified framework is proposed for the CPA and CPT of a non-Gaussian control system. The simulation results demonstrate that the proposed strategy can provide a consistent benchmark judgment criterion (threshold) for different non-Gaussian noises, and the tuned controller parameters have good performance.

Direct Position Determination of Non-Gaussian Sources for Multiple Nested Arrays: Discrete Fourier Transform and Taylor Compensation Algorithm.

This paper delves into the problem of direct position determination (DPD) for non-Gaussian sources. Existing DPD algorithms are hindered by their high computational complexity from exhaustive grid searches and a disregard for the received signal characteristics by multiple nested arrays (MNAs). To address these issues, the paper proposes a novel DPD algorithm for non-Gaussian sources with MNAs: the Discrete Fourier Transform (DFT) and Taylor compensation algorithm. Initially, the fourth-order cumulant matrix of the received signal is computed, and the vectorizing method is applied. Subsequently, a computationally efficient DPD cost function is proposed by leveraging a normalized DFT matrix to reduce complexity. Finally, first-order Taylor compensation is utilized to enhance the accuracy of the localization results. The superiority of the proposed algorithm is demonstrated through numerical simulation results.

A novel Gaussian sum quaternion constrained cubature Kalman filter algorithm for GNSS/SINS integrated attitude determination and positioning system.

The quaternion cubature Kalman filter (QCKF) algorithm has emerged as a prominent nonlinear filter algorithm and has found extensive applications in the field of GNSS/SINS integrated attitude determination and positioning system (GNSS/SINS-IADPS) data processing for unmanned aerial vehicles (UAV). However, on one hand, the QCKF algorithm is predicated on the assumption that the random model of filter algorithm, which follows a white Gaussian noise distribution. The noise in actual GNSS/SINS-IADPS is not the white Gaussian noise but rather a ubiquitous non-Gaussian noise. On the other hand, the use of quaternions as state variables is bound by normalization constraints. When applied directly in nonlinear non-Gaussian system without considering normalization constraints, the QCKF algorithm may result in a mismatch phenomenon in the filtering random model, potentially resulting in a decline in estimation accuracy. To address this issue, we propose a novel Gaussian sum quaternion constrained cubature Kalman filter (GSQCCKF) algorithm. This algorithm refines the random model of the QCKF by approximating non-Gaussian noise with a Gaussian mixture model. Meanwhile, to account for quaternion normalization in attitude determination, a two-step projection method is employed to constrain the quaternion, which consequently enhances the filtering estimation accuracy. Simulation and experimental analyses demonstrate that the proposed GSQCCKF algorithm significantly improves accuracy and adaptability in GNSS/SINS-IADPS data processing under non-Gaussian noise conditions for Unmanned Aerial Vehicles (UAVs).

K-Space Approach in Optical Coherence Tomography: Rigorous Digital Transformation of Arbitrary-Shape Beams, Aberration Elimination and Super-Refocusing beyond Conventional Phase Correction Procedures.

For the most popular method of scan formation in Optical Coherence Tomography (OCT) based on plane-parallel scanning of the illuminating beam, we present a compact but rigorous K-space description in which the spectral representation is used to describe both the axial and lateral structure of the illuminating/received OCT signals. Along with the majority of descriptions of OCT-image formation, the discussed approach relies on the basic principle of OCT operation, in which ballistic backscattering of the illuminating light is assumed. This single-scattering assumption is the main limitation, whereas in other aspects, the presented approach is rather general. In particular, it is applicable to arbitrary beam shapes without the need for paraxial approximation or the assumption of Gaussian beams. The main result of this study is the use of the proposed K-space description to analytically derive a filtering function that allows one to digitally transform the initial 3D set of complex-valued OCT data into a desired (target) dataset of a rather general form. An essential feature of the proposed filtering procedures is the utilization of both phase and amplitude transformations, unlike conventionally discussed phase-only transformations. To illustrate the efficiency and generality of the proposed filtering function, the latter is applied to the mutual transformation of non-Gaussian beams and to the digital elimination of arbitrary aberrations at the illuminating/receiving aperture. As another example, in addition to the conventionally discussed digital refocusing enabling depth-independent lateral resolution the same as in the physical focus, we use the derived filtering function to perform digital "super-refocusing." The latter does not yet overcome the diffraction limit but readily enables lateral resolution several times better than in the initial physical focus.

Matched Stochastic Resonance Enhanced Underwater Passive Sonar Detection under Non-Gaussian Impulsive Background Noise.

Remote passive sonar detection with low-frequency band spectral lines has attracted much attention, while complex low-frequency non-Gaussian impulsive noisy environments would strongly affect the detection performance. This is a challenging problem in weak signal detection, especially for the high false alarm rate caused by heavy-tailed impulsive noise. In this paper, a novel matched stochastic resonance (MSR)-based weak signal detection model is established, and two MSR-based detectors named MSR-PED and MSR-PSNR are proposed based on a theoretical analysis of the MSR output response. Comprehensive detection performance analyses in both Gasussian and non-Gaussian impulsive noise conditions are presented, which revealed the superior performance of our proposed detector under non-Gasussian impulsive noise. Numerical analysis and application verification have revealed the superior detection performance with the proposed MSR-PSNR detector compared with energy-based detection methods, which can break through the high false alarm rate problem caused by heavy-tailed impulsive noise. For a typical non-Gasussian impulsive noise assumption with α=1.5, the proposed MSR-PED and MSR-PSNR can achieve approximately 16 dB and 22 dB improvements, respectively, in the detection performance compared to the classical PED method. For stronger, non-Gaussian impulsive noise conditions corresponding to α=1, the improvement in detection performance can be more significant. Our proposed MSR-PSNR methods can overcome the challenging problem of a high false alarm rate caused by heavy-tailed impulsive noise. This work can lay a solid foundation for breaking through the challenges of underwater passive sonar detection under non-Gaussian impulsive background noise, and can provide important guidance for future research work.

Study of the response of a single-dof dynamic system under stationary non-Gaussian random loads aimed at fatigue life assessment.

Fatigue assessment of components subjected to random loads is a challenging task both due to the variability in amplitude and frequency of the loads and for the computational times required to perform classical time domain fatigue analysis. The frequency domain approach to fatigue life assessment offers a solution by utilizing the power spectral density of the random load, requiring minimal computational effort. However, frequency domain methods are limited to stationary Gaussian signals, while real-world loads often exhibit non-Gaussian characteristics. Previous research proposed formulas to extend frequency domain methods to non-Gaussian cases, but they require knowledge of the parameters related to non-Gaussianity of the component's stress (skewness and kurtosis), which would require a time domain analysis of the stress history on the component and a strong reduction of the computational advantages. This paper aims to address this gap by conducting an extensive campaign of numerical simulations to evaluate the influence of various parameters on the degree of non-Gaussianity of the response of a system. A single-dof mass-spring-damper system was subjected to non-Gaussian random loads of different natures, and the response is analyzed to determine the values of skewness and kurtosis. The study investigated the influence on non-normality indexes of the system's output of several input parameters, which include both the characteristics of the input load and the properties of the dynamic system. The findings contribute to a better understanding of non-Gaussianity in dynamic systems and pave the way for conducting efficient fatigue analyses in the frequency domain. Future work will extend the study to non-stationary random loads, further advancing the understanding of non-Gaussianity and non-stationarity in dynamic systems.

Relating pre-treatment non-Gaussian intravoxel incoherent motion diffusion-weighted imaging to human papillomavirus status and response in oropharyngeal carcinoma.

Background and purpose: Diffusion-weighted imaging (DWI) is a promising technique for response assessment in head-and-neck cancer. Recently, we optimized Non-Gaussian Intravoxel Incoherent Motion Imaging (NG-IVIM), an extension of the conventional apparent diffusion coefficient (ADC) model, for the head and neck. In the current study, we describe the first application in a group of patients with human papillomavirus (HPV)-positive and HPV-negative oropharyngeal squamous cell carcinoma. The aim of this study was to relate ADC and NG-IVIM DWI parameters to HPV status and clinical treatment response. Materials and methods: Thirty-six patients (18 HPV-positive, 18 HPV-negative) were prospectively included. Presence of progressive disease was scored within one year. The mean pre-treatment ADC and NG-IVIM parameters in the gross tumor volume were compared between HPV-positive and HPV-negative patients. In HPV-negative patients, ADC and NG-IVIM parameters were compared between patients with and without progressive disease. Results: ADC, the NG-IVIM diffusion coefficient D, and perfusion fraction f were significantly higher, while pseudo-diffusion coefficient D* and kurtosis K were significantly lower in the HPV-negative compared to HPV-positive patients. In the HPV-negative group, a significantly lower D was found for patients with progressive disease compared to complete responders. No relation with ADC was observed. Conclusion: The results of our single-center study suggest that ADC is related to HPV status, but not an independent response predictor. The NG-IVIM parameter D, however, was independently associated to response in the HPV-negative group. Noteworthy in the opposite direction as previously thought based on ADC.

Validation of causal inference data using DirectLiNGAM in an environmental small-scale model and calculation settings.

The development of data science has been needed in environmental fields such as marine, weather, and soil data. In general, the datasets are large in some cases, but they are often small because they contain observation data that the analyses themselves are limited. In such a case, the data are statistically evaluated by increasing or decreasing the levels of factors using differential analysis, resulting in the essential factors are estimated. However, there is no consistent approach to the means of assessing strong associations as a group between factors. Causal inference method has the possibility to output effective results for small data, and the results are expected to provide important information for understanding the potential highly association between factors, not necessarily the inference with big data. Here, we describe essential checkpoints and settings for the calculation by a direct method for learning a linear non-Gaussian structural equation model (DirectLiNGAM) and validation methods for the calculation results by using DirectLiNGAM with small-scale model data as an additional discussion of DirectLiNGAM portion of the related research article. Thus, this study provides the statistical validation methods for the association networks, treatments, and interventions for structural inference as a group of essential factors.•Causal inference with DirectLiNGAM•Validation of correlation coefficient and feature importance•Validation using causal effect object and propensity scores.

Skewness and Staging: Does the Floor Effect Induce Bias in Multilevel AR(1) Models?

Multilevel autoregressive models are popular choices for the analysis of intensive longitudinal data in psychology. Empirical studies have found a positive correlation between autoregressive parameters of affective time series and the between-person measures of psychopathology, a phenomenon known as the staging effect. However, it has been argued that such findings may represent a statistical artifact: Although common models assume normal error distributions, empirical data (for instance, measurements of negative affect among healthy individuals) often exhibit the floor effect, that is response distributions with high skewness, low mean, and low variability. In this paper, we investigated whether-and to what extent-the floor effect leads to erroneous conclusions by means of a simulation study. We describe three dynamic models which have meaningful substantive interpretations and can produce floor-effect data. We simulate multilevel data from these models, varying skewness independent of individuals' autoregressive parameters, while also varying the number of time points and cases. Analyzing these data with the standard multilevel AR(1) model we found that positive bias only occurs when modeling with random residual variance, whereas modeling with fixed residual variance leads to negative bias. We discuss the implications of our study for data collection and modeling choices.

Staging liver fibrosis by a continuous-time random-walk diffusion model.

PURPOSE: Noninvasive assessment of liver fibrosis holds significant clinical importance. We aimed to evaluate the clinical potential of using a continuous-time random-walk diffusion model (CTRW) for staging liver fibrosis. METHODS: This prospective study included 52 patients suspected of liver disease and scheduled for liver biopsy. All patients underwent multi-b value diffusion-weighted imaging (DWI) using a 1.5 T MR scanner to derive the anomalous diffusion coefficient (D) and temporal (α) and spatial (ß) diffusion heterogeneity indexes sourced from the CTRW. The mono-exponential DWI-derived apparent diffusion coefficient (ADC), transient elastography-derived liver stiffness measurement (LSM), aspartate aminotransferase-to-platelet ratio index (APRI), and fibrosis-4 (FIB-4) index were calculated. We assessed and compared the correlations of these parameters with fibrosis stages and their efficacy in staging liver fibrosis. RESULTS: Significant correlations with fibrosis stages were found for APRI (r = 0.336), FIB-4 (r = 0.351), LSM (r = 0.523), D (r = -0.458), and ADC (r = -0.473). Significant differences were observed between APRI, LSM, D, and ADC of different fibrosis stages. The diagnostic performance of an index that combined D, α, ß, ADC, and LSM was superior to that of ADC or LSM alone for fibrosis stage F ≥ 2 and better than the index that combined D, α, ß for fibrosis stage F ≥ 4. CONCLUSIONS: Accurate liver fibrosis staging was achieved with a model that combined CTRW-derived parameters (D, α, and ß), ADC, and LSM. The model could serve as a reliable tool for noninvasive fibrosis evaluation.

Probing Renal Microstructure and Function with Advanced Diffusion MRI: Concepts, Applications, Challenges, and Future Directions.

Diffusion measurements in the kidney are affected not only by renal microstructure but also by physiological processes (i.e., glomerular filtration, water reabsorption, and urine formation). Because of the superposition of passive tissue diffusion, blood perfusion, and tubular pre-urine flow, the limitations of the monoexponential apparent diffusion coefficient (ADC) model in assessing pathophysiological changes in renal tissue are becoming apparent and motivate the development of more advanced diffusion-weighted imaging (DWI) variants. These approaches take advantage of the fact that the length scale probed in DWI measurements can be adjusted by experimental parameters, including diffusion-weighting, diffusion gradient directions and diffusion time. This forms the basis by which advanced DWI models can be used to capture not only passive diffusion effects, but also microcirculation, compartmentalization, tissue anisotropy. In this review, we provide a comprehensive overview of the recent advancements in the field of renal DWI. Following a short introduction on renal structure and physiology, we present the key methodological approaches for the acquisition and analysis of renal DWI data, including intravoxel incoherent motion (IVIM), diffusion tensor imaging (DTI), non-Gaussian diffusion, and hybrid IVIM-DTI. We then briefly summarize the applications of these methods in chronic kidney disease and renal allograft dysfunction. Finally, we discuss the challenges and potential avenues for further development of renal DWI. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 2.

Impact of gut microbiome on serum IgG4 levels in the general population: Shika-machi super preventive health examination results.

Introduction: Immunoglobulin G4 (IgG4) is a member of the human immunoglobulin G (IgG) subclass, a protein involved in immunity to pathogens and the body's resistance system. IgG4-related diseases (IgG4-RD) are intractable diseases in which IgG4 levels in the blood are elevated, causing inflammation in organs such as the liver, pancreas, and salivary glands. IgG4-RD are known to be more prevalent in males than in females, but the etiology remains to be elucidated. This study was conducted to investigate the relationship between gut microbiota (GM) and serum IgG4 levels in the general population. Methods: In this study, the relationship between IgG4 levels and GM evaluated in male and female groups of the general population using causal inference. The study included 191 men and 207 women aged 40 years or older from Shika-machi, Ishikawa. GM DNA was analyzed for the 16S rRNA gene sequence using next-generation sequencing. Participants were bifurcated into high and low IgG4 groups, depending on median serum IgG4 levels. Results: ANCOVA, Tukey's HSD, linear discriminant analysis effect size, least absolute shrinkage and selection operator logistic regression model, and correlation analysis revealed that Anaerostipes, Lachnospiraceae, Megasphaera, and [Eubacterium] hallii group were associated with IgG4 levels in women, while Megasphaera, [Eubacterium] hallii group, Faecalibacterium, Ruminococcus.1, and Romboutsia were associated with IgG4 levels in men. Linear non-Gaussian acyclic model indicated three genera, Megasphaera, [Eubacterium] hallii group, and Anaerostipes, and showed a presumed causal association with IgG4 levels in women. Discussion: This differential impact of the GM on IgG4 levels based on sex is a novel and intriguing finding.

Regularized Maximum Correntropy Criterion Kalman Filter for Uncalibrated Visual Servoing in the Presence of Non-Gaussian Feature Tracking Noise.

Some advantages of using cameras as sensor devices on feedback systems are the flexibility of the data it represents, the possibility to extract real-time information, and the fact that it does not require contact to operate. However, in unstructured scenarios, Image-Based Visual Servoing (IBVS) robot tasks are challenging. Camera calibration and robot kinematics can approximate a jacobian that maps the image features space to the robot actuation space, but they can become error-prone or require online changes. Uncalibrated visual servoing (UVS) aims at executing visual servoing tasks without previous camera calibration or through camera model uncertainties. One way to accomplish that is through jacobian identification using environment information in an estimator, such as the Kalman filter. The Kalman filter is optimal with Gaussian noise, but unstructured environments may present target occlusion, reflection, and other characteristics that confuse feature extraction algorithms, generating outliers. This work proposes RMCKF, a correntropy-induced estimator based on the Kalman Filter and the Maximum Correntropy Criterion that can handle non-Gaussian feature extraction noise. Unlike other approaches, we designed RMCKF for particularities in UVS, to deal with independent features, the IBVS control action, and simulated annealing. We designed Monte Carlo experiments to test RMCKF with non-Gaussian Kalman Filter-based techniques. The results showed that the proposed technique could outperform its relatives, especially in impulsive noise scenarios and various starting configurations.

Nonlinear back-end optimization method for VSLAM with multi-convex combined maximum correntropy criterion.

Back-end optimization plays a key role in eliminating the accumulated error in Visual Simultaneous Localization And Mapping (VSLAM). Existing back-end optimization methods are usually premised on the Gaussian noise assumption which does not always hold true due to the non-convex nature of the image and the fact that non-Gaussian noises are often encountered in real scenes. In view of this, we propose a back-end optimization method based on Multi-Convex combined Maximum Correntropy Criterion (MCMCC). A MCMCC-based cost function is first tailored for nonlinear back-end optimization in the context of VSLAM and the optimization problem is solved through Levenberg-Marquardt algorithm iteratively. Then, the proposed method is applied to ORB-SLAM3 to test its performance on public indoor and outdoor datasets. The real time performance is also validated using a RaceBot platform in real indoor and outdoor environments. In addition, the reprojection error is statistically analyzed to demonstrate the non-Gaussian characteristics in the back-end optimization process. Finally, the suggestion parameters are also provided through experiments for further study.

Accelerating High b-Value Diffusion-Weighted MRI Using a Convolutional Recurrent Neural Network (CRNN-DWI).

PURPOSE: To develop a novel convolutional recurrent neural network (CRNN-DWI) and apply it to reconstruct a highly undersampled (up to six-fold) multi-b-value, multi-direction diffusion-weighted imaging (DWI) dataset. METHODS: A deep neural network that combines a convolutional neural network (CNN) and recurrent neural network (RNN) was first developed by using a set of diffusion images as input. The network was then used to reconstruct a DWI dataset consisting of 14 b-values, each with three diffusion directions. For comparison, the dataset was also reconstructed with zero-padding and 3D-CNN. The experiments were performed with undersampling rates (R) of 4 and 6. Standard image quality metrics (SSIM and PSNR) were employed to provide quantitative assessments of the reconstructed image quality. Additionally, an advanced non-Gaussian diffusion model was employed to fit the reconstructed images from the different approaches, thereby generating a set of diffusion parameter maps. These diffusion parameter maps from the different approaches were then compared using SSIM as a metric. RESULTS: Both the reconstructed diffusion images and diffusion parameter maps from CRNN-DWI were better than those from zero-padding or 3D-CNN. Specifically, the average SSIM and PSNR of CRNN-DWI were 0.750 ± 0.016 and 28.32 ± 0.69 (R = 4), and 0.675 ± 0.023 and 24.16 ± 0.77 (R = 6), respectively, both of which were substantially higher than those of zero-padding or 3D-CNN reconstructions. The diffusion parameter maps from CRNN-DWI also yielded higher SSIM values for R = 4 (>0.8) and for R = 6 (>0.7) than the other two approaches (for R = 4, <0.7, and for R = 6, <0.65). CONCLUSIONS: CRNN-DWI is a viable approach for reconstructing highly undersampled DWI data, providing opportunities to reduce the data acquisition burden.

Quantum Bohmian-Inspired Potential to Model Non-Gaussian Time Series and Its Application in Financial Markets.

We have implemented quantum modeling mainly based on Bohmian mechanics to study time series that contain strong coupling between their events. Compared to time series with normal densities, such time series are associated with rare events. Hence, employing Gaussian statistics drastically underestimates the occurrence of their rare events. The central objective of this study was to investigate the effects of rare events in the probability densities of time series from the point of view of quantum measurements. For this purpose, we first model the non-Gaussian behavior of time series using the multifractal random walk (MRW) approach. Then, we examine the role of the key parameter of MRW, λ, which controls the degree of non-Gaussianity, in quantum potentials derived for time series. Our Bohmian quantum analysis shows that the derived potential takes some negative values in high frequencies (its mean values), then substantially increases, and the value drops again for rare events. Thus, rare events can generate a potential barrier in the high-frequency region of the quantum potential, and the effect of such a barrier becomes prominent when the system transverses it. Finally, as an example of applying the quantum potential beyond the microscopic world, we compute quantum potentials for the S&P financial market time series to verify the presence of rare events in the non-Gaussian densities and demonstrate deviation from the Gaussian case.

Parallel Kalman filter group integrated particle filter method for the train nonlinear operational status high-precision estimation under non-Gaussian environment.

For the problem of multi-mode state estimation in actual train operation, this paper proposes a nonlinear non-gaussian high-precision parallel Kalman filter group (NN-HEKFG) integrated Particle Filter. A multi-model Gaussian decomposition of the probability density function for state equations and measurement equations is performed, and each local state model is represented by a multi-dimensional high-order polynomial to establish the expanded dimensional state model. Then, by updating the mean and variance of the local state expanded dimensional model and in turn solving the particle filtering posterior probability density distribution function, the global estimation results are obtained. In reducing the number of Gaussian terms, a new parameter reduction criterion is established, which can effectively carry out the re-identification of parameters such as weights and means, so as to avoid the problem of parameter explosion. The superiority of NN-HEKFG over particle filters and Gaussian sum filters and its effectiveness for train running state estimation are verified by simulating the multi-model running state of trains.

Preoperative prediction of VETC in hepatocellular carcinoma using non-Gaussian diffusion-weighted imaging at high b values: a pilot study.

Background: Vessels encapsulating tumor clusters (VETC) have been considered an important cause of hepatocellular carcinoma (HCC) metastasis. Purpose: To compare the potential of various diffusion parameters derived from the monoexponential model and four non-Gaussian models (DKI, SEM, FROC, and CTRW) in preoperatively predicting the VETC of HCC. Methods: 86 HCC patients (40 VETC-positive and 46 VETC-negative) were prospectively enrolled. Diffusion-weighted images were acquired using six b-values (range from 0 to 3000 s/mm2). Various diffusion parameters derived from diffusion kurtosis (DK), stretched-exponential (SE), fractional-order calculus (FROC), and continuous-time random walk (CTRW) models, together with the conventional apparent diffusion coefficient (ADC) derived from the monoexponential model were calculated. All parameters were compared between VETC-positive and VETC-negative groups using an independent sample t-test or Mann-Whitney U test, and then the parameters with significant differences between the two groups were combined to establish a predictive model by binary logistic regression. Receiver operating characteristic (ROC) analyses were used to assess diagnostic performance. Results: Among all studied diffusion parameters, only DKI_K and CTRW_α significantly differed between groups (P=0.002 and 0.004, respectively). For predicting the presence of VETC in HCC patients, the combination of DKI_K and CTRW_α had the larger area under the ROC curve (AUC) than the two parameters individually (AUC=0.747 vs. 0.678 and 0.672, respectively). Conclusion: DKI_K and CTRW_α outperformed traditional ADC for predicting the VETC of HCC.