%load_ext rmagic

import numpy as np
X = np.array([0,1,2,3,4])
Y = np.array([3,5,4,6,7])

%%R -i X,Y -o XYcoef
XYlm = lm(Y~X)
XYcoef = coef(XYlm)
print(summary(XYlm))
par(mfrow=c(2,2))
plot(XYlm)