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!