标签:GWAS GEMMA tencent 软件 txt com gemma
首先基因型文件满足plink格式,生成包含有ped.map及二进制的bed.bim和fam即可;
./gemma-0.98.1-linux-static -bfile g -gk 2 -p p.txt -o gemma 亲缘关系矩阵;
./gemma-0.98.1-linux-static -bfile g -k ./output/gemma.sXX.txt -lmm 1 -p p.txt -c c.txt -o gemma_out 关联分析;
主要参考邓飞笔记 | GWAS 操作流程5-1:根红苗正的GWAS分析软件:GEMMA - 云+社区 - 腾讯云 (tencent.com)笔记 | GWAS 操作流程5-2:利用GEMMA软件进行LMM+PCA+协变量 - 云+社区 - 腾讯云 (tencent.com)两篇文章
中间出错有,采用协变量成分时如下:
一直报错:GSL ERROR: matrix is singular in lu.c at line 266 errno 1,
出不来结果,中间无数次修改,还看到有博主说删除第一列截距就能出来结果,可我的还是不行。
再最后一遍尝试把中间的性别和世代都删除了忽然就出来结果了。
当然也不能确定结果就正确,也不太明白协变量应该留那些,目前就是留了
这样协变量就能做出来了。
如果有用也希望能帮助到你,有错误也希望能告诉我下。一起入门学习,借鉴经验。
参考文献:
使用GEMMA进行GWAS分析 - 简书 (jianshu.com)
GWAS | 使用GEMMA进行全基因组关联分析 - 云+社区 - 腾讯云 (tencent.com)
标签:GWAS,GEMMA,tencent,软件,txt,com,gemma 来源: https://blog.csdn.net/rongrongnancy/article/details/123218778
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。