#!/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
base = com.importr('base')
stats = com.importr('stats')
xtable = importr('xtable')
stargazer = importr('stargazer')
texreg = importr('texreg')

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+Error(ID)')

test=stats.aov(formula,data=dfr)

test_sum= base.summary(test)

print test_sum