set.seed(42)
n = 10000
df = cbind(ID = 1:1000, Var1 = runif(n, min=0, max=1), Var2 = ifelse(runif(n, min=0, max=1)>=0.5,1,0), Var3 = ifelse(runif(n, min=0, max=1)>0.75,1,0), Var4 = round(runif(n, min=20, max=90)) ) 
df = as.data.frame(df)

summary(df)

modelo_prueba = lm('Var2 ~ Var3 + Var4', df)
summary(modelo_prueba)

library(pwr)
### Antes de ejecutar library(pwr) DEBEN INSTALAR la librería, lo que se hace con el siguiente código: install.packages('pwr')
pwr.t.test(d=0.5, sig.level=0.1, power=0.9, type="two.sample", alternative = 'two.sided')