写R程序的人，相信都会遇到过“cannot allocate vector of size”或者“无法分配大小为…的矢量”这样的错误。原因很简单，基本都是产生一个大矩阵等对象时发生的，最干脆的解决办法...

使用R计算相关系数的函数为： cor.test(X,Y,method="") method可以为"spearman","pearson" and "kendall",分别对应三种相关系数的计算和检验。...

还在纠结用如何用R读取Excel的同学有福了。昨天逛CRAN的时候发现了一个xlsx包，它给出的介绍是可以读取、写入Excel 2007/2003文件并支持格式的设置。简单地来说，将Excel读取为数...

Metropolis Hasting（下面简称MH）是蒙特卡罗马尔科夫链中一种重要的抽样方法。本文简要介绍MH算法，并给出一个实例。 MH算法在参数空间随机取值，作为起始点。按照参数的概率分布生成随机...

有时,在推断进化树的过程中,我们希望phylip格式的文件能够将序列的名称完整得保留下来。通常情况下,我们的名称会超过10字符, 而ClustalX和MUSCLE等软件给出的phylip文件是严格的p...

在以往R语言的教材、帮助示例和各类读书笔记中，有关图形中标注公式的介绍只有指定符号和数字的介绍。如何像SPSS那样在计算结束后将带参数的模型计算公式一并绘制在图形中呢？其实提取模型参数很简单，关键是如...

在做生物信息数据挖掘、分析之后为了可视化，我们常常需要画各种图。这里给大家推荐一个很不错的网站。http://addictedtor.free.fr/graphiques/ 这个网站里面包含了大量的利...

在同一页面上绘制多张图用layout()函数，在一幅图上添加新图，得到组合图用par()函数，下面举个例子。 data <- read.table("http://www.ats.ucla.ed...

目前用于二代测序Reads中错误校正的工具比较多，如SOAP的correction、Quake等，目前quake效果比较不错。下面给大家简单介绍quake软件的安装方法 安装 Quake在使用中需要一...

此图代表我的心 [code lang="R"] n=50000; r=0.7;r_e=(1-r*r)^.5; X=rnorm(n); Y=X*r+r_e*rnorm(n); Y=ifelse(X>...