First time using the RPy2 implementation in Python. Attempting to do an one-way ANOVA with two factors. It works in R on another machine, but Python does not like the syntax. Any thoughts are appreciated!
from rpy2.robjects import aov
huanova = aov(formula = df1['human_den'] ~ df1['region']+df1['years'])
Error message points at the tilda.
huanova = aov(formula=df1['human_den'] ~ df1['region']+df1['years'])
^
SyntaxError: invalid syntax
As per the documentation about Formulae in Rpy2, you have to pass the formula as a string. This is one way of doing it:
from rpy2.robjects import aov
from rpy2.robjects import Formula
formula = Formula('human_den ~ region + years')
env = formula.environment
env['human_den'] = df1['human_den']
env['region'] = df1['region']
env['years'] = df1['years']
huanova = aov(formula = formula)
The above answer from musically_ut will work but you need load the 'stats' package that contains the aov function using importr, as shown in this question:
from rpy2.robjects import Formula
from rpy2.robjects.packages import importr
stats = importr('stats')
formula = Formula('human_den ~ region + years')
env = formula.environment
env['human_den'] = df1['human_den']
env['region'] = df1['region']
env['years'] = df1['years']
huanova = stats.aov(formula = formula)
Tilde ~ is Unary in python, but you are using it as Binary. You may want:
huanova = aov(formula = df1['human_den'] + ~ df1['region']+df1['years'])
Notice that I've added a PLUS before Tilde.
