If the list contains m elements, then lengths of each element must be m, too.
This function is intended to generate a list of x0 (the initial community)
with one missing species, to prepare the parameter simulations_compare in
estimateAFromSimulations
.