#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)