Jacobian of moment conditions for linear GMM estimation
[M,e] = gmmLinJ(b,gmmopt,Y,X,Z)
b = parameter vector fed to function
gmmopt = gmm options tlist
Y = "dependent" variables
X = "independent" variables
Z = instruments (can be same as X)
M = Jacobian (rows(m) x k)