Ecogwiki

Usage examples

Code

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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())

Math

Inline:

Equation: