Get Omega and Omega^-1 and derivatives
rxSymInv(invobj, theta, pow = 0, dTheta = 0)
rxSymInv object
Thetas to be used for calculation
Power of Omega, either 1 or -1 for inverse. If 0, return an environment with all possible deriatives and Omegas caluclated
the theta number to take the derivative. That is d(Omega)/dTheta. When dTheta = 0, just return the Omega or Omega^-1 matrix. Also can be 0.5. When 0.5 and pow is -1, return -1/2*log(det(omegaInv)).
Matrix based on parameters or environment with all the matrixes calculated in variables omega, omegaInv, dOmega, dOmegaInv.