make a genitic map from zero

46 views
Skip to first unread message

Luigi Faino

unread,
Dec 23, 2024, 7:11:26 AM12/23/24
to R/qtl discussion
Dear,

I'm new to R/qtl and i need some help. I need to build a new genetic map from scratch meaning that i have an F2 population and marker genotypes.
I see that R/qtl can be used but i can not really understand how.  From what i see, I need to start from a geno file and a pheno file and in the geno file i need to have chromosome and position of each marker. However, I do not have this information. 

which king of file do i need to have to make a linkage map starting only from the marker information?

thanks
L

Karl Broman

unread,
Dec 23, 2024, 7:12:33 AM12/23/24
to rqtl...@googlegroups.com
Assign arbitrary values. For example, put "un" for the chromosome for all markers, and give them positions 1, 2, 3, ...

karl

--
You received this message because you are subscribed to the Google Groups "R/qtl discussion" group.
To unsubscribe from this group and stop receiving emails from it, send an email to rqtl-disc+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/rqtl-disc/4f299064-7ac9-4173-803e-c8d66e6a367cn%40googlegroups.com.

Luigi Faino

unread,
Dec 23, 2024, 3:42:17 PM12/23/24
to R/qtl discussion
ok i did it and it works.

however i have another problem. I have many linkage groups and i would like to keep only the largest 11. do you have any idea on how to do it?

L

Karl Broman

unread,
Dec 23, 2024, 3:44:49 PM12/23/24
to R/qtl discussion
It’s hard to say without knowing more about the situation.
There is a tutorial on genetic map construction with R/qtl: https://rqtl.org/tutorials/geneticmaps.pdf

karl

Luigi Faino

unread,
Dec 23, 2024, 4:18:28 PM12/23/24
to R/qtl discussion
this is an example of my data:

n.mar length ave.spacing max.spacing 1 355 3540 10 10 2 351 3500 10 10 3 344 3430 10 10 4 340 3390 10 10 5 261 2600 10 10 6 243 2420 10 10 7 207 2060 10 10 8 174 1730 10 10 9 173 1720 10 10 10 152 1510 10 10 11 96 950 10 10 12 22 210 10 10 13 19 180 10 10 14 11 100 10 10 15 10 90 10 10 16 10 90 10 10 17 10 90 10 10 18 9 80 10 10 19 8 70 10 10 20 8 70 10 10 21 8 70 10 10 22 7 60 10 10 23 7 60 10 10 24 7 60 10 10 25 6 50 10 10 26 6 50 10 10


I want to keep only the 11 top linkage groups

L

Karl Broman

unread,
Dec 27, 2024, 11:03:28 AM12/27/24
to R/qtl discussion
You can use drop.markers to remove the other markers from the dataset, or movemarker to move them to some other “unlinked” linkage group.

karl

Luigi Faino

unread,
Dec 27, 2024, 5:30:57 PM12/27/24
to R/qtl discussion
Dear

 what do you think about this map? can you advice any adjustments or modification?

LRplotMAP.jpeg

Samuel Wang

unread,
Dec 27, 2024, 9:14:25 PM12/27/24
to rqtl...@googlegroups.com
Dear L:

What methods were used to generate the genotype data?

If the genotype data were generated from the next-generation sequencing method, I would suggest removing the genotype as missing with a genotype quality less than 20, and the marker missing rate should less than 0.25.

The genotype would become easier to build the reasonable linkage map.

Samuel

On Dec 28, 2024, at 06:30, Luigi Faino <luigi...@gmail.com> wrote:

Dear

 what do you think about this map? can you advice any adjustments or modification?

L<RplotMAP.jpeg>

Reply all
Reply to author
Forward
0 new messages