Deciphering Interactions in Moving Animal Groups

1Collective motion occurs across a variety of scales in nature, offering a wealth of fascinating phenomena which have attracted a lot of attention –. The self-organized motion of social animals is particularly intriguing because the behavioral rules the individuals actually follow and from which these remarkable collective phenomena emerge often remain largely unknown due to the tremendous difficulties to collect quality field data and/or perform controlled experiments in the laboratory. This situation does not prevent a thriving modeling activity, thanks to the relative ease by which numerical simulations can be conducted. However, most models of moving animal groups are built from general considerations, educated guesses following qualitative observations, or ideas developed along purely theoretical lines of thought –. Even when authors strive to build a model from data, as in the recent paper by Lukeman et al., this model building amounts to writing down a fairly complicated structure a priori, involving many implicit assumptions, and to fit collective data to determine effective parameters, yielding a best-fit model.

2On the other hand, recent studies within the physics community of simple, minimal models for collective motion have revealed an emerging picture of universality classes –: Take, for instance, the Vicsek model, arguably one of the simplest models exhibiting collective motion. In this model, point particles move at constant speed and choose, at discrete time-steps, their new heading to be the average of that of their neighbors located within unit distance. Many of these behavioral restrictions can be relaxed without changing the emerging collective properties. Fluctuations of speed can be allowed, some short-range repulsion (conferring a finite size to the particles) can be added, even explicit alignment can be replaced by inelastic collisions, etc., all these changes will still produce the remarkable nonlinear high-density high-order bands emerging near onset of collective motion, and, deeper in the ordered moving phase, the anomalously strong number fluctuations which have become a landmark of the collective motion of polarly aligning self-propelled particles –. The Vicsek model, in this context, is one of the simplest members of a large universality class defined by all models sharing the same large-scale properties. This universality class can be embodied in the continuous field equations that physicists are now able to derive. With such a viewpoint, different models in this class merely differ in the numerical values of their parameters –, very much like different fluids are commonly described by the Navier-Stokes equations and differ only in their viscosity and other constitutive parameters.

3Significant features nevertheless may be altered when a qualitatively important feature is changed, such as the symmetry of the aligning interaction, or added, as when local attraction/repulsion between individuals is also considered, In this latter case, for instance, no strong clustering and high density band appears when attraction is sufficiently strong, and finite groups may keep cohesion in open space as most natural groups do. These models yield a more complex phase diagram where collectively moving groups may assume gas-like, liquid-like or even moving crystal states as the two parameters controlling alignment and cohesion are varied.

4So, it remains important to know how individuals make behavioral choices when interacting with others, not only from a social ethology and cognitive viewpoint, but also because i) different behavioral rules may make a difference in small enough groups and ii) the analysis of local-scale data that this requires may lead to discover features eventually found to give rise to different qualitative collective properties. A recent instance can be found in the results on the structure of starling flocks gathered by Ballerini et al.: They have ignited an ongoing debate about the possibility that individuals might interact mostly with neighbors determined by topological rules and not by metric criteria as assumed in most models. While this message has intrinsic value for the study of decision-making processes in animal groups, it was also shown recently that such metric-free, topological interactions are relevant, in the sense that they give rise to collective properties that are qualitatively different from those of metric models. Thus, in this case, an individual-level ingredient suggested by data, which had been only partially and theoretically considered before,,, defines new classes of collective properties. Given that animals are likely to possess more sophisticated behavior than, say, sub-cellular filaments displaced by molecular motors, one can expect more hidden features to play an important role at the collective level. This is a central finding of the recent work by Katz et al. where a careful analysis of groups of two and three fish revealed that the mechanisms at play are, at least in the golden shiners studied there, much more subtly intertwined that in existing fish models. Indeed they concluded that alignment emerges from attraction and repulsion as opposed to being an explicit tendency among fish. Whether fish display some mechanisms of active alignment or only attraction/repulsion is likely to lead to different patterns as interactions accumulate over time. In short, extracting interaction rules from individual scale data is crucial not only for animal behavior studies, but also because heretofore overlooked features can be found decisive in governing the emergent collective properties of moving animal groups.

5Here, we assess the power of a bottom-up methodology to build models for animal group motion from data gathered at the individual scale in groups of increasing sizes. We use data obtained by recording the motion of barred flagtails ( Kuhlia mugil) in a tank. In natural conditions, the barred flagtail form schools with a few thousands individuals along the reef margin of rocky shorelines, from just below the breaking surf to a depth of a few meters. However the size of these schools is much smaller than in species like the sardine or the Atlantic herring.

6Our analysis is incremental: in a previous work we characterized the spontaneous behavior of a single fish, including wall-avoidance behavior. Here, using pairs of fish, we first characterize the response function of one fish depending on the position and orientation of the other fish. Then we calibrate multiple fish interactions, using data in larger groups. At each step, the already-determined factors and parameters are kept unchanged and the new terms introduced in the stimulus-response function and the corresponding new parameters are determined from data with nonlinear regression routines (see Statistical Analysis in Materials and Methods). The resulting model is validated by comparing extensive simulations to the original data. Often, different functional forms are tested and we determine which one is most faithful to the data. When no significant difference is found, the simplest version is retained, following a principle of parsimony.

7Experiments with 1 to 30 fish were performed in shallow circular swimming pools that let the fish form quasi 2-dimensional schools (see Fig. 1A and Video S1, S2, S3, S4). At the collective level, we observe a transition from schooling to shoaling behavior when the density of fish increases in the tank: the group polarization, which measures the degree of alignment, is high in groups of two and five fish, even if sometimes we do observe some breaks in the synchronization, while in larger groups, when, it remains low (Fig. 1B). Within each group size, we notice some variability, the most striking effect being an increase of the synchronization level with the individuals velocity in groups of two fish.

8For every group size, fish move continuously and quickly synchronize their speed to a well defined, but replicate-dependent value (Fig. S1). The fish trajectories are smooth, differentiable and the instantaneous speed has a well-defined mean and root mean square fluctuations of about 10–20% which are found to be uncorrelated to, the angular velocity of the fish orientation (Fig. S2). On this basis, fish can be modeled as self-propelled particles moving in 2D space at constant speed and the only dynamical variable retained is. Moreover, since the recorded trajectories, be they extracted from a single fish or from small groups in the tank, are always irregular/stochastic, our model takes the form of coupled stochastic differential equations for the angular velocities of each fish. Note that if noise acts on rather than the fish position or heading, trajectories are smooth and differentiable, as observed.

9We have shown elsewhere that single fish trajectories in barred flagtails are very well described by an Ornstein-Uhlenbeck process acting on the instantaneous curvature, or, equivalently, on. When the fish is away from the tank wall, the distribution of is nearly Gaussian with zero mean and variance, where is the characteristic time of the (exponentially decaying) autocorrelation function of. To avoid collisions with the tank walls, we found that a single fish adjusts its current turning speed towards a (time-dependent) target value where is a parameter, is the distance to the point of impact on the wall should the fish continue moving straight ahead, and is the angle between the current heading of the fish and the normal to the point of impact (see Fig. 2A). In short, obeys the stochastic differential equation:where is a Wiener process of variance reflecting the stochasticity of the behavioral response. Non-linear regression analysis of the above model against our experimental data yielded excellent agreement and accurate estimations of and. Note that in the present work we adopted a slightly different form for the wall avoidance term with regards to the exponentially decreasing one of Ref., since it actually prevents fish from crossing the tank boundary, while both ansatz are similar as fish moves away from tank walls (Fig. S8A).

