small cleaning...

1.  adding new code for NA check for missing data
2.  making fitting more robust by using REML and fixing SigmaR when
nonconverging