# This function computes an estimator of the covariance matrix and the square root of its inverse and permutes its rows and columns if it is necessary to make the block structure appear. ```r Sigma_estimation(E, k = NULL, nb_nn0 = NULL, big = FALSE, reorder = FALSE, inv_12 = FALSE, method_k = "Cattell", times = 10, method_0 = "Elbow", N = 10) ``` ## Arguments - `E`: the observation matrix such that each of its row has a block structure correlation matrix Sigma to estimate up to a permutation of its columns and rows. - `k`: numerical or NULL, the rank for the low rank approximation. If NULL the rank is computed using the slope_change function applied on the eigenvalues of the low rank part of Sigma. Default to NULL. - `nb_nn0`: numerical or NULL, corresponds to the number of non null values to keep in the estimation of the covariance matrix. If NULL the number of non null values is computed using the slope_change function to the Frobenius norm of the difference between the empirical correlation matrix and its estimation with nb_nn0 non null values. Default to NULL. - `big`: logical, default to FALSE. If the dataset is too big the empirical correlation is calculated by crossprod(E) * 1 / n to fasten the computation - `reorder`: logical, default to FALSE. Whether or not the columns of E are permuted. If TRUE a hierarchical clustering is first performed and the columns are permuted according to it. - `inv_12`: logical, default to FALSE Whether or not computing the square root of the inverse of the covariance matrix. - `method_k`: character if "Cattell" (the default) then the Cattell criterion if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="cattell1966",package="BlockCov",cached_env=.Rdpack.currefs) is performed on the singular values of the covariance matrix. to estimate the number of rank use in the low rank approximation, while "PA" use the parrallel analysis if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="Horn1965",package="BlockCov",cached_env=.Rdpack.currefs) wich can be more accurate if the number of rows of E is not to small but which is much slower. - `times`: numeric the number of resampling done for the "PA" method, ignored if metod_k is different from "PA". - `method_0`: character if "Elbow" (the default) then the Elbow criterion (see if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="blc;textual",package="BlockCov",cached_env=.Rdpack.currefs) for details) is performed to estimate the number of rank use in the low rank approximation, while "BL" use the approach proposed in if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_citeOnly(keys="bickel2008;textual",package="BlockCov",cached_env=.Rdpack.currefs) based on cross-validation wich can be more accurate if the number of rows of E is not to small but which is much slower. - `N`: numeric the number of fold used for the "BL" method. Ignored if method_0 is different from "BL" ## Returns A list with the elements - **Sigma_est**: estimator of the covariance matrix - **k**: rank of the low rank part of the covariance matrix - **nb_nn0**: number of non null values of the upper triangular part of the covariance matrix - **S_inv_12**: square root of the inverse of the estimated covariance matrix - **order**: permutation to apply to the rows and the columns of the covariance to make the block structure appear ## Examples ```r n <- 30 q <- 100 Sigma <- Simu_Sigma(q = q, diag = FALSE, equal = TRUE) Matrix::image(Sigma) E <- matrix(rnorm(n * q), ncol = q) %*% chol(as.matrix(Sigma)) res <- Sigma_estimation(E, inv_12 = TRUE) Matrix::image(res$Sigma_est) Matrix::image(res$S_inv_12) ``` ## References if(!exists(".Rdpack.currefs")) .Rdpack.currefs \<-new.env();Rdpack::insert_all_ref(.Rdpack.currefs)

Sigma_estimation function