library(dplyr)
library(readr)
library(tidyr)

df <- read_delim("/Users/msserpa/Downloads/microdados_censo_superior_2016/DADOS/DM_CURSO.CSV", delim="|") %>% as.data.frame

df <- df %>% filter(!is.na(QT_MATRICULA_CURSO)) %>% filter(NO_MUNICIPIO_CURSO == "Porto Alegre") %>%  group_by(NO_IES) %>% summarize(N=n(), mean=sum(QT_MATRICULA_CURSO)) %>% as.data.frame()

df