10The stimulus/response function of a single fish in the tank is directly expressed by how varies with the relative position of the fish and the wall. We now assume that this framework holds when two fish and are present in the tank by defining how, for fish, its turning speed is modulated by the combined stimuli due to the wall and to fish. Almost all existing fish behavior models, on the basis of common sense, intuition, and sometimes experimental evidence –, offer a combination of three basic ingredients: short distance repulsion (to avoid collisions), alignment for intermediate distances, and attraction up to some maximal range. Here, we dispose of repulsion not only because we want to allow for the rare experimentally observed over- and under-passings events, but mostly because we do not need to incorporate it explicitly to avoid collisions (see below and Video S1, S2, S3, S4). In contrast with most existing “zonal” models, and because there is little cognitive/physiological evidence for a sudden switch between alignment and attraction, we want to allow for continuous, distance-dependent weighting between alignment and attraction in agreement with the recent findings of Katz et al.. These two factors a priori depend on the geometrical quantities defining the location of fish from the viewpoint of fish: their distance is involved, but also, the angular position of fish with respect to, the current heading of fish, as well as their relative heading difference (Fig. 2A). The main angular variable for explicit alignment is, as usual,, whereas for attraction it is; both may also depend on. The stimulus/response function of fish thus combines a priori wall avoidance, alignment and attraction in some unknown function with parameters and (reaction to the wall),, and:.

11Next, in the spirit of an expansion around the no-interaction case, we write the expression for above as the sum of three terms:where the “main” variables have been placed first for each term. The wall avoidance term depends explicitly on to reflect a possible screening of the wall by the other fish. We have tested the influence of this by introducing a dependence in the wall avoidance term determined for the single-fish behavior. Essentially, was made smaller for. But this brought no significant improvement, so we keep as found previously.

12On general grounds, one expects that the relative importance of the positional interaction (attraction) to the velocity interaction (alignment) increases with. Given that the fish are constrained in a rather small tank, a limited range of inter-distances is effectively explored. In the spirit, again, of a small-distance expansion, a satisfactory choice is given by a linear dependence of on, while is independent of. Of course, such a functional choice cannot be correct at large distances since then would take large unrealistic values, meaning that the fish would spend enormous amounts of energy turning toward a distant “neighbor” (see the Discussion for more comments on this point).

13The attraction interaction must depend on, the relative angle with the other fish position: it is reasonable to assume that a fish is not attracted much towards a neighbor located behind, and of course this term must be zero when the other fish is right ahead, yielding. A simple, compatible, trigonometric function representing the leading term of a Fourier expansion is the sine function. We thus write where is a parameter controlling the weight of the positional information. Finally, we neglect the possible dependence on: the way a fish would turn toward the position of a neighbor does not depend on the orientation of that fish. This is especially natural when this interaction dominates, i.e. when the neighbor is far away. Moreover knowing the other fish orientation is a cognitively expensive and/or time consuming process at larger distances.

14The alignment interaction is mostly characterized by its functional dependence on. The main constraint here is that (the two fish are then already aligned). Here again, the simplest choice is as in most models –. Including higher harmonics (e.g. ) would allow to account for the few observed nematic alignment events where a fish remains anti-aligned with its neighbors. However, incorporating this term did not improve the faithfulness of the model to our dataset, so we keep only the leading sine function. In principle, the strength of alignment can also depend on: less attention may be paid to “back neighbors”. We have tested simple and reasonable choices for the dependence of on, e.g., but this did not lead to significant improvement so we kept no angular position dependence in the alignment interaction. We thus write, finally: where is a parameter controlling the weight of the orientational information.

15To summarize the case of two fish and, the stimulus/response function in the general evolution equation (1) is thus finally written:

16Using nonlinear regression analysis, the faithfulness to our data of the model consisting of Eqs. (1) and (3) was found very good for each of our two-fish recordings and the 5 parameters,,, and were estimated for each fish. We find clear dependences of the estimated parameters on, the average speed of each fish (see Fig. 3A). In particular,,, and are found proportional to, whereas and no significant -dependence appears for. Results regarding this last parameter are the least convincing, with a large dispersion of individual values. This is mostly due to the confinement of fish in the tank: the positional interaction never dominates alignment, preventing its accurate estimation. Nevertheless it is crucial to note here that without these positional interactions the model fails to match the data. Furthermore, we have tested a posteriori our ansatz by testing each contribution (either wall avoidance, neighbor position or neighbor orientation) after the other twos have been subtracted from the fish response according to Eq. (3). Results show an excellent agreement between our ansatz and the mean fish response (for more details see Fig. S8 B–D).

