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

refineqtl can give min(diff(a)) warnings #78

Closed
kbroman opened this issue Sep 18, 2019 · 0 comments
Closed

refineqtl can give min(diff(a)) warnings #78

kbroman opened this issue Sep 18, 2019 · 0 comments
Labels
bug

Comments

@kbroman
Copy link
Owner

@kbroman kbroman commented Sep 18, 2019

It was reported in the R/qtl discussion group that refineqtl() can give a bunch of warnings about min(diff(a)), like this:

In min(diff(a)) : no non-missing arguments to min; returning Inf

I suspect this is in the case of a chromosome with a single marker, in which case diff(a) would be length 0. Here's an example:

data(fake.bc)
fake.bc <- fake.bc[c(2,5),]
fake.bc <- drop.markers(fake.bc, markernames(fake.bc, chr=2)[-1])
fake.bc <- calc.genoprob(fake.bc)
qtl <- makeqtl(fake.bc, chr=c(2,5), pos=c(0, 10), what="prob")
rqtl <- refineqtl(fake.bc, qtl=qtl, method="hk")
@kbroman kbroman added the bug label Sep 18, 2019
kbroman added a commit that referenced this issue Sep 18, 2019
 - warning in `refineqtl()` that's like `min(diff(a))`
 - I think this is due to case of a single marker on a chromosome
@kbroman kbroman closed this Oct 4, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Linked pull requests

Successfully merging a pull request may close this issue.

None yet
1 participant
You can’t perform that action at this time.