The preset invention provides a method of obtaining an array of a plurality of regions on a substrate, including obtaining, using a prior distribution representing a probability distribution of parameters of a regression model used to estimate the array, a first posterior distribution representing the probability distribution of the parameters, obtaining, using the first posterior distribution as the prior distribution representing the probability distribution of the parameters, a second posterior distribution representing the probability distribution of the parameters, and updating the regression model by deciding the parameters based on the second posterior distribution and obtaining, using the updated regression model, the array of the plurality of regions on a substrate from the second position measurement data.