Pular para o conteúdo principal

agrinewbras.com.br

Portal de Materiais AcadêmicosCódigos | Santos H. B. Dias
CÓDIGOS DE PROGRAMAÇÃO

Scripts em R utilizados nos vídeos do canal

Nesta página você encontra todos os códigos utilizados nos vídeos sobre análises estatísticas. Cada vídeo possui seu próprio script comentado para facilitar o aprendizado.

📺
Vídeos disponíveis

Miniatura

Análise Estatística de Experimento DBC no R

Aprenda a importar os dados, executar a ANOVA e interpretar corretamente os resultados utilizando o software R.

📅 Julho 2026⏱ 18 min
Assistir vídeo

Importação dos dados



install.packages("readxl")
install.packages("ExpDes")

library(readxl)
require(ExpDes)

 dados = read_excel("D:/OneDrive/Pesquisa/TCC_IC/TCC_2026_Pedro/Resultados.xlsx")

str(dados)
summary(dados)

#transformar 
dados$BLOCO = as.factor(dados$BLOCO)
dados$TRAT = as.factor(dados$TRAT)
str(dados)


dados_Alt = subset(dados,!is.na(ALTURA_CM))

rbd(
  treat=dados_Alt$TRAT,
  block=dados_Alt$BLOCO,
  resp=dados_Alt$ALTURA_CM,
  quali = TRUE,
  mcomp = "tukey",
  hvar = "oneillmathews",
  sigT = 0.05,
  sigF = 0.05
)

dados_BROTOS = subset(dados, !is.na(N_BROTOS))

rbd(
  treat = dados_BROTOS$TRAT,
  block = dados_BROTOS$BLOCO,
  resp = dados_BROTOS$N_BROTOS,
  quali = TRUE,
  mcomp = "tukey",
  sigT = 0.05,
  sigF = 0.05
)

dados_PESO_BRUTO = subset(dados,!is.na(peso_bruto_kg))

rbd(
  treat = dados_PESO_BRUTO$TRAT,
  block = dados_PESO_BRUTO$BLOCO,
  resp = dados_PESO_BRUTO$peso_bruto_kg,
  quali = TRUE,
  mcomp = "tukey",
  sigT = 0.05,
  sigF = 0.05
)


dados_PESO_AMIDO = subset(dados,!is.na(peso_amido_kg))

rbd(
  treat = dados_PESO_AMIDO$TRAT,
  block = dados_PESO_AMIDO$BLOCO,
  resp = dados_PESO_AMIDO$peso_amido_kg,
  quali = TRUE,
  mcomp = "tukey",
  sigT = 0.05,
  sigF = 0.05
)

dados_RAIZ = subset(dados,!is.na(N_raiz_UN))

rbd(
  treat = dados_RAIZ$TRAT,
  block = dados_RAIZ$BLOCO,
  resp = dados_RAIZ$N_raiz_UN,
  quali = TRUE,
  mcomp = "tukey",
  sigT = 0.05,
  sigF = 0.05
)

dados_NDVI = subset(dados,!is.na(NDVI))

rbd(
  treat = dados_NDVI$TRAT,
  block = dados_NDVI$BLOCO,
  resp = dados_NDVI$NDVI,
  quali = TRUE,
  mcomp = "tukey",
  sigT = 0.05,
  sigF = 0.05
)


Miniatura

Teste de Normalidade e Homogeneidade

Aprenda a verificar os pressupostos da ANOVA utilizando Shapiro-Wilk e Bartlett.

📅 Em breve⏱ 15 min
Assistir vídeo

Teste de Shapiro-Wilk



shapiro.test(residuals(modelo))

Teste de Bartlett



bartlett.test(
Produtividade ~ Tratamento,
data = dados)

Miniatura

Gráficos da ANOVA no R

Construção de gráficos para apresentação de resultados utilizando ggplot2.

📅 Em breve⏱ 20 min
Assistir vídeo

Gráfico de médias



library(ggplot2)

ggplot(
dados,
aes(
Tratamento,
Produtividade)
)+
geom_boxplot()

Perguntas Frequentes

Posso utilizar estes códigos em pesquisas?

Sim. Todos os scripts são disponibilizados gratuitamente para fins acadêmicos e podem ser adaptados conforme a necessidade do seu experimento.

Os códigos são atualizados?

Sim. Sempre que um vídeo receber melhorias ou correções, o script desta página também será atualizado.

© 2026 Santos Henrique Brant Dias