我有Dartseq的SNP数据集(148个个体和12,123个基因座)。
我正在dartR中对其进行分析。
我已经在dartR中计算了Ho(gl.Ho(gl))和Hs(gl.Hs(gl))。
我想绘制Ho v Hs。
我该怎么办?
我认为这应该很容易,但是我没有取得任何进展。
我知道以下内容:
head(Ho)
39086092-43-A/T 39030115-36-G/T 39036131-46-C/T 39031387-28-T/A 39043148-8-C/T
0.006944444 0.000000000 0.014184397 0.076388889 0.007246377
39052543-5-C/T
0.000000000
str(Ho) # # Shows the structure of the data
Named num [1:12119] 0.00694 0 0.01418 0.07639 0.00725 ...
- attr(*, "names")= chr [1:12119] "39086092-43-A/T" "39030115-36-G/T" "39036131-46-C/T" "39031387-28-T/A" ...
class(Ho) # numeric
[1]“数字”
head(Hs)
39086092-43-A/T 39030115-36-G/T 39036131-46-C/T 39031387-28-T/A 39043148-8-C/T
0.006920332 0.013888210 0.014083799 0.192105517 0.049438143
39052543-5-C/T
0.027391975
PS。在我找到有关Markdown帮助的帖子后,今天上午进行了编辑。现在至少更容易阅读。