#ANALIZO LOS DATOS DE FLUOROMETRO #packages

#abro el df

df <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQ4Nk96Jx1n-pLC1e7lXlUgMtqAFc1xOTyNCP6_72cvgtcZXjITM2-wocUbgKWin63eOHWZP3ezwJoB/pub?gid=1563909949&single=true&output=csv")

#cambio colutorio como factor

df <- df %>% 
  mutate(Colutorio = as.factor(Colutorio))

#glimse

glimpse(df)
Rows: 427
Columns: 4
$ id        <int> 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4,…
$ Colutorio <fct> Cymenol, Cymenol, Cymenol, Cymenol, Cymenol, Cymenol, Cymenol, Cymenol, Cy…
$ Tiempo    <chr> "T0", "T1", "T2", "T3", "T4", "T5", "T6", "T0", "T1", "T2", "T3", "T4", "T…
$ Reduccion <dbl> 100.0, 60.1, 63.6, 39.2, 30.5, 48.0, 22.6, 100.0, 94.0, 51.0, 56.7, 56.5, …

#hist

df %>% 
  # pivot_longer(Colutorio:Tiempo, 
  #             names_to = "Variables", 
  #             values_to = "Values") %>% 
  ggplot(aes(x = Reduccion)) + 
  geom_histogram() + 
  facet_grid(Tiempo ~ Colutorio)

#grafica lineal

df %>% 
ggline(x = "Tiempo", 
       y = "Reduccion", 
       color = "Colutorio",
       add = c("mean_se", "dotplot")) + 
  labs(title = "Profile plot de reducción en tiempo y tipo de colutorio") + 
  theme_pubr()

#bonferroni

pairwise.t.test( df$Reduccion, 
                 df$Colutorio,
           
                 p.adj = "bonferroni")

    Pairwise comparisons using t tests with pooled SD 

data:  df$Reduccion and df$Colutorio 

        CPC     CPC+Cy  Cymenol
CPC+Cy  1.0000  -       -      
Cymenol 0.1441  0.0042  -      
Placebo 1.1e-06 4.8e-09 0.0154 

P value adjustment method: bonferroni 

#modelo de aov

aov1 <- aov(Reduccion ~ Tiempo + Colutorio + Tiempo*Colutorio , 
           data = df)
summary(aov1)
                  Df Sum Sq Mean Sq F value   Pr(>F)    
Tiempo             6  67791   11299  12.268 1.11e-12 ***
Colutorio          3  50974   16991  18.449 3.16e-11 ***
Tiempo:Colutorio  18  20039    1113   1.209     0.25    
Residuals        399 367484     921                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
df %>% 
  with(aov(Reduccion ~ Tiempo + Colutorio + Tiempo:Colutorio + Error(id))) %>% 
  summary()

Error: id
          Df Sum Sq Mean Sq
Colutorio  1   5080    5080

Error: Within
                  Df Sum Sq Mean Sq F value   Pr(>F)    
Tiempo             6  67791   11299  12.328 9.66e-13 ***
Colutorio          3  48611   16204  17.680 8.59e-11 ***
Tiempo:Colutorio  18  20039    1113   1.215    0.245    
Residuals        398 364768     917                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

#tukey

TukeyHSD(aov1)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = Reduccion ~ Tiempo + Colutorio + Tiempo * Colutorio, data = df)

$Tiempo
            diff         lwr         upr     p adj
T1-T0 -36.850820 -53.1355806 -20.5660587 0.0000000
T2-T0 -34.357377 -50.6421380 -18.0726161 0.0000000
T3-T0 -32.803279 -49.0880396 -16.5185178 0.0000001
T4-T0 -26.296721 -42.5814822 -10.0119604 0.0000492
T5-T0 -15.601639 -31.8864003   0.6831216 0.0703924
T6-T0 -12.978689 -29.2634495   3.3060724 0.2177958
T2-T1   2.493443 -13.7913183  18.7782036 0.9993384
T3-T1   4.047541 -12.2372200  20.3323019 0.9902257
T4-T1  10.554098  -5.7306626  26.8388593 0.4678243
T5-T1  21.249180   4.9644194  37.5339413 0.0024487
T6-T1  23.872131   7.5873702  40.1568921 0.0003543
T3-T2   1.554098 -14.7306626  17.8388593 0.9999579
T4-T2   8.060656  -8.2241052  24.3454167 0.7643500
T5-T2  18.755738   2.4709768  35.0404986 0.0124575
T6-T2  21.378689   5.0939276  37.6634495 0.0022373
T4-T3   6.506557  -9.7782036  22.7913183 0.8999462
T5-T3  17.201639   0.9168784  33.4864003 0.0306798
T6-T3  19.824590   3.5398292  36.1093511 0.0063684
T5-T4  10.695082  -5.5896790  26.9798429 0.4509527
T6-T4  13.318033  -2.9667281  29.6027937 0.1914359
T6-T5   2.622951 -13.6618101  18.9077118 0.9991161

$Colutorio
                     diff         lwr       upr     p adj
CPC+Cy-CPC      -5.882418 -16.9321696  5.167334 0.5167936
Cymenol-CPC      9.930357  -0.5322569 20.392971 0.0698592
Placebo-CPC     23.231250  12.7686359 33.693864 0.0000001
Cymenol-CPC+Cy  15.812775   4.7630227 26.862527 0.0014420
Placebo-CPC+Cy  29.113668  18.0639156 40.163420 0.0000000
Placebo-Cymenol 13.300893   2.8382788 23.763507 0.0061939

$`Tiempo:Colutorio`
                               diff         lwr        upr     p adj
