# import plotly.graph_objs as go
# from plotly.tools import FigureFactory as FF
import os
import sys
import numpy as np
import pandas as pd
import scipy
import scipy.stats as stats
import statsmodels
import statsmodels.api as sm
from statsmodels.formula.api import ols
import util_ler_dados as udata
data = udata.obterDados2()
arq_destino = '03-Anova-2/_ANOVA_TODAS_ANALISES_2_Fatores_com_Repeticao-[Tam-Prob-Alg].txt'
if os.path.exists(arq_destino):
os.remove(arq_destino)
arq_destino = open(arq_destino, 'w+')
lm = ols(formula='percentual_k_unordered ~ algoritmo * probabilidade_erro * size_of_array', data=data).fit()
anova = sm.stats.anova_lm(lm, typ=2)
tit = ' ANOVA para ALgoritmo X Tamanho X Probabilidade'
hr = '=' * 60 #len(tit)
anov = anova.head(10)
s = '%s\n%s\n%s\n%s\n\n' %(hr,tit,hr,anova.to_string())
arq_destino.write(s)
# import plotly.plotly as py
# import plotly.graph_objs as go
# from plotly.tools import FigureFactory as FF
import os
import sys
import numpy as np
import pandas as pd
import scipy
import scipy.stats as stats
import statsmodels
import statsmodels.api as sm
from statsmodels.formula.api import ols
import util_ler_dados as udata
df_tudo = udata.obterDados2()
arq_destino = '03-Anova-2/_ANOVA_TODAS_ANALISES_2_Fatores_com_Repeticao-[Prob-Alg].txt'
if os.path.exists(arq_destino):
os.remove(arq_destino)
arq_destino = open(arq_destino, 'w+')
for tam in udata.TAMANHOS:
data = udata.obterDfPorTam2(tam=tam, df=df_tudo)
# print(data.head())
lm = ols(formula='percentual_k_unordered ~ algoritmo * probabilidade_erro',
data=data).fit()
anova = sm.stats.anova_lm(lm, typ=2) # Type 2 ANOVA DataFrame
tit = ' ANOVA para Tamanho = %s ' % (tam)
