Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Warning INFO #14

Open
ddrinan opened this issue Jul 7, 2017 · 2 comments
Open

Warning INFO #14

ddrinan opened this issue Jul 7, 2017 · 2 comments

Comments

@ddrinan
Copy link

ddrinan commented Jul 7, 2017

Matthew,

Hi, I asked previously (in another thread) about warnings, and wanted to double back on it one more time. When running SCAT on a large SNP data set (+6000 loci), I often get the message Warning INFO=4. I've tried to track this down, and best I can tell, it is related to the dpotrf function (computes the Cholesky factorization of a real symmetric positive definite matrix). My understanding is that if INFO=4, then "the leading minor of order 4 is not positive definite, and the factorization could not be completed." (information from here). That being said, I don't fully understand what is being computed and how it would affect the results. Do you have any insight on whether this is a fatal error such that the results can't be trusted? Also, do you have any idea what would cause this error, or more importantly, how do I get around it?

Thanks,
Dan Drinan

@stephens999
Copy link
Owner

it could be caused if you have some locations very close together?
Possible solution would be to merge locations that are very close.

Difficult to predict what the impact on results would be....my best guess
would be either very little effect or huge effect... so if the results "generally make sense" I would
speculate they are OK.

@ddrinan
Copy link
Author

ddrinan commented Jul 7, 2017

Haha, I understand. We do have some sites that are relatively close. So, I'll play around with combining them. Generally, the results look good. Thanks so much for your help.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants