1. 写GWAS专辑的初衷当时没有想过写GWAS专辑,主要是因为有一些东西我需要弄明白,就写了一些博客,然后就同步更新到公众号上,然后就汇总到了一个专辑。现在写了17篇,包括数据的质控,使用plink软件进行GLM模型分析,使用GEMMA软件进行LMM模型分析,添加因子协变量和数值协变量的方法等等。 2. 只有对自己有帮助的,才是对别人有帮助的我比较大的收获是用R语言进行GLM模型的GWAS分析,包括添加因子协变量,添加数值协变量,添加PCA协变量,从编程到结果理解,发现GWAS分析“知其然也知其所以然”。我不仅仅是调包侠了,还可以编写一些东西,学习的快乐就是这么朴实无华且枯燥。 3. 真正的快乐前一段时间和朋友聊天,介绍我最近取得了一些项目的进展,那种成就感体验非常棒,和山珍美味锦衣美食是不一样的(好像我经常体验山珍美味一样,新疆之行九天胖三斤说明了吃不胖主要是因为穷这个朴实无华的真理(新疆归来之不敢发朋友圈))。马斯洛的人本主义(对的,我也对心理学有研究,什么佛洛依德梦的解析,荣格集体无意识都有涉猎),对人的需求分为了几类,有一种是自我实现的需求,而这种突破困难,历经千辛万苦找到答案后的快乐,那种狂喜,应该可以归为自我实现的快乐,是比较高级的快乐,马斯洛称之为:“高峰体验”(禁止开车警告)。 4. 所以我到底想说什么?我想说的是,我需要设立目标,不是做简单的事情,而是做困难的事情。简单的事情是什么呢,搞一个PPT分享一下,写一个基础教程分享一下,凉白开一样,没有味道。搞就搞一个大的,辛苦学习理解之后豁然开朗的那种(比如--一文讲清楚动物模型中的母体效应)。 「所以呢?」 「继续学习GWAS」 既然学习GWAS,手动用R编写了GLM模型,可以用R编写LMM模型么?二分类性状GWAS分析,plink可以做logistic,如果考虑亲缘关系矩阵的GLMM(广义线性混合模型)如何处理呢?阈值性状(多于二分类的分类性状)的GWAS如何分析呢? 5. GWAS的哔哩哔哩教程在很久以前,我录制了两个GWAS的教程,那时候我技术还比较青涩,现在已经变成老司机了(禁止开车警告),包括代码的显著提示,结果的注意事项,都可以通过视频编辑的形式做成大家喜闻乐见的形式。(试想一下,用鬼畜做GWAS教程是个什么感觉???)。而我准备尝试一下。 之前的哔哩哔哩视频: 6. 所以呢,计划是什么?继续学习GWAS,我要把相关东西理解清楚,然后把我理解的过程总结一下,毕竟听“没有结过婚的人,大谈婚姻对人的伤害”,是不地道。我相信我的问题不是我一个人的问题,而是后来者都会遇到的问题,写出来,分享一下,才是价值所在,正如我之前所言: ❝ 7. ^*^广告时间哔哩哔哩毕竟是个学习的地方,而微信的“视频号”则是另一个风格,时长一分钟,非常适合用精简的内容讲一个合适的段子,毕竟我这个“育种数据分析之放飞自我”的博主段子讲得好真不是盖的,不信扫码看一下我的主页:
|
|