17Note that these results mean also that the wall avoidance is actually governed by, the time it would take the fish to hit the wall, rather than the distance. Conversely,, the relaxation time of the angular velocity, is better expressed as the ratio between a characteristic length and the speed. These -dependences were then incorporated explicitly in the model:withwhere, and are now constants over all fish. Running again our nonlinear regressions using this form, and using data for all replicate, allows for a more accurate estimation of the parameters,,, and now the same for all fish. We find,,, and.

18To validate this experimental finding, these parameter values were used in simulations of the model which were compared directly to the data. Good agreement is found not only for statistical quantifiers of the emergent synchronization between the two fish (see Fig. 3C), but in fact also for the dynamics: see for instance Video S1, S2, S5, S6 and the time series of polarization which show the same intermittent behavior (Fig. 3B). We emphasize that the model captures the experimental observation that the orientational order is lower when the swimming speed is lower, and is better in faster groups (Fig. 3B, C).

19Can multiple-fish interactions be factorized into pairs? This is often taken for granted, following a typical physics approach where this assumption is routinely made. However, recent work has suggested that this is not valid when describing pedestrian interactions in a crowd. Even more recently, Katz et al. argued that this is also the case for groups of three golden shiners (but see for the case of birds). Here, our data set is too small to allow for an in-depth analysis of group behavior at the level of detail that was accomplished above for two fish, mostly because many more variables are involved, but the quality of the pair approximation can be evaluated a posteriori. Assuming that multiple fish interactions are indeed essentially made of the sum of the pair interactions involved, Eq. (5) is extended towhere is the (current) neighborhood of fish which contains individuals. In our observations with fish, individuals mostly stayed together, suggesting that individuals remains aware of all others. Using all-to-all, equal-weight coupling, we found good agreement between data and simulations of Eqs. (4) and (6) (see Fig. S3). This justifies a posteriori the factorization in pairs and the use of two-fish parameters for groups, but also the overall normalization factor in Eq. (6), which indicates that, in the stimulus response of a fish, wall avoidance and the averaged influence of neighbors keep, on average, the same relative importance irrespective of the group size. The raw, “force-like” un-normalized superposition would yield too strong a coupling.

20For the larger group sizes, all-to-all equal-weight coupling quickly becomes unrealistic, and one must determine the set of neighbors a fish interacts with. In principle, abundant data recorded in larger tanks would allow to discriminate between alternative choices, but our experimental recordings are too short for this. Nevertheless, many choices can be eliminated: the usual one, which consists in cutting off interactions at fixed distances (zonal models), is inconsistent with our continuous weighting of alignment and attraction with fish inter-distance. Based on an analysis of starling flocks, Ballerini et al. have argued that these birds actually pay attention to their 6–8 closest neighbors, irrespective of the density of the flock. Coming back to our observations, this non-metric choice of neighbors can, however, lead to unrealistic situations when, for instance, a fish is leading a small group, since then this fish will only pay attention to those behind, even if individuals are located at intermediate distances ahead (but see Fig. S7). A simple, reasonable, non-metric solution is that of neighbors determined by the Voronoi tessellation around each individual: this allows for continuous weighting between alignment and attraction and avoids the caveat mentioned above in the case of a fixed number of closest neighbors. Moreover, given the rather small inter-distances observed, individuals beyond the first shell of Voronoi neighbors are largely screened out, so that our final choice was that of the first shell of Voronoi neighbors (see Fig. 2B). Using this, the validation of the model simulated with fish using the parameters is again quite satisfactory (see Fig. S3).

