`df` é o data frame com os dados que irei trabalhar

    df <- data.frame(matrix(NA,nrow = 30, ncol = 6))
    colnames(df) <- c("ano","mes","x1","x2","x3","x4")
    
    df$ano <- c(2002,2002,2002,2003,2003,2003,2004,2004,2004,2004,2002,2002,2002,2003,2003,2003,2004,2004,2004,2004,2002,2002,2002,2003,2003,2003,2004,2004,2004,2004)
    df$mes <- c(5,9,7,3,8,2,4,1,6,10,2,11,12,10,12,2,8,9,10,7,6,4,5,7,3,1,9,5,2,6)
    df$x1 <- c(2,4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36,38,40,42,44,46,48,50,52,54,56,58,60)
    df$x2 <- c(3,6,9,12,15,18,21,24,27,30,33,36,39,42,45,48,51,54,57,60,63,66,69,72,75,78,81,84,87,90)
    df$x3 <- c(4,8,12,16,20,24,28,32,36,40,44,48,52,56,60,64,68,72,76,80,84,88,92,96,100,104,108,112,116,120)
    df$x4 <- c(5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100,105,110,115,120,125,130,135,140,145,150)
    
    View(df)

`resumo` é o data frame onde quero que minhas análises sejam registradas de forma a exportar em .csv ao final

    resumo <- data.frame(matrix(NA,nrow = 12, ncol = 7))
    colnames(resumo) <- c("mes","media_diaria_X1","media_mensal_X1","media_diaria_X4","media_mensal_X4","media_diaria_x3_divido_x4","media_mensal_x3_divido_x4")
    resumo$mes <- c(1,2,3,4,5,6,7,8,9,10,11,12)
    
    View(resumo)
    
    library(dplyr)

Quando eu utilizo o pacote dplyr ele me retorna valores com as colunas nomeadas iguais as colunas do data frame `df` visto que a seleção é realizada pelas colunas de `df`. 
A grande questão é que gostaria que os resultados obtidos abaixo fossem registrados no data frame `resumo` com as colunas nomeadas de acordo com o que defini acima .

Abaixo são os cálculos que consegui fazer e que pretendo fazer com o auxílio de alguém:

Calculo das medias diarias de x1 e x4

    media_diaria <- df %>%
      group_by(mes) %>%
      summarise_at(vars(x1,x4), mean, na.rm = TRUE)
    
    View(media_diaria)

[![view(media_diaria][1]][1]

Calculo das medias mensais (referente a soma mensal) de x1 e x4

    soma_mensal <- df %>% 
      group_by(ano,mes) %>%
      summarise_at(vars(x1,x4), sum, na.rm = TRUE)
    media_mensal <- soma_mensal %>% 
      group_by(mes) %>% 
      summarise_at(vars(x1,x4), mean, na.rm = TRUE)
    
    View(soma_mensal) 

   

[![view(soma_mensal)][2]][2]

    View(media_mensal)

[![view(media_mensal)][3]][3]

A partir desse ponto os cálculos não deram certo e não consigo resolver.  

Cálculos da média diária com o resultado obtido da dividão de x3 por x4. Não deu certo!

    media_diaria_divisao <- df %>%
      group_by(mes) %>%
      summarise_at((df$x3/df$x4), mean, na.rm = TRUE)
    
    Error: `.vars` must be a character/numeric vector or a `vars()` object, not a double vector.
    Run `rlang::last_error()` to see where the error occurred.

Cálculo das medias mensais (referente a soma mensal) de x1 e x4. Também não deu certo!

    soma_mensal <- df %>% 
      group_by(ano,mes) %>%
      summarise_at((df$x3/df$x4), sum, na.rm = TRUE)
    media_mensal <- soma_mensal %>% 
      group_by(mes) %>% 
      summarise_at((df$x3/df$x4), mean, na.rm = TRUE)

    Error: `.vars` must be a character/numeric vector or a `vars()` object, not a double vector.
    Run `rlang::last_error()` to see where the error occurred.

Resumindo minhas duas dúvidas:

1. Como consigo realizar as 2 análises acima em relação à divisão de x3/x4 com o `dplyr`? Caso utilizar o `dplyr` não seja uma boa opção, qual seria uma solução possível utilizando funções? Estou tentando evitar trabalhar com loops. 
 
2. Como faço para alimentar o dataframe `resumo` com os resultados das análises acima?


  [1]: https://i.sstatic.net/hau9l.png
  [2]: https://i.sstatic.net/jlkiy.png
  [3]: https://i.sstatic.net/6J06v.png