linux - 如何从 PLINK 2.0 上的剂量文件创建 MAP/PED 文件

我是 PLINK 和 GWAS 研究的新手。我不确定如何编写脚本将我的 INFO 文件或剂量数据转换为 MAP 文件，以进行基因型和 1 表型之间的关联研究。所以，第一件事：我有一个 INFO 文件，其中包含剂量信息和 SNP、ID、位置、等位基因 1 和 2 的信息以及插补信息（见下文）

rs148087467:60523:T:G 60523 TG 0.001 0.161 0.998 0 -1 -1 -1 --- 10:60684:A:C 60684 AC 0.029 0.713 0.978 0 -1 -1 -1

我想创建一个脚本，但不知道如何将其转换为地图文件。我了解 MAP 文件包含有关染色体代码、变异标识符 (rs)、位置/遗传距离 (morgans) 和碱基对坐标 (bp 单位) 的信息。

我假设也会由此生成一个 ped 文件？

有人可以帮我如何创建脚本！我对这一切都很陌生，所以我道歉。

谢谢！

标签： linuxunix