This paper estimates time-varying COVID-19 reproduction numbers worldwide solely based on the number of reported infected cases, allowing for under-reporting. Estimation is based on a moment condition that can be derived from an agent-based stochastic network model of COVID-19 transmission. The outcomes in terms of the reproduction number and the trajectory of per-capita cases through the end of 2020 are very diverse. The reproduction number depends on the transmission rate and the proportion of susceptible population, or the herd immunity effect. Changes in the transmission rate depend on changes in the behavior of the virus, reflecting mutations and vaccinations, and changes in people's behavior, reflecting voluntary or government mandated isolation. Over our sample period, neither mutation nor vaccination are major factors, so one can attribute variation in the transmission rate to variations in behavior. Evidence based on panel data models explaining transmission rates for nine European countries indicates that the diversity of outcomes results from the non-linear interaction of mandatory containment measures, voluntary precautionary isolation, and the economic incentives that governments provided to support isolation. These effects are precisely estimated and robust to various assumptions. As a result, countries with seemingly different social distancing policies achieved quite similar outcomes in terms of the reproduction number. These results imply that ignoring the voluntary component of social distancing could introduce an upward bias in the estimates of the effects of lock-downs and support policies on the transmission rates. The full set of estimation results and the replication package are available on the authors' websites.