#SETWD

setwd("~/MAESTRIA EPIDEMIOLOGIA ICESI/ensayos_clinicos/tamano_muestra")

Instalar paquetes

Llamar paquetes

library(dplyr)
## Warning: package 'dplyr' was built under R version 4.5.3
## 
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(tidyr)
## Warning: package 'tidyr' was built under R version 4.5.3
library(stats)

EJERCICIO 1 - Diferencia de promedios

“Eficacia del entrenamiento muscular respiratorio en el destete de la ventilación mecánica en pacientes con ventilación mecánica por 48 o más horas: un ensayo clínico controlado”

¿El entrenamiento muscular respiratorio afecta el tiempo del retiro de la ventilación mecánica?

#Parámetros

sd <- 52.8        
delta_grupos <- 24       
alfa <- 0.05     
poder <- 0.80     
dropout <- 0.10

Cálculo de tamaño de muestra en destete de ventilación mecánica

VentilacionM <- power.t.test(
  delta = 24,
  sd = 52.8,
  sig.level = 0.05,
  power = 0.80,
  type = "two.sample",
  alternative = "two.sided"
)

VentilacionM
## 
##      Two-sample t test power calculation 
## 
##               n = 76.9494
##           delta = 24
##              sd = 52.8
##       sig.level = 0.05
##           power = 0.8
##     alternative = two.sided
## 
## NOTE: n is number in *each* group

n por grupo redondeado

n_por_grupo_funcion <- ceiling(VentilacionM$n)

Ajuste por pérdidas

n_ajustado_funcion <- ceiling(n_por_grupo_funcion / (1 - 0.10))
cat("n por grupo según power.t.test:", n_por_grupo_funcion, "\n")
## n por grupo según power.t.test: 77
cat("n por grupo ajustado por pérdidas:", n_ajustado_funcion, "\n")
## n por grupo ajustado por pérdidas: 86
cat("Tamaño total:", n_ajustado_funcion * 2, "\n")
## Tamaño total: 172

EJERCICIO 2 - Diferencia de proporciones

solucionsln <- power.prop.test(
  p1 = 0.60,
  p2 = 0.75,
  sig.level = 0.05,
  power = 0.80,
  alternative = "two.sided"
)

solucionsln
## 
##      Two-sample comparison of proportions power calculation 
## 
##               n = 151.8689
##              p1 = 0.6
##              p2 = 0.75
##       sig.level = 0.05
##           power = 0.8
##     alternative = two.sided
## 
## NOTE: n is number in *each* group

#Sin Ajuste

n_por_grupo <- ceiling(solucionsln$n)

#con ajuste de 20% de pérdidas

n_ajustado <- ceiling(n_por_grupo / (1 - 0.20))
n_total <- n_ajustado * 2

#Resumen

cat("n por grupo sin ajuste:", n_por_grupo, "\n")
## n por grupo sin ajuste: 152
cat("n por grupo ajustado por pérdidas:", n_ajustado, "\n")
## n por grupo ajustado por pérdidas: 190
cat("Tamaño total:", n_total, "\n")
## Tamaño total: 380