21This is however not true anymore for larger groups which display too high a polarization when using the parameters (whereas distance predictions remains satisfactory, see Fig. S3). Our approach actually allows to further investigate this discrepancy. We estimate the parameters at the individual scale for each fish with our nonlinear least-square procedure using the Ito-integrated version of the Ornstein-Uhlenbeck process of Eqs (4) and (6) for each fish time series (see Statistical Analysis). Thanks to this parametric inversion strategy, we have been able to extract the parameter values for each replicate separately (Fig. 4A). The model predictions with these replicate-based parameters yield a near-perfect match with the data (Fig. 4B). The results confirm that, within the limits of statistical accuracy, the parameters and their v-dependence remain about the same up to N = 10, in agreement with the above findings; but in larger groups there is a decreased tendency of fish to react to their neighbors, which both concerns the alignment and positional interactions (Fig. 4A).

22Characterizing and modeling the interactions between individuals and their behavioral consequences is a crucial step to understand the emergence of complex collective animal behaviors. With the recent progress in tracking technologies, high precision datasets on moving animal groups are now available, thus opening the way to a fine-scale analysis of individual behavior, –. Here we adopted a bottom-up modeling strategy for deciphering interactions in fish shoaling together. This strategy is based on a step-by-step quantification of the spontaneous motion of a single fish and of the combined effects of local interactions with neighbors and obstacles on individuals motion. At each step, one model ingredient is considered and checked against experimental data. The required parameters are determined using a dedicated inversion procedure and the numerical values of these parameters are kept unchanged in the following steps, yielding, in the end, a model without any free parameter. Such an incremental procedure fosters the explicit enunciation of the rationale behind each functional choice, and differs from searching the best set of free parameters to fit large-group data,. Proceeding step by step also puts stronger constraints on matching, since the incorporation of additional behavioral features at each step assumes the stability of the previously explored behaviors and of the corresponding model parameters. Using pairs of fish, we were able to show how positional and directional stimuli combine, and the crucial role of the swimming speed in the alignment interaction. At intermediate sizes, multiple fish interactions could be faithfully factorized into pair interactions albeit in a normalized form. However we found that at even larger group sizes our incremental modeling approach fails to accurately reproduce the collective dynamics.

23We explored this point further, still considering the statistical behavior of each fish separately, but only using the data corresponding to the large-group experiments. We concluded that our model could still grasp the observed individual and collective features but with smaller positional and alignment coefficients. We believe that this decrease in reactivity to neighbors is a consequence of the high density already imposed by confinement effects. Indeed, our model predicts that large groups adopting the high neighbor reactivity found in smaller groups would remain polarized also in open space, keeping group cohesion with an average distance to neighbors of about two body lengths (Fig. S6). Since the largest groups we observed in the tank are already characterized by such a typical neighbors distance due to confinement effects, we argue that lower interaction strengths may simply indicate the fish vanishing need to actively react to neighbors position and heading in order to maintain a high density. This could be, for instance, a physiological consequence of the density per se: the physiological and behavioral consequences, for an individual, of living in dense groups, known as group effect, have been described in numerous species from insects to vertebrates,. Our results investigation suggests that this sensitivity may be represented in a quite straightforward manner, preserving the model shape of Eqs. (4) and (6) and only modifying the interaction parameters. This conjecture, of course, could only be validated by experiments on large groups conducted in open space or larger tanks. While we believe in a positive answer, namely that without too strong a confinement, individuals would react to the perceived neighbors the same way regardless of the overall group size, we leave this question for future investigations on group effect in fish schools.

