Entering edit mode
9.8 years ago
JorgeGonzalez
▴
20
Hi
I'm working in the search loci under selection from genetic data (SNPs) by bayesian approach (BayesScan 2.1) and I got the following plot, the result is very good but I'd like visualize better way, pots instead of code ( the code is the name of given locus).
I used the plot_R.r
script and the following functions
> plot_bayescan("tesis_fst.txt",0,FDR=0.05)
$outliers
[1] 149 195 251 273 434 494 500 576 592 623 635 655 795 852 853 897
[17] 903 959 1022 1088 1178 1181 1355 1414 1418 1523 1801 1829 1856 1878 1915 2063
[33] 2131 2205 2353 2362 2423 2427 2545 2637 2639 2657 2706 2717 2732 2757 2806 2835
[49] 2850 2911 2966 2993 3019 3030 3120 3190 3221 3281 3296 3345 3408 3416 3506 3523
[65] 3529 3530 3612 3643 3705 3760 3882 3940 3984 4023 4156 4201 4435 4440 4461 4578
[81] 4677 4775 4869 4908 5002 5008 5031 5088 5113 5125 5364 5472 5515 5555 5581 5584
[97] 5722 5746 5760 5782 5807
$nb_outliers
[1] 101
Can someone please help me?
It looks like the plot is not displayed here.