#Desabilite a linha abaixo caso o pacote ainda não foi instalado no sistema
#install.packages("FrF2", repos = "http://cran.rstudio.com/") # Instala o pacote

library(FrF2) # Carrega o pacote para uso


#exemplo do slide (retirado do livro do Jain)
plan.person = FrF2(nruns = 4,
                   nfactors =  2,
                   replications = 1,
                   repeat.only = FALSE,
                   factor.names = list(
                     Rede = c("crossbar", "omega"),
                     Acesso = c("Aleatorio", "Matriz")),
                   randomize = FALSE)

summary(plan.person)

# Variável de resposta: Throughput
resultados = c(0.6041, 0.4220, 0.7922, 0.4717)

# comente a linha acima e tire o comentário da linha abaixo caso a variável de resposta a ser analisada seja o número de ciclos para transmissão
resultados = c(3, 5, 2, 4)


plan.atualizado = add.response(design = plan.person, response = resultados)

summary(plan.atualizado)

MEPlot(plan.atualizado)

IAPlot(plan.atualizado)


# Modelo linear = resultados depende (~) de Rede e Acesso
plan.formula = lm(plan.atualizado$resultados~(plan.atualizado$Rede*plan.atualizado$Acesso))

# A tabela fornecida na linha abaixo diz na coluna "Estimate" os valores de q0 (média), qA (Redes1), qB (Acesso1) e qAB (Redes e Acesso)
summary(plan.formula)

# Cálculo das Somas dos Quadrados
plan.anova = anova(plan.formula)

# Valores das médias dos quadrados podem ser observados na tabela gerada pela linha abaixo
# na coluna "Mean Sq". 
summary(plan.anova$`Mean Sq`)

# Cálculo da Soma dos Quadrados Total
SST = plan.anova$"Mean Sq"[1] + plan.anova$"Mean Sq"[2] + plan.anova$"Mean Sq"[3]

#Cálculo da influência de cada fator
InfluenciaA = plan.anova$"Mean Sq"[1] / SST
InfluenciaB = plan.anova$"Mean Sq"[2] / SST
InfluenciaAB = plan.anova$"Mean Sq"[3] / SST

cat("Influencia devido ao fator Redes:", InfluenciaA)
cat("Influencia devido ao fator Acesso:", InfluenciaB)
cat("Influencia devido a interação dos fatores:", InfluenciaAB)