Get_Gy1 function

Compute the feedback matrix from a GVAR model with global factors