Tesis mediciones LAE y materiales_2019

packages

library("tidyverse")
library("ggthemes")
library("forcats")
library("ggplot2")
library("BlandAltmanLeh")

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(Control=LT.control.- LT.control.) %>% 
  mutate(LAE = LT.control. - LT.LAE.) %>% 
  mutate(Chemfil = LT.control. - LT.Chemfil.) %>% 
  mutate(Ketac = LT.control. - LT.Ketac.Molar.) %>% 
  mutate(Resina= LT.control. - LT.Resina.)

boxplot de las diferencias

bland altman para control vs lae

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

bland altman para control vs chemfil

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

Chemfil

bland.altman.stats(df2$LT.control.,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.control.,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.control.,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 

LT lae

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

lower.limit.ci.lower lower.limit.ci.upper   mean.diff.ci.lower   mean.diff.ci.upper 
                   0                    0                    0                    0 
upper.limit.ci.lower upper.limit.ci.upper 
                   0                    0 

Visualizo en otra grafica para LAE vs control

df2 %>% 
ggplot(mapping = aes(x = LT.control., y = LT.LAE.)) +
  geom_point(color = "black", size = 1) +
  labs(title = "Diagrama de dispersión", x = "método A", y = "método B") +
  geom_smooth(method = "lm", se = TRUE, color = "red", lwd = 0.5) +
  geom_abline(intercept = 0, slope = 1, lwd = 0.7, col = "blue") +
  theme(axis.line = element_line(colour = "black"), 
        panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(), 
        panel.border = element_blank(), 
        panel.background = element_blank()) +
  theme_bw() + theme(plot.title = element_text(hjust = 0.5)) 

Para control vs chemfil

df2 %>% 
ggplot(mapping = aes(x = LT.control., y = LT.Chemfil.)) +
  geom_point(color = "black", size = 1) +
  labs(title = "Diagrama de dispersión", x = "método A", y = "método B") +
  geom_smooth(method = "lm", se = TRUE, color = "red", lwd = 0.5) +
  geom_abline(intercept = 0, slope = 1, lwd = 0.7, col = "blue") +
  theme(axis.line = element_line(colour = "black"), 
        panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(), 
        panel.border = element_blank(), 
        panel.background = element_blank()) +
  theme_bw() + theme(plot.title = element_text(hjust = 0.5)) 

Para control vs Ketac-molar

df2 %>% 
ggplot(mapping = aes(x = LT.control., y = LT.Ketac.Molar.)) +
  geom_point(color = "black", size = 1) +
  labs(title = "Diagrama de dispersión", x = "método A", y = "método B") +
  geom_smooth(method = "lm", se = TRUE, color = "red", lwd = 0.5) +
  geom_abline(intercept = 0, slope = 1, lwd = 0.7, col = "blue") +
  theme(axis.line = element_line(colour = "black"), 
        panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(), 
        panel.border = element_blank(), 
        panel.background = element_blank()) +
  theme_bw() + theme(plot.title = element_text(hjust = 0.5)) 

Para control vs resina compuest

df2 %>% 
ggplot(mapping = aes(x = LT.control., y = LT.Resina.)) +
  geom_point(color = "black", size = 1) +
  labs(title = "Diagrama de dispersión", x = "método A", y = "método B") +
  geom_smooth(method = "lm", se = TRUE, color = "red", lwd = 0.5) +
  geom_abline(intercept = 0, slope = 1, lwd = 0.7, col = "blue") +
  theme(axis.line = element_line(colour = "black"), 
        panel.grid.major = element_blank(), 
        panel.grid.minor = element_blank(), 
        panel.border = element_blank(), 
        panel.background = element_blank()) +
  theme_bw() + theme(plot.title = element_text(hjust = 0.5)) 

reformateo y abro el nuevo df

bland antmamm