Re: [R] R SIG mailing lists

2019-04-05 Thread H
On 04/05/2019 12:39 AM, Bert Gunter wrote: > https://www.r-project.org/mail.html > > Found immediately by a web search on "R Mailing lists" . Please make a > minimal effort yourself before posting, or let us know if you have already > done so but came up empty. > > You might try here for genomics

Re: [R] isSingular for lm?

2019-04-05 Thread Fox, John
Hi Peter, > On Apr 5, 2019, at 7:22 AM, peter dalgaard wrote: > > 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 recolle

Re: [R] R SIG mailing lists

2019-04-05 Thread Powers, Randall - BLS via R-help
I think you mean "Healthcare" not "Healthrcare." Please take minimal time to check your spelling before replying to a post. :) -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Bert Gunter Sent: Friday, April 5, 2019 12:39 AM To: H Cc: R Mailing List Sub

Re: [R] R SIG mailing lists

2019-04-05 Thread Marc Schwartz via R-help
Hi, To supplement Bert's reply, there is not a specific SIG list for R in healthcare, however, there are a number of us, both in the pre-clinical and clinical realm, that participate here on R-Help. As Bert noted, for genomic applications, the Bioconductor folks have their own lists. If you h

Re: [R] isSingular for lm?

2019-04-05 Thread peter dalgaard
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 sil

[R] isSingular for lm?

2019-04-05 Thread Witold E Wolski
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.res