VP.gllvm function

Calculate variance partitioning