在用GEMMA软件进行GWAS分析的时候,需要加入Q矩阵,有一个疑问

在用GEMMA软件进行GWAS分析的时候,需要加入Q矩阵,有一个疑问。为什么这个命令中,cat ../03.STRUCTURE/admixture.4.Q|cut -d " " -f 1-3 >4.Q。用的是admixture.4.Q的结果,是跟群体结构的CV error有关吗?而且后边为什么只取了前3列,生成4.Q。是所有的都只取前3列吗?

请先 登录 后评论

1 个回答

Ti Amo

和CV error有关系。群体分析选择的k是几这里就选择 admixture.k.Q得结果。 选择前三列不是必须的,应当是选择前 k-1 列,也就是不要最后一列。

请先 登录 后评论