T1:CPC-T0:CPC         -5.775000e+01 -97.8998446 -17.600155 0.0000454
T2:CPC-T0:CPC         -4.685000e+01 -86.9998446  -6.700155 0.0049628
T3:CPC-T0:CPC         -4.151250e+01 -81.6623446  -1.362655 0.0326020
T4:CPC-T0:CPC         -3.625000e+01 -76.3998446   3.899845 0.1482983
T5:CPC-T0:CPC         -1.706250e+01 -57.2123446  23.087345 0.9988243
T6:CPC-T0:CPC         -1.157500e+01 -51.7248446  28.574845 0.9999993
T0:CPC+Cy-T0:CPC      -7.105427e-14 -42.4029619  42.402962 1.0000000
T1:CPC+Cy-T0:CPC      -5.279231e+01 -95.1952696 -10.389346 0.0014307
T2:CPC+Cy-T0:CPC      -5.166154e+01 -94.0645004  -9.258577 0.0022105
T3:CPC+Cy-T0:CPC      -5.211538e+01 -94.5183465  -9.712423 0.0018587
T4:CPC+Cy-T0:CPC      -3.870000e+01 -81.1029619   3.702962 0.1341981
T5:CPC+Cy-T0:CPC      -2.580000e+01 -68.2029619  16.602962 0.8811843
T6:CPC+Cy-T0:CPC      -3.110769e+01 -73.5106542  11.295270 0.5578368
T0:Cymenol-T0:CPC      8.526513e-14 -40.1498446  40.149845 1.0000000
T1:Cymenol-T0:CPC     -2.991250e+01 -70.0623446  10.237345 0.5233366
T2:Cymenol-T0:CPC     -2.816250e+01 -68.3123446  11.987345 0.6543733
T3:Cymenol-T0:CPC     -2.827500e+01 -68.4248446  11.874845 0.6461221
T4:Cymenol-T0:CPC     -2.918750e+01 -69.3373446  10.962345 0.5779992
T5:Cymenol-T0:CPC     -1.861250e+01 -58.7623446  21.537345 0.9954516
T6:Cymenol-T0:CPC     -7.337500e+00 -47.4873446  32.812345 1.0000000
T0:Placebo-T0:CPC     -2.842171e-14 -40.1498446  40.149845 1.0000000
T1:Placebo-T0:CPC     -9.937500e+00 -50.0873446  30.212345 1.0000000
T2:Placebo-T0:CPC     -1.400000e+01 -54.1498446  26.149845 0.9999659
T3:Placebo-T0:CPC     -1.293125e+01 -53.0810946  27.218595 0.9999930
T4:Placebo-T0:CPC     -3.375000e+00 -43.5248446  36.774845 1.0000000
T5:Placebo-T0:CPC     -2.843750e+00 -42.9935946  37.306095 1.0000000
T6:Placebo-T0:CPC     -5.293750e+00 -45.4435946  34.856095 1.0000000
T2:CPC-T1:CPC          1.090000e+01 -29.2498446  51.049845 0.9999998
T3:CPC-T1:CPC          1.623750e+01 -23.9123446  56.387345 0.9994874
T4:CPC-T1:CPC          2.150000e+01 -18.6498446  61.649845 0.9688502
T5:CPC-T1:CPC          4.068750e+01   0.5376554  80.837345 0.0423541
T6:CPC-T1:CPC          4.617500e+01   6.0251554  86.324845 0.0064043
T0:CPC+Cy-T1:CPC       5.775000e+01  15.3470381 100.152962 0.0001883
T1:CPC+Cy-T1:CPC       4.957692e+00 -37.4452696  47.360654 1.0000000
T2:CPC+Cy-T1:CPC       6.088462e+00 -36.3145004  48.491423 1.0000000
T3:CPC+Cy-T1:CPC       5.634615e+00 -36.7683465  48.037577 1.0000000
T4:CPC+Cy-T1:CPC       1.905000e+01 -23.3529619  61.452962 0.9971636
T5:CPC+Cy-T1:CPC       3.195000e+01 -10.4529619  74.352962 0.4978809
T6:CPC+Cy-T1:CPC       2.664231e+01 -15.7606542  69.045270 0.8427585
T0:Cymenol-T1:CPC      5.775000e+01  17.6001554  97.899845 0.0000454
T1:Cymenol-T1:CPC      2.783750e+01 -12.3123446  67.987345 0.6779376
T2:Cymenol-T1:CPC      2.958750e+01 -10.5623446  69.737345 0.5478199
T3:Cymenol-T1:CPC      2.947500e+01 -10.6748446  69.624845 0.5563091
T4:Cymenol-T1:CPC      2.856250e+01 -11.5873446  68.712345 0.6248545
T5:Cymenol-T1:CPC      3.913750e+01  -1.0123446  79.287345 0.0676505
T6:Cymenol-T1:CPC      5.041250e+01  10.2626554  90.562345 0.0012000
T0:Placebo-T1:CPC      5.775000e+01  17.6001554  97.899845 0.0000454
T1:Placebo-T1:CPC      4.781250e+01   7.6626554  87.962345 0.0034227
T2:Placebo-T1:CPC      4.375000e+01   3.6001554  83.899845 0.0153896
T3:Placebo-T1:CPC      4.481875e+01   4.6689054  84.968595 0.0105389
T4:Placebo-T1:CPC      5.437500e+01  14.2251554  94.524845 0.0002162
T5:Placebo-T1:CPC      5.490625e+01  14.7564054  95.056095 0.0001701
T6:Placebo-T1:CPC      5.245625e+01  12.3064054  92.606095 0.0005042
T3:CPC-T2:CPC          5.337500e+00 -34.8123446  45.487345 1.0000000
T4:CPC-T2:CPC          1.060000e+01 -29.5498446  50.749845 0.9999999
T5:CPC-T2:CPC          2.978750e+01 -10.3623446  69.937345 0.5327425
T6:CPC-T2:CPC          3.527500e+01  -4.8748446  75.424845 0.1880155
T0:CPC+Cy-T2:CPC       4.685000e+01   4.4470381  89.252962 0.0124152
T1:CPC+Cy-T2:CPC      -5.942308e+00 -48.3452696  36.460654 1.0000000
T2:CPC+Cy-T2:CPC      -4.811538e+00 -47.2145004  37.591423 1.0000000
T3:CPC+Cy-T2:CPC      -5.265385e+00 -47.6683465  37.137577 1.0000000
T4:CPC+Cy-T2:CPC       8.150000e+00 -34.2529619  50.552962 1.0000000
T5:CPC+Cy-T2:CPC       2.105000e+01 -21.3529619  63.452962 0.9880547
T6:CPC+Cy-T2:CPC       1.574231e+01 -26.6606542  58.145270 0.9998885
T0:Cymenol-T2:CPC      4.685000e+01   6.7001554  86.999845 0.0049628
T1:Cymenol-T2:CPC      1.693750e+01 -23.2123446  57.087345 0.9989578
T2:Cymenol-T2:CPC      1.868750e+01 -21.4623446  58.837345 0.9951741
T3:Cymenol-T2:CPC      1.857500e+01 -21.5748446  58.724845 0.9955851
T4:Cymenol-T2:CPC      1.766250e+01 -22.4873446  57.812345 0.9979539
T5:Cymenol-T2:CPC      2.823750e+01 -11.9123446  68.387345 0.6488774
T6:Cymenol-T2:CPC      3.951250e+01  -0.6373446  79.662345 0.0605781
T0:Placebo-T2:CPC      4.685000e+01   6.7001554  86.999845 0.0049628
T1:Placebo-T2:CPC      3.691250e+01  -3.2373446  77.062345 0.1251800
T2:Placebo-T2:CPC      3.285000e+01  -7.2998446  72.999845 0.3175733
T3:Placebo-T2:CPC      3.391875e+01  -6.2310946  74.068595 0.2550958
T4:Placebo-T2:CPC      4.347500e+01   3.3251554  83.624845 0.0169296
T5:Placebo-T2:CPC      4.400625e+01   3.8564054  84.156095 0.0140702
T6:Placebo-T2:CPC      4.155625e+01   1.4064054  81.706095 0.0321452
T4:CPC-T3:CPC          5.262500e+00 -34.8873446  45.412345 1.0000000
T5:CPC-T3:CPC          2.445000e+01 -15.6998446  64.599845 0.8802590
T6:CPC-T3:CPC          2.993750e+01 -10.2123446  70.087345 0.5214576
T0:CPC+Cy-T3:CPC       4.151250e+01  -0.8904619  83.915462 0.0643775
T1:CPC+Cy-T3:CPC      -1.127981e+01 -53.6827696  31.123154 0.9999999
T2:CPC+Cy-T3:CPC      -1.014904e+01 -52.5520004  32.253923 1.0000000
T3:CPC+Cy-T3:CPC      -1.060288e+01 -53.0058465  31.800077 1.0000000
T4:CPC+Cy-T3:CPC       2.812500e+00 -39.5904619  45.215462 1.0000000
T5:CPC+Cy-T3:CPC       1.571250e+01 -26.6904619  58.115462 0.9998924
T6:CPC+Cy-T3:CPC       1.040481e+01 -31.9981542  52.807770 1.0000000
T0:Cymenol-T3:CPC      4.151250e+01   1.3626554  81.662345 0.0326020
T1:Cymenol-T3:CPC      1.160000e+01 -28.5498446  51.749845 0.9999993
T2:Cymenol-T3:CPC      1.335000e+01 -26.7998446  53.499845 0.9999866
T3:Cymenol-T3:CPC      1.323750e+01 -26.9123446  53.387345 0.9999887
T4:Cymenol-T3:CPC      1.232500e+01 -27.8248446  52.474845 0.9999974
T5:Cymenol-T3:CPC      2.290000e+01 -17.2498446  63.049845 0.9369108
T6:Cymenol-T3:CPC      3.417500e+01  -5.9748446  74.324845 0.2413571
T0:Placebo-T3:CPC      4.151250e+01   1.3626554  81.662345 0.0326020
T1:Placebo-T3:CPC      3.157500e+01  -8.5748446  71.724845 0.4019703
T2:Placebo-T3:CPC      2.751250e+01 -12.6373446  67.662345 0.7010227
T3:Placebo-T3:CPC      2.858125e+01 -11.5685946  68.731095 0.6234596
T4:Placebo-T3:CPC      3.813750e+01  -2.0123446  78.287345 0.0899740
T5:Placebo-T3:CPC      3.866875e+01  -1.4810946  78.818595 0.0774574
T6:Placebo-T3:CPC      3.621875e+01  -3.9310946  76.368595 0.1494643
T5:CPC-T4:CPC          1.918750e+01 -20.9623446  59.337345 0.9929352
T6:CPC-T4:CPC          2.467500e+01 -15.4748446  64.824845 0.8700186
T0:CPC+Cy-T4:CPC       3.625000e+01  -6.1529619  78.652962 0.2336045
T1:CPC+Cy-T4:CPC      -1.654231e+01 -58.9452696  25.860654 0.9997265
T2:CPC+Cy-T4:CPC      -1.541154e+01 -57.8145004  26.991423 0.9999249
T3:CPC+Cy-T4:CPC      -1.586538e+01 -58.2683465  26.537577 0.9998713
T4:CPC+Cy-T4:CPC      -2.450000e+00 -44.8529619  39.952962 1.0000000
T5:CPC+Cy-T4:CPC       1.045000e+01 -31.9529619  52.852962 1.0000000
T6:CPC+Cy-T4:CPC       5.142308e+00 -37.2606542  47.545270 1.0000000
T0:Cymenol-T4:CPC      3.625000e+01  -3.8998446  76.399845 0.1482983
T1:Cymenol-T4:CPC      6.337500e+00 -33.8123446  46.487345 1.0000000
T2:Cymenol-T4:CPC      8.087500e+00 -32.0623446  48.237345 1.0000000
T3:Cymenol-T4:CPC      7.975000e+00 -32.1748446  48.124845 1.0000000
T4:Cymenol-T4:CPC      7.062500e+00 -33.0873446  47.212345 1.0000000
T5:Cymenol-T4:CPC      1.763750e+01 -22.5123446  57.787345 0.9979990
T6:Cymenol-T4:CPC      2.891250e+01 -11.2373446  69.062345 0.5986912
T0:Placebo-T4:CPC      3.625000e+01  -3.8998446  76.399845 0.1482983
T1:Placebo-T4:CPC      2.631250e+01 -13.8373446  66.462345 0.7806525
T2:Placebo-T4:CPC      2.225000e+01 -17.8998446  62.399845 0.9537965
T3:Placebo-T4:CPC      2.331875e+01 -16.8310946  63.468595 0.9239568
T4:Placebo-T4:CPC      3.287500e+01  -7.2748446  73.024845 0.3160200
T5:Placebo-T4:CPC      3.340625e+01  -6.7435946  73.556095 0.2840336
T6:Placebo-T4:CPC      3.095625e+01  -9.1935946  71.106095 0.4460352
T6:CPC-T5:CPC          5.487500e+00 -34.6623446  45.637345 1.0000000
T0:CPC+Cy-T5:CPC       1.706250e+01 -25.3404619  59.465462 0.9995300
T1:CPC+Cy-T5:CPC      -3.572981e+01 -78.1327696   6.673154 0.2598981
T2:CPC+Cy-T5:CPC      -3.459904e+01 -77.0020004   7.803923 0.3231630
T3:CPC+Cy-T5:CPC      -3.505288e+01 -77.4558465   7.350077 0.2967960
T4:CPC+Cy-T5:CPC      -2.163750e+01 -64.0404619  20.765462 0.9828446
T5:CPC+Cy-T5:CPC      -8.737500e+00 -51.1404619  33.665462 1.0000000
T6:CPC+Cy-T5:CPC      -1.404519e+01 -56.4481542  28.357770 0.9999876
T0:Cymenol-T5:CPC      1.706250e+01 -23.0873446  57.212345 0.9988243
T1:Cymenol-T5:CPC     -1.285000e+01 -52.9998446  27.299845 0.9999938
T2:Cymenol-T5:CPC     -1.110000e+01 -51.2498446  29.049845 0.9999997
T3:Cymenol-T5:CPC     -1.121250e+01 -51.3623446  28.937345 0.9999997
T4:Cymenol-T5:CPC     -1.212500e+01 -52.2748446  28.024845 0.9999981
T5:Cymenol-T5:CPC     -1.550000e+00 -41.6998446  38.599845 1.0000000
T6:Cymenol-T5:CPC      9.725000e+00 -30.4248446  49.874845 1.0000000
T0:Placebo-T5:CPC      1.706250e+01 -23.0873446  57.212345 0.9988243
T1:Placebo-T5:CPC      7.125000e+00 -33.0248446  47.274845 1.0000000
T2:Placebo-T5:CPC      3.062500e+00 -37.0873446  43.212345 1.0000000
T3:Placebo-T5:CPC      4.131250e+00 -36.0185946  44.281095 1.0000000
T4:Placebo-T5:CPC      1.368750e+01 -26.4623446  53.837345 0.9999781
T5:Placebo-T5:CPC      1.421875e+01 -25.9310946  54.368595 0.9999541
T6:Placebo-T5:CPC      1.176875e+01 -28.3810946  51.918595 0.9999990
T0:CPC+Cy-T6:CPC       1.157500e+01 -30.8279619  53.977962 0.9999998
T1:CPC+Cy-T6:CPC      -4.121731e+01 -83.6202696   1.185654 0.0698590
T2:CPC+Cy-T6:CPC      -4.008654e+01 -82.4895004   2.316423 0.0946001
T3:CPC+Cy-T6:CPC      -4.054038e+01 -82.9433465   1.862577 0.0839207
T4:CPC+Cy-T6:CPC      -2.712500e+01 -69.5279619  15.277962 0.8179261
T5:CPC+Cy-T6:CPC      -1.422500e+01 -56.6279619  28.177962 0.9999840
T6:CPC+Cy-T6:CPC      -1.953269e+01 -61.9356542  22.870270 0.9958597
T0:Cymenol-T6:CPC      1.157500e+01 -28.5748446  51.724845 0.9999993
T1:Cymenol-T6:CPC     -1.833750e+01 -58.4873446  21.812345 0.9963564
T2:Cymenol-T6:CPC     -1.658750e+01 -56.7373446  23.562345 0.9992636
T3:Cymenol-T6:CPC     -1.670000e+01 -56.8498446  23.449845 0.9991753
T4:Cymenol-T6:CPC     -1.761250e+01 -57.7623446  22.537345 0.9980433
T5:Cymenol-T6:CPC     -7.037500e+00 -47.1873446  33.112345 1.0000000
T6:Cymenol-T6:CPC      4.237500e+00 -35.9123446  44.387345 1.0000000
T0:Placebo-T6:CPC      1.157500e+01 -28.5748446  51.724845 0.9999993
T1:Placebo-T6:CPC      1.637500e+00 -38.5123446  41.787345 1.0000000
T2:Placebo-T6:CPC     -2.425000e+00 -42.5748446  37.724845 1.0000000
T3:Placebo-T6:CPC     -1.356250e+00 -41.5060946  38.793595 1.0000000
T4:Placebo-T6:CPC      8.200000e+00 -31.9498446  48.349845 1.0000000
T5:Placebo-T6:CPC      8.731250e+00 -31.4185946  48.881095 1.0000000
T6:Placebo-T6:CPC      6.281250e+00 -33.8685946  46.431095 1.0000000
T1:CPC+Cy-T0:CPC+Cy   -5.279231e+01 -97.3345611  -8.250054 0.0037380
T2:CPC+Cy-T0:CPC+Cy   -5.166154e+01 -96.2037918  -7.119285 0.0055257
T3:CPC+Cy-T0:CPC+Cy   -5.211538e+01 -96.6576380  -7.573131 0.0047293
T4:CPC+Cy-T0:CPC+Cy   -3.870000e+01 -83.2422534   5.842253 0.2059408
T5:CPC+Cy-T0:CPC+Cy   -2.580000e+01 -70.3422534  18.742253 0.9260120
T6:CPC+Cy-T0:CPC+Cy   -3.110769e+01 -75.6499457  13.434561 0.6632986
T0:Cymenol-T0:CPC+Cy   1.563194e-13 -42.4029619  42.402962 1.0000000
T1:Cymenol-T0:CPC+Cy  -2.991250e+01 -72.3154619  12.490462 0.6425831
T2:Cymenol-T0:CPC+Cy  -2.816250e+01 -70.5654619  14.240462 0.7582485
T3:Cymenol-T0:CPC+Cy  -2.827500e+01 -70.6779619  14.127962 0.7513117
T4:Cymenol-T0:CPC+Cy  -2.918750e+01 -71.5904619  13.215462 0.6922704
T5:Cymenol-T0:CPC+Cy  -1.861250e+01 -61.0154619  23.790462 0.9980239
T6:Cymenol-T0:CPC+Cy  -7.337500e+00 -49.7404619  35.065462 1.0000000
T0:Placebo-T0:CPC+Cy   4.263256e-14 -42.4029619  42.402962 1.0000000
T1:Placebo-T0:CPC+Cy  -9.937500e+00 -52.3404619  32.465462 1.0000000
T2:Placebo-T0:CPC+Cy  -1.400000e+01 -56.4029619  28.402962 0.9999884
T3:Placebo-T0:CPC+Cy  -1.293125e+01 -55.3342119  29.471712 0.9999977
T4:Placebo-T0:CPC+Cy  -3.375000e+00 -45.7779619  39.027962 1.0000000
T5:Placebo-T0:CPC+Cy  -2.843750e+00 -45.2467119  39.559212 1.0000000
T6:Placebo-T0:CPC+Cy  -5.293750e+00 -47.6967119  37.109212 1.0000000
T2:CPC+Cy-T1:CPC+Cy    1.130769e+00 -43.4114841  45.673023 1.0000000
T3:CPC+Cy-T1:CPC+Cy    6.769231e-01 -43.8653303  45.219176 1.0000000
T4:CPC+Cy-T1:CPC+Cy    1.409231e+01 -30.4499457  58.634561 0.9999951
T5:CPC+Cy-T1:CPC+Cy    2.699231e+01 -17.5499457  71.534561 0.8854868
T6:CPC+Cy-T1:CPC+Cy    2.168462e+01 -22.8576380  66.226869 0.9908437
T0:Cymenol-T1:CPC+Cy   5.279231e+01  10.3893458  95.195270 0.0014307
T1:Cymenol-T1:CPC+Cy   2.287981e+01 -19.5231542  65.282770 0.9659312
T2:Cymenol-T1:CPC+Cy   2.462981e+01 -17.7731542  67.032770 0.9238799
T3:Cymenol-T1:CPC+Cy   2.451731e+01 -17.8856542  66.920270 0.9273378
T4:Cymenol-T1:CPC+Cy   2.360481e+01 -18.7981542  66.007770 0.9514295
T5:Cymenol-T1:CPC+Cy   3.417981e+01  -8.2231542  76.582770 0.3486105
T6:Cymenol-T1:CPC+Cy   4.545481e+01   3.0518458  87.857770 0.0196553
T0:Placebo-T1:CPC+Cy   5.279231e+01  10.3893458  95.195270 0.0014307
T1:Placebo-T1:CPC+Cy   4.285481e+01   0.4518458  85.257770 0.0438275
T2:Placebo-T1:CPC+Cy   3.879231e+01  -3.6106542  81.195270 0.1312127
T3:Placebo-T1:CPC+Cy   3.986106e+01  -2.5419042  82.264020 0.1003042
T4:Placebo-T1:CPC+Cy   4.941731e+01   7.0143458  91.820270 0.0050746
T5:Placebo-T1:CPC+Cy   4.994856e+01   7.5455958  92.351520 0.0041851
T6:Placebo-T1:CPC+Cy   4.749856e+01   5.0955958  89.901520 0.0099622
T3:CPC+Cy-T2:CPC+Cy   -4.538462e-01 -44.9960995  44.088407 1.0000000
T4:CPC+Cy-T2:CPC+Cy    1.296154e+01 -31.5807149  57.503792 0.9999992
T5:CPC+Cy-T2:CPC+Cy    2.586154e+01 -18.6807149  70.403792 0.9242029
T6:CPC+Cy-T2:CPC+Cy    2.055385e+01 -23.9884072  65.096100 0.9957513
T0:Cymenol-T2:CPC+Cy   5.166154e+01   9.2585766  94.064500 0.0022105
T1:Cymenol-T2:CPC+Cy   2.174904e+01 -20.6539234  64.152000 0.9816760
T2:Cymenol-T2:CPC+Cy   2.349904e+01 -18.9039234  65.902000 0.9537873
T3:Cymenol-T2:CPC+Cy   2.338654e+01 -19.0164234  65.789500 0.9562016
T4:Cymenol-T2:CPC+Cy   2.247404e+01 -19.9289234  64.877000 0.9724613
T5:Cymenol-T2:CPC+Cy   3.304904e+01  -9.3539234  75.452000 0.4218474
T6:Cymenol-T2:CPC+Cy   4.432404e+01   1.9210766  86.727000 0.0281078
T0:Placebo-T2:CPC+Cy   5.166154e+01   9.2585766  94.064500 0.0022105
T1:Placebo-T2:CPC+Cy   4.172404e+01  -0.6789234  84.127000 0.0606772
T2:Placebo-T2:CPC+Cy   3.766154e+01  -4.7414234  80.064500 0.1714837
T3:Placebo-T2:CPC+Cy   3.873029e+01  -3.6726734  81.133250 0.1332127
T4:Placebo-T2:CPC+Cy   4.828654e+01   5.8835766  90.689500 0.0075832
T5:Placebo-T2:CPC+Cy   4.881779e+01   6.4148266  91.220750 0.0062883
T6:Placebo-T2:CPC+Cy   4.636779e+01   3.9648266  88.770750 0.0145838
T4:CPC+Cy-T3:CPC+Cy    1.341538e+01 -31.1268687  57.957638 0.9999982
T5:CPC+Cy-T3:CPC+Cy    2.631538e+01 -18.2268687  70.857638 0.9099214
T6:CPC+Cy-T3:CPC+Cy    2.100769e+01 -23.5345611  65.549946 0.9941502
T0:Cymenol-T3:CPC+Cy   5.211538e+01   9.7124227  94.518347 0.0018587
T1:Cymenol-T3:CPC+Cy   2.220288e+01 -20.2000773  64.605847 0.9762542
T2:Cymenol-T3:CPC+Cy   2.395288e+01 -18.4500773  66.355847 0.9430490
T3:Cymenol-T3:CPC+Cy   2.384038e+01 -18.5625773  66.243347 0.9458636
T4:Cymenol-T3:CPC+Cy   2.292788e+01 -19.4750773  65.330847 0.9650856
T5:Cymenol-T3:CPC+Cy   3.350288e+01  -8.9000773  75.905847 0.3917181
T6:Cymenol-T3:CPC+Cy   4.477788e+01   2.3749227  87.180847 0.0243881
T0:Placebo-T3:CPC+Cy   5.211538e+01   9.7124227  94.518347 0.0018587
T1:Placebo-T3:CPC+Cy   4.217788e+01  -0.2250773  84.580847 0.0533451
T2:Placebo-T3:CPC+Cy   3.811538e+01  -4.2875773  80.518347 0.1543380
T3:Placebo-T3:CPC+Cy   3.918413e+01  -3.2188273  81.587097 0.1191125
T4:Placebo-T3:CPC+Cy   4.874038e+01   6.3374227  91.143347 0.0064632
T5:Placebo-T3:CPC+Cy   4.927163e+01   6.8686727  91.674597 0.0053477
T6:Placebo-T3:CPC+Cy   4.682163e+01   4.4186727  89.224597 0.0125341
T5:CPC+Cy-T4:CPC+Cy    1.290000e+01 -31.6422534  57.442253 0.9999992
T6:CPC+Cy-T4:CPC+Cy    7.592308e+00 -36.9499457  52.134561 1.0000000
T0:Cymenol-T4:CPC+Cy   3.870000e+01  -3.7029619  81.102962 0.1341981
T1:Cymenol-T4:CPC+Cy   8.787500e+00 -33.6154619  51.190462 1.0000000
T2:Cymenol-T4:CPC+Cy   1.053750e+01 -31.8654619  52.940462 1.0000000
T3:Cymenol-T4:CPC+Cy   1.042500e+01 -31.9779619  52.827962 1.0000000
T4:Cymenol-T4:CPC+Cy   9.512500e+00 -32.8904619  51.915462 1.0000000
T5:Cymenol-T4:CPC+Cy   2.008750e+01 -22.3154619  62.490462 0.9937652
 [ reached getOption("max.print") -- omitted 128 rows ]
