Finds initial cliques using a sparse PCA on bootstraps sub-samples

boot_FitSparsePCA(M, B, r, min.size = 1, cores = 1, unique = TRUE)

Arguments

M

Gaussian proxy for the original dataset. Possibly obtained using PLNmodels.

B

Number of bootstrap samples.

r

Number of missing actors.

min.size

Minimum number of neighbors of missing actors.

cores

Number of cores for parallel computation (uses mclapply, not available for Windows).

unique

Boolean for keeping only unique results.

Value

  • cliqueList: list of all cliques found. Each element is itself a list of size r.

  • nb_occ: vector of the number of times each cliques has been found by sPCA.

Examples

data=generate_missing_data(n=100,p=10,r=1,type="scale-free", plot=TRUE)
PLNfit=norm_PLN(data$Y) boot_FitSparsePCA(PLNfit$MO, B=100, r=1)
#> $cliqueList #> $cliqueList[[1]] #> $cliqueList[[1]][[1]] #> [1] 3 4 5 7 8 9 10 #> #> #> $cliqueList[[2]] #> $cliqueList[[2]][[1]] #> [1] 3 7 8 9 10 #> #> #> $cliqueList[[3]] #> $cliqueList[[3]][[1]] #> [1] 8 9 10 #> #> #> $cliqueList[[4]] #> $cliqueList[[4]][[1]] #> [1] 3 5 8 9 10 #> #> #> $cliqueList[[5]] #> $cliqueList[[5]][[1]] #> [1] 5 8 9 10 #> #> #> $cliqueList[[6]] #> $cliqueList[[6]][[1]] #> [1] 2 7 8 9 10 #> #> #> $cliqueList[[7]] #> $cliqueList[[7]][[1]] #> [1] 3 4 5 8 9 10 #> #> #> $cliqueList[[8]] #> $cliqueList[[8]][[1]] #> [1] 2 3 5 7 8 9 10 #> #> #> $cliqueList[[9]] #> $cliqueList[[9]][[1]] #> [1] 3 8 9 10 #> #> #> $cliqueList[[10]] #> $cliqueList[[10]][[1]] #> [1] 3 5 7 8 9 10 #> #> #> $cliqueList[[11]] #> $cliqueList[[11]][[1]] #> [1] 3 4 5 6 7 8 9 10 #> #> #> $cliqueList[[12]] #> $cliqueList[[12]][[1]] #> [1] 2 8 9 10 #> #> #> $cliqueList[[13]] #> $cliqueList[[13]][[1]] #> [1] 2 3 5 8 9 10 #> #> #> $cliqueList[[14]] #> $cliqueList[[14]][[1]] #> [1] 2 3 4 5 8 9 10 #> #> #> $cliqueList[[15]] #> $cliqueList[[15]][[1]] #> [1] 2 3 4 5 7 8 9 10 #> #> #> $cliqueList[[16]] #> $cliqueList[[16]][[1]] #> [1] 2 3 7 8 9 10 #> #> #> $cliqueList[[17]] #> $cliqueList[[17]][[1]] #> [1] 7 8 9 10 #> #> #> $cliqueList[[18]] #> $cliqueList[[18]][[1]] #> [1] 1 3 4 5 6 7 8 9 10 #> #> #> #> $nb_occ #> [1] 12 3 21 9 3 1 3 8 6 17 4 1 2 1 6 1 1 1 #>