python3调用R语言干货

倾然丶 夕夏残阳落幕 提交于 2020-12-06 19:14:23

 

R语言知识:https://www.w3cschool.cn/r/r_lists.html

 

1. 安装库rpy2

1. 下载与本地对应python版本模块,pip install rpy2是安装不上的

下载地址是:http://www.lfd.uci.edu/~gohlke/pythonlibs/#rpy2  这是python下包的专用地址
需要下载版本和平台都相对应的whl包,我下的是rpy2-2.9.4-cp36-cp36m-win32.whl

pip install rpy2-2.9.4-cp36-cp36m-win32.whl安装即可。

 如果还不行,参考:https://www.cnblogs.com/caiyishuai/p/9520214.html

2. 安装broom --》R语言的一个库--》与R脚本有关,可以忽略

install.packages('broom')

 

3. 写R脚本

library(broom)

test <- function() {
  # x <- c(1:1200000)
  # y <- c(1:1200000)
  x <- c(151, 174, 138, 186, 128, 136, 179, 163, 152, 131)
  y <- c(63, 81, 56, 91, 47, 57, 76, 72, 62, 48)
  
  relation <- lm(y ~ x)
  data <- summary(relation)
  
  data_dict <- c()
  
  newData <- c(data)
  data_dict["residuals"] <- newData["residuals"]
  data_dict["coefficients"] = newData["coefficients"]
  data_dict["aliased"] = newData["aliased"]
  data_dict["sigma"] = newData["sigma"]
  data_dict["df"] = newData["df"]
  data_dict["r.squared"] = newData["r.squared"]
  data_dict["adj.r.squared"] = newData["adj.r.squared"]
  data_dict["fstatistic"] = newData["fstatistic"]
  data_dict["cov.unscaled"] = newData["cov.unscaled"]
  data_dict["p.value"] = c(broom::glance(data))["p.value"]
  
  return(data_dict)
}


# result <- test()
# print(result)

 

4. 写python脚本

报错: RuntimeError: R_USER not defined.

解决方案,各种搜索都是环境变量的问题,于是我各种加

 还tm不行..........................................又懒得重启

stackflow找到答案

os模块的运用,直接看脚本

 

import os
os.environ['R_HOME'] = r'C:\Program Files\R\R-3.6.0'
os.environ['R_USER'] = r'C:\python3.6.3\Lib\site-packages\rpy2' #path depe

import rpy2.robjects as robjects   # ----------------------------------------------> 一定要注意这句,不能放到最上面,因为要先添加环境变量,才能找到这个rpy2。一定要注意
robjects.r.source(r'C:\code\r_test\test_one\test.R')
a = robjects.r('test()')
print(type(a))
# print(list(a))
from pandas import DataFrame
print(a[0])
print(a[0][0])

 

 

打印结果,以及转换数据类型,参考:http://rpy.sourceforge.net/rpy2/doc-2.2/html/vector.html#creating-vectors                  https://blog.csdn.net/suzyu12345/article/details/50587267

 

5. python传值给R脚本,如何实现, 形参方法1

R脚本: 这个脚本的关键在于如何将list转换为c

library(broom)

test <- function(list_data) {
  # print(list_data)
  # print(class(list_data))
  # r语言list 转换成 vector: v = as.vector(unlist(你的list))
  x = c(as.vector(unlist(list_data['x'])))
  y = c(as.vector(unlist(list_data['y'])))

  
  relation <- lm(y ~ x)
  data <- summary(relation)
  print(data)

  return(0)
}

 

python脚本

import os

os.environ['R_HOME'] = r'C:\Program Files\R\R-3.6.0'
os.environ['R_USER'] = r'C:\python3.6.3\Lib\site-packages\rpy2' #path depe

from pandas import DataFrame as df
import rpy2.robjects as robjects
import time
robjects.r.source(r'C:\code\r_test\test_one\test.R')

time1 = time.time()

y =  robjects.ListVector({
    "x":[1, 2, 3],
    "y":[1, 2, 3], # 这里可以给float

    })
a = robjects.r["test"](y)

 

6. python传值给R脚本,如何实现, 形参方法2:类似python的args

R语言脚本

library(broom)

test <- function(...) {
  list_data <- list(...) # 类似python的args,可以传递多个参数
  print(list_data)
  print(class(list_data))
  x = c(as.vector(unlist(list_data[1]))) # 注意R是从1开始的
  y = c(as.vector(unlist(list_data[2])))
  print(x)
  print(y)


  
  relation <- lm(y ~ x)
  data <- summary(relation)
  print(data)

  return(0)
}

python语言

import os

os.environ['R_HOME'] = r'C:\Program Files\R\R-3.6.0'
os.environ['R_USER'] = r'C:\python3.6.3\Lib\site-packages\rpy2' #path depe

from pandas import DataFrame as df
import rpy2.robjects as robjects
import time
robjects.r.source(r'C:\code\r_test\test_one\test.R')

x =  robjects.IntVector([151, 174, 138, 186, 128, 136, 179, 163, 152, 131])
y =  robjects.IntVector([63, 81, 56, 91, 47, 57, 76, 72, 62, 48])

a = robjects.r["test"](x, y)


 

 

 

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!