pleio.glm.fit function

Compute Intermediate Statisics for Testing Pleiotropy