纸上得来终觉浅,绝知此事要躬行
学习要明确输入输出,多举例子,多找到应用场景,多进行检索练习——费曼学习法
之前总感觉WGCNA已经学习的差不多了,已经能够实际使用了,结果时间一长,就会把这个技能给生疏了,所以,记录下学习的过程,有助于自己的认识深入。


WGCNA第一次实际上是把众多的基因进行了一个归类,将具有相似的生物学功能的基因放置在一起,这个过程可以理解为相关性,也可以是一种聚类的思维。许多的WGCNA教程都只是把代码跑了一遍,而这幅图非常好,是果子学生信课程图片,可视化展示了这个过程,摘录到此处。
将相似的基因归类为一类,那么也就是一个模块,可以继续分析。紧接着面对的问题就是,后续分析怎样利用到模块,而不是单纯继续对基因进行分析。模块思维就在这里使用了,将这些类似的基因,统一成一个值,就叫做epigene,我的理解,这个指标更像是一个聚类之后的评分。

后续的分析,就以新的epigene为基础,进行相应的分析。
当得到了上面的epigene的矩阵文件后,紧接着就是怎么和我们的性状联系,下面的图也是从果子学生信中摘录。

这就是熟悉的两个矩阵求相关性的问题,在之前的批量相关性分析中,其实也涉及到这个操作,批量处理——基因与免疫细胞的相关性。
明白原理后,每一个步骤都会变得愉悦起来
通过上面两个矩阵之间的相关性分析,就可以得到模块epigene和性状之间的关系

