knitr rmarkdown

#Asignatura: Disenos Experimentales 
#Facultad  : Ingeniería Estadistica e Informatica
#Universidad Nacional del Altiplano
library(knitr)
## Warning: package 'knitr' was built under R version 3.6.3
library(readxl)
## Warning: package 'readxl' was built under R version 3.6.3
DCAM <- read_excel("DCAM.xlsx")
View(DCAM)
head(DCAM)
## # A tibble: 6 x 2
##   Dietas Peso   
##   <chr>  <chr>  
## 1 Dieta1 1281.25
## 2 Dieta1 1211.91
## 3 Dieta1 1280.01
## 4 Dieta1 1214.01
## 5 Dieta1 1222.20
## 6 Dieta1 1272.00
str(DCAM)
## tibble [40 x 2] (S3: tbl_df/tbl/data.frame)
##  $ Dietas: chr [1:40] "Dieta1" "Dieta1" "Dieta1" "Dieta1" ...
##  $ Peso  : chr [1:40] "1281.25" "1211.91" "1280.01" "1214.01" ...
DCAM$Dietas <- factor(DCAM$Dietas)
DCAM$Peso <- as.numeric(DCAM$Peso)
head(DCAM)
## # A tibble: 6 x 2
##   Dietas  Peso
##   <fct>  <dbl>
## 1 Dieta1 1281.
## 2 Dieta1 1212.
## 3 Dieta1 1280.
## 4 Dieta1 1214.
## 5 Dieta1 1222.
## 6 Dieta1 1272
str(DCAM)
## tibble [40 x 2] (S3: tbl_df/tbl/data.frame)
##  $ Dietas: Factor w/ 4 levels "Dieta1","Dieta2",..: 1 1 1 1 1 1 1 1 1 1 ...
##  $ Peso  : num [1:40] 1281 1212 1280 1214 1222 ...
#Graficos-Homogeneidad
library(ggplot2)
ggplot(DCAM) + geom_bar(aes(Dietas, fill = Dietas))

#Tabla ANOVA
fit<-aov(Peso~Dietas, data=DCAM)
summary(fit)
##             Df Sum Sq Mean Sq F value Pr(>F)  
## Dietas       3  31487   10496    3.13 0.0375 *
## Residuals   36 120712    3353                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Prueba de Shaphiro-Wilk
shapiro.test(residuals(fit))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(fit)
## W = 0.9631, p-value = 0.2135
#Grafico-Normalidad
library(car)
## Warning: package 'car' was built under R version 3.6.3
## Loading required package: carData
## Warning: package 'carData' was built under R version 3.6.3
qqPlot(fit)
## [1] 15 18
mtext("By:Amarillas", side = 3, adj = 1, family = "mono")

#Prueba de Bartlett
bartlett.test(Peso~Dietas, data=DCAM)
## 
##  Bartlett test of homogeneity of variances
## 
## data:  Peso by Dietas
## Bartlett's K-squared = 14.635, df = 3, p-value = 0.002157
#Grafico-Cajas
boxplot(Peso~Dietas,data=DCAM, col=c("red","green","blue","orange"), main="Diagrama de Cajas de los Tratamientos")
mtext("By:Amarillas", side = 3, adj = 1, family = "mono")