24Our approach yielded a novel type of fish school model whose main features are its built-in balancing mechanism between positional and orientational information, a topological interaction neighborhood, and explicit dependencies on fish speed. Note that similar features were recently uncovered for another species thanks to a novel data analysis procedure. The smooth transition from a dominant alignment reaction when a neighbor is close to attraction when it is far away is in line with a simple additive physiological integration of both information. The linear dependence of the positional interaction strength on fish inter-distance obviously cannot hold for sparse groups, and will have to be modified by introducing a long-distance saturation when dealing with situations where confinement effects are weaker. Even if we claim that a Voronoi neighborhood was the best choice to account for our data thus extending the relevance of topological interactions, we also checked that our conclusions were robust against this choice, by testing a simple K-Nearest Neighbors network of interactions (which remains topological ). We computed the model predictions with the parameters estimated for groups of N = 2 fish, but considering only the K nearest neighbors for increasing values of K (K = 1 to 7, and 10). The results are reported in Fig. S7; the main impact of a lower level of connectivity is a decrease of polarization, but it does not lead to better predictions at the collective scale. Interestingly, the best predictions were found with a number of nearest neighbors that corresponds to the average number of neighbors belonging to the first shell in a Voronoi neighborhood (, Fig S7–B). This number of influential nearest neighbors is remarkably similar to the one found in starlings and in contrast with recent results found by Herbert-Read et al. in mosquito fish. Further dedicated experiments will be required to discriminate between alternative choices of the relevant neighborhood.

25The speed dependence of the parameters, directly derived from our data, is in contrast with most previous fish school models. It leads to an increase of group polarization with swimming speed, a direct consequence of the predominance of alignment at high speed (see Video S7). In natural conditions, this mechanism could be involved in the transitions from shoaling at low speed often associated with feeding behavior to polarized schooling at high speed associated with searching for food. Such speed change could also be elicited by the detection of a threat and abrupt transitions can occur when fish suddenly increase their speed, for instance generating a flash expansion (see Video S8). The question of whether the propagation of such an excitation wave within large schools can generate an efficient collective evasion call for further experimental tests.

26The reason why our approach was fruitful in spite of the limited amount of data available lies largely in the suitable properties of the behavior of the fish studied: the smooth fluctuations of tangential speed and their de-correlation from angular velocity variations were essential in limiting the number of variables at play but also allowed for a faithful account of single fish behavior by a simple Ornstein-Uhlenbeck process. Clearly it is likely that more complicated solutions will be needed for other species where tangential and angular accelerations are intimately coupled and/or the underlying stochastic process is not as transparent. Nevertheless, we expect that, pending sufficient amounts of data, our approach could be successfully applied to more complex situations occurring in various biological systems at different scales of organization.

27Our experiments were all carried out in full accordance with the ethical guidelines of our research institutions and comply with the European legislation for animal welfare. The welfare of fishes in the tanks was optimized with a continuous seawater flow, a suitable temperature, and oxygen content. The maximum density in the holding tank was lower than. During the experiments, low mortality occurred (five individuals). At the end of the experiment, the fish were released at their capture site.

28The experiments were performed from April to June 2001 at the Sea Turtle Survey and Discovery Centre of Reunion Island. Barred flagtail Kuhlia mugil (Forster) were caught in March 2001 in the coastal area around Reunion Island. 80–100 fishes were conveyed to the marine station and housed in a holding tank of 4 m diameter and 1.2 m depth. Fishes were fed daily ad libitum with a mixture of aquaria flake-food and pieces of fish flesh. Fishes were considered acclimatized when all of them feed on the aquaria flake-food. This weaning period lasted 15 days. Experiments were performed in a circular tank similar to the holding tank. Opaque curtains were placed around and above the tank to obtain diffuse lighting and to reduce external disturbances from the environment. The tank was supplied with a continuous flow of seawater. Since currents may influence fish behavior, the seawater inlet pipe was placed vertically and the water flow was stopped throughout the observation periods. A digital video camera (Sony model CDR-TRV 900E) was fixed at 5 meters above the tank and tilted at to observe the whole tank. The remotely operated video camera was fitted with a polarizing filter and a wide-angle lens. Groups of N = 1 to 30 fish were introduced in the experimental tank and acclimatized to their new environment for a period of 20 min. Their behavior was then recorded at 24 fps for 2 mins. Prior to each trial, the fish were deprived of food for 12 hours to standardize the hunger level and were transferred to the experimental tank. The relative shallowness of the water ensured quasi two-dimensional motion. Five replicates per group size using different individuals were performed. Eighty per cent of the trials were performed in the morning to avoid possible conditions of strong wind that may disturb the fish, and sunshine that may render light inside the tank unsuitable for video recording. A first data processing consisted in sampling 12 images per second out of the 24 images recorded by the video camera. A custom-made tracking software was then used to extract high-quality, smooth trajectories from the video recordings, with crossing ambiguities resolved by eyes (see Video S3, S4). In order to get even higher precision data, the head position and the orientation of each fish in groups of N = 2 were acquired with a manual tracking software (Video S1, S2).

