<?xml version="1.0" encoding="UTF-8"?>
<!-- generator="bbPress/1.0.3" -->
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom">
	<channel>
		<title>大杂烩 &#187; Recent Posts</title>
		<link>http://www.yelinsky.com/notes/</link>
		<description>记录生物和计算机领域的点点滴滴</description>
		<language>en-US</language>
		<pubDate>Sat, 19 May 2012 04:26:56 +0000</pubDate>
		<generator>http://bbpress.org/?v=1.0.3</generator>
		<textInput>
			<title><![CDATA[Search]]></title>
			<description><![CDATA[Search all topics from these forums.]]></description>
			<name>q</name>
			<link>http://www.yelinsky.com/notes/search.php</link>
		</textInput>
		<atom:link href="http://www.yelinsky.com/notes/rss/" rel="self" type="application/rss+xml" />

		<item>
			<title>yelin on "C/C++比较字符串"</title>
			<link>http://www.yelinsky.com/notes/topic/39#post-39</link>
			<pubDate>Mon, 13 Feb 2012 10:27:35 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">39@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;C/C++字符串比较&#60;br /&#62;
1. 区分大小写&#60;br /&#62;
stricmp&#60;/p&#62;
&#60;p&#62;2. 不区分大小写&#60;br /&#62;
在Windows中用stricmp，在Linux中用strcasecmp
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "C/C++程序Segmentation fault"</title>
			<link>http://www.yelinsky.com/notes/topic/38#post-38</link>
			<pubDate>Thu, 09 Feb 2012 21:28:05 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">38@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;今天在写一个C++程序的时候，编译通过，但运行的时候出现Segmentation fault错误提示，在网上搜索了一下发现Segmentation fault往往是由于最内存操作不当引起的。经测试发现我的程序中的Segmentation fault原因是声明的数组太大，在&#60;a href=&#34;http://www.diybl.com/course/3_program/c++/cppjs/2007925/73592.html&#34;&#62;这里&#60;/a&#62;找到了解决方法：&#60;/p&#62;
&#60;p&#62;方法一:&#60;br /&#62;
  在VC的Project   setting里的link选项卡里把栈开大一点(windows里默认是4M)   &#60;/p&#62;
&#60;p&#62;方法二:&#60;br /&#62;
  声明成全局或static的,这两种变量不压栈,想开多大都可以   &#60;/p&#62;
&#60;p&#62;方法三:    &#60;/p&#62;
&#60;p&#62;	int   *A   =   new   int[90000];&#60;br /&#62;
	.....&#60;br /&#62;
	delete   A;&#60;/p&#62;
&#60;p&#62;方法四:&#60;br /&#62;
  用vector&#60;br /&#62;
	#include   &#38;lt;vector&#38;gt;   &#60;/p&#62;
&#60;p&#62;	using   namespace   std;   &#60;/p&#62;
&#60;p&#62;	void   main()&#60;br /&#62;
	{&#60;br /&#62;
	      vector&#38;lt;int&#38;gt;   A(90000);&#60;br /&#62;
	      A[0]   =   1;&#60;br /&#62;
	}&#60;/p&#62;
&#60;p&#62;我用了第二种方法。
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "awk字符串函数"</title>
			<link>http://www.yelinsky.com/notes/topic/37#post-37</link>
			<pubDate>Tue, 07 Feb 2012 10:14:40 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">37@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;awk是linux下的一种强大的文本处理工具，下面是awk主要的字符串处理函数：&#60;/p&#62;
&#60;p&#62;length(s)          返回s长度&#60;br /&#62;
gsub(r,s)           在整个$0中用s替代r&#60;br /&#62;
gsub(r,s,t)         在整个t中用s替代r&#60;br /&#62;
index(s,t)          返回s中字符串t的第一位置&#60;br /&#62;
match(s,r)         检查s是否包含匹配r的字符串&#60;br /&#62;
split(s,a,fs)        在fs上将s分成序列a&#60;br /&#62;
sprint(fmt,exp) 返回经fmt格式化后的exp&#60;br /&#62;
sub(r,s)             用$0中最左边最长的子串代替s&#60;br /&#62;
substr(s,p)        返回字符串s中从p开始的后缀部分&#60;br /&#62;
substr(s,p,n)     返回字符串s中从p开始长度为n的后缀部分
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Linux 加载、卸载ntfs分区方法"</title>
			<link>http://www.yelinsky.com/notes/topic/36#post-36</link>
			<pubDate>Thu, 02 Feb 2012 15:48:47 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">36@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Linux上加载ntfs分区：&#60;/p&#62;
