dadi

这是我用的命令 cd $dir/sfs

awk '{ print $2}' $data/africa.group.txt | sort | uniq -c

# 无外群，生成foldfs

easySFS.py -i$data/Africadadi.vcf.gz -p $data/africa.group.txt -a --proj=74,56,60 -o./

## 生成了3个群体的、以及三个群体两两之间的sfs文件

cd dadi/

mkdir 2D

mkdir 3D

python3 $script/sfs_draw.py -i ./NorthAfrica-Yourawildgoat.sfs -o 2D/2D

绘制出的图是空的

默认排序时间排序

2 个回答

Ti Amo 13小时前

检查输入文件./NorthAfrica-Yourawildgoat.sfs是不是有一行全是0。easySFS.py里面投影值设置根据自己样本数目进行更改

随风 9小时前

这个我明白了，样本数修改后，运行结果如下，这个正常吗？