library("tidyverse")

#df

df <- read_csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vSiNayKcnFBJyT4W3uUwJiVhQjCzkigIa0ONek8aQQGLhw6_9kStkMtMffUjFwOoyeDW7gU33me8vYN/pub?gid=119729383&single=true&output=csv")

#agrupo

#ANOVA test

#summary

summary(aov1)
            Df Sum Sq Mean Sq F value Pr(>F)
df$Marca     4  107.2   26.81   2.004  0.109
Residuals   49  655.6   13.38               

#Intervalos de confianza

library("Publish")
library("prodlim")

#IC

ci.mean(df$Halos~df$Marca, data=df)
 df$Marca    mean  CI-95%       
 Control(+)  24.00 [23.06;24.94]
 Corega Tabs 24.67 [22.54;26.79]
 Efferdent   25.33 [22.89;27.78]
 Polident    24.00 [21.83;26.17]
 Walgreens   27.75 [24.79;30.71]
LS0tCnRpdGxlOiAiMjAyM19UYWJsZXRhcyBlZmVydmVzY2VudGVzIEMuYWxiaWNhbnMiCm91dHB1dDogaHRtbF9ub3RlYm9vawotLS0KCmBgYHtyfQpsaWJyYXJ5KCJ0aWR5dmVyc2UiKQpgYGAKI2RmCmBgYHtyfQpkZiA8LSByZWFkX2NzdigiaHR0cHM6Ly9kb2NzLmdvb2dsZS5jb20vc3ByZWFkc2hlZXRzL2QvZS8yUEFDWC0xdlNpTmF5S2NuRkJKeVQ0VzN1VXdKaVZoUWpDemtpZ0lhME9OZWs4YVFRR0xodzZfOWtTdGtNdE1mZlVqRndPb3llRFc3Z1UzM21lOHZZTi9wdWI/Z2lkPTExOTcyOTM4MyZzaW5nbGU9dHJ1ZSZvdXRwdXQ9Y3N2IikKYGBgCiNhZ3J1cG8KYGBge3J9CmRmICU+JSAKICBncm91cF9ieShNYXJjYSkgJT4lIAogIHN1bW1hcmlzZShwcm9tZWRpbyA9IG1lYW4oSGFsb3MpLCBkZXN2aWFjaW9uID0gc2QoSGFsb3MpKQpgYGAKI0FOT1ZBIHRlc3QgCmBgYHtyfQphb3YxIDwtIGFvdihkZiRIYWxvc35kZiRNYXJjYSkKYGBgCiNzdW1tYXJ5IApgYGB7cn0Kc3VtbWFyeShhb3YxKQpgYGAKI0ludGVydmFsb3MgZGUgY29uZmlhbnphCmBgYHtyfQpsaWJyYXJ5KCJQdWJsaXNoIikKbGlicmFyeSgicHJvZGxpbSIpCmBgYAojSUMgCmBgYHtyfQpjaS5tZWFuKGRmJEhhbG9zfmRmJE1hcmNhLCBkYXRhPWRmKQpgYGAKCg==