Cuadrado latino

library(readxl)
## Warning: package 'readxl' was built under R version 4.0.4
Cuadrado_latino <- read_excel("C:/Users/ASUS/OneDrive/Escritorio/Cuadrado latino.xlsx")

library(collapsibleTree)
## Warning: package 'collapsibleTree' was built under R version 4.0.5
collapsibleTreeSummary(Cuadrado_latino,
                       hierarchy = c('noches','paciente','farmaco','calidad_sue'))

##Analisis descriptivo

library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.0.5
ggplot(Cuadrado_latino)+aes(farmaco, calidad_sue)+geom_point()+facet_grid(paciente ~ noches)+coord_flip()

aggregate(Cuadrado_latino$calidad_sue, list(Cuadrado_latino$farmaco, Cuadrado_latino$noches, Cuadrado_latino$paciente),mean)
##                Group.1 Group.2   Group.3    x
## 1          Etclorvinol    dia1 paciente1 2.92
## 2  Secobarbitol sodico    dia2 paciente1 2.43
## 3              Placebo    dia3 paciente1 2.19
## 4          Glutetimida    dia4 paciente1 2.71
## 5    hidrato de cloral    dia5 paciente1 2.71
## 6    hidrato de cloral    dia1 paciente2 2.86
## 7              Placebo    dia2 paciente2 1.64
## 8  Secobarbitol sodico    dia3 paciente2 3.02
## 9          Etclorvinol    dia4 paciente2 3.03
## 10         Glutetimida    dia5 paciente2 3.03
## 11 Secobarbitol sodico    dia1 paciente3 1.97
## 12         Etclorvinol    dia2 paciente3 2.50
## 13         Glutetimida    dia3 paciente3 2.47
## 14   hidrato de cloral    dia4 paciente3 2.65
## 15             Placebo    dia5 paciente3 1.89
## 16             Placebo    dia1 paciente4 1.99
## 17         Glutetimida    dia2 paciente4 2.39
## 18   hidrato de cloral    dia3 paciente4 2.37
## 19 Secobarbitol sodico    dia4 paciente4 2.33
## 20         Etclorvinol    dia5 paciente4 2.71
## 21         Glutetimida    dia1 paciente5 2.64
## 22   hidrato de cloral    dia2 paciente5 2.31
## 23         Etclorvinol    dia3 paciente5 2.44
## 24             Placebo    dia4 paciente5 1.86
## 25 Secobarbitol sodico    dia5 paciente5 2.78
tapply(Cuadrado_latino$calidad_sue, list(Cuadrado_latino$farmaco,Cuadrado_latino$noches, Cuadrado_latino$paciente), mean)
## , , paciente1
## 
##                     dia1 dia2 dia3 dia4 dia5
## Etclorvinol         2.92   NA   NA   NA   NA
## Glutetimida           NA   NA   NA 2.71   NA
## hidrato de cloral     NA   NA   NA   NA 2.71
## Placebo               NA   NA 2.19   NA   NA
## Secobarbitol sodico   NA 2.43   NA   NA   NA
## 
## , , paciente2
## 
##                     dia1 dia2 dia3 dia4 dia5
## Etclorvinol           NA   NA   NA 3.03   NA
## Glutetimida           NA   NA   NA   NA 3.03
## hidrato de cloral   2.86   NA   NA   NA   NA
## Placebo               NA 1.64   NA   NA   NA
## Secobarbitol sodico   NA   NA 3.02   NA   NA
## 
## , , paciente3
## 
##                     dia1 dia2 dia3 dia4 dia5
## Etclorvinol           NA  2.5   NA   NA   NA
## Glutetimida           NA   NA 2.47   NA   NA
## hidrato de cloral     NA   NA   NA 2.65   NA
## Placebo               NA   NA   NA   NA 1.89
## Secobarbitol sodico 1.97   NA   NA   NA   NA
## 
## , , paciente4
## 
##                     dia1 dia2 dia3 dia4 dia5
## Etclorvinol           NA   NA   NA   NA 2.71
## Glutetimida           NA 2.39   NA   NA   NA
## hidrato de cloral     NA   NA 2.37   NA   NA
## Placebo             1.99   NA   NA   NA   NA
## Secobarbitol sodico   NA   NA   NA 2.33   NA
## 
## , , paciente5
## 
##                     dia1 dia2 dia3 dia4 dia5
## Etclorvinol           NA   NA 2.44   NA   NA
## Glutetimida         2.64   NA   NA   NA   NA
## hidrato de cloral     NA 2.31   NA   NA   NA
## Placebo               NA   NA   NA 1.86   NA
## Secobarbitol sodico   NA   NA   NA   NA 2.78

##Hipotesis## \[H_0: \mu_A = \mu_B = \mu_C\] Modelo \[y_{ijk} = \mu + \tau_i + \beta_j + \delta_k + \epsilon_{ijk}\\i=j=k=1,\dots,3\]

mod_sue = lm(calidad_sue ~ noches + paciente + farmaco, data = Cuadrado_latino)
anova(mod_sue)
## Analysis of Variance Table
## 
## Response: calidad_sue
##           Df  Sum Sq Mean Sq F value    Pr(>F)    
## noches     4 0.36622 0.09155  1.9119 0.1731404    
## paciente   4 0.61126 0.15281  3.1911 0.0529803 .  
## farmaco    4 2.08326 0.52081 10.8758 0.0005825 ***
## Residuals 12 0.57465 0.04789                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Se visualiza el p-valor del farmaco el cual fue el tratamiento usado, se evidencia que el p-valor es menor a 5% por lo tanto se rechaza la H_o y se determina que los farmacos que afectan la calidad del sueño de los ancianos presentan diferencias en la accion de los mismos, el más efectivo podria ser el farmaco Etclorvinol y Glutetimida.