Low-level interface

Overview

The package rinterface is provided as a lower-level interface, for situations where either the use-cases addressed by robjects are not covered, or for the cases where the layer in robjects has an excessive cost in terms of performance.

The package can be imported with:

>>> import rpy2.rinterface as rinterface

Initialization

One has to initialize R before much can be done. The function initr() lets one initialize the embedded R.

This is done with the function initr().

rpy2.rinterface.initr()

Initialize the embedded R.

Parameters:

interactive – Should R run in interactive or non-interactive mode?

if None the value in _DEFAULT_R_INTERACTIVE will be used. :param _want_setcallbacks: Should custom rpy2 callbacks for R frontends be set?. :param _c_stack_limit: Limit for the C Stack. if None the value in _DEFAULT_C_STACK_LIMIT will be used.

>>> rinterface.initr()

Initialization should only be performed once. To avoid unpredictable results when using the embedded R, subsequent calls to initr() will not have any effect.

The functions rpy2.rinterface_lib.embedded.get_initoptions() and rpy2.rinterface_lib.embedded.set_initoptions() can be used to modify the options. Default parameters for the initialization are otherwise in the module variable _options.

Note

If calling initr() returns an error stating that R_HOME is not defined, you should either have the R executable in your path (PATH on unix-alikes, or Path on Microsoft Windows) or have the environment variable R_HOME defined.

Should the initialization fail, a mismatch between the version of the R rpy2 was compiled against and the R rpy2 is run with should be investigated. The variable rpy2.rinterface.R_BUILD_VERSION contains information about the R version rpy2 was built against. rpy2 is relatively independent of R versions, but changes in the R C API might cause problems.

Ending R

Ending the R process is possible, but starting it again with initr() does appear to lead to an R process that is hardly usable. For that reason, the use of endr() should be considered carefully, if at all.

rpy2.rinterface.endr()[source]

Note

When writing a GUI for R, a developper may want to either prevent a user to call R quit(), or ensure that specific code is executed before terminating R (for example a confirmation dialog window “do you really want to terminate ?”). This can be done by replacing the callback cleanup with an appropriate function (see Clean up).

R space and Python space

When using the RPy2 package, two realms are co-existing: R and Python.

rpy2.rinterface_lib.sexp.Sexp objects can be considered as Python envelopes pointing to data stored and administered in the R space.

R variables exist within an embedded R workspace, and can be accessed from Python through their python object representations.

We distinguish two kinds of R objects: named objects and anonymous objects. Named objects have an associated symbol in the R workspace (a “variable name”) while anonymous objects don’t, but are protected from garbage collection on the R side for as long as they are used on the Python side.

Named objects

For example, the following R code is creating two objects, named x and hyp respectively, in the global environment. Those two objects could be accessed from Python using their names.

x <- c(1,2,3)

hyp <- function(x, y) sqrt(x^2 + y^2)

By default R starts with two environments: baseenv and globalenv. Both are instances of class rpy2.rinterface.SexpEnvironment in rpy2.

globalenv

The global environment (globalenv) can be seen as the root (or topmost) environment, and is in fact a list, that is a sequence, of environments.

When an R library (package in R’s terminology) is loaded, is it added to the existing sequence of environments. Unless specified, it is inserted in second position. The first position being always attributed to the global environment. The library is said to be attached to the current search path.

baseenv

The base package has a namespace (baseenv), that can be accessed as an environment.

Note

Depending on what is in globalenv and on the attached packages, base objects can be masked when starting the search from globalenv. Use baseenv when you want to be sure to access a function you know to be in the base namespace.

Anonymous objects

Anonymous R objects do not have an associated symbol, yet are protected from garbage collection.

Such objects can be created when using the constructor for an Sexp* class.

For example:

>>> x = rinterface.IntVector((1,2,3))

creates a fully usable R vector, but it does not have an associtated R symbol (it is in memory, but cannot be called by name fomr R). It is also protected from garbage collection until, until x is deleted and the Python garbage collector destroys x.

