Tesis mediciones LAE y materiales_2019

packages

cargo df desde la web y creo el objeto

df <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQgBuSiktQ2-rhwyRVO3s0HK3b0p-Wq2v3owxNGcg-To57mHrEfmH6AhuChavjOCMvhjYvIHuNL1liP/pub?gid=474470130&single=true&output=csv")

grafica

df %>% 
 ggplot(aes(x = fct_reorder(Material, `Medida`), y = `Medida`, color=Material)) +
  geom_boxplot()+
  theme_minimal()+
  coord_flip()

cargo plantilla

df2 <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQgBuSiktQ2-rhwyRVO3s0HK3b0p-Wq2v3owxNGcg-To57mHrEfmH6AhuChavjOCMvhjYvIHuNL1liP/pub?gid=1845685926&single=true&output=csv")

Nuevas variables de diferencias entre mediciones

Lt_diferencias <- df2 %>% 
  mutate(Chemfil = LT.LAE. - LT.Chemfil.) %>% 
  mutate(Ketac = LT.LAE. - LT.Ketac.Molar.) %>% 
  mutate(Resina= LT.LAE. - LT.Resina.)

grafica box plot

grafica

bland altman para LAE y chemfil

bland.altman.plot(df2$LT.LAE., df2$LT.Chemfil. , 
                  graph.sys = "ggplot2")

bland altman para LAE y ketac

bland.altman.plot(df2$LT.LAE.,df2$LT.Ketac.Molar. , 
                  graph.sys = "ggplot2")

bland altman para LAE y Resina

bland.altman.plot(df2$LT.LAE.,df2$LT.Resina. , 
                  graph.sys = "ggplot2")

Chemfil

bland.altman.stats(df2$LT.LAE.,df2$LT.Chemfil.)$CI.lines

lower.limit.ci.lower lower.limit.ci.upper   mean.diff.ci.lower   mean.diff.ci.upper 
          -0.7741487            0.1665095            0.3053784            0.8484677 
upper.limit.ci.lower upper.limit.ci.upper 
           0.9873367            1.9279949 

ketac

bland.altman.stats(df2$LT.LAE.,df2$LT.Ketac.Molar.)$CI.lines

lower.limit.ci.lower lower.limit.ci.upper   mean.diff.ci.lower   mean.diff.ci.upper 
         -0.67014167           0.01065676           0.11116273           0.50422189 
upper.limit.ci.lower upper.limit.ci.upper 
          0.60472786           1.28552628 

resina

bland.altman.stats(df2$LT.LAE.,df2$LT.Resina.)$CI.lines

lower.limit.ci.lower lower.limit.ci.upper   mean.diff.ci.lower   mean.diff.ci.upper 
         -1.41553123          -0.23183628          -0.05708791           0.62631868 
upper.limit.ci.lower upper.limit.ci.upper 
          0.80106705           1.98476199 

otro grafico para chemfil

Ketac-molar

resina compuesta

reformateo y abro el nuevo df

df3 <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQgBuSiktQ2-rhwyRVO3s0HK3b0p-Wq2v3owxNGcg-To57mHrEfmH6AhuChavjOCMvhjYvIHuNL1liP/pub?gid=602806611&single=true&output=csv")

bland.altman segun materiales

df3 %>% 
ggplot(aes(x = LT_promedio, y = LT_dif, shape=Longitud)) +
  geom_point(alpha = 0.5) +
  geom_hline(yintercept = mean(df3$LT_dif), colour = "blue", size = 0.5) +
  geom_hline(yintercept = mean(df3$LT_dif) - (1.96 * sd(df3$LT_dif)), colour = "red", size = 0.5) +
  geom_hline(yintercept = mean(df3$LT_dif) + (1.96 * sd(df3$LT_dif)), colour = "red", size = 0.5) +
  labs(title = "Bland antmann", 
             x = "Mean of measurement", 
             y = "Difference", 
             shape = "Length") + 
  theme_minimal()