file.choose()
[1] "C:\\Users\\LENOVO\\Documents\\R Studio work\\trabajo encargado +\\Grado_alcoholico.xlsx"
ruta_Emsamble <- "C:\\Users\\LENOVO\\Documents\\R Studio work\\trabajo encargado +\\Grado_alcoholico.xlsx"
excel_sheets(ruta_Emsamble)
[1] "VINOS POR PRUEBA FQ."
casoDBCA1 <- read_excel(ruta_Emsamble)
print(casoDBCA1)
GRADO <- factor(casoDBCA1$Grado_alcohólico_volumétrico_porcentajeVolumen)
VARi <- factor(casoDBCA1$Variedad)
GRADO1 <- as.numeric(GRADO)
par(mfrow=c(1,1))
boxplot(split(GRADO1,VARi),xlab="variedad", ylab="grado de %alcohol")

resaov<-aov(GRADO1 ~ VARi)
anova(resaov)
Analysis of Variance Table

Response: GRADO1
          Df  Sum Sq Mean Sq F value Pr(>F)
VARi       6  339.71  56.619  0.7115 0.6438
Residuals 24 1909.96  79.582               
cv.model(resaov)
[1] 57.37486
euc.lm <- lm(GRADO1 ~ VARi)
anova(euc.lm , test="F")
Analysis of Variance Table

Response: GRADO1
          Df  Sum Sq Mean Sq F value Pr(>F)
VARi       6  339.71  56.619  0.7115 0.6438
Residuals 24 1909.96  79.582               
shapiro.test(euc.lm$res)

    Shapiro-Wilk normality test

data:  euc.lm$res
W = 0.96325, p-value = 0.3548
fitb <- fitted(resaov)
res_stb <- rstandard(resaov)
plot(fitb,res_stb,xlab="Valores predichos", ylab="valores estandarizados",abline(h=0))

#Prueba de Levene

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

#Método de la diferencia mínima significativa, Least Significant Difference (LSD)

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

Study: resaov ~ "VARi"

LSD t Test for GRADO1 

Mean Square Error:  79.58183 

VARi,  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 GRADO1 

Mean Square Error:  79.58183 