Note

To finalize the recovery of the memory used, the R garbage collector must be also be called. This should happen automatically while running R code when a threshold of memory usage is reached, but it is also possible to call explicitly both garbage collectors. See Memory management and garbage collection for more details.

Pass-by-value paradigm

The root of the R language is functional, with arguments passed by value. R is actually using tricks to lower memory usage, such as only copying an object when needed (that is when the object is modified in a local block), but copies of objects are nevertheless frequent. This can remain unnoticed by a user until large objects are in use or a large number of modification of objects are performed, in which case performance issues may appear. An infamous example is when the column names for a matrix are changed, bringing a system to its knees when the matrix is very large, as the whole matrix ends up being copied.

On the contrary, Python is using pointer objects passed around through function calls, and since rpy2, is a Python-to-R interface the Python approach was conserved.

Although being contrived, the example below will illustrate the point. With R, renaming a column is like:

# create a matrix
m <- matrix(1:10, nrow = 2,
            dimnames = list(c("1", "2"),
                            c("a", "b", "c", "d", "e")))
# rename the third column
i <- 3
colnames(m)[i] <- "foo"

With rpy2.rinterface:

# import and initialize
import rpy2.rinterface as ri
ri.initr()

# make a function to rename column i
def rename_col_i(m, i, name):
    m.do_slot("dimnames")[1][i] = name

# create a matrix
matrix = ri.baseenv["matrix"]
rlist = ri.baseenv["list"]
m = matrix(ri.baseenv[":"](1, 10),
           nrow = 2,
           dimnames = rlist(ri.StrSexpVector(("1", "2")),
                            ri.StrSexpVector(("a", "b", "c", "d", "e"))))

Now we can check that the column names

>>> tuple(m.do_slot("dimnames")[1])
('a', 'b', 'c', 'd', 'e')

And rename the third column (remembering that R vectors are 1-indexed while Python sequences are 0-indexed).

>>> i = 3-1
>>> rename_col_i(m, i, ri.StrSexpVector(("foo", )))
>>> tuple(m.do_slot("dimnames")[1])
('a', 'b', 'foo', 'd', 'e')

Unlike with the R code, neither the matrix or the vector with the column names are copied. Whenever this is not a good thing, R objects can be copied the way Python objects are usually copied (using copy.deepcopy(), Sexp implements Sexp.__deepcopy__()).

Parsing and evaluating R code

The R C-level function for parsing an arbitrary string a R code is exposed as the function parse()

>>> expression = ri.parse('1 + 2')

The resulting expression is a nested list of R statements.

>>> len(expression)
1
>>> len(expression[0])
3

The R code 1 + 2 translates to an expression of length 3: +(1, 2), that is a call to the function + (or rather the symbol associated with the function) with the arguments 1 and 2.

>>> expression[0][0].typeof
<RTYPES.SYMSXP: 1>
>>> tuple(expression[0][1])
(1.0,)
>>> tuple(expression[0][2])
(2.0,)

Note

The expression must be evaluated if the result from its execution is wanted.

rpy2.rinterface.parse()

Calling Python functions from R

As could be expected from R’s functional roots, functions are first-class objects. This means that the use of callback functions as passed as parameters is not seldom, and this also means that the Python programmer has to either be able write R code for functions as arguments, or have a way to pass Python functions to R as genuine R functions. That last option is becoming possible, in other words one can write a Python function and expose it to R in such a way that the embedded R engine can use as a regular R function.

As an example, let’s consider the R function optim() that looks for optimal parameters for a given cost function. The cost function should be passed in the call to optim() as it will be repeatedly called as the parameter space is explored, and only Python coding skills are necessary as the code below demonstrates it.

from rpy2.robjects.vectors import FloatVector
from rpy2.robjects.packages import importr
import rpy2.rinterface as ri
stats = importr('stats')

