stephens999 / ashr

An R package for adaptive shrinkage
GNU General Public License v3.0
79 stars 35 forks source link

Many warnings of singular system in R 4.2 #131

Closed frederikziebell closed 2 years ago

frederikziebell commented 2 years ago

As of R4.2, the following code

beta = c(rep(0,100),rnorm(100))
sebetahat = abs(rnorm(200,0,1))
betahat = rnorm(200,beta,sebetahat)
beta.ash = ash(betahat, sebetahat)

results in many warnings like these: warning: solve(): system is singular (rcond: 2.12669e-17); attempting approx solution.

pcarbo commented 2 years ago

@frederikziebell Thanks for reporting this. This is actually a problem with mixsqp. Can you please try installing the new mixsqp and let me know if the warnings go away?

frederikziebell commented 2 years ago

This is working without warnings, thanks!