R(又称R语言)是一款开源的跨平台的数值统计和数值图形化展现工具。通俗点说,R是用来做统计和画图的。R拥有自己的脚本语言和大量的统计、图形库(得益于开源社区),这让她看起来既美又实用。与其他同类软件(如 SPSS)相比,R的特点是纯命令行的,这倒也好,我们更应该把注意力放在数据本身,而非统计工具的UI。
R虽说有一套自己的语言,还挺完备,但它最专业的还是做统计和画图,而像连接数据库、文本处理、文件操作等这些脏活可不能委屈R来做,这些得有其他语言来负责,我的选择是咱最熟悉的、做这些脏活最棒的Python。那接下来的问题很清楚了,R和Python如何一起工作?我总结了2个方法来进行操作。
01. R和Python只共享文件
Python把源数据处理干净,生成格式化的文件放在预定的目录下,做个定时器让R去读文件,最终输出统计结果和图表。
这种做法一定程度上可行,除了做定时器外,还可以让Python即时执行”rscript”命令调用R脚本来工作,只是这种办法限制太大,只能够交换文件,Python不能对R进行精确的控制。
02. Python直接调用R的函数
R是开源项目,肯定会有一些第三方库实现Python与R互通。
果然,我找到了rpy2,可以实现使用python读取R的对象、调用R的方法以及Python与R数据结构转换等。实际上除了Python,其他语言与R互通的第三方包也大大的有。
最后我选择第2种方法,来让R与Python一起工作。下面开始进行操作讲解。
关于rpy2.robjects是rpy2对R的一个高级封装,该模块里包含了一个R对象和一系列的R数据结构。使用rpy2的大多数情况,只需要跟这个模块打交道即可。rpy2的安装在此不多讲了,直接体验一下R如何与Python无缝整合吧。
了解R实例
R实例是指rpy2.robjects.R,它是在Python中的嵌入式R进程,把R当作从python走向R的通道来看就可以了。通过R实例,我们可以读取R的内置变量、调用R的函数、甚至,直接把它当作R的解析器来用。
访问R的对象
在R的命令行中,我们直接输入对象名来访问R的内置对象,如pi、letters:
在R控制台中访问R对象
>pi
[1] 3.141593
> letters
[1]"a"b""c""d""e""f""g""h""i""j""k""l""m""n""o""p""q"
[18]"r""s""t""u""v""w""x""y""z"
而使用R实例,python访问R对象也很简单,而且方法很多:
在python中访问R对象:
>>> from rpy2 import robjects
>>> r. robjects.r
>>> print r['pi']
[1] 3.141593
>>> print r('pi')
[1] 3.141593
>>> print r.pi
[1] 3.141593
在这段代码中,我们用了三种方式来访问R对象,把R实例当作字典、方法及一个类对象。在实际中,使用哪一种方式要因习惯而异,我喜欢的方法是使用第三种,把R实例当作自己人,直接使用”.”来访问R对象。但这种方法有一个缺陷,就是不能访问带名字空间的R对象或函数,而其他两种方式是可以的,这点将在随后说明。
调用R函数
通过R实例,我们可以轻易地实现用Python调用R的函数。下面我们分别在R控制台和python命令行下读一个数据文件并画一张点图。
R控制台读取文件画点图:
> data = read. table( 'data.csv')
> data
total sms gprs
A 10 4 6
B 9 4 5
C 9 4 5
D 10 5 5
> mtx . data. matrix(data)
> dotchart (mtx)
代码解读:
结果如下:
接下来用python来做一遍同样的事情,我们之前了解到,使用R实例可以直接访问R对象,还可以直接调用R的函数,其实在Python看来,对象和函数是相同的东西,函数也是一种对象罢了。现在来试一下调用”read.table()”函数读入一 个数据文件data.csv:
>>> data r. read. table('data.csv')
Traceback most recent call last):
File "<stdin»", line 1, in cmodule>
File "/Library/python/2.5/site-packages/Il-l.1.7-py2.-macosx-10.6-i386.egg/ init_ ·Py"5.
AttributeError: 'R' object has no attribute ' read'
出错了!怎么回事?在上面我提到过了,使用“.”引用的方式不能访问带有名字空间的R对象和函数,read.table是表示在read包下面的table函数,通过”.”的形式调用失败,必须要用字典的方式或参数的方式来获得:
>>> data = r['read.table']('data.csv')
>>> data
<RDataFrame - Python: 0x4417601 R: 0x2955dc8>
>>> print data
total sms gprs
A 10 4 6
B 9 4 5
C 9 4 5
D 10 5 5
>>> mtx r['data matrix'](data)
>>> r. dotchart(mtx)
这段代码得到的结果与在R控制台下画点图的效果是一样的。最后一行r.dotchart(mtx)是直接通过”.”来调用R的函数dotchart的,在没有名实空间的情况下,是正常的。如果你为了避免太多不可控制的出 错机会,你可以统一地使用字典的方式来访问R对象和方法,这是最保险的方法,虽然我个人认为看起来有点别扭。
R实例就是一R控制台。
其实R实例就是一个可交互的R控制台,只不过交互对象是Python与R罢了,为了证明R实例具有R控制台的特性,来做个实验,写一串R脚本,作为Python一个字符串变量的内容,把该字符串传给R实例,然后把R实例当作方法来调用:
>>> rscript = """
... x = c(1,3,2,5,4,6)
... plot(x)
... """
>>> r(rscript)
<Robject- Python:0x44f990/R:0x8c63b0>
出来的结果这样:
注意:把r实例当作控制台,只能够通过r(r代码)的方式来使用r实例,字典的方式行不通。
加载自定义函数
在实际应用中,使用R语言来编写自己的函数同样是不可避免的,在R控制台中,可以使用source(‘script_path’)的方法来加载自定义R脚本。
在Python中使用自己R脚本中的函数也同样方便:
使用r.source(‘script_path’)即可把自定义函数加载到全局环境 中,再使用R自定义方法名就可以实现调用,我就是这样做的,在此不再详述,同学们自己动手玩一下。
R Vector与Python list
向量(Vector)是R的一个最重要的也是最常用的数据类型,可以理解为一个二维数据,对应Python的list。在R控制台中,声明一个变 量:“x <- 1”,X会被声明成一个向量,而其第一个值是1。
R常常用c()函数来创建一个由多个值组成的向量,例如c(1,2,3,4)。Python要与R打交道,除了访问R对象和调用R函数,还有就是要学会如何转换常见的数据类型。
rpy2提供了几个类,供我们把Python的list转换成R的Vector。分别是 robjects.IntVector,robjects.Boolvector,robjects.StringVector,robjects.FloatVector。
以IntVector为例,将Python的list转换成R的Vector:robjects.IntVector([1,2,3,4,5])。
下面来使用刚学到的类型转换知识画上一个例子的散点图来结束此次体验:
>>> x = [1,3,2,5,4,6]
>>> r.plot(robjeats . Intvector(x))
<RObject - Python:0x44f788/R:0x8c63b0>
写在最后
rpy2提供的不仅仅是上面这些,上面的知识只是rpy2所提供的20%,但是已经足以解决80%的问题。rpy2还提供了更低级的API,你可以做更多的事情,例如你可以实现另一个robjects对象来支持使用”.”来访问带名字空间的对象和函数。
End.
作者:Golonger
来源:简书