A stochastic individual-based model (IBM) of mosquitofish population dynamics in experimental ponds was constructed in order to increase, virtually, the number of replicates of control populations in an ecotoxicology trial, and thus to increase the statistical power of the experiments. In this context, great importance had to be paid to model calibration as this conditions the use of the model as a reference for statistical comparisons. Accordingly, model calibration required that both mean behaviour and variability behaviour of the model were in accordance with real data. Currently, identifying parameter values from observed data is still an open issue for IBMs, especially when the parameter space is large. Our model included 41 parameters: 30 driving the model expectancy and 11 driving the model variability. Under these conditions, the use of “Latin hypercube” sampling would most probably have “missed” some important combinations of parameter values. Therefore, complete factorial design was preferred. Unfortunately, due to the constraints of the computational capacity, cost-acceptable “complete designs” were limited to no more than nine parameters, the calibration question becoming a parameter selection question. In this study, successive “complete designs” were conducted with different sets of parameters and different parameter values, in order to progressively narrow the parameter space. For each “complete design”, the selection of a maximum of nine parameters and their respective n values was carefully guided by sensitivity analysis. Sensitivity analysis was decisive in selecting parameters that were both influential and likely to have strong interactions. According to this strategy, the model of mosquitofish population dynamics was calibrated on real data from two different years of experiments, and validated on real data from another independent year. This model includes two categories of agents; fish and their living environment. Fish agents have four main processes: growth, survival, puberty and reproduction. The outputs of the model are the length frequency distribution of the population and the 16 scalar variables describing the fish populations. In this study, the length frequency distribution was parameterized by 10 scalars in order to be able to perform calibration. The recently suggested notion of “probabilistic distribution of the distributions” was also applied to our case study, and was shown to be very promising for comparing length frequency distributions (as such).