Question: How to get the beginning and end-positions of a QTL-interval after doing QTL-analysis with R/QTL?
gravatar for nx8
4.5 years ago by
nx810 wrote:


I've done QTL-Analysis using R/QTL. I would now like to find all the genes lying in the QTL-regions. Is there a function that will give me the positions of the flanking markers of the QTL-interval in my genetic map? So that I can say the QTL interval is from position x to position y and check the genes within this region.

Many thanks!


ADD COMMENTlink modified 2.7 years ago by Matthias20 • written 4.5 years ago by nx810
gravatar for Matthias
2.7 years ago by
Matthias20 wrote:

Supposing you run scan1 (name of return object is data.scan1) on your data and want to get the 95% confidence interval (upper and lower marker) for the peak on chromosome 1, type

CI = bayesint(data.scan1, chr=1, prob=0.95)

lower_marker = find.marker(data, chr=CI[1,1], pos=CI[1,2])
lower_lod = CI[1,3]
print(paste("Lower marker:", lower_marker, "LOD:", lower_lod, sep=" ")) 

upper_marker = find.marker(data, chr=CI[3,1], pos=CI[3,2]) 
upper_lod = CI[3,3]
print(paste("Upper marker:", upper_marker, "LOD:", upper_lod, sep=" "))
ADD COMMENTlink written 2.7 years ago by Matthias20

Hi. Question: would this work also for composite-interval-mapping (CIM) ? Or how it can be calculated for CIM? Thanks

ADD REPLYlink written 16 months ago by boludopublico0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1371 users visited in the last hour