In this version I have:
replaced mclapply with future_lapply,
allowed the adjacency matrix to be set to NULL,
changed the order of the arguments in scov,
implemented the option joint_estimation=TRUE.