# # Exemplo Componentes de variância
# Planejamento de Experimentos I
# Profa. Cibele Russo

# Exemplo teares

# Carregando pacotes necessários
library(lme4)

dados<-read.csv('https://raw.githubusercontent.com/cibelerusso/PlanejamentodeExperimentos/main/Dados/teares.csv', sep=',')
attach(dados)

# Carregando a biblioteca necessária
library(ggplot2)

# Gráfico de dispersão
ggplot(dados, aes(x = Tear, y = Resistência)) +
  geom_point() +
  labs(title = "Gráfico de Dispersão", x = "Tear", y = "Resistência")

# Boxplot por Tear
ggplot(dados, aes(x = factor(Tear), y = Resistência)) +
  geom_boxplot() +
  labs(title = "Boxplot por Tear", x = "Tear", y = "Resistência")

# Definindo o tamanho da amostra

n_teares <- 4      # Número de peças por operador
n_pecas <- 4 # Número de observações por tear
n_total <- n_teares * n_pecas

# Criando o dataframe
dados <- data.frame(Tear = factor(Tear),
                    Medicao = Resistência)

# Ajustando o modelo com efeitos aleatórios
modelo <- lmer(Medicao ~ (1 | Tear), data = dados)

# Resumo do modelo
summary(modelo)

# Ajustando o modelo sem efeitos 
modelo2 <- lm(Medicao ~ 1, data=dados)

anova(modelo, modelo2)