A robust-design formulation of the incidence function model of metapopulation dynamics applied to two species of rails