In this paper we study the impact of cosmic variance and observational uncertainties in constraining the mass and occupation fraction, , of dark matter (DM) halos hosting Lyα-emitting galaxies (LAEs) at high redshift. To this end, we construct mock catalogs from an N-body simulation to match the typical size of observed fields at z = 3.1 (). In our model a DM halo with mass in the range can only host one detectable LAE at most. We proceed to explore the parameter space determined by , , and with a Markov Chain Monte Carlo algorithm using the angular correlation function and the LAEs' number density as observational constraints. We find that the preferred minimum and maximum masses in our model span a wide range , followed by a wide range in the occupation fraction . As a consequence, the median mass, M50, of all the consistent models has a large uncertainty . However, we find that the same individual models have a relatively tight scatter around the median mass dex. We are also able to show that is uniquely determined by , regardless of . We argue that upcoming large surveys covering at least 25 deg2 should be able to put tighter constraints on and through the LAE number density distribution width constructed over several fields of ∼1 deg2.