&#60;p&#62;sudo mkdir temp&#60;br /&#62;
sudo mount /dev/sda1 ./temp -t ntfs -o umask=0222&#60;/p&#62;
&#60;p&#62;将umask=0222改为umask=0000，可以把读写权限开放给所有用户&#60;/p&#62;
&#60;p&#62;卸载分区：&#60;br /&#62;
sudo umount ./temp&#60;/p&#62;
&#60;p&#62;----&#60;br /&#62;
加载U盘：&#60;br /&#62;
sudo mount /dev/sdx1 ./temp -t vfat -o iocharset=utf8
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "C/C++统计一个字符在另一个字符中出现的次数"</title>
			<link>http://www.yelinsky.com/notes/topic/35#post-35</link>
			<pubDate>Tue, 31 Jan 2012 21:15:53 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">35@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;C/C++计算字符出现次数代码：&#60;/p&#62;
&#60;p&#62;----------------------------------1----------------------------------&#60;br /&#62;
#include&#38;lt;stdio.h&#38;gt;&#60;br /&#62;
#include&#38;lt;string.h&#38;gt;&#60;br /&#62;
void main()&#60;br /&#62;
{char str1[2000],str2[2000],*p1,*p2;&#60;br /&#62;
 int sum=0,i;&#60;br /&#62;
 printf(&#34;请输入二个字符串!\n&#34;);&#60;br /&#62;
 scanf(&#34;%s%s&#34;,str1,str2);&#60;br /&#62;
 p1=str1;&#60;br /&#62;
 p2=str2;&#60;br /&#62;
 while(strstr(p1,p2)!=NULL)  //判断有无出现&#60;br /&#62;
 {&#60;br /&#62;
  sum++;                  //有的话，次数加1&#60;br /&#62;
  p1=strstr(p1,p2)+1;     //取出现位置首地址加1形成新的串，再判断&#60;br /&#62;
 }&#60;br /&#62;
 printf(&#34;第二个字符串在第一个字符串中出现的次数为%d!\n&#34;,sum);&#60;br /&#62;
}&#60;/p&#62;
&#60;p&#62;----------------------------------2----------------------------------&#60;br /&#62;
#include &#38;lt;iostream&#38;gt;&#60;/p&#62;
&#60;p&#62;#include &#38;lt;string&#38;gt;&#60;/p&#62;
&#60;p&#62;using namespace std;&#60;/p&#62;
&#60;p&#62;int main()&#60;/p&#62;
&#60;p&#62;{&#60;/p&#62;
&#60;p&#62;            int n=0,sum=0;&#60;/p&#62;
&#60;p&#62;            string s1(&#34;qwertyqwertyqwerty&#34;);&#60;/p&#62;
&#60;p&#62;            string s2(&#34;ert&#34;);&#60;/p&#62;
&#60;p&#62;            string temp=s1;&#60;/p&#62;
&#60;p&#62;            while (1){                                 //循环：查找到子串就把子串删去&#60;/p&#62;
&#60;p&#62;                        n=temp.find(s2);                     //返回子串的位置&#60;/p&#62;
&#60;p&#62;                        if (n!=-1){        //n=-1表示未找到子串&#60;/p&#62;
&#60;p&#62;                                    temp=temp.substr(n+s2.length(),temp.length()-s2.length());   //开一个新的字符串变量temp存储删去子串后的字符串&#60;/p&#62;
&#60;p&#62;                                    sum++;           //出现次数&#60;/p&#62;
&#60;p&#62;                        }&#60;/p&#62;
&#60;p&#62;                        else break;&#60;/p&#62;
&#60;p&#62;            }&#60;/p&#62;
&#60;p&#62;            cout&#38;lt;&#38;lt;sum;&#60;/p&#62;
&#60;p&#62;}&#60;/p&#62;
&#60;p&#62;转自：&#60;a href=&#34;http://blog.csdn.net/qucooln/article/details/3280617&#34;&#62;http://blog.csdn.net/qucooln/article/details/3280617&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Linux结束进程释放内存"</title>
			<link>http://www.yelinsky.com/notes/topic/34#post-34</link>
			<pubDate>Tue, 31 Jan 2012 15:44:45 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">34@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Ubuntu系统结束进程释放内存命令&#60;/p&#62;
&#60;p&#62;kill -l PID   或者&#60;br /&#62;
kill -9 PID&#60;/p&#62;
&#60;p&#62;PID用ps或top命令查看
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Python正则表达式"</title>
			<link>http://www.yelinsky.com/notes/topic/33#post-33</link>
			<pubDate>Mon, 16 Jan 2012 10:20:58 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">33@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;&#60;strong&#62;1. 字符串查找 (match, search)&#60;/strong&#62;&#60;br /&#62;
match checks for a match only at the beginning of the string, while search checks for a match anywhere in the string&#60;/p&#62;
&#60;p&#62;#!/usr/bin/python&#60;br /&#62;
import re&#60;br /&#62;
line = &#34;Cats are smarter than dogs&#34;;&#60;br /&#62;
matchObj = re.match( r'dogs', line, re.M&#124;re.I)&#60;br /&#62;
if matchObj:&#60;br /&#62;
   print &#34;match --&#38;gt; matchObj.group() : &#34;, matchObj.group()&#60;br /&#62;
else:&#60;br /&#62;
   print &#34;No match!!&#34;&#60;/p&#62;
&#60;p&#62;matchObj = re.search( r'dogs', line, re.M&#124;re.I)&#60;br /&#62;
if matchObj:&#60;br /&#62;
   print &#34;search --&#38;gt; matchObj.group() : &#34;, matchObj.group()&#60;br /&#62;
else:&#60;br /&#62;
   print &#34;No match!!&#34;&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;2. 字符串替换&#60;/strong&#62;&#60;br /&#62;
#!/usr/bin/python&#60;br /&#62;
import re&#60;br /&#62;
a = &#34;Hello world！&#34;&#60;br /&#62;
b = re.sub(r'world', &#34;python&#34;, a)&#60;br /&#62;
print b&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.tutorialspoint.com/python/python_reg_expressions.htm&#34;&#62;http://www.tutorialspoint.com/python/python_reg_expressions.htm&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "二项分布，泊松分布，正态分布"</title>
			<link>http://www.yelinsky.com/notes/topic/32#post-32</link>
			<pubDate>Sun, 15 Jan 2012 11:58:52 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">32@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;二项分布有两个参数，一个 n 表示试验次数，一个 p 表示一次试验成功概率。&#60;/p&#62;
&#60;p&#62;现在考虑一列二项分布，其中试验次数 n 无限增加，而 p 是 n 的函数。&#60;/p&#62;
&#60;p&#62;如果 np 存在有限极限 λ，则这列二项分布就趋于参数为  λ 的 泊松分布。反之，如果 np 趋于无限大（如 p 是一个定值），则根据德莫佛-拉普拉斯(De'Moivre-Laplace)中心极限定理，这列二项分布将趋近于正态分布。&#60;/p&#62;
&#60;p&#62;实际运用中当 n 很大时一般都用正态分布来近似计算二项分布，但是如果同时 np 又比较小（比起 n来说很小），那么用泊松分布近似计算更简单些，毕竟泊松分布跟二项分布一样都是离散型分布。&#60;/p&#62;
&#60;p&#62;转载自：&#60;a href=&#34;http://emuch.net&#34;&#62;小木虫&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Kmer分析工具"</title>
			<link>http://www.yelinsky.com/notes/topic/31#post-31</link>
			<pubDate>Tue, 03 Jan 2012 11:00:42 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">31@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;jellyfish是一个用于对DNA序列中的kmer进行分析和统计的工具，支持多核并行处理。&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.cbcb.umd.edu/software/jellyfish/&#34;&#62;http://www.cbcb.umd.edu/software/jellyfish/&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "并行压缩与解压缩"</title>
			<link>http://www.yelinsky.com/notes/topic/30#post-30</link>
			<pubDate>Fri, 30 Dec 2011 16:42:16 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">30@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;piz是一个可以运行在多核机器上进行并行压缩与解压缩到软件，可以完全替代gzip，大大提高压缩与解压缩到速度。&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.zlib.net/pigz/&#34;&#62;http://www.zlib.net/pigz/&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Wordpress留言链接在新窗口中打开"</title>
			<link>http://www.yelinsky.com/notes/topic/29#post-29</link>
			<pubDate>Wed, 07 Dec 2011 17:01:04 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">29@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;修改文件的方法让&#60;a href=&#34;http://www.yelinsky.com/blog/archives/tag/wordpress&#34;&#62;Wordpress&#60;/a&#62;留言链接在新窗口中打开：&#60;/p&#62;
&#60;p&#62;打开wp-includes目录下的comment-template.php文件，找到第get_comment_author_link()这个函数，在第else $return 这里的a标签里面加入target='_blank' 属性，即可。
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "GNU parallel 教程"</title>
			<link>http://www.yelinsky.com/notes/topic/28#post-28</link>
			<pubDate>Thu, 10 Nov 2011 12:52:50 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">28@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;GNU parallel使用说明及示例。&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.gnu.org/software/parallel/man.html&#34;&#62;http://www.gnu.org/software/parallel/man.html&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "什么是fragment recruitment和recruitment plots？"</title>
			<link>http://www.yelinsky.com/notes/topic/26#post-26</link>
			<pubDate>Tue, 08 Nov 2011 22:26:19 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">26@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;&#60;strong&#62;What is fragment recruitment?&#60;/strong&#62;&#60;br /&#62;
Fragment recruitment is the process of aligning short sequences to a reference sequence and filtering out low quality matches. The short sequences are typically sequencing reads generated via Sanger or Roche 454 sequencing machines. Reference sequences are typically finished or draft genomes of microbes and viruses. However it is worth noting that reference sequences can be other things including assemblies derived from metagenomic reads.&#60;br /&#62;
In the fragment recruitment process sequencing reads are aligned using NCBI BLAST to all the available reference sequences. The BLAST parameters are tuned to allow large gaps and very low percent identity matches. Aligned sequencing reads are then stringently filtered so that only reads that align over ~90% of their length are retained.&#60;/p&#62;
&#60;p&#62;From a biological perspective, reads have to belong to microbes of the same genus as the reference to be recruited. However this is not always the case as there are numerous examples of highly conserved genes that will recruit reads from even distantly related organisms.&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;What are fragment recruitment plots?&#60;/strong&#62;&#60;br /&#62;
Fragment recruitment plots are a visual representation of the fragment recruitment data. Effectively these are fancy percent identity plots modified to display metagenomic data. In these plots the position on the reference sequence is shown on the x-axis and the percent identity of the alignments is shown on the y-axis.&#60;/p&#62;
&#60;p&#62;from &#60;a href=&#34;http://gos.jcvi.org&#34;&#62;http://gos.jcvi.org&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "去除文件中的重复行"</title>
			<link>http://www.yelinsky.com/notes/topic/25#post-25</link>
			<pubDate>Fri, 04 Nov 2011 21:27:09 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">25@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;去除文本文件中的重复行，用awk命令，非常简单：&#60;/p&#62;
&#60;p&#62;awk '!a[$0]++' input.txt &#38;gt; output.txt&#60;/p&#62;
&#60;p&#62;或者用sort:&#60;/p&#62;
&#60;p&#62;sort -u input.txt&#38;gt; output.txt
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "一些高通量测序问题"</title>
			<link>http://www.yelinsky.com/notes/topic/24#post-24</link>
			<pubDate>Mon, 31 Oct 2011 11:25:57 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">24@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;叶林，你好！&#60;/p&#62;
&#60;p&#62;我是xxx所的博士生。&#60;br /&#62;
老板近期申请了一个项目是想针对生物信息这一领域从计算机系统结构的角度做一些优化工作，打算是从基因测序开始着手。目前我在做一些前期调研的工作，主要是看了一些测序的算法，但有些实际应用中的问题不是太清楚。&#60;br /&#62;
看到了你的博客，非常佩服，获得了不少有用的信息。还有一些疑问，想请教一下：&#60;/p&#62;
&#60;p&#62;1. 博文中有提到微生物领域的数据量不是太大，不知这个数据量是多少？微生物测序主要是resequencing还是de novo sequencing?&#60;br /&#62;
2. 是否mapping在工作站上可以完成，而assembly需要大型机？需要大型机的原因是CPU受限还是内存容量受限？&#60;br /&#62;
3. 现在mapping和assembly的软件非常多，结果相差似乎也很大，不知实际中常用的一般是什么？耗时的情况如何？比如100M读长100bp的reads的mapping和assembly时间分别是多少？&#60;br /&#62;
4. 目前考虑到了在mapping过程中加入indel增加对gap的支持，不知道现在如100bp的短序列比对对Gap是否有要求？&#60;/p&#62;
&#60;p&#62;呵呵，问题较多，数字有个大致的就好，还望不吝赐教！&#60;br /&#62;
多谢多谢！&#60;/p&#62;
&#60;p&#62;----------------------------------&#60;/p&#62;
&#60;p&#62;**你好&#60;/p&#62;
&#60;p&#62;1.&#60;a href=&#34;http://www.yelinsky.com/blog/archives/tag/%e9%ab%98%e9%80%9a%e9%87%8f%e6%b5%8b%e5%ba%8f&#34;&#62; 高通量测序&#60;/a&#62;数据量是多少，要根据样品的复杂度确定，纯菌可能几十M就够了吧，复杂度环境样品几百G也可能不够。&#60;br /&#62;
2. assembly主要是内存受限制，我现在用的工作站96G内存，经常不够。有些单位用几百G甚至1T内存的电脑去做组装。mapping在工作站上进行应该问题不大。&#60;br /&#62;
3. 具体哪个软件好，我现在也不是非常清楚。组装过程时间都不是问题，主要是占内存。&#60;br /&#62;
4. 我不太懂。&#60;/p&#62;
&#60;p&#62;叶林
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Velvet使用说明"</title>
			<link>http://www.yelinsky.com/notes/topic/23#post-23</link>
			<pubDate>Thu, 20 Oct 2011 22:14:19 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">23@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;序列组装软件Velvet使用说明&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://wiki.bioinformatics.ucdavis.edu/index.php/Illumina_assembly_using_velvet&#34;&#62;http://wiki.bioinformatics.ucdavis.edu/index.php/Illumina_assembly_using_velvet&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Ubuntu安装easy_install"</title>
			<link>http://www.yelinsky.com/notes/topic/22#post-22</link>
			<pubDate>Mon, 17 Oct 2011 20:01:34 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">22@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;使用easy_install可是使ubuntu中的很多软件安装变得非常简单，安装easy_install的命令如下：&#60;/p&#62;
&#60;p&#62;sudo apt-get install python-setuptools
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "在Ubuntu上安装Velvet"</title>
			<link>http://www.yelinsky.com/notes/topic/21#post-21</link>
			<pubDate>Mon, 17 Oct 2011 19:44:01 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">21@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;&#60;a href=&#34;http://www.ebi.ac.uk/~zerbino/velvet/&#34;&#62;Velvet&#60;/a&#62;是著名的短序列组装软件，只能在Linux系统中运行。在Ubuntu中的安装方法如下：&#60;/p&#62;
&#60;p&#62;1. 安装zlib1g-dev和pdflatex&#60;br /&#62;
sudo apt-get install zlib1g-dev&#60;br /&#62;
sudo aptitude install texlive&#60;br /&#62;
2. 下载并安装velvet_1.1.06.tgz&#60;br /&#62;
tar xvfz velvet_1.1.06.tgz&#60;br /&#62;
cd velvet_1.1.06&#60;br /&#62;
make MAXKMERLENGTH=61&#60;/p&#62;
&#60;p&#62;参考：&#60;br /&#62;
（1） &#60;a href=&#34;http://molecularevolution.org/resources/activities/velvetbowtie_activity#compilingvelvet&#34;&#62;http://molecularevolution.org/resources/activities/velvetbowtie_activity#compilingvelvet&#60;/a&#62;&#60;br /&#62;
（2） &#60;a href=&#34;http://socwiki.wordpress.com/2011/06/24/velvet-and-oases-for-transcriptome-assembly/&#34;&#62;http://socwiki.wordpress.com/2011/06/24/velvet-and-oases-for-transcriptome-assembly/&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "从Metagenome数据中筛选rRNA序列的工具：rRNASelector"</title>
			<link>http://www.yelinsky.com/notes/topic/20#post-20</link>
			<pubDate>Fri, 14 Oct 2011 11:52:12 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">20@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;rRNASelector是一个用于从Metagenome数据中筛选rRNA序列的工具，基于Hidden Markov模型，使用Java语言开发。&#60;/p&#62;
&#60;p&#62;下载地址：&#60;a href=&#34;http://sw.ezbiocloud.net/rrnaselector&#34;&#62;http://sw.ezbiocloud.net/rrnaselector&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Ubuntu无法上网networking disabled解决方法"</title>
			<link>http://www.yelinsky.com/notes/topic/19#post-19</link>
			<pubDate>Thu, 13 Oct 2011 13:24:37 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">19@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Ubuntu休眠之后（suspend）就无法上网了，显示networking disabled，这是Ubuntu一个严重的bug，解决方法如下：&#60;/p&#62;
&#60;p&#62;sudo service network-manager stop&#60;br /&#62;
sudo rm /var/lib/NetworkManager/NetworkManager.state&#60;br /&#62;
sudo service network-manager start
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Ubuntu安装SSH Server"</title>
			<link>http://www.yelinsky.com/notes/topic/18#post-18</link>
			<pubDate>Thu, 13 Oct 2011 10:51:55 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">18@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Ubuntu似乎没有自带SSH Server程序，在Ubuntu上安装SSH Server的命令如下：&#60;br /&#62;
sudo apt-get install openssh-server&#60;/p&#62;
&#60;p&#62;启动与重启：&#60;br /&#62;
sudo /etc/init.d/ssh start&#60;br /&#62;
sudo /etc/init.d/ssh restart
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Windows设置环境变量"</title>
			<link>http://www.yelinsky.com/notes/topic/17#post-17</link>
			<pubDate>Thu, 13 Oct 2011 10:28:48 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">17@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;通过命令行设置Windows（Win7, Vista, XP) 设置环境变量以Python为例，方法如下：&#60;/p&#62;
&#60;p&#62;开始-运行-cmd&#60;br /&#62;
set PATH=%PATH%; C:\Python27&#60;/p&#62;
&#60;p&#62;查看所有环境变量：&#60;br /&#62;
set PATH&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.yelinsky.com/notes/topic/16&#34;&#62;Linux设置环境变量的方法&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Linux设置环境变量"</title>
			<link>http://www.yelinsky.com/notes/topic/16#post-16</link>
			<pubDate>Wed, 12 Oct 2011 23:49:09 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">16@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Linux设置环境变量，方法很多，下面是其中一种方法：&#60;/p&#62;
&#60;p&#62;1. 编辑.bashrc&#60;br /&#62;
vim ~/.bashrc&#60;br /&#62;
添加：&#60;br /&#62;
PATH=$PATH:XXXXXX&#60;br /&#62;
export PATH&#60;br /&#62;
2. 执行&#60;br /&#62;
source ~/.bashrc&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://www.yelinsky.com/notes/topic/16&#34;&#62;设置Windows环境变量的方法&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Ubuntu设置固定IP"</title>
			<link>http://www.yelinsky.com/notes/topic/15#post-15</link>
			<pubDate>Wed, 12 Oct 2011 23:42:36 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">15@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Ubuntu Desktop设置固定IP：&#60;/p&#62;
&#60;p&#62;System-Preferences-Network Connections.&#60;br /&#62;
IPV4 Settings-Add
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "awk字符串截取"</title>
			<link>http://www.yelinsky.com/notes/topic/14#post-14</link>
			<pubDate>Mon, 10 Oct 2011 11:53:07 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">14@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;awk字符串截取函数substr的基本用法：&#60;br /&#62;
substr(s,p) 返回字符串s中从p开始的后缀部分&#60;br /&#62;
substr(s,p,n) 返回字符串s中从p开始长度为n的后缀部分&#60;/p&#62;
&#60;p&#62;去掉字符串最后两个字符：&#60;br /&#62;
substr(s,1,length(s)-2)
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "用R绘制热图(heatmap)的方法"</title>
			<link>http://www.yelinsky.com/notes/topic/13#post-13</link>
			<pubDate>Sat, 17 Sep 2011 23:54:43 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">13@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;用R绘制热图(heatmap)并生成PDF文件：&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;1. Install gplots and RColorBrewer for more color options&#60;/strong&#62;&#60;/p&#62;
&#60;p&#62;$ R&#60;br /&#62;
&#38;gt; install.packages(&#34;gplots&#34;)&#60;br /&#62;
&#38;gt; install.packages(&#34;RColorBrewer&#34;)&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;2. Copy the script below and save it as “heatmap.r”&#60;/strong&#62;&#60;/p&#62;
&#60;p&#62;library(RColorBrewer)&#60;br /&#62;
library(gplots)&#60;br /&#62;
x=read.table(&#34;matrix.dat&#34;, header=TRUE)&#60;br /&#62;
mat=data.matrix(x)&#60;/p&#62;
&#60;p&#62;pdf(&#34;heatmap.pdf&#34;, height=10, width=10)&#60;br /&#62;
heatmap.2(mat,&#60;br /&#62;
Rowv=TRUE,&#60;br /&#62;
Colv=TRUE,&#60;br /&#62;
#    dendrogram= c(&#34;none&#34;),&#60;br /&#62;
distfun = dist,&#60;br /&#62;
hclustfun = hclust,&#60;br /&#62;
xlab = &#34;X data&#34;, ylab = &#34;Y data&#34;,&#60;br /&#62;
key=TRUE,&#60;br /&#62;
keysize=1,&#60;br /&#62;
trace=&#34;none&#34;,&#60;br /&#62;
density.info=c(&#34;none&#34;),&#60;br /&#62;
margins=c(10, 8),&#60;br /&#62;
col=brewer.pal(10,&#34;PiYG&#34;)&#60;br /&#62;
#    col=redgreen(75),&#60;br /&#62;
)&#60;br /&#62;
dev.off()&#60;/p&#62;
&#60;p&#62;The commented line dendrogram=c(“none”) will suppress clustering of data in both axes. Another option for coloring the heatmap is given by col=redgreen(75), commonly used in log data in microarray chips.&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;3. Generate a data matrix as&#60;/strong&#62;&#60;/p&#62;
&#60;p&#62;        a       b       c       d       e&#60;br /&#62;
k       0       0       1       1       5&#60;br /&#62;
l       1       2       1       6       7&#60;br /&#62;
m       1       1       2       9       1&#60;br /&#62;
n       0       0       1       3       2&#60;br /&#62;
o       9       8       0       6       5&#60;/p&#62;
&#60;p&#62;and save it as “matrix.dat”&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;4. Run as&#60;/strong&#62;&#60;/p&#62;
&#60;p&#62;    $ R –vanilla &#38;lt; heatmap.r&#60;/p&#62;
&#60;p&#62;5. See the heat map in heatmap.pdf&#60;/p&#62;
&#60;p&#62;原文链接：&#60;a href=&#34;http://enotacoes.wordpress.com/2007/11/16/easy-guide-to-drawing-heat-maps-to-pdf-with-r-with-color-key/&#34;&#62;Easy guide to drawing heat maps to PDF with R (with color key)&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "awk处理两个文件的方法"</title>
			<link>http://www.yelinsky.com/notes/topic/12#post-12</link>
			<pubDate>Wed, 31 Aug 2011 20:23:36 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">12@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;&#60;strong&#62;awk处理两个文件的方法：&#60;/strong&#62;&#60;/p&#62;
&#60;p&#62;awk 'NR==FNR{...}NR&#38;gt;FNR{...}' file1 file2&#60;br /&#62;
# 读入file1的时候，已读入file1的记录数FNR一定等于awk已读入的总记录数NR，因为file1是awk读入的首个文件，故读入file1时执行前一个命令块{...}&#60;br /&#62;
# 读入file2的时候，已读入的总记录数NR一定&#38;gt;读入file2的记录数FNR，故读入file2时执行后一个命令块{...}&#60;/p&#62;
&#60;p&#62;awk 'NR==FNR{...;next}{...}' file1 file2&#60;br /&#62;
# 读入file1时，满足NR==FNR，先执行前一个命令块，但因为其中有next命令，故后一个命令块{...}是不会执行的&#60;br /&#62;
# 读入file2时，不满足NR==FNR，前一个命令块{..}不会执行，只执行后一个命令块{...}
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Linux中cut命令的用法"</title>
			<link>http://www.yelinsky.com/notes/topic/11#post-11</link>
			<pubDate>Mon, 29 Aug 2011 12:21:38 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">11@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;cut -c &#60;em&#62;m-n&#60;/em&#62; filename&#60;br /&#62;
显示&#60;strong&#62;每行&#60;/strong&#62;从开头算起&#60;em&#62;m-n&#60;/em&#62;的字符&#60;/p&#62;
&#60;p&#62;cut -f m-n&#60;br /&#62;
显示第m栏到第n栏(使用tab分隔)&#60;/p&#62;
&#60;p&#62;-b、-c、-f分别表示字节、字符、字段&#60;/p&#62;
&#60;p&#62;cut -d: -f 3 filename&#60;br /&#62;
-d用来定义分隔符，默认为tab键。
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "标准氨基酸代码（符号）表"</title>
			<link>http://www.yelinsky.com/notes/topic/10#post-10</link>
			<pubDate>Sat, 27 Aug 2011 13:05:12 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">10@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;标准氨基酸代码（符号）表&#60;br /&#62;
Standard IUB/IUPAC amino acid codes (symbols) table&#60;/p&#62;
&#60;p&#62;One-letter code	Amino acid	Three-letter code&#60;br /&#62;
A----	Alanine	Ala&#60;br /&#62;
B----	Aspartic acid or Asparagine	Asx&#60;br /&#62;
C----	Cysteine	Cys&#60;br /&#62;
D----	Aspartic acid	Asp&#60;br /&#62;
E----	Glutamic acid	Glu&#60;br /&#62;
F----	Phenylalanine	Phe&#60;br /&#62;
G----	Glycine	Gly&#60;br /&#62;
H----	Histidine	His&#60;br /&#62;
I----	Isoleucine	Ile&#60;br /&#62;
K----	Lysine	Lys&#60;br /&#62;
L----	Leucine	Leu&#60;br /&#62;
M----	Methionine	Met&#60;br /&#62;
N----	Asparagine	Asn&#60;br /&#62;
P----	Proline	Pro&#60;br /&#62;
Q----	Glutamine	Gln&#60;br /&#62;
R----	Arginine	Arg&#60;br /&#62;
S----	Serine	Ser&#60;br /&#62;
T----	Threonine	Thr&#60;br /&#62;
V----	Valine	Val&#60;br /&#62;
W----	Tryptophan	Trp&#60;br /&#62;
X----	Unknown	Xxx&#60;br /&#62;
Y----	Tyrosine	Tyr&#60;br /&#62;
Z----	Glutamic acid or Glutamine	Glx
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "标准核酸代码（符号）表"</title>
			<link>http://www.yelinsky.com/notes/topic/9#post-9</link>
			<pubDate>Sat, 27 Aug 2011 13:02:03 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">9@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;标准核酸（DNA、RNA）代码（符号）表&#60;br /&#62;
Standard IUB/IUPAC nucleic acid codes (symbols) table&#60;/p&#62;
&#60;p&#62;Code	Nucleic Acid(s)&#60;br /&#62;
A----	Adenine&#60;br /&#62;
C----	Cytosine&#60;br /&#62;
G----	Guanine&#60;br /&#62;
T----	Thymine&#60;br /&#62;
U----	Uracil&#60;br /&#62;
M----	A or C (amino)&#60;br /&#62;
R----	A or G (purine)&#60;br /&#62;
W----	A or T (weak)&#60;br /&#62;
S----	C or G (strong)&#60;br /&#62;
Y----	C or T (pyrimidine)&#60;br /&#62;
K----	G or T (keto)&#60;br /&#62;
V----	A or C or G&#60;br /&#62;
H----	A or C or T&#60;br /&#62;
D----	A or G or T&#60;br /&#62;
B----	C or G or T&#60;br /&#62;
N----	A or G or C or T (any)
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "短序列搜索匹配程序"</title>
			<link>http://www.yelinsky.com/notes/topic/8#post-8</link>
			<pubDate>Sat, 27 Aug 2011 00:14:26 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">8@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;PatMaN是一个C++编写的从大量序列中搜索查找某个短序列片段的程序，允许mismatch。可用于在大量序列中搜索Primer或Probe。&#60;/p&#62;
&#60;p&#62;&#60;a href=&#34;http://bioinf.eva.mpg.de/patman/&#34;&#62;http://bioinf.eva.mpg.de/patman/&#60;/a&#62;&#60;/p&#62;
&#60;p&#62;只能在Linux下运行。
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Ubuntu中查看CPU及内存信息的方法"</title>
			<link>http://www.yelinsky.com/notes/topic/7#post-7</link>
			<pubDate>Fri, 26 Aug 2011 22:29:55 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">7@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;&#60;strong&#62;1. 查看CPU信息&#60;/strong&#62;&#60;br /&#62;
cat /proc/cpuinfo&#60;/p&#62;
&#60;p&#62;&#60;strong&#62;2. 查看内存信息&#60;/strong&#62;&#60;br /&#62;
free -m&#60;br /&#62;
显示的数据的单位是M&#60;/p&#62;
&#60;p&#62;以上命令，除了Ubuntu，其它Linux系统也基本都可以用。
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Smith–Waterman算法实现的序列对齐程序"</title>
			<link>http://www.yelinsky.com/notes/topic/6#post-6</link>
			<pubDate>Thu, 25 Aug 2011 17:34:44 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">6@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;Smith Waterman算法是一种与Blast类似的寻找相似序列的算法。Blast算法的运行速度要比Smith Waterman算法快，但是Smith Waterman算法比BLAST算法更精确。&#60;br /&#62;
&#60;a href=&#34;http://en.wikipedia.org/wiki/Smith%E2%80%93Waterman_algorithm&#34;&#62;http://en.wikipedia.org/wiki/Smith%E2%80%93Waterman_algorithm&#60;/a&#62;&#60;/p&#62;
&#60;p&#62;JAligner是用java写的一个Smith–Waterman算法实现的序列对齐程序&#60;br /&#62;
&#60;a href=&#34;http://jaligner.sourceforge.net/&#34;&#62;http://jaligner.sourceforge.net/&#60;/a&#62;
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "Blast结果m8格式意义"</title>
			<link>http://www.yelinsky.com/notes/topic/5#post-5</link>
			<pubDate>Mon, 15 Aug 2011 12:11:39 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">5@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;进行Blast比对，用参数-m 8可以以列表的方式输出结果，结果中从左到右每一列的意义分别是：&#60;br /&#62;
[00] Query id&#60;br /&#62;
[01] Subject id&#60;br /&#62;
[02] % identity&#60;br /&#62;
[03] alignment length&#60;br /&#62;
[04] mismatches&#60;br /&#62;
[05] gap openings&#60;br /&#62;
[06] q. start&#60;br /&#62;
[07] q. end&#60;br /&#62;
[08] s. start&#60;br /&#62;
[09] s. end&#60;br /&#62;
[10] e-value&#60;br /&#62;
[11] bit score
&#60;/p&#62;</description>
		</item>
		<item>
			<title>yelin on "纯CSS图形"</title>
			<link>http://www.yelinsky.com/notes/topic/4#post-4</link>
			<pubDate>Mon, 15 Aug 2011 10:59:42 +0000</pubDate>
			<dc:creator>yelin</dc:creator>
			<guid isPermaLink="false">4@http://www.yelinsky.com/notes/</guid>
			<description>&#60;p&#62;只用CSS，在浏览器里绘制各种图形：&#60;br /&#62;
&#60;em&#62;正方形&#60;br /&#62;
长方形&#60;br /&#62;
圆形&#60;br /&#62;
椭圆形&#60;br /&#62;
三角形&#60;br /&#62;
平行四边行&#60;br /&#62;
梯形&#60;br /&#62;
六角星形&#60;br /&#62;
五角星形&#60;br /&#62;
五边形&#60;br /&#62;
六边形&#60;br /&#62;
八边形&#60;/em&#62;&#60;br /&#62;
链接地址：&#60;a href=&#34;http://css-tricks.com/examples/ShapesOfCSS/?=derp&#34;&#62;http://css-tricks.com/examples/ShapesOfCSS/?=derp&#60;/a&#62;
&#60;/p&#62;</description>
		</item>

	</channel>
</rss>

