CRD with unequal
Tue Apr 16 2024 10:26:29 GMT+0000 (Coordinated Universal Time)
Saved by
@mrv
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')
content_copyCOPY
Comments