VARi,  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.
```r
SNK.test(resaov, \VARi\,console=TRUE)

<!-- rnb-source-end -->

<!-- rnb-output-begin eyJkYXRhIjoiXG5TdHVkeTogcmVzYW92IH4gXFxWQVJpXFxcblxuU3R1ZGVudCBOZXdtYW4gS2V1bHMgVGVzdFxuZm9yIEdSQURPMSBcblxuTWVhbiBTcXVhcmUgRXJyb3I6ICA3OS41ODE4MyBcblxuVkFSaSwgIG1lYW5zXG5cbkdyb3VwcyBhY2NvcmRpbmcgdG8gcHJvYmFiaWxpdHkgb2YgbWVhbnMgZGlmZmVyZW5jZXMgYW5kIGFscGhhIGxldmVsKCAwLjA1IClcblxuTWVhbnMgd2l0aCB0aGUgc2FtZSBsZXR0ZXIgYXJlIG5vdCBzaWduaWZpY2FudGx5IGRpZmZlcmVudC5cbiJ9 -->

Study: resaov ~

Student Newman Keuls Test for GRADO1

Mean Square Error: 79.58183

VARi, means

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

Means with the same letter are not significantly different.




<!-- rnb-output-end -->

<!-- rnb-frame-begin eyJtZXRhZGF0YSI6eyJjbGFzc2VzIjoiZGF0YS5mcmFtZSIsIm5yb3ciOjcsIm5jb2wiOjksInN1bW1hcnkiOnsiRGVzY3JpcHRpb24iOiJkZiBbNyDDlyA5XSJ9fSwicmRmIjoiSDRzSUFBQUFBQUFBQmwxVFRXZ1RRUlNlNW0rYlFHb3hJaGFLSnFBTkdoczJLVnNqOFdja0phSlN4V3FoSHFmWlNiSjBPeHRuTjdXbGh5Z2VGQlE4aUlnOWlIaFVFVS9ldk9paDlPQzVWdzg5aUtlQ2VQRlFuY20rV1pwOU1IenZmVFBmdnArWm5adFptRW90cEJCQ1VSU0xEYUZvWExnb1BuKzdQbGxCS0JZUndSQ0tvYVRFVlhFb0l4d1pqSWlsWWYyVnRBMThIUFVORndIejJFZmpuYlQzV044Y3Y3czV6bkY1dm0rQmZzeDg4VE4vOGlPZWVQTzYyYXQrNysxS2xmWVc1M3BmanUwOGVvNExlWE43UFYzQjJhMzZqNmQzREp4YnUzck4vdk5WSEVyTGc3SlNzZndpRVJydWQrSFhGdzF5SFBqMTdOdVQvRU04VnFYMGQvc3p6bTJ2dnl3MEgrQWt4SWMyTW5uejhDZXNiZjA5dXZQNGc0b0RmUlo2UW9NOVh0eUYrQWpnTU9BRVV1YnJLOENYUXpNNkMzZ0JzQXA0THFRL0ZjcWo5Qm5BMG4wZjFSMFVRdm9TR3R4WCtoT0E1d0duQVBXUTNnRCtkRWl2WTlCRC9tazBnS0hIRTJka21ickNHVVh5Z253eWNYbnUwc3lORWtSUjF6UFZVK1BnUkZ5cWRtY3RGcmhrVmJrM3kwYmdHbnJnbmpGQytaUGN1VmZjWDRNR0d5TTFza2c1bzE2MnpnbHJBSHN3WUYzU1hiRmF6RkhKa3pYQ2x5bWpYQldXdU9JUjJ5SVFhYlBFWHFTTkpRalQxMm1MazJ5Tlc0NXRxelB4VzJ1Y3RNTURhdGpFVmNVcE1tVVNqeFNiWE5RdG9yMlFSSE02bnVVd0lZckl4NUFJaWNVUUI0blJMcE1UTUNjYjdTNWJtcHpXUS92UkRtdkpwTWovbmVSS0EycjcvS1JmUnVRZnlCTnFFSlMxTEtiR0VyZkZCRzMxWlpPdXFIbUxpK2pmUTdIRExlYXBSZ1hyRmoxSGpGSXhEY2RXVEw5MXRQY2ZlNEl0bTZ3RUFBQT0ifQ== -->

<div data-pagedtable="false">
  <script data-pagedtable-source type="application/json">
{"columns":[{"label":[""],"name":["_rn_"],"type":[""],"align":["left"]},{"label":["GRADO1"],"name":[1],"type":["dbl"],"align":["right"]},{"label":["std"],"name":[2],"type":["dbl"],"align":["right"]},{"label":["r"],"name":[3],"type":["int"],"align":["right"]},{"label":["se"],"name":[4],"type":["dbl"],"align":["right"]},{"label":["Min"],"name":[5],"type":["dbl"],"align":["right"]},{"label":["Max"],"name":[6],"type":["dbl"],"align":["right"]},{"label":["Q25"],"name":[7],"type":["dbl"],"align":["right"]},{"label":["Q50"],"name":[8],"type":["dbl"],"align":["right"]},{"label":["Q75"],"name":[9],"type":["dbl"],"align":["right"]}],"data":[{"1":"16.60000","2":"6.348228","3":"5","4":"3.989532","5":"8","6":"24","7":"13.00","8":"17.0","9":"21.00","_rn_":"Cabernet Franc"},{"1":"10.00000","2":"11.313708","3":"2","4":"6.308004","5":"2","6":"18","7":"6.00","8":"10.0","9":"14.00","_rn_":"Cabernet sauvignon"},{"1":"15.00000","2":"NA","3":"1","4":"8.920865","5":"15","6":"15","7":"15.00","8":"15.0","9":"15.00","_rn_":"Carmenere"},{"1":"11.62500","2":"8.749490","3":"8","4":"3.154002","5":"1","6":"25","7":"4.75","8":"10.5","9":"16.25","_rn_":"Italia"},{"1":"21.66667","2":"13.576941","3":"3","4":"5.150464","5":"6","6":"30","7":"17.50","8":"29.0","9":"29.50","_rn_":"Malbeck"},{"1":"16.77778","2":"8.393119","3":"9","4":"2.973622","5":"3","6":"27","7":"10.00","8":"19.0","9":"22.00","_rn_":"Negra Criolla"},{"1":"18.33333","2":"8.736895","3":"3","4":"5.150464","5":"11","6":"28","7":"13.50","8":"16.0","9":"22.00","_rn_":"Syrah"}],"options":{"columns":{"min":{},"max":[10],"total":[9]},"rows":{"min":[10],"max":[10],"total":[7]},"pages":{}}}
  </script>
</div>

<!-- rnb-frame-end -->

<!-- rnb-frame-begin eyJtZXRhZGF0YSI6eyJjbGFzc2VzIjoiZGF0YS5mcmFtZSIsIm5yb3ciOjcsIm5jb2wiOjIsInN1bW1hcnkiOnsiRGVzY3JpcHRpb24iOiJkZiBbNyDDlyAyXSJ9fSwicmRmIjoiSDRzSUFBQUFBQUFBQnAxUndVN0NRQkJkU2tzRENZVG9sUU1IRTA4MFZhUHh1QWFDOGFBbUtBblhwVjFMd3pKYnR5M3FwK2duNklmd1YraHU2UmpzMFVrMk8vUGFOL3ZlekdRME8ydk5Xb1NRT3JIdEdxazdPaVhPOUhFOHVDVEV0blJSSXpacG12dFYvM1NvRXdOMjlISHArYWVKTDNvNkxZTDZtOTd6cHFlby8yN2lnM3FrQ0hwTWQvY1J3ZWdhUGpabS8weis2bk9BclhoYTlyWktzSEU5dVJyZG4yQVZLWmtuYVlYWVZQTEYyeWVqTVBlV2lUa1BsdmpBdzV0aWk3Sm8zL0ZJc2Y1UXhWSUkxTlFac2psWHdMUCtXREVJc1ArUXFSVUhyampLdU1tWWlKRjA4RXRLV2I2T0k1QlFkUllJbHFJNEJGc2h5NWozcExSdVhXMHJGRmNtV1N4Qmt5eXpza2FGWEZNVm9KdURtVUE0Q0JZNUxBY1hmdVY3UFlISVBGcEF1OU11YjNjdnQzWXlyTytTM2tESEhLSVkwTDhqdEdPQm5VTyt4dkhwUlJSNzhCSVZRNFpHTlpwNm1kUXpReVNRQXBIQ090bitBRlY0SXZmSUFnQUEifQ== -->

<div data-pagedtable="false">
  <script data-pagedtable-source type="application/json">
{"columns":[{"label":[""],"name":["_rn_"],"type":[""],"align":["left"]},{"label":["GRADO1"],"name":[1],"type":["dbl"],"align":["right"]},{"label":["groups"],"name":[2],"type":["chr"],"align":["left"]}],"data":[{"1":"21.66667","2":"a","_rn_":"Malbeck"},{"1":"18.33333","2":"a","_rn_":"Syrah"},{"1":"16.77778","2":"a","_rn_":"Negra Criolla"},{"1":"16.60000","2":"a","_rn_":"Cabernet Franc"},{"1":"15.00000","2":"a","_rn_":"Carmenere"},{"1":"11.62500","2":"a","_rn_":"Italia"},{"1":"10.00000","2":"a","_rn_":"Cabernet sauvignon"}],"options":{"columns":{"min":{},"max":[10],"total":[2]},"rows":{"min":[10],"max":[10],"total":[7]},"pages":{}}}
  </script>
</div>

<!-- rnb-frame-end -->

<!-- rnb-chunk-end -->


<!-- rnb-chunk-begin -->


<!-- rnb-source-begin eyJkYXRhIjoiYGBgclxuc2NoZWZmZS50ZXN0KHJlc2FvdiwgXCJWQVJpXCIsY29uc29sZT1UUlVFKVxuYGBgIn0= -->

```r
scheffe.test(resaov, "VARi",console=TRUE)