#Grupos
library(agricolae)
## Warning: package 'agricolae' was built under R version 3.6.3
Grupos1 <- LSD.test(fit, "Dietas", group = T, console = T); 
## 
## Study: fit ~ "Dietas"
## 
## LSD t Test for Peso 
## 
## Mean Square Error:  3353.124 
## 
## Dietas,  means and individual ( 95 %) CI
## 
##            Peso      std  r      LCL      UCL     Min     Max
## Dieta1 1245.105 31.95019 10 1207.967 1282.243 1211.91 1290.96
## Dieta2 1290.734 94.04522 10 1253.596 1327.872 1165.13 1445.09
## Dieta3 1293.067 31.69722 10 1255.929 1330.205 1255.58 1355.47
## Dieta4 1323.716 50.42288 10 1286.578 1360.854 1233.42 1380.61
## 
## Alpha: 0.05 ; DF Error: 36
## Critical Value of t: 2.028094 
## 
## least Significant Difference: 52.52039 
## 
## Treatments with the same letter are not significantly different.
## 
##            Peso groups
## Dieta4 1323.716      a
## Dieta3 1293.067     ab
## Dieta2 1290.734     ab
## Dieta1 1245.105      b
Grupos2 <- LSD.test(fit, "Dietas", group = F, console = T)
## 
## Study: fit ~ "Dietas"
## 
## LSD t Test for Peso 
## 
## Mean Square Error:  3353.124 
## 
## Dietas,  means and individual ( 95 %) CI
## 
##            Peso      std  r      LCL      UCL     Min     Max
## Dieta1 1245.105 31.95019 10 1207.967 1282.243 1211.91 1290.96
## Dieta2 1290.734 94.04522 10 1253.596 1327.872 1165.13 1445.09
## Dieta3 1293.067 31.69722 10 1255.929 1330.205 1255.58 1355.47
## Dieta4 1323.716 50.42288 10 1286.578 1360.854 1233.42 1380.61
## 
## Alpha: 0.05 ; DF Error: 36
## Critical Value of t: 2.028094 
## 
## Comparison between treatments means
## 
##                 difference pvalue signif.        LCL        UCL
## Dieta1 - Dieta2    -45.629 0.0866       .  -98.14939   6.891386
## Dieta1 - Dieta3    -47.962 0.0722       . -100.48239   4.558386
## Dieta1 - Dieta4    -78.611 0.0044      ** -131.13139 -26.090614
## Dieta2 - Dieta3     -2.333 0.9287          -54.85339  50.187386
## Dieta2 - Dieta4    -32.982 0.2110          -85.50239  19.538386
## Dieta3 - Dieta4    -30.649 0.2444          -83.16939  21.871386
#Prueba-Bonferroni-Ducan
LSD.test(fit, "Dietas", p.adj= "bonferroni",console=TRUE) -> bon; 
## 
## Study: fit ~ "Dietas"
## 
## LSD t Test for Peso 
## P value adjustment method: bonferroni 
## 
## Mean Square Error:  3353.124 
## 
## Dietas,  means and individual ( 95 %) CI
## 
##            Peso      std  r      LCL      UCL     Min     Max
## Dieta1 1245.105 31.95019 10 1207.967 1282.243 1211.91 1290.96
## Dieta2 1290.734 94.04522 10 1253.596 1327.872 1165.13 1445.09
## Dieta3 1293.067 31.69722 10 1255.929 1330.205 1255.58 1355.47
## Dieta4 1323.716 50.42288 10 1286.578 1360.854 1233.42 1380.61
## 
## Alpha: 0.05 ; DF Error: 36
## Critical Value of t: 2.791972 
## 
## Minimum Significant Difference: 72.3021 
## 
## Treatments with the same letter are not significantly different.
## 
##            Peso groups
## Dieta4 1323.716      a
## Dieta3 1293.067     ab
## Dieta2 1290.734     ab
## Dieta1 1245.105      b
duncan.test(fit, "Dietas",console=TRUE)
## 
## Study: fit ~ "Dietas"
## 
## Duncan's new multiple range test
## for Peso 
## 
## Mean Square Error:  3353.124 
## 
## Dietas,  means
## 
##            Peso      std  r     Min     Max
## Dieta1 1245.105 31.95019 10 1211.91 1290.96
## Dieta2 1290.734 94.04522 10 1165.13 1445.09
## Dieta3 1293.067 31.69722 10 1255.58 1355.47
## Dieta4 1323.716 50.42288 10 1233.42 1380.61
## 
## Alpha: 0.05 ; DF Error: 36 
## 
## Critical Range
##        2        3        4 
## 52.52039 55.21329 56.96961 
## 
## Means with the same letter are not significantly different.
## 
##            Peso groups
## Dieta4 1323.716      a
## Dieta3 1293.067     ab
## Dieta2 1290.734     ab
## Dieta1 1245.105      b