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

HGTfinder errror #2

Open
gsun2unl opened this issue May 3, 2019 · 1 comment
Open

HGTfinder errror #2

gsun2unl opened this issue May 3, 2019 · 1 comment

Comments

@gsun2unl
Copy link

gsun2unl commented May 3, 2019

Hi, I am currently trying to use HGT finder following the readme file, however, it gave me this error:
Invalid R value for subject: 5*******

And eventually, the program was halted with this error:

Error in smooth.spline(lambda, pi0, df = smooth.df) : missing or infinite values in inputs are not allowed Calls: qvalue -> pi0est -> smooth.spline

Can you give me some clue on fixing them?
Thank you!

@Tinyman392
Copy link
Collaborator

The error stems from the Invalid R value for your subject. More specifically, this message is triggered when the bit score of the query to itself is lower than the bit score between the query and the subject (which theoretically shouldn't happen, but can in practice). We've seen this occur with sequences that contain N (wildcard) values.

Off the top of my head I do not remember if the message is more likely to be triggered if the query has many N characters or if it's the subject. You'd need to check both and possibly exclude that (those) sequence(s) from the analysis.

The Error in the qvalue (smooth.spline) call is triggered most likely due to having an invalid P-value due to the invalid R value.

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