r - 如何将 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 格式的输出文件吗?
解决方案
推荐阅读
- reactjs - 我无法让我的钩子正常工作
- mysql - 我需要通过多个分隔符将 mysql 字段(字符串)拆分为 json 对象
- python - 使用 python,我需要每天使用 MongoDB api 将文档从 Azure CosmosDB 传输到 Postrgres,可能使用 Azure 函数
- javascript - 跟随光标位置将图像添加到 DOM
- python - 在通过 gerrit-python-api 包采摘樱桃时确定其是合并冲突还是相同的树?
- excel - Excel Power Query 合并和转换列
- python - 根据另一个数组的元素从 numpy 数组中删除元素
- javascript - 指定replicationRegions时的Javascript NestedStack CREATE_FAILED
- react-native - UI Kitten 中的 Top Navigation 组件如何在 React Native Stack Navigator 中用作自定义标题
- python - 字典更新序列元素#0的长度为0;2 是必需的