Errors are caused for certain models when running `fit_mgx_model()`

statdivlab / enviromtx

A package to analyze environmental metatranscriptomics data

Other

1 stars 1 forks source link

Errors are caused for certain models when running `fit_mgx_model()` #9

Open svteichman opened 5 months ago

svteichman commented 5 months ago

Specifically, there are examples (see some cases here) in which the call to raoBust::gee_test calls geeasy::geelm() which creates the error "Error in .solve.checkCond(a, tol) : 'a' is computationally singular, rcond(a)=0".

@svteichman to try to figure out why these errors are happening.

svteichman commented 5 months ago

Notes on this issue:

when the model does converge for the examples in the linked issue, the beta values are not unusually large
when the model fails to converge, it runs into case in geeasy:::update_beta where the magnitude of the beta values gets very large, and eventually the hessian has a very large magnitude and cannot be inverted, leading to the error message Error in .solve.checkCond(a, tol): 'a' is computationally singular, ...