Study: resaov ~ "VARi"

Scheffe Test for GRADO1 

Mean Square Error  : 79.58183 

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

Study: resaov ~ "VARi"

Duncan's new multiple range test
for GRADO1 

Mean Square Error:  79.58183 

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

Study: resaov ~ "VARi"

LSD t Test for GRADO1 
P value adjustment method: bonferroni 

Mean Square Error:  79.58183 

VARi,  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= "VARi",  dispersion="se", sig.level=0.05)
summary(sk)
Goups of means at sig.level = 0.05 
tukey_result <- TukeyHSD(resaov, "VARi", conf.level = 0.95)
print(tukey_result)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = GRADO1 ~ VARi)

$VARi
                                        diff        lwr      upr     p adj
Cabernet sauvignon-Cabernet Franc -6.6000000 -30.567485 17.36749 0.9715681
Carmenere-Cabernet Franc          -1.6000000 -32.980804 29.78080 0.9999980
Italia-Cabernet Franc             -4.9750000 -21.306088 11.35609 0.9540023
Malbeck-Cabernet Franc             5.0666667 -15.853870 25.98720 0.9850261
Negra Criolla-Cabernet Franc       0.1777778 -15.800546 16.15610 1.0000000
Syrah-Cabernet Franc               1.7333333 -19.187203 22.65387 0.9999639
Carmenere-Cabernet sauvignon       5.0000000 -30.084806 40.08481 0.9991607
Italia-Cabernet sauvignon          1.6250000 -21.022145 24.27214 0.9999845
Malbeck-Cabernet sauvignon        11.6666667 -14.484004 37.81734 0.7793182
Negra Criolla-Cabernet sauvignon   6.7777778 -15.616318 29.17187 0.9553747
Syrah-Cabernet sauvignon           8.3333333 -17.817337 34.48400 0.9434243
Italia-Carmenere                  -3.3750000 -33.759333 27.00933 0.9997993
Malbeck-Carmenere                  6.6666667 -26.411606 39.74494 0.9942844
Negra Criolla-Carmenere            1.7777778 -28.418415 31.97397 0.9999952
Syrah-Carmenere                    3.3333333 -29.744939 36.41161 0.9998861
Malbeck-Italia                    10.0416667  -9.352190 29.43552 0.6454740
Negra Criolla-Italia               5.1527778  -8.766979 19.07254 0.8916290
Syrah-Italia                       6.7083333 -12.685523 26.10219 0.9186204
Negra Criolla-Malbeck             -4.8888889 -23.986638 14.20886 0.9801776
Syrah-Malbeck                     -3.3333333 -26.723204 20.05654 0.9991607
Syrah-Negra Criolla                1.5555556 -17.542194 20.65330 0.9999673
plot(tukey_result)

