blup_indexes function

Stability indexes based on a mixed-effect model