R Notebook

file.choose()

[1] "C:\\Users\\LENOVO\\Downloads\\Copia de Grado_alcohólico(1).xlsx"

ruta_Ensamble<-"C:\\Users\\LENOVO\\Downloads\\Copia de Grado_alcohólico(1).xlsx"

excel_sheets(ruta_Ensamble)

[1] "VINOS POR PRUEBA FQ."

casoDBCA1<-read_excel(ruta_Ensamble)

print(casoDBCA1)

VARD<-factor(casoDBCA1$Variedad)

GRADVOL<-as.vector(casoDBCA1$`Grado_alcohólico_volumétrico_%Volumen`)

GRADVOL1<-as.numeric(GRADVOL)

par(mfrow=c(1,1))

boxplot(split(GRADVOL1,VARD), xlab="Variedad",ylab="`Grado_alcohólico_volumétrico_%Volumen`")

resaov<-aov(GRADVOL1~VARD)

anova(resaov)

Analysis of Variance Table

Response: GRADVOL1
          Df  Sum Sq Mean Sq F value Pr(>F)
VARD       6  3.3901 0.56502  0.7298 0.6302
Residuals 24 18.5813 0.77422               

cv.model(resaov)

[1] 6.887564

euc.lm <- lm(GRADVOL1~VARD)

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

Analysis of Variance Table

Response: GRADVOL1
          Df  Sum Sq Mean Sq F value Pr(>F)
VARD       6  3.3901 0.56502  0.7298 0.6302
Residuals 24 18.5813 0.77422               

shapiro.test(euc.lm$res)

    Shapiro-Wilk normality test

data:  euc.lm$res
W = 0.94655, p-value = 0.1253

fitb <- fitted(resaov)

plot(fitb,res_stb,xlab = "`Grado_alcohólico_volumétrico_%Volumen`",ylab = "variedad")

leveneTest(GRADVOL1~VARD , center ="median")

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

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

Study: resaov ~ "VARD"

LSD t Test for GRADVOL1 

Mean Square Error:  0.7742194 

VARD,  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, "VARD",console=TRUE)

Study: resaov ~ "VARD"

HSD Test for GRADVOL1 

Mean Square Error:  0.7742194 

VARD,  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, "VARD",console=TRUE)

Study: resaov ~ "VARD"

Student Newman Keuls Test
for GRADVOL1 

Mean Square Error:  0.7742194 

VARD,  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, "VARD",console=TRUE)

Study: resaov ~ "VARD"

Scheffe Test for GRADVOL1 

Mean Square Error  : 0.7742194 

VARD,  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, "VARD",console=TRUE)

Study: resaov ~ "VARD"

Duncan's new multiple range test
for GRADVOL1 

Mean Square Error:  0.7742194 

VARD,  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, "VARD", p.adj= "bon",console=TRUE)

Study: resaov ~ "VARD"

LSD t Test for GRADVOL1 
P value adjustment method: bonferroni 

Mean Square Error:  0.7742194 

VARD,  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= "VARD",  dispersion="se", sig.level=0.05)

summary(sk)

Goups of means at sig.level = 0.05