addtoqtl()

23 views
Skip to first unread message

nata...@gmail.com

unread,
Aug 8, 2021, 7:58:10 AM8/8/21
to R/qtl discussion
Hi,
I encountered an error with addtoqtl() after identifing an additional QTL with addqtl()

>rqtl<-addtoqtl(mycross, rqtl, s.aq[,1], s.aq[,2], drop.lod.profile=F)
> rqtl
QTL object containing genotype probabilities. 

name  chr     pos n.gen
Q1 LG6@27.0  LG6 27.0000     3
Q2 LG19@8.4 LG19  8.3968     3
Q3 LG6@10.1  LG6 10.1451     3
>     (bay <-bayesint(rqtl, qtl.index = 3))
Error in bayesint(rqtl, qtl.index = 3) : 
  qtl.index should be between 1 and 2

 What should I do?
Thanks!

Karl Broman

unread,
Aug 8, 2021, 1:29:29 PM8/8/21
to R/qtl discussion
You need to run refineqtl again before bayesint.

karl

nata...@gmail.com

unread,
Aug 8, 2021, 1:56:11 PM8/8/21
to R/qtl discussion
Ok thank you

ב-יום ראשון, 8 באוגוסט 2021 בשעה 20:29:29 UTC+3, ‪Karl Broman‬‏ כתב/ה:
Reply all
Reply to author
Forward
0 new messages