ggplot(casoDBCA1, aes(x = Variedad, y = Grado_alcohólico_volumétrico_porcentajeVolumen)) +
  geom_bar(stat = "identity", fill = "skyblue", color = "black") +
  labs(title = "Gráfico de Barras",
       x = "VARIEDAD",
       y = "Grado_alcohólico_volumétrico_porcentajeVolumen")

LS0tDQp0aXRsZTogIkdyYWRvIGFsY29ob2xpY28gdmlub3MiDQphdXRob3I6ICJKZXNzaWNhIFRoYWxpYSBNaXJhbmRhIE9saXZhIg0KZGF0ZTogIjE1LzA3LzIwMjQiDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KYGBge3J9DQpmaWxlLmNob29zZSgpDQpgYGANCmBgYHtyfQ0KcnV0YV9FbXNhbWJsZSA8LSAiQzpcXFVzZXJzXFxMRU5PVk9cXERvY3VtZW50c1xcUiBTdHVkaW8gd29ya1xcdHJhYmFqbyBlbmNhcmdhZG8gK1xcR3JhZG9fYWxjb2hvbGljby54bHN4Ig0KYGBgDQpgYGB7cn0NCmV4Y2VsX3NoZWV0cyhydXRhX0Vtc2FtYmxlKQ0KYGBgDQpgYGB7cn0NCmNhc29EQkNBMSA8LSByZWFkX2V4Y2VsKHJ1dGFfRW1zYW1ibGUpDQpgYGANCmBgYHtyfQ0KcHJpbnQoY2Fzb0RCQ0ExKQ0KYGBgDQpgYGB7cn0NCkdSQURPIDwtIGZhY3RvcihjYXNvREJDQTEkR3JhZG9fYWxjb2jDs2xpY29fdm9sdW3DqXRyaWNvX3BvcmNlbnRhamVWb2x1bWVuKQ0KYGBgDQpgYGB7cn0NClZBUmkgPC0gZmFjdG9yKGNhc29EQkNBMSRWYXJpZWRhZCkNCmBgYA0KYGBge3J9DQpHUkFETzEgPC0gYXMubnVtZXJpYyhHUkFETykNCmBgYA0KYGBge3J9DQpwYXIobWZyb3c9YygxLDEpKQ0KYGBgDQpgYGB7cn0NCmJveHBsb3Qoc3BsaXQoR1JBRE8xLFZBUmkpLHhsYWI9InZhcmllZGFkIiwgeWxhYj0iZ3JhZG8gZGUgJWFsY29ob2wiKQ0KYGBgDQpgYGB7cn0NCnJlc2FvdjwtYW92KEdSQURPMSB+IFZBUmkpDQpgYGANCmBgYHtyfQ0KYW5vdmEocmVzYW92KQ0KYGBgDQpgYGB7cn0NCmN2Lm1vZGVsKHJlc2FvdikNCmBgYA0KYGBge3J9DQpldWMubG0gPC0gbG0oR1JBRE8xIH4gVkFSaSkNCmBgYA0KYGBge3J9DQphbm92YShldWMubG0gLCB0ZXN0PSJGIikNCmBgYA0KYGBge3J9DQpzaGFwaXJvLnRlc3QoZXVjLmxtJHJlcykNCmBgYA0KYGBge3J9DQpmaXRiIDwtIGZpdHRlZChyZXNhb3YpDQpgYGANCmBgYHtyfQ0KcmVzX3N0YiA8LSByc3RhbmRhcmQocmVzYW92KQ0KYGBgDQpgYGB7cn0NCnBsb3QoZml0YixyZXNfc3RiLHhsYWI9IlZhbG9yZXMgcHJlZGljaG9zIiwgeWxhYj0idmFsb3JlcyBlc3RhbmRhcml6YWRvcyIsYWJsaW5lKGg9MCkpDQpgYGANCiNQcnVlYmEgZGUgTGV2ZW5lDQpgYGB7cn0NCmxldmVuZVRlc3QoR1JBRE8xIH4gVkFSaSwgY2VudGVyID0gIm1lZGlhbiIpDQpgYGANCiNNw6l0b2RvIGRlIGxhIGRpZmVyZW5jaWEgbcOtbmltYSBzaWduaWZpY2F0aXZhLCBMZWFzdCBTaWduaWZpY2FudCBEaWZmZXJlbmNlIChMU0QpDQpgYGB7cn0NCm91dExTRCA8LUxTRC50ZXN0KHJlc2FvdiwgIlZBUmkiLGNvbnNvbGU9VFJVRSkNCmBgYA0KYGBge3J9DQpvdXRIU0Q8LUhTRC50ZXN0KHJlc2FvdiwgIlZBUiIsY29uc29sZT1UUlVFKQ0KYGBgDQpgYGB7cn0NClNOSy50ZXN0KHJlc2FvdiwgIlZBUmkiLGNvbnNvbGU9VFJVRSkNCmBgYA0KYGBge3J9DQpzY2hlZmZlLnRlc3QocmVzYW92LCAiVkFSaSIsY29uc29sZT1UUlVFKQ0KYGBgDQpgYGB7cn0NCmR1bmNhbi50ZXN0KHJlc2FvdiwgIlZBUmkiLGNvbnNvbGU9VFJVRSkNCmBgYA0KYGBge3J9DQpMU0QudGVzdChyZXNhb3YsICJWQVJpIiwgcC5hZGo9ICJib24iLGNvbnNvbGU9VFJVRSkNCmBgYA0KYGBge3J9DQpzayA8LSBTSyhyZXNhb3YsIHdoaWNoPSAiVkFSaSIsICBkaXNwZXJzaW9uPSJzZSIsIHNpZy5sZXZlbD0wLjA1KQ0KYGBgDQpgYGB7cn0NCnN1bW1hcnkoc2spDQpgYGANCmBgYHtyfQ0KdHVrZXlfcmVzdWx0IDwtIFR1a2V5SFNEKHJlc2FvdiwgIlZBUmkiLCBjb25mLmxldmVsID0gMC45NSkNCmBgYA0KYGBge3J9DQpwcmludCh0dWtleV9yZXN1bHQpDQpgYGANCmBgYHtyfQ0KcGxvdCh0dWtleV9yZXN1bHQpDQpgYGANCmBgYHtyfQ0KZ2dwbG90KGNhc29EQkNBMSwgYWVzKHggPSBWYXJpZWRhZCwgeSA9IEdyYWRvX2FsY29ow7NsaWNvX3ZvbHVtw6l0cmljb19wb3JjZW50YWplVm9sdW1lbikpICsNCiAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIsIGZpbGwgPSAic2t5Ymx1ZSIsIGNvbG9yID0gImJsYWNrIikgKw0KICBsYWJzKHRpdGxlID0gIkdyw6FmaWNvIGRlIEJhcnJhcyIsDQogICAgICAgeCA9ICJWQVJJRURBRCIsDQogICAgICAgeSA9ICJHcmFkb19hbGNvaMOzbGljb192b2x1bcOpdHJpY29fcG9yY2VudGFqZVZvbHVtZW4iKQ0KYGBgDQoNCg==