An analytic approach for interpretable predictive models in high‐dimensional data in the presence of interactions with exposures

Abstract

Predicting a phenotype and understanding which variables improve that prediction are two very challenging and overlapping problems in the analysis of high‐dimensional (HD) data such as those arising from genomic and brain imaging studies. It is often believed that the number of truly important predictors is small relative to the total number of variables, making computational approaches to variable selection and dimension reduction extremely important. To reduce dimensionality, commonly used two‐step methods first cluster the data in some way, and build models using cluster summaries to predict the phenotype. It is known that important exposure variables can alter correlation patterns between clusters of HD variables, that is, alter network properties of the variables. However, it is not well understood whether such altered clustering is informative in prediction. Here, assuming there is a binary exposure with such network‐altering effects, we explore whether the use of exposure‐dependent clustering relationships in dimension reduction can improve predictive modeling in a two‐step framework. Hence, we propose a modeling framework called ECLUST to test this hypothesis, and evaluate its performance through extensive simulations. With ECLUST, we found improved prediction and variable selection performance compared to methods that do not consider the environment in the clustering step, or to methods that use the original data as features. We further illustrate this modeling framework through the analysis of three data sets from very different fields, each with HD data, a binary exposure, and a phenotype of interest. Our method is available in the eclust CRAN package.

Parallel and nonparallel genomic responses contribute to herbicide resistance in Ipomoea purpurea, a common agricultural weed

“…gene command from the ape package (vers 5.0; [100]) in R [101]. This genetic distance matrix was reduced to 2 dimensions by multidimensional scaling using the cmdscale and eclust commands [102]. These two dimensions were then used to hierarchically cluster the populations into 2 groups using kmeans clustering.…”

Blood transcriptome differentiates clinical clusters for asthma

“…In the subsequent phase, we employ K-means clustering using ECLUST, 21 considering the cluster counts and centroids established in the initial phase. The resulting cluster are shown in Fig.…”

Operation of a TCA cycle subnetwork in the mammalian nucleus

“…The input values for each putative interactor corresponded to the average (across all biological replicates) log 2 fold changes in the abundance of the putative interactor in cells expressing a bait enzyme relative to the parental HeLa cells. Hierarchical clustering was performed on the centered and scaled input values with the “eclust” package ( 65 ) using the “ward.D2” linkage and the “Euclidean” distance metric. For visualization, the “ComplexHeatmap” package ( 66 ) was used in R.…”