From what I can tell, these numbers are not that unreasonable. For 2016 1-year ACS data, table S1901 shows mean incomes by household type. (See:
factfinder.census.gov/.../0400000US25) From table S1901:
There are 2,579,398 households with a mean household income of $101,911.
There are 1,633,661 family households with a mean HHI of $122,310.
There are 1,194,726 married-couple family households with a mean HHI of $143,966.
There are 945,737 non-family households with a mean HHI of $62,500.
With weighted 2016 1-year ACS PUMS for MA I get:
There are 2,579,453 households with a mean household income of $101,593.
There are 1,634,746 family households with a mean HHI of $124,355.
There are 1,195,011 married-couple family households with a mean HHI of $144,204.
There are 944,707 non-family households with a mean HHI of $62,207.
-And for comparison, with PUMS-
There are 17,396 same-sex couple households with a mean HHI of $181,920.
Stata code:
...load data...
gen hincp_adj = hincp * adjinc/1000000
gen hh = 1
gen fam_hh = inrange(hht, 1, 3)
gen mcfam_hh = hht == 1
gen nonfam_hh = inrange(hht, 4, 7)
gen ssmc_hh = inlist(ssmc, 1, 2)
tab1 hh fam_hh mcfam_hh nonfam_hh ssmc_hh [fw=wgtp], sum(hincp_adj)