# cost function, callable from R
@ri.rternalize
def cost_f(x):
    # Rosenbrock Banana function as a cost function
    # (as in the R man page for optim())
    x1, x2 = x
    return 100 * (x2 - x1 * x1)**2 + (1 - x1)**2

# starting parameters
start_params = FloatVector((-1.2, 1))

# call R's optim() with our cost funtion
res = stats.optim(start_params, cost_f)

For convenience, the code example uses the higher-level interface robjects whenever possible.

The lower-level function rternalize() will take an arbitray Python function and return an rinterface.SexpClosure instance, that is a object that can be used by R as a function.

rpy2.rinterface.rternalize()[source]

Make a Python function callable from R.

Takes an arbitrary Python function and wrap it in such a way that it can be called from the R side.

This factory can be used as a decorator, and has an optional named argument “signature” that can be True or False (default is False). When True, the R function wrapping the Python one will have a matching signature or a close one, as detailed below.

The Python ellipsis arguments`*args` and **kwargs are translated to the R ellipsis .

For example:

will be visible in R with the signature:

The only limitation is that whenever *args and **kwargs are both present in the Python declaration they must be consecutive. For example:

is a valid Python declaration. However, it cannot be “rternalized” because the R ellipsis can only appear at most once in the signature of a given function. Trying to apply the decorator rternalize would raise an exception.

The following Python function can be “rternalized”:

It is visible to R with the signature

Python function definitions can allow the optional naming of required arguments. The mapping of signatures between Python and R is then quasi-indentical since R does it for unnamed arguments. The check that all arguments are present is still performed on the Python side.

Example:

>>> _ = foo(1, 2, 3)
x: 1, y: 2, z: 3
ValueErro: None
>>> _ = foo(1)
TypeError: rternalized foo() missing 2 required keyword-only arguments: 'y' and 'z'
>>> _ = foo(1, z=2, y=3)
x: 1, y: 3, z: 2
>>> _ = foo(1, z=2, y=3)
x: 1, y: 3, z: 2

Note that whenever the Python function has an ellipsis (either *args or **kwargs) earlier parameters in the signature that are positional-or-keyword are considered to be positional arguments in a function call.

Parameters:

function – A Python callable object. This is a positional

argument with a default value None to allow the decorator function without parentheses when optional argument is not wanted.

Returns:

A wrapped R object that can be use like any other rpy2

object.

Interactive features

The embedded R started from rpy2 is interactive by default, which means that a number of interactive features present when working in an interactive R console will be available for use.

Such features can be called explicitly by the rpy2 user, but can also be triggered indirectly, as some on the R functions will behave differently when run interactively compared to when run in the so-called BATCH mode.

Note

However, interactive use may mean the ability to periodically check and process events. This is for example the case with interactive graphics devices or with the HTML-based help system (see Processing interactive events).

I/O with the R console

See Console I/O.

Processing interactive events

An interactive R session is can start interactive activities that require a continuous monitoring for events. A typical example is the interactive graphical devices (plotting windows), as they can be resized and the information they display is refreshed.

However, to do so the R process must be instructed to process pending interactive events. This is done by the R console for example, but rpy2 is designed as a library rather than as a threaded R process running within Python (yet this can be done as shown below).

The way to restore interactivity is to simply call the function rinterface_lib.callbacks.process_revents() at regular intervals.

A higher-level interface is available, running the processing of R events in a thread (see Section R event loop).

Multithreading

R is quite not friendly to multithreading, and trying to play with threads at the C level can quickly result in an embedded R crashing. Since we are using R’s C API and Python can do multithreading, getting to such software failure will not be too hard. However, multithreading should not be considered impossible, or even very difficult to achieve when applying the one guideline below.

rpy2 has a lock that can be used as a context manager. Interactions with R that a code author knows should never be interrupted by thread switching can simply be wrapped in a thread-locked block as follows:

from rpy2.rinterface_lib import openrlib

with openrlib.rlock:
    # (put interactions with R that should not be interrupted by
    # thread switching here).
    pass

