A model that allows for two hidden variables and spatial autocorrelation is proposed.
Accuracy in predicting species biomass varies within the spatially resolved areas.
The specific hidden variable models spatial unmeasured effect.
We found temporally and spatially differentiated functional networks.
Combining structure learning from data and experts' knowledge in the model architecture is optimum.