主要是用同一批数据, 测试不同软件和方法对结果的影响, 不同方法有:
BLUPF90构建G矩阵, 计算PCA
BLUPF90构建H矩阵, 计算PCA
PLINK构建G矩阵, 计算PCA
GCTA构建G矩阵...模拟一套数据, 5个世代, 最后三代有基因型数据, 每个世代400个个体, SNP为50K.
1. blupf90构建G矩阵的PCA
blupf90如果想要进行GBLUP分析, 不写系谱信息即可, 示例...sol se
#OPTION use_yams
#OPTION missing -999
OPTION plotpca
运行preGSf90后, 会生成pc1vspc2文件, 里面包括PC1和PC2两列,...level = 0.95, size = 1)
p <- p + geom_hline(yintercept = 0)
p <- p + geom_vline(xintercept = 0)
p <...sol se
#OPTION use_yams
#OPTION missing -999
OPTION plotpca
运行preGSf90后, 会生成pc1vspc2文件, 里面包括PC1和PC2两列,