aov2 <- lm(Reduccion ~ Tiempo:Colutorio, 
                 data = df)
summary(aov2)

Call:
lm(formula = Reduccion ~ Tiempo:Colutorio, data = df)

Residuals:
    Min      1Q  Median      3Q     Max 
-132.19  -13.69    0.00   12.05  144.41 

Coefficients: (1 not defined because of singularities)
                          Estimate Std. Error t value Pr(>|t|)    
(Intercept)                 94.706      7.587  12.483  < 2e-16 ***
TiempoT0:ColutorioCPC        5.294     10.730   0.493 0.622020    
TiempoT1:ColutorioCPC      -52.456     10.730  -4.889 1.47e-06 ***
TiempoT2:ColutorioCPC      -41.556     10.730  -3.873 0.000126 ***
TiempoT3:ColutorioCPC      -36.219     10.730  -3.376 0.000809 ***
TiempoT4:ColutorioCPC      -30.956     10.730  -2.885 0.004125 ** 
TiempoT5:ColutorioCPC      -11.769     10.730  -1.097 0.273374    
TiempoT6:ColutorioCPC       -6.281     10.730  -0.585 0.558605    
TiempoT0:ColutorioCPC+Cy     5.294     11.332   0.467 0.640642    
TiempoT1:ColutorioCPC+Cy   -47.499     11.332  -4.192 3.41e-05 ***
TiempoT2:ColutorioCPC+Cy   -46.368     11.332  -4.092 5.18e-05 ***
TiempoT3:ColutorioCPC+Cy   -46.822     11.332  -4.132 4.39e-05 ***
TiempoT4:ColutorioCPC+Cy   -33.406     11.332  -2.948 0.003386 ** 
TiempoT5:ColutorioCPC+Cy   -20.506     11.332  -1.810 0.071108 .  
TiempoT6:ColutorioCPC+Cy   -25.814     11.332  -2.278 0.023255 *  
TiempoT0:ColutorioCymenol    5.294     10.730   0.493 0.622020    
TiempoT1:ColutorioCymenol  -24.619     10.730  -2.294 0.022283 *  
TiempoT2:ColutorioCymenol  -22.869     10.730  -2.131 0.033671 *  
TiempoT3:ColutorioCymenol  -22.981     10.730  -2.142 0.032812 *  
TiempoT4:ColutorioCymenol  -23.894     10.730  -2.227 0.026513 *  
TiempoT5:ColutorioCymenol  -13.319     10.730  -1.241 0.215226    
TiempoT6:ColutorioCymenol   -2.044     10.730  -0.190 0.849033    
TiempoT0:ColutorioPlacebo    5.294     10.730   0.493 0.622020    
TiempoT1:ColutorioPlacebo   -4.644     10.730  -0.433 0.665399    
TiempoT2:ColutorioPlacebo   -8.706     10.730  -0.811 0.417611    
TiempoT3:ColutorioPlacebo   -7.637     10.730  -0.712 0.476999    
TiempoT4:ColutorioPlacebo    1.919     10.730   0.179 0.858165    
TiempoT5:ColutorioPlacebo    2.450     10.730   0.228 0.819501    
TiempoT6:ColutorioPlacebo       NA         NA      NA       NA    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 30.35 on 399 degrees of freedom
Multiple R-squared:  0.2742,    Adjusted R-squared:  0.225 
F-statistic: 5.582 on 27 and 399 DF,  p-value: 9.674e-16

