ResampleEMtree: Resampling procedure for edges probability

In Rmomal/EMtree: Infers Direct Species Association Networks using Tree Averaging

Resampling procedure for edges probability

Description

Resampling procedure for edges probability

Usage

ResampleEMtree(
  counts,
  covar_matrix = NULL,
  unlinked = NULL,
  O = NULL,
  user_covariance_estimation = NULL,
  v = 0.8,
  S = 100,
  maxIter = 30,
  cond.tol = 1e-10,
  eps = 0.001,
  cores = 3,
  init = FALSE
)

Arguments

counts	Data of observed counts with dimensions n x p, either a matrix, data.frame or tibble.
covar_matrix	matrix of covariates, should have the same number of rows as the count matrix.
unlinked	An optional vector of nodes which are not linked with each other
O	Matrix of offsets, with dimension n x p
user_covariance_estimation	A user-provided function for the estimation of a covariance
v	The proportion of observed data to be taken in each sub-sample. It is the ratio (sub-sample size)/n
S	Total number of wanted sub-samples.
maxIter	Maximum number of EMtree iterations at each sub-sampling.
cond.tol	Tolerance for the psi matrix.
eps	Precision parameter controlling the convergence of weights beta
cores	Number of cores, can be greater than 1 if data involves less than about 32 species.
init	boolean: should the resampling be carried out with different initial points (TRUE), or with different initial data (FALSE)

Value

Returns a list which contains the Pmat data.frame, and vectors of EMtree maximum iterations and running times in each resampling.

• Pmat: S x p(p-1)/2 matrix with edge probabilities for each resample

• maxIter: EMtree maximum iterations in each resampling.

• times: EMtree running times in each resampling.

Examples

n=100
p=12
S=5
set.seed(2021)
simu=data_from_scratch("erdos",p=p,n=n)
G=1*(simu$omega!=0) ; diag(G) = 0
# With default evaluation, using the PLNmodel paradigm:
default_resample=ResampleEMtree(simu$data, S=S,cores = 1)

# With provided correlation estimation function:
estimSigma<-function(counts, covar_matrix, sample){
Dum_Sigma = cov2cor(cov(counts[sample,]))
}
custom_resample=ResampleEMtree(simu$data,S=S,cores = 1,user_covariance_estimation=estimSigma)

# We then run the stability selection to find the optimal selection frequencies,
# for a stability of 85%:
stab_default=StATS(default_resample$Pmat, nlambda=50, stab.thresh=0.8,plot=TRUE)
stab_custom=StATS(custom_resample$Pmat, nlambda=50, stab.thresh=0.8,plot=TRUE)

#Check quality of result
table(pred=1*(stab_default$freqs_opt>0.9), truth=ToVec(G))
table(pred=1*(stab_custom$freqs_opt>0.9), truth=ToVec(G))

Rmomal/EMtree documentation built on Dec. 14, 2024, 8:16 a.m.