Stata - descriptive statistic for Panel data

Question

I am trying to create a table showing the number of neighborhoods and metropolitan areas for each year of my (unbalanced) panel data set. I would like something like this :

enter image description here

Also, I would like to display the number of coastal neighborhoods (I have a dummy indicating the coastal status of the neighborhood) on another columns.

I have tried this code,

tabout year if tr_pop>0 & pri!=. using "$results/table1/table1.tex", stats(count summ i_shore) style(tex) replace f(0c ) clab(Neighborhoods)

but I only got one columns about the number of neighborhoods.