This R package implements the Random Forest with Canonical Correlation Analysis (RFCCA) method described in Alakus et al. (2021) <doi:10.1093/bioinformatics/btab158>. The theoretical details are given in the next section. The last section presents a data analysis using this method.

Multi-view data denotes many kinds of data that provide information about a subject from multiple sources. Integration of multiple feature sets and investigating the relationships between them may help to understand their interactions and obtain more meaningful interpretations. Canonical correlation analysis (CCA) is a multivariate statistical method that analyzes the relationship between two multivariate data sets, \(X\) and \(Y\). CCA searches for linear combinations of each of the two data sets, \(Xa\) and \(Yb\), having maximum correlation.

The multi-view data collection may also include some subject-related covariates such as age, gender and IQ as one of the views. We want to identify the linear association between two sets of variables, \(X\) and \(Y\), but we assume that this association can depend on a third set of subject-related covariates \(Z\).

We consider the following setting: let \(X\) and \(Y\) be two multivariate data sets of dimension \(p\) and \(q\), respectively, and let \(Z\) be a \(r\)-dimensional vector of subject-related covariates. We assume that the canonical correlation between \(X\) and \(Y\) depends on \(Z\). RFCCA uses an unsupervised random forest based on the set of covariates \(Z\) to find subgroups of observations with similar canonical correlations between \(X\) and \(Y\). This random forest consists of many unsupervised decision trees with a specialized splitting criterion \[\sqrt{n_Ln_R}*|\rho_L - \rho_R|\] where \(\rho_L\) and \(\rho_R\) are the canonical correlation estimations of left and right nodes, and \(n_L\) and \(n_R\) are the left and right node sizes, respectively. This splitting rule seeks to increase the canonical correlation heterogeneity. Following the random forest growing, RFCCA builds the Bag of Observations for Prediction (BOP), which is the set of training observations that are in the same terminal nodes as the observation of interest, for a new observation. Then, it applies CCA to the observations in BOP to estimate the canonical correlation of the new observation.

We perform a hypothesis test to evaluate the global effect of the
subject-related covariates on distinguishing between canonical
correlations. Define the unconditional canonical correlation between
\(X\) and \(Y\) as \(\rho_{\tiny CCA}(X,Y)\) which is found by
computing CCA with all \(X\) and \(Y\), and the conditional canonical
correlation between \(X\) and \(Y\) given \(Z\) as \(\rho(X,Y
| Z)\) which is found by RFCCA. If there is a global effect of
\(Z\) on correlations between \(X\) and \(Y\), \(\rho(X,Y |
Z)\) should be significantly different from \(\rho_{\tiny CCA}(X,Y)\). We conduct a
permutation test for the null hypothesis \[H_0 : \rho(X,Y | Z) = \rho_{\tiny
CCA}(X,Y)\] We estimate a *p*-value with the permutation
test. If the \(p\)-value is less than
the pre-specified significance level \(\alpha\), we reject the null
hypothesis.

We will show how to use RFCCA package on a generated data set. The
data set has three sets of variables: \(X \in
R^{n \times 2}\), \(Y \in R^{n \times
2}\) and \(Z \in R^{n \times
7}\). The sample size (\(n\)) is
300. The correlation between \(X\) and
\(Y\) depends only on \(Z_1,...,Z_5\) (*i.e.* \(Z_6\) and \(Z_7\) are noise variables). Firstly, we
load the data and split the data into train and test sets:

```
library(RFCCA)
data(data)
set.seed(2345)
<- sample(dim(data$X)[1], round(dim(data$X)[1]*0.7))
smp <- data$X[smp,]
train.X <- data$Y[smp,]
train.Y <- data$Z[smp,]
train.Z <- data$Z[-smp,] test.Z
```

Then, we train the random forest with `rfcca()`

.

`<- rfcca(X = train.X, Y = train.Y, Z = train.Z, ntree = 100) rf.obj `

Before analyzing the prediction results, we have to check first if the the global effect of \(Z\) is significant by applying the global significance test.

```
<- global.significance(X = train.X, Y = train.Y, Z = train.Z,
test.obj ntree = 100, nperm = 10)
$pvalue
test.obj#> [1] 0
```