这幅图是运行WGCNA作者给出的操作手册中的示例文件,得到的相关性图。颜色只是作者为了能够更好的可视化,其中这个给某一部分特征性上色的技巧可以迁移。
而这幅图就是我们最需要的一幅图,我们能够从这幅图中找到与性状有关的模块,进而能够将这些模块基因提取出来,再进一步分析,一定程度上,这已经缩小了一部分研究范围。虽然离聚焦到一个基因上而言还有一段距离,但是能够有这样的范围缩小,也会极大地有利于研究。
引用一下
当你想用纯生信发文章的时候,对一个技术的要求是苛刻的,恨不能用某个技术直接代替Western blot,直接给我定位到最重要的单个基因上去。
而,对于科研人员来说,一个技术能准确的把范围缩小到一群基因,已经很满意了。接下来通过多个技术和实验来确定最重要的基因就行了。
WGCNA: an R package for weighted correlation network analysis | BMC Bioinformatics
Tutorials for WGCNA R package (ucla.edu)
我正在尝试从Postgresql表(table1)中获取数据,该表由另一个相关表(property)的字段(table2)过滤。在纯SQL中,我会这样编写查询:SELECT*FROMtable1JOINtable2USING(table2_id)WHEREtable2.propertyLIKE'query%'这工作正常:scope:my_scope,->(query){includes(:table2).where("table2.property":query)}但我真正需要的是使用LIKE运算符进行过滤,而不是严格相等。然而,这是行不通的:scope:my_scope,->(que
目录前言滤波电路科普主要分类实际情况单位的概念常用评价参数函数型滤波器简单分析滤波电路构成低通滤波器RC低通滤波器RL低通滤波器高通滤波器RC高通滤波器RL高通滤波器部分摘自《LC滤波器设计与制作》,侵权删。前言最近需要学习放大电路和滤波电路,但是由于只在之前做音乐频谱分析仪的时候简单了解过一点点运放,所以也是相当从零开始学习了。滤波电路科普主要分类滤波器:主要是从不同频率的成分中提取出特定频率的信号。有源滤波器:由RC元件与运算放大器组成的滤波器。可滤除某一次或多次谐波,最普通易于采用的无源滤波器结构是将电感与电容串联,可对主要次谐波(3、5、7)构成低阻抗旁路。无源滤波器:无源滤波器,又称
最近在学习CAN,记录一下,也供大家参考交流。推荐几个我觉得很好的CAN学习,本文也是在看了他们的好文之后做的笔记首先是瑞萨的CAN入门,真的通透;秀!靠这篇我竟然2天理解了CAN协议!实战STM32F4CAN!原文链接:https://blog.csdn.net/XiaoXiaoPengBo/article/details/116206252CAN详解(小白教程)原文链接:https://blog.csdn.net/xwwwj/article/details/105372234一篇易懂的CAN通讯协议指南1一篇易懂的CAN通讯协议指南1-知乎(zhihu.com)视频推荐CAN总线个人知识总
深度学习部署:Windows安装pycocotools报错解决方法1.pycocotools库的简介2.pycocotools安装的坑3.解决办法更多Ai资讯:公主号AiCharm本系列是作者在跑一些深度学习实例时,遇到的各种各样的问题及解决办法,希望能够帮助到大家。ERROR:Commanderroredoutwithexitstatus1:'D:\Anaconda3\python.exe'-u-c'importsys,setuptools,tokenize;sys.argv[0]='"'"'C:\\Users\\46653\\AppData\\Local\\Temp\\pip-instal
我写了一个非常简单的rake任务来尝试找到这个问题的根源。namespace:foodotaskbar::environmentdoputs'RUNNING'endend当在控制台中执行rakefoo:bar时,输出为:RUNNINGRUNNING当我执行任何rake任务时会发生这种情况。有没有人遇到过这样的事情?编辑上面的rake任务就是写在那个.rake文件中的所有内容。这是当前正在使用的Rakefile。requireFile.expand_path('../config/application',__FILE__)OurApp::Application.load_tasks这里
我完全不是程序员,正在学习使用Ruby和Rails框架进行编程。我目前正在使用Ruby1.8.7和Rails3.0.3,但我想知道我是否应该升级到Ruby1.9,因为我真的没有任何升级的“遗留”成本。缺点是什么?我是否会遇到与普通gem的兼容性问题,或者甚至其他我不太了解甚至无法预料的问题? 最佳答案 你应该升级。不要坚持从1.8.7开始。如果您发现不支持1.9.2的gem,请避免使用它们(因为它们很可能不被维护)。如果您对gem是否兼容1.9.2有任何疑问,您可以在以下位置查看:http://www.railsplugins.or
我正在尝试按Rails相关模型中的字段进行排序。我研究的所有解决方案都没有解决如果相关模型被另一个参数过滤?元素模型classItem相关模型:classPriority我正在使用where子句检索项目:@items=Item.where('company_id=?andapproved=?',@company.id,true).all我需要按相关表格中的“位置”列进行排序。问题在于,在优先级模型中,一个项目可能会被多家公司列出。因此,这些职位取决于他们拥有的company_id。当我显示项目时,它是针对一个公司的,按公司内的职位排序。完成此任务的正确方法是什么?感谢您的帮助。PS-我
如何学习ruby的正则表达式?(对于假人) 最佳答案 http://www.rubular.com/在Ruby中使用正则表达式时是一个很棒的工具,因为它可以立即将结果可视化。 关于ruby-我如何学习ruby的正则表达式?,我们在StackOverflow上找到一个类似的问题: https://stackoverflow.com/questions/1881231/
这个问题在这里已经有了答案:关闭10年前。PossibleDuplicate:HowdoIgeneratealistofnuniquerandomnumbersinRuby?我想做的事:Random.rand(0..10).timesdoputsRandom.rand(0..10)end但如果随机数已经显示过,则无法再次显示。如何最轻松地做到这一点?
假设您有一个可执行文件foo.rb,其库bar.rb的布局如下:/bin/foo.rb/lib/bar.rb在foo.rb的header中放置以下要求以在bar.rb中引入功能:requireFile.dirname(__FILE__)+"../lib/bar.rb"只要对foo.rb的所有调用都是直接的,这就可以正常工作。如果你把$HOME/project和符号链接(symboliclink)foo.rb放入$HOME/usr/bin,然后__FILE__解析为$HOME/usr/bin/foo.rb,因此无法找到bar.rb关于foo.rb的目录名.我意识到像rubygems这