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

mydata = pd.DataFrame.from_csv('http://chymera.eu/data/test/ER_nlme.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

formula = robjects.Formula('ER~COI+Error(ID)')
test = stats.aov(formula, data=dfr)
test_sum= base.summary(test)
print test_sum