Hola Victor:
Me permito hacerte alguna sugerencia sobre tu código:
- Codificar los valores de sexo numéricamente puede llevar a confusión, ¿cuál es 0 y cuál 1?, no hay ningún inconveniente en usar 'H' y 'M'.
- Dentro de la función utilizas la variable todos que no se ha pasado a la función como parámetro, eso hace que la función sea poco reusable, pues depende del nombre del conjunto de datos a usar, yo la definiría como
subgr <- function(todos, a, s, e)
Un saludo