That lock is already used in a handful of critical low-level accesses to the R API in the rpy2 code base (e.g., when a protection/unprotection stack is used for R objects transiently protected from garbage collection, or when the embedded R is initialized) but can be safely reused and nested in higher level code.

Note

Web Server Gateway Interfaces (WSGIs) for Python scripts can use multithreading to optimize resources, allowing one process to handle several connections as they are presumably of higher latency than what happens on the server side.

When using rpy2 to build services that run R code, attention should be paid to whether threads are used, and if the case the lock mentioned in this section should be used to ensure that coherent results results are computed by R even in the presence of multithreading.

Classes

Sexp

The class Sexp is the base class for all R objects.

class rpy2.rinterface.Sexp[source]
__sexp__

Python C capsule wrapping the pointer to the underlying R object (SEXP)

named

R does not count references for its object. This method returns the NAMED value (an integer). See the R-extensions manual for further details.

typeof

Internal R type for the underlying R object

>>> letters.typeof
<RTYPES.STRSXP: 16>
__deepcopy__(self)

Make a deep copy of the object, calling the R-API C function c:function::Rf_duplicate() for copying the R object wrapped.

New in version 2.0.3.

do_slot(name)

R objects can be given attributes. In R, the function attr lets one access an object’s attribute; it is called do_slot() in the C interface to R.

Parameters:

name – string

Return type:

instance of Sexp

>>> matrix = rinterface.globalenv.find("matrix")
>>> letters = rinterface.globalenv.find("letters")
>>> m = matrix(letters, ncol = 2)
>>> [x for x in m.do_slot("dim")]
[13, 2]
>>>
do_slot_assign(name, value)[source]

Assign value to the slot with the given name, creating the slot whenver not already existing.

Parameters:
  • name – string

  • value – instance of Sexp

rsame(sexp_obj)[source]

Tell whether the underlying R object for sexp_obj is the same or not.

Return type:

boolean

Underlying R object

The underlying R object is a pointer to a c:type::SEXPREC as defined in R’s Rinternals.h. That object is wrapped in a c:type::SexpObj and placed in a Python capsule.

The capsule is providing a relatively safe mechanism to exchange underlying R objects between rpy2 objects.

from rpy2.rinterface import SexpIntVector, SexpFloatVector
vector1=SexpIntVector((1, 2, 3))
vector2=SexpFloatVector((4.0, 5.0, 6.0))

vector1.__sexp__ = vector2.__sexp_

R arrays (vectors) inherit from SexpVector

Overview

In R there are no scalars, only arrays (called “vectors” when unidimensional). Anything like a one-value variable is in fact a vector of length 1.

To use again the constant pi:

>>> pi = rinterface.globalenv.find('pi')
>>> len(pi)
1
>>> pi
<rinterface.FloatSexpVector - Python:0x2b20325d2660 / R:0x16d5248>
>>> pi[0]
3.1415926535897931
>>>

The letters of the (western) alphabet are:

>>> letters = rinterface.globalenv.find('letters')
>>> len(letters)
26
>>> LETTERS = rinterface.globalenv.find('LETTERS')

R types

R vectors all have a type, sometimes referred to in R as a mode. Rpy2 has chosen to map R types to child classes of rpy2.rinterface_lib.sexp.SexpVector, as shown in the inheritance diagram below, with the numpy array interface implemented for some of them.

Inheritance diagram of rpy2.rinterface.SexpVector, rpy2.rinterface.IntSexpVector, rpy2.rinterface.FloatSexpVector, rpy2.rinterface.ByteSexpVector, rpy2.rinterface.ComplexSexpVector, rpy2.rinterface.StrSexpVector, rpy2.rinterface.ListSexpVector, rpy2.rinterface.PairlistSexpVector, rpy2.rinterface.ExprSexpVector, rpy2.rinterface.LangSexpVector, rpy2.rinterface.BoolSexpVector

C-level R array objects

Indexing

The indexing is working like it would on regular Python tuples or lists. The indexing starts at 0 (zero), which differs from R, where indexing start at 1 (one).

