import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.formula.api import ols

A=[10,12,13,11,10,14,15,13]
B=[9,11,10,12,13]
C=[11,12,15,14,12,13]
All_scores=A+B+C
print(All_scores)

company_name=(['A']*len(A)+['B']*len(B)+['C']*len(C))
print(company_name)

df=pd.DataFrame({'company':company_name, 'output':All_scores})
print(df)

model=ols('output~company',data=df).fit()
anova_table1=sm.stats.anova_lm(model,typ=1)
print(anova_table1)

p_val= 0.212944

if p_val>0.05:
    print('Accept H0')
else:
    print('Reject H0')