#plot los residuos

plot(aov1, 1)

#Normal Q-Q

plot(aov1, 2)

#prueba de shapiro-wilk normality test

# Extract the residuals
aov_residuals <- residuals(object = aov1)
# Run Shapiro-Wilk test
shapiro.test(x = aov_residuals )

    Shapiro-Wilk normality test

data:  aov_residuals
W = 0.93659, p-value = 1.645e-12
performance::check_model(aov1)

#ANALIZO PARA MÉTODO DE CONFOCAL #abro el df

df1 <- read_csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vQ4Nk96Jx1n-pLC1e7lXlUgMtqAFc1xOTyNCP6_72cvgtcZXjITM2-wocUbgKWin63eOHWZP3ezwJoB/pub?gid=1791723434&single=true&output=csv")

── Column specification ──────────────────────────────────────────────────────────────────────
cols(
  id = col_double(),
  Colutorio = col_character(),
  Tiempo = col_character(),
  Reduccion = col_double()
)

#cambio colutorio como factor

df1 <- df1 %>% 
  mutate(Colutorio = as.factor(Colutorio))

#glimse

glimpse(df1)
Rows: 378
Columns: 4
$ id        <dbl> 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4,…
$ Colutorio <fct> CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC, CPC,…
$ Tiempo    <chr> "T0", "T1", "T2", "T3", "T4", "T5", "T6", "T0", "T1", "T2", "T3", "T4", "T…
$ Reduccion <dbl> 100.0, 43.9, 63.4, 71.4, 85.6, 71.8, 73.9, 100.0, 55.4, 50.0, 54.9, 76.0, …

#head df1

head(df1)

#hist confocal

df1 %>% 
  # pivot_longer(Colutorio:Tiempo, 
  #             names_to = "Variables", 
  #             values_to = "Values") %>% 
  ggplot(aes(x = Reduccion)) + 
  geom_histogram() + 
  facet_grid(Tiempo ~ Colutorio)

#grafica lineal confocal

df1 %>% 
ggline(x = "Tiempo", 
       y = "Reduccion", 
       color = "Colutorio",
       add = c("mean_se", "dotplot")) + 
  labs(title = "Profile plot de reducción en tiempo y tipo de colutorio") + 
  theme_pubr()

#bonferroni confocal

pairwise.t.test( df1$Reduccion, 
                 df1$Colutorio,
           
                 p.adj = "bonferroni")

    Pairwise comparisons using t tests with pooled SD 

data:  df1$Reduccion and df1$Colutorio 

            CPC     CPC+Cymenol Cymenol
CPC+Cymenol 0.16    -           -      
Cymenol     2.9e-07 2.0e-13     -      
Placebo     2.1e-12 < 2e-16     0.13   

P value adjustment method: bonferroni 

#modelo de aov para confocal

aov3 <- aov(Reduccion ~ Tiempo + Colutorio + Tiempo*Colutorio , 
           data = df1)
summary(aov3)
                  Df Sum Sq Mean Sq F value   Pr(>F)    
Tiempo             6  49689    8281  26.584  < 2e-16 ***
Colutorio          3  59639   19880  63.815  < 2e-16 ***
Tiempo:Colutorio  18  17402     967   3.103 2.63e-05 ***
Residuals        350 109031     312                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
df1 %>% 
  with(aov(Reduccion ~ Tiempo + Colutorio + Tiempo:Colutorio + Error(id))) %>% 
  summary()

Error: id
          Df Sum Sq Mean Sq
Colutorio  1  510.5   510.5

Error: Within
                  Df Sum Sq Mean Sq F value  Pr(>F)    
Tiempo             6  49689    8281  26.548 < 2e-16 ***
Colutorio          3  59292   19764  63.359 < 2e-16 ***
Tiempo:Colutorio  18  17402     967   3.099 2.7e-05 ***
Residuals        349 108867     312                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

#tukey

TukeyHSD(aov3)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = Reduccion ~ Tiempo + Colutorio + Tiempo * Colutorio, data = df1)

$Tiempo
            diff         lwr        upr     p adj
T1-T0 -35.062963 -45.1362001 -24.989726 0.0000000
T2-T0 -33.362963 -43.4362001 -23.289726 0.0000000
T3-T0 -29.192593 -39.2658297 -19.119355 0.0000000
T4-T0 -26.116667 -36.1899038 -16.043430 0.0000000
T5-T0 -18.846296 -28.9195334  -8.773059 0.0000012
T6-T0 -14.320370 -24.3936075  -4.247133 0.0006254
T2-T1   1.700000  -8.3732371  11.773237 0.9988403
T3-T1   5.870370  -4.2028668  15.943607 0.5975137
T4-T1   8.946296  -1.1269408  19.019533 0.1190995
T5-T1  16.216667   6.1434295  26.289904 0.0000542
T6-T1  20.742593  10.6693555  30.815830 0.0000001
T3-T2   4.170370  -5.9028668  14.243607 0.8830146
T4-T2   7.246296  -2.8269408  17.319533 0.3354214
T5-T2  14.516667   4.4434295  24.589904 0.0004919
T6-T2  19.042593   8.9693555  29.115830 0.0000009
T4-T3   3.075926  -6.9973112  13.149163 0.9715558
T5-T3  10.346296   0.2730592  20.419533 0.0397034
T6-T3  14.872222   4.7989851  24.945459 0.0003159
T5-T4   7.270370  -2.8028668  17.343607 0.3313455
T6-T4  11.796296   1.7230592  21.869533 0.0103079
T6-T5   4.525926  -5.5473112  14.599163 0.8363046

$Colutorio
                         diff         lwr       upr     p adj
CPC+Cymenol-CPC     -6.780000 -13.1794040 -0.380596 0.0330540
Cymenol-CPC         17.273469  10.7646660 23.782273 0.0000000
Placebo-CPC         24.919481  17.9810727 31.857888 0.0000000
Cymenol-CPC+Cymenol 24.053469  17.6540654 30.452873 0.0000000
Placebo-CPC+Cymenol 31.699481  24.8635935 38.535368 0.0000000
Placebo-Cymenol      7.646011   0.7076034 14.584419 0.0242203

$`Tiempo:Colutorio`
                                       diff         lwr        upr     p adj
T1:CPC-T0:CPC                 -4.592857e+01 -70.9167112 -20.940432 0.0000000
T2:CPC-T0:CPC                 -4.377857e+01 -68.7667112 -18.790432 0.0000001
T3:CPC-T0:CPC                 -3.803571e+01 -63.0238540 -13.047575 0.0000093
T4:CPC-T0:CPC                 -3.788571e+01 -62.8738540 -12.897575 0.0000104
T5:CPC-T0:CPC                 -2.895714e+01 -53.9452826  -3.969003 0.0056623
T6:CPC-T0:CPC                 -1.601429e+01 -41.0024255   8.973854 0.8133308
T0:CPC+Cymenol-T0:CPC         -2.131628e-13 -24.5681411  24.568141 1.0000000
T1:CPC+Cymenol-T0:CPC         -4.609333e+01 -70.6614744 -21.525192 0.0000000
T2:CPC+Cymenol-T0:CPC         -4.764000e+01 -72.2081411 -23.071859 0.0000000
T3:CPC+Cymenol-T0:CPC         -5.181333e+01 -76.3814744 -27.245192 0.0000000
T4:CPC+Cymenol-T0:CPC         -4.867333e+01 -73.2414744 -24.105192 0.0000000
T5:CPC+Cymenol-T0:CPC         -3.404000e+01 -58.6081411  -9.471859 0.0001264
T6:CPC+Cymenol-T0:CPC         -2.980000e+01 -54.3681411  -5.231859 0.0024441
T0:Cymenol-T0:CPC             -2.415845e-13 -24.9881397  24.988140 1.0000000
T1:Cymenol-T0:CPC             -3.335000e+01 -58.3381397  -8.361860 0.0003158
T2:Cymenol-T0:CPC             -2.682857e+01 -51.8167112  -1.840432 0.0192781
T3:Cymenol-T0:CPC             -1.347143e+01 -38.4595683  11.516711 0.9655943
T4:Cymenol-T0:CPC             -7.114286e+00 -32.1024255  17.873854 0.9999994
T5:Cymenol-T0:CPC             -4.735714e+00 -29.7238540  20.252425 1.0000000
T6:Cymenol-T0:CPC             -4.185714e+00 -29.1738540  20.802425 1.0000000
T0:Placebo-T0:CPC             -2.700062e-13 -26.6374465  26.637446 1.0000000
T1:Placebo-T0:CPC             -8.372727e+00 -35.0101737  18.264719 0.9999955
T2:Placebo-T0:CPC             -8.954545e+00 -35.5919919  17.682901 0.9999825
T3:Placebo-T0:CPC             -7.100000e+00 -33.7374465  19.537446 0.9999999
T4:Placebo-T0:CPC             -4.563636e+00 -31.2010828  22.073810 1.0000000
T5:Placebo-T0:CPC             -3.218182e+00 -29.8556283  23.419265 1.0000000
T6:Placebo-T0:CPC             -3.954545e+00 -30.5919919  22.682901 1.0000000
T2:CPC-T1:CPC                  2.150000e+00 -22.8381397  27.138140 1.0000000
T3:CPC-T1:CPC                  7.892857e+00 -17.0952826  32.880997 0.9999950
T4:CPC-T1:CPC                  8.042857e+00 -16.9452826  33.030997 0.9999927
T5:CPC-T1:CPC                  1.697143e+01  -8.0167112  41.959568 0.7157926
T6:CPC-T1:CPC                  2.991429e+01   4.9261460  54.902425 0.0031389
T0:CPC+Cymenol-T1:CPC          4.592857e+01  21.3604303  70.496713 0.0000000
T1:CPC+Cymenol-T1:CPC         -1.647619e-01 -24.7329030  24.403379 1.0000000
T2:CPC+Cymenol-T1:CPC         -1.711429e+00 -26.2795697  22.856713 1.0000000
T3:CPC+Cymenol-T1:CPC         -5.884762e+00 -30.4529030  18.683379 1.0000000
T4:CPC+Cymenol-T1:CPC         -2.744762e+00 -27.3129030  21.823379 1.0000000
T5:CPC+Cymenol-T1:CPC          1.188857e+01 -12.6795697  36.456713 0.9913503
T6:CPC+Cymenol-T1:CPC          1.612857e+01  -8.4395697  40.696713 0.7758037
T0:Cymenol-T1:CPC              4.592857e+01  20.9404317  70.916711 0.0000000
T1:Cymenol-T1:CPC              1.257857e+01 -12.4095683  37.566711 0.9852757
T2:Cymenol-T1:CPC              1.910000e+01  -5.8881397  44.088140 0.4636794
T3:Cymenol-T1:CPC              3.245714e+01   7.4690031  57.445283 0.0005884
T4:Cymenol-T1:CPC              3.881429e+01  13.8261460  63.802425 0.0000049
T5:Cymenol-T1:CPC              4.119286e+01  16.2047174  66.180997 0.0000007
T6:Cymenol-T1:CPC              4.174286e+01  16.7547174  66.730997 0.0000004
T0:Placebo-T1:CPC              4.592857e+01  19.2911250  72.566018 0.0000001
T1:Placebo-T1:CPC              3.755584e+01  10.9183977  64.193291 0.0000807
T2:Placebo-T1:CPC              3.697403e+01  10.3365795  63.611472 0.0001207
T3:Placebo-T1:CPC              3.882857e+01  12.1911250  65.466018 0.0000327
T4:Placebo-T1:CPC              4.136494e+01  14.7274886  68.002382 0.0000050
T5:Placebo-T1:CPC              4.271039e+01  16.0729431  69.347836 0.0000018
T6:Placebo-T1:CPC              4.197403e+01  15.3365795  68.611472 0.0000031
T3:CPC-T2:CPC                  5.742857e+00 -19.2452826  30.730997 1.0000000
T4:CPC-T2:CPC                  5.892857e+00 -19.0952826  30.880997 1.0000000
T5:CPC-T2:CPC                  1.482143e+01 -10.1667112  39.809568 0.9050486
T6:CPC-T2:CPC                  2.776429e+01   2.7761460  52.752425 0.0114244
T0:CPC+Cymenol-T2:CPC          4.377857e+01  19.2104303  68.346713 0.0000000
T1:CPC+Cymenol-T2:CPC         -2.314762e+00 -26.8829030  22.253379 1.0000000
T2:CPC+Cymenol-T2:CPC         -3.861429e+00 -28.4295697  20.706713 1.0000000
T3:CPC+Cymenol-T2:CPC         -8.034762e+00 -32.6029030  16.533379 0.9999899
T4:CPC+Cymenol-T2:CPC         -4.894762e+00 -29.4629030  19.673379 1.0000000
T5:CPC+Cymenol-T2:CPC          9.738571e+00 -14.8295697  34.306713 0.9996232
T6:CPC+Cymenol-T2:CPC          1.397857e+01 -10.5895697  38.546713 0.9374757
T0:Cymenol-T2:CPC              4.377857e+01  18.7904317  68.766711 0.0000001
T1:Cymenol-T2:CPC              1.042857e+01 -14.5595683  35.416711 0.9990955
T2:Cymenol-T2:CPC              1.695000e+01  -8.0381397  41.938140 0.7181609
T3:Cymenol-T2:CPC              3.030714e+01   5.3190031  55.295283 0.0024477
T4:Cymenol-T2:CPC              3.666429e+01  11.6761460  61.652425 0.0000272
T5:Cymenol-T2:CPC              3.904286e+01  14.0547174  64.030997 0.0000041
T6:Cymenol-T2:CPC              3.959286e+01  14.6047174  64.580997 0.0000026
T0:Placebo-T2:CPC              4.377857e+01  17.1411250  70.416018 0.0000008
T1:Placebo-T2:CPC              3.540584e+01   8.7683977  62.043291 0.0003477
T2:Placebo-T2:CPC              3.482403e+01   8.1865795  61.461472 0.0005088
T3:Placebo-T2:CPC              3.667857e+01  10.0411250  63.316018 0.0001479
T4:Placebo-T2:CPC              3.921494e+01  12.5774886  65.852382 0.0000247
T5:Placebo-T2:CPC              4.056039e+01  13.9229431  67.197836 0.0000092
T6:Placebo-T2:CPC              3.982403e+01  13.1865795  66.461472 0.0000158
T4:CPC-T3:CPC                  1.500000e-01 -24.8381397  25.138140 1.0000000
T5:CPC-T3:CPC                  9.078571e+00 -15.9095683  34.066711 0.9999221
T6:CPC-T3:CPC                  2.202143e+01  -2.9667112  47.009568 0.1826458
T0:CPC+Cymenol-T3:CPC          3.803571e+01  13.4675732  62.603855 0.0000055
T1:CPC+Cymenol-T3:CPC         -8.057619e+00 -32.6257601  16.510522 0.9999893
T2:CPC+Cymenol-T3:CPC         -9.604286e+00 -34.1724268  14.963855 0.9997046
T3:CPC+Cymenol-T3:CPC         -1.377762e+01 -38.3457601  10.790522 0.9464158
T4:CPC+Cymenol-T3:CPC         -1.063762e+01 -35.2057601  13.930522 0.9983545
T5:CPC+Cymenol-T3:CPC          3.995714e+00 -20.5724268  28.563855 1.0000000
T6:CPC+Cymenol-T3:CPC          8.235714e+00 -16.3324268  32.803855 0.9999835
T0:Cymenol-T3:CPC              3.803571e+01  13.0475745  63.023854 0.0000093
T1:Cymenol-T3:CPC              4.685714e+00 -20.3024255  29.673854 1.0000000
T2:Cymenol-T3:CPC              1.120714e+01 -13.7809969  36.195283 0.9971474
T3:Cymenol-T3:CPC              2.456429e+01  -0.4238540  49.552425 0.0613274
T4:Cymenol-T3:CPC              3.092143e+01   5.9332888  55.909568 0.0016467
T5:Cymenol-T3:CPC              3.330000e+01   8.3118603  58.288140 0.0003272
T6:Cymenol-T3:CPC              3.385000e+01   8.8618603  58.838140 0.0002212
T0:Placebo-T3:CPC              3.803571e+01  11.3982678  64.673161 0.0000576
T1:Placebo-T3:CPC              2.966299e+01   3.0255405  56.300433 0.0110251
T2:Placebo-T3:CPC              2.908117e+01   2.4437224  55.718615 0.0150194
T3:Placebo-T3:CPC              3.093571e+01   4.2982678  57.573161 0.0054501
T4:Placebo-T3:CPC              3.347208e+01   6.8346315  60.109524 0.0012021
T5:Placebo-T3:CPC              3.481753e+01   8.1800860  61.454979 0.0005110
T6:Placebo-T3:CPC              3.408117e+01   7.4437224  60.718615 0.0008197
T5:CPC-T4:CPC                  8.928571e+00 -16.0595683  33.916711 0.9999431
T6:CPC-T4:CPC                  2.187143e+01  -3.1167112  46.859568 0.1932997
T0:CPC+Cymenol-T4:CPC          3.788571e+01  13.3175732  62.453855 0.0000062
T1:CPC+Cymenol-T4:CPC         -8.207619e+00 -32.7757601  16.360522 0.9999846
T2:CPC+Cymenol-T4:CPC         -9.754286e+00 -34.3224268  14.813855 0.9996125
T3:CPC+Cymenol-T4:CPC         -1.392762e+01 -38.4957601  10.640522 0.9398361
T4:CPC+Cymenol-T4:CPC         -1.078762e+01 -35.3557601  13.780522 0.9979465
T5:CPC+Cymenol-T4:CPC          3.845714e+00 -20.7224268  28.413855 1.0000000
T6:CPC+Cymenol-T4:CPC          8.085714e+00 -16.4824268  32.653855 0.9999885
T0:Cymenol-T4:CPC              3.788571e+01  12.8975745  62.873854 0.0000104
T1:Cymenol-T4:CPC              4.535714e+00 -20.4524255  29.523854 1.0000000
T2:Cymenol-T4:CPC              1.105714e+01 -13.9309969  36.045283 0.9976830
T3:Cymenol-T4:CPC              2.441429e+01  -0.5738540  49.402425 0.0658283
T4:Cymenol-T4:CPC              3.077143e+01   5.7832888  55.759568 0.0018155
T5:Cymenol-T4:CPC              3.315000e+01   8.1618603  58.138140 0.0003636
T6:Cymenol-T4:CPC              3.370000e+01   8.7118603  58.688140 0.0002463
T0:Placebo-T4:CPC              3.788571e+01  11.2482678  64.523161 0.0000640
T1:Placebo-T4:CPC              2.951299e+01   2.8755405  56.150433 0.0119493
T2:Placebo-T4:CPC              2.893117e+01   2.2937224  55.568615 0.0162434
T3:Placebo-T4:CPC              3.078571e+01   4.1482678  57.423161 0.0059335
T4:Placebo-T4:CPC              3.332208e+01   6.6846315  59.959524 0.0013195
T5:Placebo-T4:CPC              3.466753e+01   8.0300860  61.304979 0.0005631
T6:Placebo-T4:CPC              3.393117e+01   7.2937224  60.568615 0.0009014
T6:CPC-T5:CPC                  1.294286e+01 -12.0452826  37.930997 0.9787853
T0:CPC+Cymenol-T5:CPC          2.895714e+01   4.3890018  53.525284 0.0041853
T1:CPC+Cymenol-T5:CPC         -1.713619e+01 -41.7043316   7.431951 0.6636988
T2:CPC+Cymenol-T5:CPC         -1.868286e+01 -43.2509982   5.885284 0.4752165
T3:CPC+Cymenol-T5:CPC         -2.285619e+01 -47.4243316   1.711951 0.1112043
T4:CPC+Cymenol-T5:CPC         -1.971619e+01 -44.2843316   4.851951 0.3565582
T5:CPC+Cymenol-T5:CPC         -5.082857e+00 -29.6509982  19.485284 1.0000000
T6:CPC+Cymenol-T5:CPC         -8.428571e-01 -25.4109982  23.725284 1.0000000
T0:Cymenol-T5:CPC              2.895714e+01   3.9690031  53.945283 0.0056623
T1:Cymenol-T5:CPC             -4.392857e+00 -29.3809969  20.595283 1.0000000
T2:Cymenol-T5:CPC              2.128571e+00 -22.8595683  27.116711 1.0000000
T3:Cymenol-T5:CPC              1.548571e+01  -9.5024255  40.473854 0.8584713
T4:Cymenol-T5:CPC              2.184286e+01  -3.1452826  46.830997 0.1953782
T5:Cymenol-T5:CPC              2.422143e+01  -0.7667112  49.209568 0.0720225
T6:Cymenol-T5:CPC              2.477143e+01  -0.2167112  49.759568 0.0555441
T0:Placebo-T5:CPC              2.895714e+01   2.3196964  55.594589 0.0160251
T1:Placebo-T5:CPC              2.058442e+01  -6.0530309  47.221862 0.4391875
T2:Placebo-T5:CPC              2.000260e+01  -6.6348491  46.640044 0.5036001
T3:Placebo-T5:CPC              2.185714e+01  -4.7803036  48.494589 0.3104605
T4:Placebo-T5:CPC              2.439351e+01  -2.2439400  51.030953 0.1296269
T5:Placebo-T5:CPC              2.573896e+01  -0.8984854  52.376408 0.0745992
T6:Placebo-T5:CPC              2.500260e+01  -1.6348491  51.640044 0.1016664
T0:CPC+Cymenol-T6:CPC          1.601429e+01  -8.5538554  40.582427 0.7873370
T1:CPC+Cymenol-T6:CPC         -3.007905e+01 -54.6471887  -5.510907 0.0020373
T2:CPC+Cymenol-T6:CPC         -3.162571e+01 -56.1938554  -7.057573 0.0007179
T3:CPC+Cymenol-T6:CPC         -3.579905e+01 -60.3671887 -11.230907 0.0000330
T4:CPC+Cymenol-T6:CPC         -3.265905e+01 -57.2271887  -8.090907 0.0003468
T5:CPC+Cymenol-T6:CPC         -1.802571e+01 -42.5938554   6.542427 0.5554975
T6:CPC+Cymenol-T6:CPC         -1.378571e+01 -38.3538554  10.782427 0.9460746
T0:Cymenol-T6:CPC              1.601429e+01  -8.9738540  41.002425 0.8133308
T1:Cymenol-T6:CPC             -1.733571e+01 -42.3238540   7.652425 0.6745414
T2:Cymenol-T6:CPC             -1.081429e+01 -35.8024255  14.173854 0.9983671
T3:Cymenol-T6:CPC              2.542857e+00 -22.4452826  27.530997 1.0000000
T4:Cymenol-T6:CPC              8.900000e+00 -16.0881397  33.888140 0.9999464
T5:Cymenol-T6:CPC              1.127857e+01 -13.7095683  36.266711 0.9968571
T6:Cymenol-T6:CPC              1.182857e+01 -13.1595683  36.816711 0.9936562
T0:Placebo-T6:CPC              1.601429e+01 -10.6231608  42.651732 0.8921546
T1:Placebo-T6:CPC              7.641558e+00 -18.9958880  34.279005 0.9999993
T2:Placebo-T6:CPC              7.059740e+00 -19.5777062  33.697187 0.9999999
T3:Placebo-T6:CPC              8.914286e+00 -17.7231608  35.551732 0.9999840
T4:Placebo-T6:CPC              1.145065e+01 -15.1867971  38.088096 0.9985341
T5:Placebo-T6:CPC              1.279610e+01 -13.8413426  39.433550 0.9921879
T6:Placebo-T6:CPC              1.205974e+01 -14.5777062  38.697187 0.9967082
T1:CPC+Cymenol-T0:CPC+Cymenol -4.609333e+01 -70.2341698 -21.952497 0.0000000
T2:CPC+Cymenol-T0:CPC+Cymenol -4.764000e+01 -71.7808365 -23.499164 0.0000000
T3:CPC+Cymenol-T0:CPC+Cymenol -5.181333e+01 -75.9541698 -27.672497 0.0000000
T4:CPC+Cymenol-T0:CPC+Cymenol -4.867333e+01 -72.8141698 -24.532497 0.0000000
T5:CPC+Cymenol-T0:CPC+Cymenol -3.404000e+01 -58.1808365  -9.899164 0.0000804
T6:CPC+Cymenol-T0:CPC+Cymenol -2.980000e+01 -53.9408365  -5.659164 0.0017297
T0:Cymenol-T0:CPC+Cymenol     -2.842171e-14 -24.5681411  24.568141 1.0000000
T1:Cymenol-T0:CPC+Cymenol     -3.335000e+01 -57.9181411  -8.781859 0.0002104
T2:Cymenol-T0:CPC+Cymenol     -2.682857e+01 -51.3967125  -2.260430 0.0149636
T3:Cymenol-T0:CPC+Cymenol     -1.347143e+01 -38.0395697  11.096713 0.9582047
T4:Cymenol-T0:CPC+Cymenol     -7.114286e+00 -31.6824268  17.453855 0.9999992
T5:Cymenol-T0:CPC+Cymenol     -4.735714e+00 -29.3038554  19.832427 1.0000000
T6:Cymenol-T0:CPC+Cymenol     -4.185714e+00 -28.7538554  20.382427 1.0000000
T0:Placebo-T0:CPC+Cymenol     -5.684342e-14 -26.2438561  26.243856 1.0000000
T1:Placebo-T0:CPC+Cymenol     -8.372727e+00 -34.6165834  17.871129 0.9999939
T2:Placebo-T0:CPC+Cymenol     -8.954545e+00 -35.1984016  17.289311 0.9999766
T3:Placebo-T0:CPC+Cymenol     -7.100000e+00 -33.3438561  19.143856 0.9999998
T4:Placebo-T0:CPC+Cymenol     -4.563636e+00 -30.8074925  21.680220 1.0000000
T5:Placebo-T0:CPC+Cymenol     -3.218182e+00 -29.4620379  23.025674 1.0000000
T6:Placebo-T0:CPC+Cymenol     -3.954545e+00 -30.1984016  22.289311 1.0000000
T2:CPC+Cymenol-T1:CPC+Cymenol -1.546667e+00 -25.6875031  22.594170 1.0000000
T3:CPC+Cymenol-T1:CPC+Cymenol -5.720000e+00 -29.8608365  18.420836 1.0000000
T4:CPC+Cymenol-T1:CPC+Cymenol -2.580000e+00 -26.7208365  21.560836 1.0000000
T5:CPC+Cymenol-T1:CPC+Cymenol  1.205333e+01 -12.0875031  36.194170 0.9867750
T6:CPC+Cymenol-T1:CPC+Cymenol  1.629333e+01  -7.8475031  40.434170 0.7274621
T0:Cymenol-T1:CPC+Cymenol      4.609333e+01  21.5251922  70.661474 0.0000000
T1:Cymenol-T1:CPC+Cymenol      1.274333e+01 -11.8248078  37.311474 0.9784073
T2:Cymenol-T1:CPC+Cymenol      1.926476e+01  -5.3033792  43.832903 0.4067551
T3:Cymenol-T1:CPC+Cymenol      3.262190e+01   8.0537637  57.190046 0.0003561
T4:Cymenol-T1:CPC+Cymenol      3.897905e+01  14.4109065  63.547189 0.0000025
T5:Cymenol-T1:CPC+Cymenol      4.135762e+01  16.7894780  65.925760 0.0000003
T6:Cymenol-T1:CPC+Cymenol      4.190762e+01  17.3394780  66.475760 0.0000002
T0:Placebo-T1:CPC+Cymenol      4.609333e+01  19.8494772  72.337189 0.0000001
T1:Placebo-T1:CPC+Cymenol      3.772061e+01  11.4767500  63.964462 0.0000482
T2:Placebo-T1:CPC+Cymenol      3.713879e+01  10.8949318  63.382644 0.0000731
T3:Placebo-T1:CPC+Cymenol      3.899333e+01  12.7494772  65.237189 0.0000190
T4:Placebo-T1:CPC+Cymenol      4.152970e+01  15.2858409  67.773553 0.0000027
T5:Placebo-T1:CPC+Cymenol      4.287515e+01  16.6312954  69.119008 0.0000009
T6:Placebo-T1:CPC+Cymenol      4.213879e+01  15.8949318  68.382644 0.0000017
T3:CPC+Cymenol-T2:CPC+Cymenol -4.173333e+00 -28.3141698  19.967503 1.0000000
T4:CPC+Cymenol-T2:CPC+Cymenol -1.033333e+00 -25.1741698  23.107503 1.0000000
T5:CPC+Cymenol-T2:CPC+Cymenol  1.360000e+01 -10.5408365  37.740836 0.9437141
T6:CPC+Cymenol-T2:CPC+Cymenol  1.784000e+01  -6.3008365  41.980836 0.5394930
T0:Cymenol-T2:CPC+Cymenol      4.764000e+01  23.0718589  72.208141 0.0000000
T1:Cymenol-T2:CPC+Cymenol      1.429000e+01 -10.2781411  38.858141 0.9216180
T2:Cymenol-T2:CPC+Cymenol      2.081143e+01  -3.7567125  45.379570 0.2491579
T3:Cymenol-T2:CPC+Cymenol      3.416857e+01   9.6004303  58.736713 0.0001149
T4:Cymenol-T2:CPC+Cymenol      4.052571e+01  15.9575732  65.093855 0.0000007
T5:Cymenol-T2:CPC+Cymenol      4.290429e+01  18.3361446  67.472427 0.0000001
T6:Cymenol-T2:CPC+Cymenol      4.345429e+01  18.8861446  68.022427 0.0000000
T0:Placebo-T2:CPC+Cymenol      4.764000e+01  21.3961439  73.883856 0.0000000
T1:Placebo-T2:CPC+Cymenol      3.926727e+01  13.0234166  65.511129 0.0000155
T2:Placebo-T2:CPC+Cymenol      3.868545e+01  12.4415984  64.929311 0.0000238
T3:Placebo-T2:CPC+Cymenol      4.054000e+01  14.2961439  66.783856 0.0000059
T4:Placebo-T2:CPC+Cymenol      4.307636e+01  16.8325075  69.320220 0.0000008
T5:Placebo-T2:CPC+Cymenol      4.442182e+01  18.1779621  70.665674 0.0000003
T6:Placebo-T2:CPC+Cymenol      4.368545e+01  17.4415984  69.929311 0.0000005
T4:CPC+Cymenol-T3:CPC+Cymenol  3.140000e+00 -21.0008365  27.280836 1.0000000
T5:CPC+Cymenol-T3:CPC+Cymenol  1.777333e+01  -6.3675031  41.914170 0.5478394
T6:CPC+Cymenol-T3:CPC+Cymenol  2.201333e+01  -2.1275031  46.154170 0.1349310
T0:Cymenol-T3:CPC+Cymenol      5.181333e+01  27.2451922  76.381474 0.0000000
T1:Cymenol-T3:CPC+Cymenol      1.846333e+01  -6.1048078  43.031474 0.5018182
T2:Cymenol-T3:CPC+Cymenol      2.498476e+01   0.4166208  49.552903 0.0405273
T3:Cymenol-T3:CPC+Cymenol      3.834190e+01  13.7737637  62.910046 0.0000043
T4:Cymenol-T3:CPC+Cymenol      4.469905e+01  20.1309065  69.267189 0.0000000
T5:Cymenol-T3:CPC+Cymenol      4.707762e+01  22.5094780  71.645760 0.0000000
T6:Cymenol-T3:CPC+Cymenol      4.762762e+01  23.0594780  72.195760 0.0000000
T0:Placebo-T3:CPC+Cymenol      5.181333e+01  25.5694772  78.057189 0.0000000
T1:Placebo-T3:CPC+Cymenol      4.344061e+01  17.1967500  69.684462 0.0000006
T2:Placebo-T3:CPC+Cymenol      4.285879e+01  16.6149318  69.102644 0.0000009
T3:Placebo-T3:CPC+Cymenol      4.471333e+01  18.4694772  70.957189 0.0000002
T4:Placebo-T3:CPC+Cymenol      4.724970e+01  21.0058409  73.493553 0.0000000
T5:Placebo-T3:CPC+Cymenol      4.859515e+01  22.3512954  74.839008 0.0000000
T6:Placebo-T3:CPC+Cymenol      4.785879e+01  21.6149318  74.102644 0.0000000
T5:CPC+Cymenol-T4:CPC+Cymenol  1.463333e+01  -9.5075031  38.774170 0.8836685
T6:CPC+Cymenol-T4:CPC+Cymenol  1.887333e+01  -5.2675031  43.014170 0.4133363
T0:Cymenol-T4:CPC+Cymenol      4.867333e+01  24.1051922  73.241474 0.0000000
T1:Cymenol-T4:CPC+Cymenol      1.532333e+01  -9.2448078  39.891474 0.8505005
T2:Cymenol-T4:CPC+Cymenol      2.184476e+01  -2.7233792  46.412903 0.1693235
T3:Cymenol-T4:CPC+Cymenol      3.520190e+01  10.6337637  59.770046 0.0000525
T4:Cymenol-T4:CPC+Cymenol      4.155905e+01  16.9909065  66.127189 0.0000003
T5:Cymenol-T4:CPC+Cymenol      4.393762e+01  19.3694780  68.505760 0.0000000
 [ reached getOption("max.print") -- omitted 128 rows ]