Note

The __getitem__ operator [ is returning a Python scalar. Casting an SexpVector into a list is only a matter of either iterating through it, or simply calling the constructor list().

Common attributes

Names

In R, vectors can be named, that is each value in the vector can be given a name (that is be associated a string). The names are added to the other as an attribute (conveniently called names), and can be accessed as such:

>>> options = rinterface.globalenv.find("options")()
>>> option_names = options.do_slot("names")
>>> [x for x in options_names]

Note

Elements in a name vector do not have to be unique. A Python counterpart is provided as rpy2.rlike.container.TaggedList.

Dim and dimnames

In the case of an array, the names across the respective dimensions of the object are accessible through the slot named dimnames.

Missing values

Note

R also has the notion of missing parameters in function calls. This is a separate concept, and more information about are given in Section Functions.

In R missing the symbol NA represents a missing value. The general rule that R scalars are in fact vectors applies here again, and the following R code is creating a vector of length 1.

x <- NA

The type of NA is logical (boolean), and one can specify a different type with the symbols NA_character_, NA_integer_, NA_real_, and NA_complex_.

In rpy2.rinterface_lib.na_values, the symbols can be accessed by through NACharacter, NAInteger, NAReal.

Those are singleton instance from respective NA<something>Type classes.

>>> my_naint = rinterface_lib.na_values.NAIntegerType()
>>> my_naint is rinterface_lib.na_values.NA_Integer
True
>>> my_naint == rinterface_lib.na_values.NA_Integer
True

NA values can be present in vectors returned by R functions.

>>> rinterface.baseenv['as.integer'](rinterface.StrSexpVector(("foo",)))[0]
NA_integer_

NA values can have operators implemented, but the results will then be missing values.

>>> rinterface.NA_Integer + 1
NA_integer_
>>> rinterface.NA_Integer * 10
NA_integer_

Warning

Python functions relying on C-level implementations might not be following the same rule for NAs.

>>> x = rinterface.IntSexpVector((1, rinterface.NA_Integer, 2))
>>> sum(x)
3
>>> max(x)
2
>>> min(x)
NA_integer_

This should be preferred way to use R’s NA as those symbol are little peculiar and cannot be retrieved with SexpEnvironment.find().

Those missing values can also be used with the rpy2.robjects layer and more documentation about their usage can be found there (see Missing values).

class rpy2.rinterface_lib.sexp.NAIntegerType[source]

Bases: int

class rpy2.rinterface_lib.sexp.NARealType[source]

Bases: float

class rpy2.rinterface_lib.sexp.NALogicalType[source]

Bases: int

class rpy2.rinterface_lib.sexp.NACharacterType[source]

Bases: CharSexp

class rpy2.rinterface_lib.sexp.NAComplexType[source]

Bases: complex

Constructors

class rpy2.rinterface.SexpVector(obj, sexptype, copy)[source]

Bases: Sexp, SexpVectorAbstract

Base abstract class for R vector objects.

R vector objects are, at the C level, essentially C arrays wrapped in the general structure for R objects.

Convenience classes are provided to create vectors of a given type:

class rpy2.rinterface.StrSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVector

R vector of strings.

get_charsxp(i: int) CharSexp[source]

Get the R CharSexp objects for the index i.

class rpy2.rinterface.IntSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVectorWithNumpyInterface

class rpy2.rinterface.ByteSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVectorWithNumpyInterface

Array of bytes.

This is the R equivalent to a Python bytesarray.

class rpy2.rinterface.FloatSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVectorWithNumpyInterface

class rpy2.rinterface.BoolSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVectorWithNumpyInterface

Array of booleans.

Note that R is internally storing booleans as integers to allow an additional “NA” value to represent missingness.

class rpy2.rinterface.ListSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVector

R list.

An R list an R vector (array) that is similar to a Python list in the sense that items in the list can be of any type, whereas most other R vectors are homogeneous (all items are of the same type).

class rpy2.rinterface.PairlistSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVector

R pairlist.

A R pairlist is rarely used outside of R’s internal libraries and a relatively small number of use cases. It is essentially a LISP-like list of (name, value) pairs.

class rpy2.rinterface.ComplexSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVector

class rpy2.rinterface.LangSexpVector(obj: SupportsSEXP | SexpCapsule | Sized)[source]

Bases: SexpVector

An R language object.

To create from a (Python) string containing R code use the classmethod from_string.

classmethod from_string(s: str) LangSexpVector_VT[source]

Create an R language object from a string.

This creates an unevaluated R language object.

Args:

s: R source code in a string.

Returns:

An instance of the class the method is called from.

SexpEnvironment

__getitem__() / __setitem__() / __delitem__()

The [ operator will only look for a symbol in the environment without looking further in the path of enclosing environments.

The following will return an exception LookupError:

>>> rinterface.globalenv["pi"]

The constant pi is defined in R’s base package, and therefore cannot be found in the Global Environment.

The assignment of a value to a symbol in an environment is as simple as assigning a value to a key in a Python dictionary:

>>> x = rinterface.IntSexpVector([123, ])
>>> rinterface.globalenv["x"] = x
>>> len(x)
1
>>> tuple(rinterface.globalenv)
('x', )

Removing an element can be done like one would do it for a Python dict:

>>> del(rinterface.globalenv['x'])
>>> len(x)
0

Note

Not all R environment are hash tables, and this may influence performances when doing repeated lookups.

Note

A copy of the R object is made in the R space.

__iter__()

The object is made iter-able.

For example, we take the base name space (that is the environment that contains R’s base objects:

>>> base = rinterface.baseenv
>>> basetypes = [x.typeof for x in base]

Warning

In the current implementation the content of the environment is evaluated only once, when the iterator is created. Adding or removing elements to the environment will not update the iterator (this is a problem, that will be solved in the near future).

find()

Whenever a search for a symbol is performed, the whole search path is considered: the environments in the list are inspected in sequence and the value for the first symbol found matching is returned.

Let’s start with an example:

>>> rinterface.globalenv.find("pi")[0]
3.1415926535897931

The constant pi is defined in the package base, that is always in the search path (and in the last position, as it is attached first). The call to get() will look for pi first in globalenv, then in the next environment in the search path and repeat this until an object is found or the sequence of environments to explore is exhausted.

We know that pi is in the base namespace and we could have gotten here directly from there:

>>> ri.baseenv.find('pi')[0]
3.1415926535897931
>>> ri.baseenv['pi'][0]
3.1415926535897931
>>> ri.globalenv["pi"][0]
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
LookupError: 'pi' not found

R can look specifically for functions, which is happening when a parsed function call is evaluated. The following example of an R interactive session should demonstrate it:

> mydate <- "hohoho"
> mydate()
Error: could not find function "mydate"
>
> date <- "hohoho"
> date()
[1] "Sat Aug  9 15:27:40 2008"

The base function date is still found, although a non-function object is present earlier on the search path.

The same behavior can be obtained from rpy2 with the optional parameter wantfun (specify that get() should return an R function).

>>> ri.globalenv['date'] = ri.StrSexpVector(['hohoho', ])
>>> ri.globalenv.find('date')[0]
'hohoho'
>>> ri.globalenv.find('date', wantfun=True)
<rinterface.SexpClosure - Python:0x7f142aa96198 / R:0x16e9500>
>>> date = ri.globalenv.find('date', wantfun=True)
>>> date()[0]
'Sat Aug  9 15:48:42 2008'

R packages as environments

In a Python programmer’s perspective, it would be nice to map loaded R packages as modules and provide access to R objects in packages the same way than Python object in modules are accessed.

This is unfortunately not possible in a completely robust way: the dot character . can be used for symbol names in R (like pretty much any character), and this can make an exact correspondance between R and Python names rather difficult. rpy uses transformation functions that translates ‘.’ to ‘_’ and back, but this can lead to complications since ‘_’ can also be used for R symbols (although this is the approach taken for the high-level interface, see Section R packages).

There is a way to provide explict access to object in R packages, since loaded packages can be considered as environments. To make it convenient to use, one can consider making a function such as the one below:

def rimport(packname):
    """ import an R package and return its environment """
    as_environment = rinterface.baseenv['as.environment']
    require = rinterface.baseenv['require']
    require(rinterface.StrSexpVector(packname),
            quiet = rinterface.BoolSexpVector((True, )))
    packname = rinterface.StrSexpVector(('package:' + str(packname)))
    pack_env = as_environment(packname)
    return pack_env
>>> class_env = rimport("class")
>>> class_env['knn']

For example, we can reimplement in Python the R function returning the search path (search).

def rsearch():
    """ Return a list of package environments corresponding to the
    R search path. """
    spath = [ri.globalenv, ]
    item = ri.globalenv.enclos()
    while not item.rsame(ri.emptyenv):
        spath.append(item)
        item = item.enclos()
    spath.append(ri.baseenv)
    return spath

As an other example, one can implement simply a function that returns from which environment an object called by get() comes from.

def wherefrom(name, startenv=ri.globalenv):
    """ when calling 'get', where the R object is coming from. """
    env = startenv
    obj = None
    retry = True
    while retry:
        try:
            obj = env[name]
            retry = False
        except LookupError, knf:
            env = env.enclos()
            if env.rsame(ri.emptyenv):
                retry = False
            else:
                retry = True
    return env
>>> wherefrom('plot').do_slot('name')[0]
'package:graphics'
>>> wherefrom('help').do_slot('name')[0]
'package:utils'

Note

Unfortunately this does not generalize to all cases: the base package does not have a name.

>>> wherefrom('get').do_slot('name')[0]
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
LookupError: The object has no such attribute.

Functions

A function with a context

In R terminology, a closure is a function (with its enclosing environment). That enclosing environment can be thought of as a context to the function.

Note

Technically, the class SexpClosure corresponds to the R types CLOSXP, BUILTINSXP, and SPECIALSXP, with only the first one (CLOSXP) being a closure.

>>> sum = rinterface.globalenv.find('sum')
>>> x = rinterface.IntSexpVector([1,2,3])
>>> s = sum(x)
>>> s[0]
6

Named arguments

Named arguments to an R function can be specified just the way they would be with any other regular Python function.

>>> rnorm = rinterface.globalenv.find('rnorm')
>>> rnorm(rinterface.IntSexpVector([1, ]),
          mean = rinterface.IntSexpVector([2, ]))[0]
0.32796768001636134

There are however frequent names for R parameters causing problems: all the names with a dot. using such parameters for an R function will either require to:

  • use the special syntax **kwargs on a dictionary with the named parameters

  • use the method rcall().

Order for named parameters

One point where function calls in R can differ from the ones in Python is that all parameters in R are passed in the order they are in the call (no matter whether the parameter is named or not), while in Python only parameters without a name are passed in order. Using the class OrdDict in the module rpy2.rlike.container, together with the method rcall(), permits calling a function the same way it would in R. For example:

import rpy2.rlike.container as rpc
args = rpc.OrdDict()
args['x'] = rinterface.IntSexpVector([1,2,3])
args[None] = rinterface.IntSexpVector([4,5])
args['y'] = rinterface.IntSexpVector([6, ])
rlist = rinterface.baseenv['list']
rl = rlist.rcall(tuple(args.items()), rinterface.globalenv)
>>> [x for x in rl.do_slot("names")]
['x', '', 'y']

closureEnv

In the example below, we inspect the environment for the function plot, that is the namespace for the package graphics.

>>> plot = rinterface.globalenv.find('plot')
>>> ls = rinterface.globalenv.find('ls')
>>> envplot_list = ls(plot.closureEnv())
>>> [x for x in envplot_ls]
>>>

Missing parameters

In R function calls can contain explicitely missing parameters.

> sum(1,,3)
Error: element 2 is empty;
   the part of the args list of 'sum' being evaluated was:
   (1, , 3)

This is used when extracting a subset of an array, with a missing parameter interpreted by the extract function [ like all elements across that dimension must be taken.

m <- matrix(1:10, nrow = 5, ncol = 2)

# extract the second column
n <- m[, 2]

# can also be written
n <- "["(m, , 2)

rinterface.MissingArg is a pointer to the singleton rinterface.MissingArgType, allowing to explicitly pass missing parameters to a function call.

For example, the extraction of the second column of a matrix with R shown above, will write almost identically in rpy2.

import rpy2.rinterface as ri
ri.initr()

matrix = ri.baseenv['matrix']
extract = ri.baseenv['[']

m = matrix(ri.IntSexpVector(range(1, 11)), nrow = 5, ncol = 2)

n = extract(m, ri.MissingArg, 2)

SexpS4

Object-Oriented programming in R exists in several flavours, and one of those is called S4. It has its own type at R’s C-API level, and because of that specificity we defined a class. Beside that, the class does not provide much specific features (see the pydoc for the class below).

An instance’s attributes can be accessed through the parent class Sexp method do_slot().

class rpy2.rinterface.SexpS4(obj)[source]

Bases: Sexp

R “S4” object.

S4 objects as one of the available forms of Object-Oriented Programming in R.

SexpExtPtr

External pointers are intended to facilitate the handling of C or C++ data structures from R. In few words they are pointers to structures external to R. They have been used to implement vectors and arrays in shared memory, or storage-based vectors and arrays.

External pointers also do not obey the pass-by-value rule and can represent a way to implement pointers in R.

Let us consider the following simple example:

ep = rinterface.SexpExtPtr.from_pyobject('hohoho')

The Python string is now encapsulated into an R external pointer, and visible as such by the embedded R process.

When thinking of sharing C-level structures between R and Python more involved examples can be considered (here still a simple example):

import ctypes
class Point2D(ctypes.Structure):
    _fields_ = [('x', ctypes.c_int),
                ('y', ctypes.c_int)]

pt = Point2D()

ep = rinterface.SexpExtPtr.from_pyobject(pt)

However, this remains a rather academic exercise unless there exists a way to access the data from R; when used in R packages, external pointers have companion functions to manipulate the C-level data structures.

In the case of external pointers and their companion functions and methods defined by R packages, the rpy2 interface lets a programmer create such external pointers directly from Python, using ctypes for example.

However, the rpy2 interface allows more than that since a programmer is able to make a Python function accessible to R has is was a function of its own. It is possible to define arbitrary Python data structures as well as functions or methods to operate on them, pass the data structure to R as an external pointer, and expose the functions and methods to R.

class rpy2.rinterface.SexpExtPtr(obj)[source]

Bases: Sexp

R “External Pointer” object.

Class diagram

Inheritance diagram of rpy2.rinterface, rpy2.rinterface_lib.sexp, rpy2.rinterface_lib.na_values

Misc. variables

R_HOME

R HOME

R_LEN_T_MAX

largest usable integer for indexing R vectors

TRUE/FALSE

R’s TRUE and FALSE

R types

Vector types

CPLXSXP

Complex

INTSXP

Integer.

LGLSXP

Boolean (logical in the R terminology)

RAWSXP

Raw (bytes) value

REALSXP

Numerical value (float / double)

STRSXP

String

VECSXP

List

LISTSXP

Paired list

LANGSXP

Language object.

EXPRSXP

Unevaluated expression.

Function types

CLOSXP

Function with an enclosure. Represented by rpy2.rinterface.SexpClosure.

BUILTINSXP

Base function

SPECIALSXP

Some other kind of function

Other types

ENVSXP

Environment. Represented by rpy2.rinterface.SexpEnvironment.

S4SXP

Instance of class S4. Represented by rpy2.rinterface.SexpS4.

Types one should not meet

PROMSXP

Promise.