import pandas as pd
import numpy as np

movies = pd.read_csv('http://www.rossmanchance.com/iscam2/data/movies03RT.txt', sep='\t')
movies.columns = ['X', 'score', 'rating', 'genre', 'box_office', 'running_time']
movies.head()

from statsmodels.stats.anova import anova_lm
from statsmodels.formula.api import ols

lm = ols('score ~ rating', movies).fit()

aovObject = anova_lm(lm)
aovObject

lm.params

lm2 = ols('score ~ rating + genre', movies).fit()

aovObject2 = anova_lm(lm2)
aovObject2

lm3 = ols('score ~ genre + rating', movies).fit()

aovObject3 = anova_lm(lm3)
aovObject3

lm4 = ols('score ~ genre + rating + box_office', movies).fit()

aovObject4 = anova_lm(lm4)
aovObject4