#intra colutorios por tiempo mediante confocal

aov4 <- lm(Reduccion ~ Tiempo:Colutorio, 
                 data = df1)
summary(aov4)

Call:
lm(formula = Reduccion ~ Tiempo:Colutorio, data = df1)

Residuals:
    Min      1Q  Median      3Q     Max 
-59.429 -10.206   0.000   8.405  77.800 

Coefficients: (1 not defined because of singularities)
                              Estimate Std. Error t value Pr(>|t|)    
(Intercept)                    96.0455     5.3216  18.048  < 2e-16 ***
TiempoT0:ColutorioCPC           3.9545     7.1113   0.556 0.578504    
TiempoT1:ColutorioCPC         -41.9740     7.1113  -5.902 8.47e-09 ***
TiempoT2:ColutorioCPC         -39.8240     7.1113  -5.600 4.34e-08 ***
TiempoT3:ColutorioCPC         -34.0812     7.1113  -4.793 2.44e-06 ***
TiempoT4:ColutorioCPC         -33.9312     7.1113  -4.771 2.69e-06 ***
TiempoT5:ColutorioCPC         -25.0026     7.1113  -3.516 0.000496 ***
TiempoT6:ColutorioCPC         -12.0597     7.1113  -1.696 0.090804 .  
TiempoT0:ColutorioCPC+Cymenol   3.9545     7.0063   0.564 0.572823    
TiempoT1:ColutorioCPC+Cymenol -42.1388     7.0063  -6.014 4.54e-09 ***
TiempoT2:ColutorioCPC+Cymenol -43.6855     7.0063  -6.235 1.30e-09 ***
TiempoT3:ColutorioCPC+Cymenol -47.8588     7.0063  -6.831 3.75e-11 ***
TiempoT4:ColutorioCPC+Cymenol -44.7188     7.0063  -6.383 5.52e-10 ***
TiempoT5:ColutorioCPC+Cymenol -30.0855     7.0063  -4.294 2.27e-05 ***
TiempoT6:ColutorioCPC+Cymenol -25.8455     7.0063  -3.689 0.000261 ***
TiempoT0:ColutorioCymenol       3.9545     7.1113   0.556 0.578504    
TiempoT1:ColutorioCymenol     -29.3955     7.1113  -4.134 4.47e-05 ***
TiempoT2:ColutorioCymenol     -22.8740     7.1113  -3.217 0.001418 ** 
TiempoT3:ColutorioCymenol      -9.5169     7.1113  -1.338 0.181677    
TiempoT4:ColutorioCymenol      -3.1597     7.1113  -0.444 0.657082    
TiempoT5:ColutorioCymenol      -0.7812     7.1113  -0.110 0.912593    
TiempoT6:ColutorioCymenol      -0.2312     7.1113  -0.033 0.974086    
TiempoT0:ColutorioPlacebo       3.9545     7.5259   0.525 0.599599    
TiempoT1:ColutorioPlacebo      -4.4182     7.5259  -0.587 0.557541    
TiempoT2:ColutorioPlacebo      -5.0000     7.5259  -0.664 0.506891    
TiempoT3:ColutorioPlacebo      -3.1455     7.5259  -0.418 0.676241    
TiempoT4:ColutorioPlacebo      -0.6091     7.5259  -0.081 0.935542    
TiempoT5:ColutorioPlacebo       0.7364     7.5259   0.098 0.922112    
TiempoT6:ColutorioPlacebo           NA         NA      NA       NA    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 17.65 on 350 degrees of freedom
Multiple R-squared:  0.5375,    Adjusted R-squared:  0.5019 
F-statistic: 15.07 on 27 and 350 DF,  p-value: < 2.2e-16

