主要是用同一批数据, 测试不同软件和方法对结果的影响, 不同方法有:
BLUPF90构建G矩阵, 计算PCA
BLUPF90构建H矩阵, 计算PCA
PLINK构建G矩阵, 计算PCA
GCTA构建G矩阵...GCTA画的PCA有出入, 怀疑是默认参数的问题....模拟一套数据, 5个世代, 最后三代有基因型数据, 每个世代400个个体, SNP为50K.
1. blupf90构建G矩阵的PCA
blupf90如果想要进行GBLUP分析, 不写系谱信息即可, 示例...2. blupf90构建H矩阵的PCA
需要定义系谱和基因型, 示例par文件:
DATAFILE
dat_f90.txt
TRAITS
10 # This is column 10 (...结论
blupf90的G矩阵, H矩阵, plink的PCA结果一致.
GCTA构建的PCA结果不太一致, 怀疑是参数默认的有问题, 回头查看一下.