rmagic.py

import numpy as np
x = np.random.randn(100)
y = 2*x + np.random.randn(100)

%load_ext rpy2.ipython

%%R -i x,y -w 500 -h 300
df <- data.frame(x,y)
m <- lm(y~x)
inter <- m$coef[1]
slop <- m$coef[2]

library(ggplot2)
p <- ggplot(df, aes(x,y))
p <- p + 
    geom_point() +
    geom_abline(intercept = inter, slope = slop)
print(p)