Using 10 permutations, the estimated *p*-value < 0.05
(where the significance level(\(\alpha\)) is 0.05). Hence, we reject the
null hypothesis. 10 permutations to compute the *p*-value is
indeed very low. Using 500 permutations, we got a *p*-value of 0.
Higher number of permutations is required to estimate the
*p*-value accurately. However, it should also be noted that
higher number of permutations increases the computational time.

Since the global effect of \(Z\) is significant, now we can analyze the predictions. We can get the out-of-bag (OOB) predictions from the grow object.

```
<- rf.obj$predicted.oob
pred.oob head(pred.oob)
#> [1] 0.5124002 0.7959757 0.8662681 0.8149739 0.6588582 0.7188775
```

Those values represent the OOB canonical correlation predictions for the training observations. The predicted canonical correlations range between 0 and 1 where the closer the values to 1, the stronger the correlation between the canonical variates, \(Xa\) and \(Yb\).

We can get the variable importance (VIMP) measures for \(Z\). VIMP measures reflect the predictive power of subject-related covariates on the estimated correlations.

```
<- vimp(rf.obj)
vimp.obj <- vimp.obj$importance
vimp
vimp#> z1 z2 z3 z4 z5 z6 z7
#> 6.461857e-03 7.730602e-03 6.496649e-04 4.488815e-03 2.326590e-04 7.559608e-05 1.650719e-05
```

Moreover, we can plot the variable importance measures. We know that the correlation between \(X\) and \(Y\) depends on \(Z_1,...,Z_5\). Therefore, it is expected to get higher VIMP measures for \(Z_1,...,Z_5\) than \(Z_6\) and \(Z_7\).

`plot.vimp(vimp.obj)`

Finally, we can get canonical correlation predictions for test \(Z\).

```
<- predict(rf.obj, test.Z)
pred.obj <- pred.obj$predicted
pred head(pred)
#> [1] 0.7681228 0.8383421 0.8466399 0.7319380 0.7363355 0.8437405
```

Those values represent the canonical correlation predictions for the test observations.

Besides the canonical correlation estimations, we can get the estimated canonical coefficients for the \(X\) and \(Y\), \(a\) and \(b\), respectively.

```
head(pred.obj$predicted.coef$coefx)
#> x1 x2
#> [1,] -0.01856821 -0.005700326
#> [2,] 0.02100048 0.007258287
#> [3,] 0.01967372 0.007064432
#> [4,] -0.02207643 -0.010359656
#> [5,] 0.01890530 0.007314632
#> [6,] -0.01913989 -0.006072593
head(pred.obj$predicted.coef$coefy)
#> y1 y2
#> [1,] -0.010205081 -0.01649033
#> [2,] 0.016810998 0.01342039
#> [3,] 0.015730490 0.01412247
#> [4,] -0.019185912 -0.01341132
#> [5,] 0.009152462 0.01813504
#> [6,] -0.014482360 -0.01430064
```

We can analyse the linear relations between the variables with the entries of the estimated canonical coefficient vectors, \(a\) and \(b\). The values under the \(X_1\) and \(X_2\) columns represent the canonical coefficients for \(X_1\) and \(X_2\), \(a_1\) and \(a_2\), respectively. Similarly, \(b_1\) and \(b_2\) for the canonical coefficients of \(Y_1\) and \(Y_2\) under the corresponding columns.

Alternatively, we can use variants of CCA for the final canonical correlation estimation. Here, we use Sparse CCA (SCCA) for the final canonical correlation estimations. The idea of SCCA is to improve the interpretability of results by forcing some of the canonical coefficients to be zero which results in a set of selected features that has the highest importance in the relationship between \(X\) and \(Y\).

```
<- predict(rf.obj, test.Z, finalcca = "scca")
pred.obj2 <- pred.obj2$predicted
pred2 head(pred2)
#> [1] 0.6867026 0.7213607 0.6790587 0.6223460 0.6776450 0.6902146
```

Again, we can look at the estimated canonical coefficients for the \(X\) and \(Y\).

```
head(pred.obj2$predicted.coef$coefx)
#> x1 x2
#> [1,] -1 0
#> [2,] -1 0
#> [3,] -1 0
#> [4,] -1 0
#> [5,] -1 0
#> [6,] -1 0
head(pred.obj2$predicted.coef$coefy)
#> y1 y2
#> [1,] 0 -1
#> [2,] -1 0
#> [3,] -1 0
#> [4,] -1 0
#> [5,] 0 -1
#> [6,] 0 -1
```

We can plot the canonical correlation estimations using classical CCA and Sparse CCA.

```
plot(pred, pred2, xlab="Classical CCA", ylab="Sparse CCA",
xlim=c(min(pred,pred2),max(pred,pred2)), ylim=c(min(pred,pred2),max(pred,pred2)),
pch = 20)
abline(0,1,col = "red")
```

For further data analysis, one can use SHAP values (Lundberg et al., 2017) to gain additional insights. SHAP values show how much each variable contributes, either positively or negatively, to the individual predictions. For an example of application to the problem in question, see Alakus et al. (2020).

Alakus, C., Larocque, D., Jacquemont, S., Barlaam, F., Martin, C.-O., Agbogba, K., Lippe, S., and Labbe, A. (2021). Conditional canonical correlation estimation based on covariates with random forests. Bioinformatics, 37(17), 2714-2721.

Lundberg, S. M., & Lee, S. I. (2017). A unified approach to interpreting model predictions. In Advances in neural information processing systems (pp. 4765-4774).

```
sessionInfo()
#> R version 4.2.0 (2022-04-22)
#> Platform: x86_64-apple-darwin17.0 (64-bit)
#> Running under: macOS Monterey 12.6.3
#>
#> Matrix products: default
#> LAPACK: /Library/Frameworks/R.framework/Versions/4.2/Resources/lib/libRlapack.dylib
#>
#> locale:
#> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
#>
#> attached base packages:
#> [1] stats graphics grDevices utils datasets methods base
#>
#> other attached packages:
#> [1] RFCCA_1.0.11 testthat_3.1.4
#>
#> loaded via a namespace (and not attached):
#> [1] viridis_0.6.2 pkgload_1.2.4 maps_3.4.0 viridisLite_0.4.0 splines_4.2.0 dotCall64_1.0-1
#> [7] rainbow_3.6 brio_1.1.3 hdrcde_3.4 highr_0.9 remotes_2.4.2 sessioninfo_1.2.2
#> [13] pillar_1.8.1 lattice_0.20-45 glue_1.6.2 colorspace_2.1-0 Matrix_1.4-1 pcaPP_2.0-1
#> [19] pkgconfig_2.0.3 devtools_2.4.3 fda_6.0.5 purrr_1.0.1 mvtnorm_1.1-3 scales_1.2.0
#> [25] processx_3.5.3 pracma_2.3.8 tibble_3.1.8 generics_0.1.3 ggplot2_3.3.6 usethis_2.1.5
#> [31] ellipsis_0.3.2 cachem_1.0.6 withr_2.5.0 cli_3.6.0 magrittr_2.0.3 crayon_1.5.1
#> [37] mclust_5.4.10 memoise_2.0.1 fds_1.8 evaluate_0.15 ps_1.7.0 ks_1.13.5
#> [43] fs_1.5.2 fansi_1.0.3 MASS_7.3-57 pkgbuild_1.3.1 tools_4.2.0 prettyunits_1.1.1
#> [49] lifecycle_1.0.3 stringr_1.4.0 munsell_0.5.0 cluster_2.1.3 callr_3.7.0 compiler_4.2.0
#> [55] rlang_1.0.6 grid_4.2.0 RCurl_1.98-1.6 rstudioapi_0.13 CCA_1.2.1 spam_2.8-0
#> [61] bitops_1.0-7 gtable_0.3.0 deSolve_1.32 PMA_1.2.1 R6_2.5.1 gridExtra_2.3
#> [67] knitr_1.39 dplyr_1.1.0 fastmap_1.1.0 utf8_1.2.2 rprojroot_2.0.3 KernSmooth_2.23-20
#> [73] desc_1.4.1 stringi_1.7.8 parallel_4.2.0 fields_14.0 vctrs_0.5.2 tidyselect_1.2.0
#> [79] xfun_0.31
```