tassel 跑PCA的时候报错了:
检查了一下发现是SNP超过了允许的数量范围,总的SNP:
tassel PCA允许的SNP数量:
作者推荐使用MDS替代PCA分析,MDS和PCA的作用比较类似,只不过MDS是基于最小线性距离进行聚类,而PCA是基于最大线性相关性进行聚类:
MDS分析需要距离矩阵作为输入,先使用tassel生成距离矩阵:
run_pipeline.pl -Xmx200G -fork1 -h hapmap.ordered.hmp.txt -DistanceMatrixPlugin -endPlugin -export DistanceMatrix -runfork1
然后运行MDS:
run_pipeline.pl -Xmx200G -fork1 -importGuess DistanceMatrix.txt -MultiDimensionalScalingPlugin -endPlugin -export MDS -runfork1
报错解决:
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!