Can't you just check for NA coefficients? > y <- rnorm(10) ; x <- rep(0,10) > coef(lm(y~x)) (Intercept) x -0.0962404 NA
so > any(is.na(coef(lm(y~x)))) [1] TRUE I have a vague recollection that at some point there might have been dragons lurking in there (? - NA coefs silently removed), but I can't see a problem with it presently. -pd > On 5 Apr 2019, at 12:14 , Witold E Wolski <wewol...@gmail.com> wrote: > > lme4 has a function isSingular to check if the fitted model is Singular, > > Although lm has the parameter singular.ok = TRUE by defualt, I could > not find a function to check if the fitted model is singular. > > What would be the correct way to implement such a function for and lm object? > Check if df.residuals == 0 > > Thanks > Witek > > > > > > -- > Witold Eryk Wolski > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd....@cbs.dk Priv: pda...@gmail.com ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.