We perform a statistical standard siren analysis of GW170817. Our analysis does not utilize knowledge of NGC 4993 as the unique host galaxy of the optical counterpart to GW170817. Instead, we consider each galaxy within the GW170817 localization region as a potential host; combining the redshifts from all of the galaxies with the distance estimate from GW170817 provides an estimate of the Hubble constant, H-0. Considering all galaxies brighter than 0.626L(B)(star) as equally likely to host a binary neutron star merger, we find H-0 = 77(-18)(+37) km s(-1) Mpc(-1) (maximum a posteriori and 68.3% highest density posterior interval; assuming a flat H-0 prior in the range [10, 220] km s(-1) Mpc(-1)). We explore the dependence of our results on the thresholds by which galaxies are included in our sample, and we show that weighting the host galaxies by stellar mass or star formation rate provides entirely consistent results with potentially tighter constraints. By applying the method to simulated gravitational-wave events and a realistic galaxy catalog we show that, because of the small localization volume, this statistical standard siren analysis of GW170817 provides an unusually informative (top 10%) constraint. Under optimistic assumptions for galaxy completeness and redshift uncertainty, we find that dark binary neutron star measurements of H-0 will converge as 40%/root(N), where N is the number of sources. While these statistical estimates are inferior to the value from the counterpart standard siren measurement utilizing NGC 4993 as the unique host, H-0 = 76(-13)(+19) km s(-1) Mpc(-1) (determined from the same publicly available data), our analysis is a proof-of-principle demonstration of the statistical approach first proposed by Bernard Schutz over 30 yr ago.