29Model parameters were estimated from each fish time series separately (typical series are shown on Fig. S4). In order to perform the estimation of the parameters,,, and in the stochastic differential equation (1), (3) and (5), we considered its discrete-time version using Ito integration over, assuming is small enough so that is constant:where i = 1,2 and is given by Eq. (3) or (5). Estimates for the parameters were obtained using a standard non-linear least squares procedure (we employed the nls package of the statistical environment R ) either separately for each fish using Eq. (3) or for all fish together using Eq. (5). Residuals given by were checked to be Gaussian-distributed (see Fig. S5) and their variance yielded.

30The model was simulated within a virtual tank, using the estimates of behavioral parameters extracted by statistical analysis from time-series in groups of fish. The fish heading (direction of motion) and position were updated by Euler integration, following:where.

31For each value, numerical simulations were performed over 120 seconds (a time corresponding to the duration of individual experiments with real fish) with a time step. A transient time of was discarded before measuring statistical averages. We computed the mean value and the variance over time of the global polarizationand of the neighbor inter-distance

32This yielded an estimation of the expected measures distribution under model hypothesis and over the typical observation time of experiments. We then computed the mean and confidence interval of such distributions, to obtain the expected mean and variance (with their confidence intervals) of alignment and of neighbor inter-distance. This provided the check of the model against experimental data. The above procedure was repeated varying the mean speed over the range covered by the experimental data, with the results plotted in Fig. 3C. The same procedure was adopted to make predictions for higher group sizes, using the stimulus/response function as determined by equation (5) with interacting neighbors defined by first neighbors in a Voronoi tessellation (For a set of points, Voronoi tessellation divides the space in different cells, each the locus of space closer to its center than to any other points in: at each time step space is divided in Voronoi cells centered around the fish position, with Voronoi neighbors being the fish lying in neighboring cells (Fig. 2B). For each experimental replicate, the same measures were repeated with the parameters extracted from the replicate, and the corresponding initial conditions (Fig. 4B).

33By construction, our method does not “learn the parameters to make the model fit”, contrasting with a more usual procedure which consists in stating an a priori model and searching a best set of free parameters that optimizes its collective patterns towards the observed collective properties (namely, make the model fit at the collective scale). In such cases, it is known that several models can adjust the data at the collective scale (because the search for best match is unconstrained and can be performed for each model, so that the collective level underdetermines the individual level).

34In the present study, once the model has been formulated, that is, once we identified in the experiments with pairs of fish the nature of stimuli (the orientation and relative position of neighboring fish, and how they combine to determine the response of a focal fish), we estimated the values of 5 parameters at the individual scale. So for each fish, we measured its behavioral response (i.e. the change of its turning speed) for each configuration of stimuli encountered in its path.

35Only then, we tested whether these parameters measured at the individual level can explain the observations at the collective scale with no free parameters. For each group independently, we thus checked that the model allows a quantitative matching concurrently at individual and collective scales. This confirmed that our model calibrated with the parameters estimated from the third derivative of the fish position (i.e. the change in the turning speed) was able to reproduce quantitatively the statistics resulting from the time integration of the coupling between fish (polarization, inter-distance). Moreover the same procedure applied separately on each group size revealed, on the one hand, the dependences of the estimated parameters on the swimming speed (using groups of N = 2 fish), and on the other hand, the modulation of interactions' strength with group size (in the largest groups).