library(rio)
p2=import("d2.xlsx")
str(p2)

## 'data.frame':    26574 obs. of  7 variables:
##  $ sex            : chr  "male" "male" "male" "male" ...
##  $ padreSector    : chr  "Agricultura" "Agricultura" "NoFijo" "Agricultura" ...
##  $ fechaNacimiento: chr  "1853-05-23" "1853-07-19" "1861-11-17" "1872-11-16" ...
##  $ edadDejaEstudio: num  15 15 15 15 0.559 0.315 15 15 15 15 ...
##  $ muere          : num  0 0 0 0 1 1 0 0 0 0 ...
##  $ naceFueraMatri : chr  "no" "no" "no" "no" ...
##  $ madreEdad      : num  35 30.6 29.3 41.2 42.1 ...

p2=na.omit(p2)

mylogit <- glm(p2$muere~p2$sex + p2$madreEdad + p2$naceFueraMatri, 
               data = p2, family = "binomial")

 sdVIs=apply(p2[,c("sex","madreEdad", "naceFueraMatri")],2,sd)

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introducidos por coerción

## Warning in var(if (is.vector(x) || is.factor(x)) x else as.double(x), na.rm =
## na.rm): NAs introducidos por coerción

list(LogitSt=sdVIs*coef(mylogit)[c(2,3,4)])

## $LogitSt
##            sex      madreEdad naceFueraMatri 
##             NA     0.06364546             NA

    data.frame()

## data frame with 0 columns and 0 rows