Update: implemented NA threshold in function 'aggMonths.R'.