Performance of alternative spatial models in empirical Douglas-fir and simulated datasets

Eduardo Pablo Cappa (1) , Facundo Muñoz (2) , Leopoldo Sanchez (2)

1 CONICET - Consejo Nacional de Investigaciones Científicas y Técnicas [Buenos Aires]
2 BioForA - Biologie intégrée pour la valorisation de la diversité des arbres et de la forêt

AbstractKey messageBased on an empirical dataset originating from the French Douglas-fir breeding program, we showed that the bidimensional autoregressive and the two-dimensional P-spline regression spatial models clearly outperformed the classical block model, in terms of both goodness of fit and predicting ability. In contrast, the differences between both spatial models were relatively small. In general, results from simulated data were well in agreement with those from empirical data.ContextEnvironmental (and/or non-environmental) global and local spatial trends can lead to biases in the estimation of genetic parameters and the prediction of individual additive genetic effects.AimsThe goal of the present research is to compare the performances of the classical a priori block design (block) and two different a posteriori spatial models: a bidimensional first-order autoregressive process (AR) and a bidimensional P-spline regression (splines).MethodsData from eight trials of the French Douglas-fir breeding program were analyzed using the block, AR, and splines models, and data from 8640 simulated datasets corresponding to 180 different scenarios were also analyzed using the two a posteriori spatial models. For each real and simulated dataset, we compared the fitted models using several performance metrics.ResultsThere is a substantial gain in accuracy and precision in switching from classical a priori blocks design to any of the two alternative a posteriori spatial methodologies. However, the differences between AR and splines were relatively small. Simulations, covering a larger though oversimplified hypothetical setting, seemed to support previous empirical findings. Both spatial approaches yielded unbiased estimations of the variance components when they match with the respective simulation data.ConclusionIn practice, both spatial models (i.e., AR and splines) suitably capture spatial variation. It is usually safe to use any of them. The final choice could be driven solely by operational reasons.

Mots clés: Forest genetics trials; Autoregressive residual; Two-dimensional P-splines; Global and local spatial trends

Domaines: Sciences du Vivant [q-bio]

Eduardo Pablo Cappa, Facundo Muñoz, Leopoldo Sanchez. Performance of alternative spatial models in empirical Douglas-fir and simulated datasets. Annals of Forest Science, 2019, 76 (2), pp.53. ⟨10.1007/s13595-019-0836-9⟩. ⟨hal-02610771⟩