df %>% group_by(c,d) %>% summarise_all(sum, na.rm = T)
df %>% group_by(c,d) %>% summarise_all(function(x) sum(!is.na(x)))

谢谢老师，但是能把代码合到一起吗？我试了下没成功

df %>%
  group_by(c,d) %>%
  summarise_all(sum, na.rm = T)  %>%
  left_join(df %>%
                 group_by(c,d) %>%
                 summarise_all(function(x) sum(!is.na(x))),
               by = c("c", "d"))

再提供一种方法，参考下

library(tidyverse)
library(rlist)

df %>%
  group_by(c,d) %>%
  nest() %>%
  mutate(dt = list.apply(data,
    function(x) {
      data.frame(
      a_non= sum(!is.na(x$a)),
      a_sum = sum(x$a, na.rm = T),
      b_non= sum(!is.na(x$b)),
      b_sum = sum(x$b, na.rm = T)
    )}
  ))%>%
  select(c,d,dt)%>%
  unnest(cols = c(dt))

嗯嗯，棒极了，谢谢两位老师，我也找到了一种：with(df,aggregate(df[,c('a','b')], by=list(Vc=c,Vd=d),
                  FUN=function(x){c(sum=sum(x,na.rm=T),N=sum(!is.na(x)))}))