Dear all, I am running lmrob() under the robustbase package, and I have a few error/warning messages I would greatly appreciate help with - I am unfamiliar with robust linear models, so my apologies in advance if the answers here are obvious. In addition, any discussion will help me get a handle on what we're doing, which always helps. (1) I get the repeated error seemingly on random species when running a loop of lmrob() over all our ecological time series (different ones almost every time): Error in eigen(ret, symmetric = TRUE) : infinite or missing values in 'x' I have checked my data many times, and if I just repeatedly run lmrob() it eventually works. Any idea what this error is in reference to? Could it be because I am using a previous version of R (3.4.1)? (2) I have several two error messages that come up with some time series; I have some idea their meaning, but again, as a newbie with RLMS, I am wondering if those time series that garner these warnings should be excluded from our meta-analysis. Also, any additional insight on the process or something I might be missing greatly welcomed! Warning in lmrob.S(x, y, control = control) : find_scale() did not converge in 'maxit.scale' (= 200) iterations with tol=1e-10, last rel.diff=0 Warning in lmrob.fit(x, y, control, init = init) : M-step did NOT converge. Returning unconverged SM-estimate Warning in lmrob.S(x, y, control = control) : S-estimated scale == 0: Probably exact fit; check your data Thank you very much in advance, Emily
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Emily S. Klein, Senior Postdoctoral Associate (she / her / hers) The Frederick S. Pardee Center for the Study of the Longer-Range Future | Boston University *Co-Chair*, ICES Working Group on the History of Fish & Fisheries (WGHIST) esklein04 at gmail.com http://www.bu.edu/pardee/ http://www.ices.dk/community/groups/Pages/WGHIST.aspx [[alternative HTML version deleted]]