首页 > 解决方案 > 如何将 STRUCTURE 结果导出到 StrataG 中的 STRUCTURE 输出文件?

问题描述

我正在使用strataG为我的微卫星数据运行 STRUCTURE。我的 Mac 正在运行 Catalina,但我不再能够运行 32 位软件 CLUMPP。为了解决这个问题,我希望能够将我的 STRUCTURE 结果导出到一个 STRUCTURE 输出文件中,以尝试在Pophelper包中运行,是否有在strataG中执行此操作的函数?

如果需要,我使用 RStudio 版本 1.4.1103 和 strataG 版本 2.4.905,

mydata <- structure(list(ID = c("Test1", "Test10", "Test11", "Test12"), pop = c("Test", 
"Test", "Test", "Test"), SalD39.1 = c("242", "244", "242", 
"242"), SalD39.2 = c("252", "244", "242", "242"), SalF56.1 = c("204", 
"204", "204", "204"), SalF56.2 = c("204", "204", "204", "204"
), SalJ81.1 = c("118", "118", "118", "118"), SalJ81.2 = c("118", 
"118", "118", "118"), SalP61.1 = c("128", "128", "128", "128"
), SalP61.2 = c("130", "128", "130", "130"), Ssa406.1 = c("411", 
"411", "411", "411"), Ssa406.2 = c("415", "411", "411", "411"
), Ssa85.1 = c("164", "168", "166", "166"), Ssa85.2 = c("164", 
"168", "166", "166"), SsaD48.1 = c("230", "238", "262", "234"
), SsaD48.2 = c("238", "238", "294", "250")), class = c("data.table", 
"data.frame"), row.names = c(NA, -4L))

library(strataG)
mydata_gtypes <- df2gtypes(mydata, ploidy = 2, id.col = 1, strata.col = 2, loc.col = 3)

安装外部程序 STRUCTURE 并在strataG中运行以下包装函数:

sr_noAdmix <- structureRun(mydata_gtypes, k.range = 2:9, burnin = 100, numreps = 100, num.k.rep = 10, noadmix = T) 

现在可以将 sr_noAdmix 导出到 STRUCTURE 格式的输出文件吗?

标签: rexportgenetics

解决方案


推荐阅读