#!/usr/bin/env python
import pandas.rpy.common as com
import pandas as pd	
import numpy as np
import rpy2.robjects as robjects
from rpy2.robjects.packages import importr
nlme = importr('nlme')
base = com.importr('base')
stats = com.importr('stats')

mydata = pd.DataFrame.from_csv('http://chymera.eu/data/test/ER_aov.csv', parse_dates=False, index_col=False)
dfr = com.convert_to_r_dataframe(mydata, True)  # convert from pandas to R and make string columns factors

formula = robjects.Formula('ER~COI')
random = robjects.Formula('~1|ID')

test=nlme.lme(formula,random,data=dfr)
test_sum= base.summary(test)

print test_sum