Python/patsy
R style statistical model description. Perfect for pandas + statsmodels.
Example:
from patsy import dmatrices
import statsmodels.api as sm
import pandas as pd
df = pd.read_csv('data.csv')
formula = 'np.log(count) ~ var1 + var2 + var3'
y, X = dmatrices(formula, data=df, return_type='dataframe')
mod = sm.GLM(y, X, family=sm.families.NegativeBinomial())
res = mod.fit()
print(res.summary())