Grado Alcohólico

ruta_Ensamble <-"C:/Users/shami/Downloads/Grado_alcohólico.xlsx"

excel_sheets(ruta_Ensamble)

[1] "VINOS POR PRUEBA FQ."

casoDBCA1 <- read_excel(ruta_Ensamble)

print(casoDBCA1)

names(Grado_alcohólico)

[1] "Grado_alcohólico_volumétrico_%Volumen" "Variedad"                             

VOL<-factor(casoDBCA1$`Grado_alcohólico_volumétrico_%Volumen`)

VAR<-factor(casoDBCA1$Variedad)

VOL1<-as.numeric(VOL)

par(mfrow=c(1,1))

boxplot(split(VOL,VAR),xlab="Grado_alcohólico_volumétrico_%Volumen", ylab="Variedad")

resaov<-aov(VOL1 ~ VAR)

anova(resaov)

Analysis of Variance Table

Response: VOL1
          Df  Sum Sq Mean Sq F value Pr(>F)
VAR        6  339.71  56.619  0.7115 0.6438
Residuals 24 1909.96  79.582               

cv.model(resaov)

[1] 57.37486

euc.lm <- lm(VOL1 ~ VAR)

anova(euc.lm , test="F")

Analysis of Variance Table

Response: VOL1
          Df  Sum Sq Mean Sq F value Pr(>F)
VAR        6  339.71  56.619  0.7115 0.6438
Residuals 24 1909.96  79.582               

boxplot(anova(euc.lm , test="F"))

shapiro.test(euc.lm$res)

    Shapiro-Wilk normality test

data:  euc.lm$res
W = 0.96325, p-value = 0.3548

plot(anova(euc.lm , test="F"))

fitb <- fitted(resaov)

res_stb <- rstandard(resaov)

plot(fitb,res_stb,xlab="Grado_alcohólico_volumétrico_%Volumen", ylab="Variedad",abline(h=0))

leveneTest(GRADO1 ~ VAR, center = "median")

Levene's Test for Homogeneity of Variance (center = "median")
      Df F value Pr(>F)
group  6  0.3228 0.9185
      24               

outLSD <-LSD.test(resaov, "VAR",console=TRUE)

Study: resaov ~ "VAR"

LSD t Test for VOL1 

Mean Square Error:  79.58183 

VAR,  means and individual ( 95 %) CI

Alpha: 0.05 ; DF Error: 24
Critical Value of t: 2.063899 

Groups according to probability of means differences and alpha level( 0.05 )

Treatments with the same letter are not significantly different.

outHSD<-HSD.test(resaov, "VAR",console=TRUE)

Study: resaov ~ "VAR"

HSD Test for VOL1 

Mean Square Error:  79.58183 

VAR,  means

Alpha: 0.05 ; DF Error: 24 
Critical Value of Studentized Range: 4.541314 

Groups according to probability of means differences and alpha level( 0.05 )

Treatments with the same letter are not significantly different.

SNK.test(resaov, "VAR",console=TRUE)

Study: resaov ~ "VAR"

Student Newman Keuls Test
for VOL1 

Mean Square Error:  79.58183 

VAR,  means

Groups according to probability of means differences and alpha level( 0.05 )

Means with the same letter are not significantly different.

scheffe.test(resaov, "VAR",console=TRUE)

Study: resaov ~ "VAR"

Scheffe Test for VOL1 

Mean Square Error  : 79.58183 

VAR,  means

Alpha: 0.05 ; DF Error: 24 
Critical Value of F: 2.508189 

Groups according to probability of means differences and alpha level( 0.05 )

Means with the same letter are not significantly different.

duncan.test(resaov, "VAR",console=TRUE)

Study: resaov ~ "VAR"

Duncan's new multiple range test
for VOL1 

Mean Square Error:  79.58183 

VAR,  means

Groups according to probability of means differences and alpha level( 0.05 )

Means with the same letter are not significantly different.

LSD.test(resaov, "VAR", p.adj= "bon",console=TRUE)

Study: resaov ~ "VAR"

LSD t Test for VOL1 
P value adjustment method: bonferroni 

Mean Square Error:  79.58183 

VAR,  means and individual ( 95 %) CI

Alpha: 0.05 ; DF Error: 24
Critical Value of t: 3.395988 

Groups according to probability of means differences and alpha level( 0.05 )

Treatments with the same letter are not significantly different.

sk <- SK(resaov, which= "VAR",  dispersion="se", sig.level=0.05)

summary(sk)

Goups of means at sig.level = 0.05