<?xml version="1.0" encoding="UTF-8"?><rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>生信菜鸟团 &#187; 总结</title>
	<atom:link href="http://www.bio-info-trainee.com/tag/%e6%80%bb%e7%bb%93/feed" rel="self" type="application/rss+xml" />
	<link>http://www.bio-info-trainee.com</link>
	<description>欢迎去论坛biotrainee.com留言参与讨论，或者关注同名微信公众号biotrainee</description>
	<lastBuildDate>Sat, 28 Jun 2025 14:30:13 +0000</lastBuildDate>
	<language>zh-CN</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>https://wordpress.org/?v=4.1.33</generator>
	<item>
		<title>R语言中的排序，集合运算，reshape，以及merge总结</title>
		<link>http://www.bio-info-trainee.com/1071.html</link>
		<comments>http://www.bio-info-trainee.com/1071.html#comments</comments>
		<pubDate>Thu, 22 Oct 2015 02:34:46 +0000</pubDate>
		<dc:creator><![CDATA[ulwvfje]]></dc:creator>
				<category><![CDATA[R]]></category>
		<category><![CDATA[总结]]></category>

		<guid isPermaLink="false">http://www.bio-info-trainee.com/?p=1071</guid>
		<description><![CDATA[一直以来我都是随便看了点R的编程教程，因为我学了一点点C，所以还算有基础，现在基 &#8230; <a href="http://www.bio-info-trainee.com/1071.html">Continue reading <span class="meta-nav">&#8594;</span></a>]]></description>
				<content:encoded><![CDATA[<p>一直以来我都是随便看了点R的编程教程，因为我学了一点点C，所以还算有基础，现在基本上简单看看教程就能懂一门语言了，区别只是熟练度而已。R用得比较多，所以还算擅长，但是很多快捷应用的地方，我总是寄希望于到时候再查资料，所以没能用心的记住，这次花了点时间好好整理了一下R里面关于数据操作的重点，我想，以后再碰到类似的数据处理要求，应该很快能解决了把。</p>
<h1>首先看看排序：</h1>
<p>在R中，和排序相关的函数主要有三个：sort()，rank()，order()。</p>
<p><strong>sort(x)</strong><strong>是对向量x</strong><strong>进行排序，返回值排序后的数值向量。</strong></p>
<p><strong>rank()</strong><strong>是求秩的函数，它的返回值是这个向量中对应元素的“排名”。</strong></p>
<p><strong>order()</strong><strong>的返回值是对应“排名”的元素所在向量中的位置。</strong></p>
<p>其中sort(x)等同于x[order(x)]</p>
<p>下面以一小段R代码来举例说明：</p>
<p>&gt; x&lt;-c(97,93,85,74,32,100,99,67)</p>
<p>&gt; sort(x)</p>
<p>[1]  32  67  74  85  93  97  99 100</p>
<p>&gt; order(x)</p>
<p>[1] 5 8 4 3 2 1 7 6</p>
<p>&gt; rank(x)</p>
<p>[1] 6 5 4 3 1 8 7 2</p>
<p>&gt; x[order(x)]</p>
<p>[1]  32  67  74  85  93  97  99 100</p>
<p>其中比较有用的order，它可以用来给数据框进行排序</p>
<p>dat[order(dat[,1]),] 以该数据框的第一列进行排序</p>
<p><strong>dat[order(dat[,1],dat[,2]),] </strong><strong>以该数据框的第一列为主要次序，第二列为次要序列进行排序</strong></p>
<p>&nbsp;</p>
<h1>然后我们看看集合运算：</h1>
<p>在R里面除了简单的对两个向量求交集并集补集之外，比较重要的就是match和 %in% 了，需要重点讲讲。</p>
<p>#首先对集合A,B,C赋值</p>
<p>&gt; A&lt;-1:10</p>
<p>&gt; B&lt;-seq(5,15,2)</p>
<p>&gt; C&lt;-1:5</p>
<p>&gt; #求A和B的并集</p>
<p>&gt; union(A,B)</p>
<p>[1]  1  2  3  4  5  6  7  8  9 10 11 13 15</p>
<p>&gt; #求A和B的交集</p>
<p>&gt; intersect(A,B)</p>
<p>[1] 5 7 9</p>
<p>&gt; #求A-B</p>
<p>&gt; setdiff(A,B)</p>
<p>[1]  1  2  3  4  6  8 10</p>
<p>&gt; #求B-A</p>
<p>&gt; setdiff(B,A)</p>
<p>[1] 11 13 15</p>
<p>&gt; #检验集合A,B是否相同</p>
<p>&gt; setequal(A,B)</p>
<p>[1] FALSE</p>
<p>&gt; #检验元素12是否属于集合C</p>
<p>&gt; is.element(12,C)</p>
<p>[1] FALSE</p>
<p>&gt; #检验集合A是否包含C</p>
<p>&gt; all(C%in%A)</p>
<p>[1] TRUE</p>
<p>&gt; all(C%in%B)</p>
<p><strong>从上面可以看到%in%</strong><strong>这个操作符只返回逻辑向量TRUE </strong><strong>或者FALSE</strong><strong>，而且返回值应该与%in%</strong><strong>这个操作符前面的向量程度相等。也就是说它相当于遍历了C</strong><strong>里面的一个个元素，判断它们是否在B</strong><strong>中出现过，然后返回是或者否即可。</strong></p>
<p><strong>而match(C,B)</strong><strong>的结果就很不一样了，它的返回结果同样与前面的向量等长，但是它并非返回逻辑向量，而是遍历了C</strong><strong>里面的一个个元素，判断它们是否在B</strong><strong>中出现过，如果出现就返回在B</strong><strong>中的索引号，如果没有出现，就返回NA</strong><strong>。 </strong></p>
<p>&gt;B</p>
<p>[1]  5  7  9 11 13 15</p>
<p>&gt;C</p>
<p>[1] 1 2 3 4 5</p>
<p>&gt;match(C,B)</p>
<p>[1] NA NA NA NA  1</p>
<p>&gt;C%in%B</p>
<p>[1] FALSE FALSE FALSE FALSE  TRUE</p>
<h1> 接下来我们看看reshape：</h1>
<p>这是一个需要安装的包，起得就是R里面最经典的把长型数据变宽，和把宽数据拉长的作用。</p>
<p>其中melt函数是把很宽的数据拉长，它就是需要指定几列数据是保证不被融合的， 其余每一列数据都必须被融合到一列了，融合后的这一列数据每个元素旁边就用列名来标记该数据来自于哪一列。</p>
<p># example of melt function<br />
library(reshape)<br />
mdata &lt;- melt(mydata, id=c("id","time"))</p>
<p>融合后的数据如下：</p>
<table>
<tbody>
<tr>
<td><strong>id</strong></td>
<td><strong>time</strong></td>
<td><strong>variable</strong></td>
<td><strong>value</strong></td>
</tr>
<tr>
<td>1</td>
<td>1</td>
<td>x1</td>
<td>5</td>
</tr>
<tr>
<td>1</td>
<td>2</td>
<td>x1</td>
<td>3</td>
</tr>
<tr>
<td>2</td>
<td>1</td>
<td>x1</td>
<td>6</td>
</tr>
<tr>
<td>2</td>
<td>2</td>
<td>x1</td>
<td>2</td>
</tr>
<tr>
<td>1</td>
<td>1</td>
<td>x2</td>
<td>6</td>
</tr>
<tr>
<td>1</td>
<td>2</td>
<td>x2</td>
<td>5</td>
</tr>
<tr>
<td>2</td>
<td>1</td>
<td>x2</td>
<td>1</td>
</tr>
<tr>
<td>2</td>
<td>2</td>
<td>x2</td>
<td>4</td>
</tr>
</tbody>
</table>
<p>可以看到<strong>variable</strong><strong>列里面有两个不同的元素，说明是把旧数据中的两列给融合了，融合后的一个很长的列就是</strong><strong>value</strong></p>
<p>而cast函数的功能就是把刚才融合好的数据给还原。</p>
<p># cast the melted data<br />
# cast(data, formula, function)<br />
subjmeans &lt;- cast(mdata, id~variable, mean)<br />
timemeans &lt;- cast(mdata, time~variable, mean)</p>
<p>可以看到它们返回的结果是：</p>
<p><strong>subjmeans</strong></p>
<table>
<tbody>
<tr>
<td><strong>id</strong></td>
<td><strong>x1</strong></td>
<td><strong>x2</strong></td>
</tr>
<tr>
<td>1</td>
<td>4</td>
<td>5.5</td>
</tr>
<tr>
<td>2</td>
<td>4</td>
<td>2.5</td>
</tr>
</tbody>
</table>
<p><strong>timemeans</strong></p>
<table>
<tbody>
<tr>
<td><strong>time</strong></td>
<td><strong>x1</strong></td>
<td><strong>x2</strong></td>
</tr>
<tr>
<td>1</td>
<td>5.5</td>
<td>3.5</td>
</tr>
<tr>
<td>2</td>
<td>2.5</td>
<td>4.5</td>
</tr>
</tbody>
</table>
<p>可以看到cast函数比较复杂一点，formula公式右边的变量是需要拆开的variable，这一列变量有多少不重复元素，就新增多少列，左边的变量是行标记了，其余所有数据都会被计算一下再放在合适的位置。</p>
<p>在reshape2这个包里面还进化出来dcast和acast函数，功能大同小异。</p>
<p>&nbsp;</p>
<h1>最后我们来看看merge函数：</h1>
<p>这个函数的功能非常强大，类似于SQL语句里面的join系列函数</p>
<p>测试数据如下，它们这两个表的连接是作者名</p>
<p><a href="http://www.bio-info-trainee.com/wp-content/uploads/2015/10/Capture1.png"><img class="alignnone  wp-image-1072" src="http://www.bio-info-trainee.com/wp-content/uploads/2015/10/Capture1.png" alt="Capture1" width="588" height="297" /></a></p>
<p>&nbsp;</p>
<p>如果要实现类似sql里面的inner join 功能，则用代码</p>
<p>m1 &lt;- merge(authors, books, by.x = "surname", by.y = "name")</p>
<p>如果要实现left join功能则用代码</p>
<p>m1 &lt;- merge(authors, books, by.x = "surname", by.y = "name",all.x=TRUE)</p>
<p>right join功能代码</p>
<p>m1 &lt;- merge(authors, books, by.x = "surname", by.y = "name",all.y=TRUE)</p>
<p>all join功能代码</p>
<p>m1 &lt;- merge(authors, books, by.x = "surname", by.y = "name",all=TRUE)</p>
<p>关于单变量匹配的总结就是这些，但对于多变量匹配呢，例如下面两个表，需要对k1,k2两个变量都相等的情况下匹配</p>
<p>x &lt;- data.frame(k1 = c(NA,NA,3,4,5), k2 = c(1,NA,NA,4,5), data = 1:5)</p>
<p>y &lt;- data.frame(k1 = c(NA,2,NA,4,5), k2 = c(NA,NA,3,4,5), data = 1:5)</p>
<p>匹配代码如下merge(x, y, by = c("k1","k2"))  #inner join</p>
<p>另外一个多行匹配的例子如下：</p>
<p><a href="http://www.bio-info-trainee.com/wp-content/uploads/2015/10/Capture2.png"><img class="alignnone  wp-image-1073" src="http://www.bio-info-trainee.com/wp-content/uploads/2015/10/Capture2.png" alt="Capture2" width="550" height="323" /></a></p>
<p>我们的测试数据如上，这两个表的连接在于作者名。下面这两个语句等价</p>
<p>merge(authors, books, by=1:2)</p>
<p>merge(authors, books, by.x=c("FirstName", "LastName"),</p>
<p>by.y=c("AuthorFirstName", "AuthorLastName"),</p>
<p>all.x=TRUE)</p>
<p>都可以把两张表关联起来。</p>
<p>当然，在我搜索资料的时候，发现了另外一个解决问题的方法：</p>
<p>A[with(A, paste(C1, C2, sep = "\r")) %in% with(B, paste(C1, C2, sep="\r")), ]</p>
<p>&nbsp;</p>
<p>&nbsp;</p>
<p>参考：<a href="http://blog.sina.com.cn/s/blog_6caea8bf0100spe9.html">http://blog.sina.com.cn/s/blog_6caea8bf0100spe9.html</a></p>
<p><a href="http://blog.sina.com.cn/s/blog_6caea8bf010159dt.html">http://blog.sina.com.cn/s/blog_6caea8bf010159dt.html</a></p>
<p><a href="http://blog.csdn.net/u014801157/article/details/24372441">http://blog.csdn.net/u014801157/article/details/24372441</a></p>
<p><a href="http://www.statmethods.net/management/reshape.html">http://www.statmethods.net/management/reshape.html</a></p>
<p><a href="https://docs.tibco.com/pub/enterprise-runtime-for-R/1.5.0_may_2013/TERR_1.5.0_LanguageRef/base/merge.html">https://docs.tibco.com/pub/enterprise-runtime-for-R/1.5.0_may_2013/TERR_1.5.0_LanguageRef/base/merge.html</a></p>
<p><a href="http://r.789695.n4.nabble.com/Matching-multiple-columns-in-a-data-frame-td890832.html">http://r.789695.n4.nabble.com/Matching-multiple-columns-in-a-data-frame-td890832.html</a></p>
<p><a href="http://bbs.pinggu.org/thread-3234639-1-1.html">http://bbs.pinggu.org/thread-3234639-1-1.html</a></p>
<p><a href="http://www.360doc.com/content/14/0821/13/18951041_403561544.shtml">http://www.360doc.com/content/14/0821/13/18951041_403561544.shtml</a></p>
<p><a href="http://developer.51cto.com/art/201312/423612.htm">http://developer.51cto.com/art/201312/423612.htm</a></p>
]]></content:encoded>
			<wfw:commentRss>http://www.bio-info-trainee.com/1071.html/feed</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
	</channel>
</rss>