#plot los residuos

plot(aov3, 1)

#Normal Q-Q

plot(aov3, 2)

#prueba de shapiro-wilk normality test

# Extract the residuals
aov_residuals1 <- residuals(object = aov3)
# Run Shapiro-Wilk test
shapiro.test(x = aov_residuals1 )

    Shapiro-Wilk normality test

data:  aov_residuals1
W = 0.95769, p-value = 5.762e-09
performance::check_model(aov3)

LS0tCnRpdGxlOiAiQW5hbGlzaXMgZGUgY29sdXRvcmlvcyBJU0RJTiAoUGxhY2Viby1DeW1lbm9sLUNQQy1DUEMrQ3ltZW5vbCIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQojQU5BTElaTyBMT1MgREFUT1MgREUgRkxVT1JPTUVUUk8KI3BhY2thZ2VzIApgYGB7cn0KcGFjbWFuOjpwX2xvYWQodGlkeXZlcnNlLCAKICAgICAgICAgICAgICAgcGVyZm9ybWFuY2UsIAogICAgICAgICAgICAgICBnZ3B1YnIpCmxpYnJhcnkoInBhdGNod29yayIpCmBgYAojYWJybyBlbCBkZgpgYGB7cn0KZGYgPC0gcmVhZC5jc3YoImh0dHBzOi8vZG9jcy5nb29nbGUuY29tL3NwcmVhZHNoZWV0cy9kL2UvMlBBQ1gtMXZRNE5rOTZKeDFuLXBMQzFlN2xYbFVnTXRxQUZjMXhPVHlOQ1A2XzcyY3ZndGNaWGpJVE0yLXdvY1ViZ0tXaW42M2VPSFdaUDNlendKb0IvcHViP2dpZD0xNTYzOTA5OTQ5JnNpbmdsZT10cnVlJm91dHB1dD1jc3YiKQpgYGAKI2NhbWJpbyBjb2x1dG9yaW8gY29tbyBmYWN0b3IKYGBge3J9CmRmIDwtIGRmICU+JSAKICBtdXRhdGUoQ29sdXRvcmlvID0gYXMuZmFjdG9yKENvbHV0b3JpbykpCmBgYAojZ2xpbXNlCmBgYHtyfQpnbGltcHNlKGRmKQpgYGAKI2hpc3QKYGBge3J9CmRmICU+JSAKICAjIHBpdm90X2xvbmdlcihDb2x1dG9yaW86VGllbXBvLCAKICAjICAgICAgICAgICAgIG5hbWVzX3RvID0gIlZhcmlhYmxlcyIsIAogICMgICAgICAgICAgICAgdmFsdWVzX3RvID0gIlZhbHVlcyIpICU+JSAKICBnZ3Bsb3QoYWVzKHggPSBSZWR1Y2Npb24pKSArIAogIGdlb21faGlzdG9ncmFtKCkgKyAKICBmYWNldF9ncmlkKFRpZW1wbyB+IENvbHV0b3JpbykKYGBgCiNncmFmaWNhIGxpbmVhbApgYGB7cn0KZGYgJT4lIApnZ2xpbmUoeCA9ICJUaWVtcG8iLCAKICAgICAgIHkgPSAiUmVkdWNjaW9uIiwgCiAgICAgICBjb2xvciA9ICJDb2x1dG9yaW8iLAogICAgICAgYWRkID0gYygibWVhbl9zZSIsICJkb3RwbG90IikpICsgCiAgbGFicyh0aXRsZSA9ICJQcm9maWxlIHBsb3QgZGUgcmVkdWNjacOzbiBlbiB0aWVtcG8geSB0aXBvIGRlIGNvbHV0b3JpbyIpICsgCiAgdGhlbWVfcHVicigpCmBgYAojYm9uZmVycm9uaQpgYGB7cn0KcGFpcndpc2UudC50ZXN0KCBkZiRSZWR1Y2Npb24sIAogICAgICAgICAgICAgICAgIGRmJENvbHV0b3JpbywKICAgICAgICAgICAKICAgICAgICAgICAgICAgICBwLmFkaiA9ICJib25mZXJyb25pIikKYGBgCiNtb2RlbG8gZGUgYW92CmBgYHtyfQphb3YxIDwtIGFvdihSZWR1Y2Npb24gfiBUaWVtcG8gKyBDb2x1dG9yaW8gKyBUaWVtcG8qQ29sdXRvcmlvICwgCiAgICAgICAgICAgZGF0YSA9IGRmKQpzdW1tYXJ5KGFvdjEpCmBgYApgYGB7cn0KZGYgJT4lIAogIHdpdGgoYW92KFJlZHVjY2lvbiB+IFRpZW1wbyArIENvbHV0b3JpbyArIFRpZW1wbzpDb2x1dG9yaW8gKyBFcnJvcihpZCkpKSAlPiUgCiAgc3VtbWFyeSgpCmBgYAojdHVrZXkKYGBge3J9ClR1a2V5SFNEKGFvdjEpCmBgYApgYGB7cn0KYW92MiA8LSBsbShSZWR1Y2Npb24gfiBUaWVtcG86Q29sdXRvcmlvLCAKICAgICAgICAgICAgICAgICBkYXRhID0gZGYpCnN1bW1hcnkoYW92MikKYGBgCiNwbG90IGxvcyByZXNpZHVvcyAKYGBge3J9CnBsb3QoYW92MSwgMSkKYGBgCiNOb3JtYWwgUS1RCmBgYHtyfQpwbG90KGFvdjEsIDIpCmBgYAojcHJ1ZWJhIGRlIHNoYXBpcm8td2lsayBub3JtYWxpdHkgdGVzdCAKYGBge3J9CiMgRXh0cmFjdCB0aGUgcmVzaWR1YWxzCmFvdl9yZXNpZHVhbHMgPC0gcmVzaWR1YWxzKG9iamVjdCA9IGFvdjEpCiMgUnVuIFNoYXBpcm8tV2lsayB0ZXN0CnNoYXBpcm8udGVzdCh4ID0gYW92X3Jlc2lkdWFscyApCmBgYApgYGB7cn0KcGVyZm9ybWFuY2U6OmNoZWNrX21vZGVsKGFvdjEpCmBgYAojQU5BTElaTyBQQVJBIE3DiVRPRE8gREUgQ09ORk9DQUwKI2Ficm8gZWwgZGYKYGBge3J9CmRmMSA8LSByZWFkX2NzdigiaHR0cHM6Ly9kb2NzLmdvb2dsZS5jb20vc3ByZWFkc2hlZXRzL2QvZS8yUEFDWC0xdlE0Tms5Nkp4MW4tcExDMWU3bFhsVWdNdHFBRmMxeE9UeU5DUDZfNzJjdmd0Y1pYaklUTTItd29jVWJnS1dpbjYzZU9IV1pQM2V6d0pvQi9wdWI/Z2lkPTE3OTE3MjM0MzQmc2luZ2xlPXRydWUmb3V0cHV0PWNzdiIpCmBgYAojY2FtYmlvIGNvbHV0b3JpbyBjb21vIGZhY3RvcgpgYGB7cn0KZGYxIDwtIGRmMSAlPiUgCiAgbXV0YXRlKENvbHV0b3JpbyA9IGFzLmZhY3RvcihDb2x1dG9yaW8pKQpgYGAKI2dsaW1zZQpgYGB7cn0KZ2xpbXBzZShkZjEpCmBgYAojaGVhZCBkZjEKYGBge3J9CmhlYWQoZGYxKQpgYGAKI2hpc3QgY29uZm9jYWwgCmBgYHtyfQpkZjEgJT4lIAogICMgcGl2b3RfbG9uZ2VyKENvbHV0b3JpbzpUaWVtcG8sIAogICMgICAgICAgICAgICAgbmFtZXNfdG8gPSAiVmFyaWFibGVzIiwgCiAgIyAgICAgICAgICAgICB2YWx1ZXNfdG8gPSAiVmFsdWVzIikgJT4lIAogIGdncGxvdChhZXMoeCA9IFJlZHVjY2lvbikpICsgCiAgZ2VvbV9oaXN0b2dyYW0oKSArIAogIGZhY2V0X2dyaWQoVGllbXBvIH4gQ29sdXRvcmlvKQpgYGAKI2dyYWZpY2EgbGluZWFsIGNvbmZvY2FsCmBgYHtyfQpkZjEgJT4lIApnZ2xpbmUoeCA9ICJUaWVtcG8iLCAKICAgICAgIHkgPSAiUmVkdWNjaW9uIiwgCiAgICAgICBjb2xvciA9ICJDb2x1dG9yaW8iLAogICAgICAgYWRkID0gYygibWVhbl9zZSIsICJkb3RwbG90IikpICsgCiAgbGFicyh0aXRsZSA9ICJQcm9maWxlIHBsb3QgZGUgcmVkdWNjacOzbiBlbiB0aWVtcG8geSB0aXBvIGRlIGNvbHV0b3JpbyIpICsgCiAgdGhlbWVfcHVicigpCmBgYAojYm9uZmVycm9uaSBjb25mb2NhbApgYGB7cn0KcGFpcndpc2UudC50ZXN0KCBkZjEkUmVkdWNjaW9uLCAKICAgICAgICAgICAgICAgICBkZjEkQ29sdXRvcmlvLAogICAgICAgICAgIAogICAgICAgICAgICAgICAgIHAuYWRqID0gImJvbmZlcnJvbmkiKQpgYGAKI21vZGVsbyBkZSBhb3YgcGFyYSBjb25mb2NhbCAKYGBge3J9CmFvdjMgPC0gYW92KFJlZHVjY2lvbiB+IFRpZW1wbyArIENvbHV0b3JpbyArIFRpZW1wbypDb2x1dG9yaW8gLCAKICAgICAgICAgICBkYXRhID0gZGYxKQpzdW1tYXJ5KGFvdjMpCmBgYApgYGB7cn0KZGYxICU+JSAKICB3aXRoKGFvdihSZWR1Y2Npb24gfiBUaWVtcG8gKyBDb2x1dG9yaW8gKyBUaWVtcG86Q29sdXRvcmlvICsgRXJyb3IoaWQpKSkgJT4lIAogIHN1bW1hcnkoKQpgYGAKI3R1a2V5CmBgYHtyfQpUdWtleUhTRChhb3YzKQpgYGAKI2ludHJhIGNvbHV0b3Jpb3MgcG9yIHRpZW1wbyBtZWRpYW50ZSBjb25mb2NhbCAKYGBge3J9CmFvdjQgPC0gbG0oUmVkdWNjaW9uIH4gVGllbXBvOkNvbHV0b3JpbywgCiAgICAgICAgICAgICAgICAgZGF0YSA9IGRmMSkKc3VtbWFyeShhb3Y0KQpgYGAKI3Bsb3QgbG9zIHJlc2lkdW9zIApgYGB7cn0KcGxvdChhb3YzLCAxKQpgYGAKI05vcm1hbCBRLVEKYGBge3J9CnBsb3QoYW92MywgMikKYGBgCiNwcnVlYmEgZGUgc2hhcGlyby13aWxrIG5vcm1hbGl0eSB0ZXN0IApgYGB7cn0KIyBFeHRyYWN0IHRoZSByZXNpZHVhbHMKYW92X3Jlc2lkdWFsczEgPC0gcmVzaWR1YWxzKG9iamVjdCA9IGFvdjMpCiMgUnVuIFNoYXBpcm8tV2lsayB0ZXN0CnNoYXBpcm8udGVzdCh4ID0gYW92X3Jlc2lkdWFsczEgKQpgYGAKYGBge3J9CnBlcmZvcm1hbmNlOjpjaGVja19tb2RlbChhb3YzKQpgYGAKCg==