Need help converting an R aggregate function into a Python groupby

Question

df <- aggregate( .~Age+Subs,df[,c(\'Age\',\'Subs\')],sum, na.rm=True)

This is the piece of code and I cant figure out what to place inside my gro