#!/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
Error: ID Df Sum Sq Mean Sq F value Pr(>F) Residuals 6 0.022 0.003666 Error: Within Df Sum Sq Mean Sq F value Pr(>F) COI 6 0.02628 0.004379 3.468 0.0083 ** Residuals 36 0.04547 0.001263 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1