The following example illustrates how to use R package with Python by using the Rpy2 package.

This problem solve in this example is to fit a normal distribution on some data samples.

from rpy2.robjects import r from rpy2.robjects.packages import importr rnorm = r('rnorm') # create an alias to the R function rnorm x = rnorm(n=1000, mean=0, sd=1) # generates the data sample (normal distribution) # load the MASS library for distribution fitting MASS = importr('MASS') params = MASS.fitdistr(x, 'normal') # perform the fitting here print params |

The result is stored in the variable **params**. The above code print the following results:

mean sd -0.07962739 1.00467570 ( 0.03177064) ( 0.02246523) |

Telling you that the distribution has mean of -0.08 and a standard deviation of 0. The bracketted numbers indicates how trustful are the results. Numbers in parentheses are confidence intervals around the parameter.

Please follow and like us: