cox reid likelihood in R

Dear all

I am trying to compute the Cox reid adjusted likelihood ratio test for mixed models in r. The adjustment needs the information matrix. Is there a package/function in R that runs the test or is there a way to extract the information matrix from an lme object?

Thanks
Akshita