PMC:7210464 / 11801-17898 JSONTXT

Annnotations TAB JSON ListView MergeView

    LitCovid-PD-FMA-UBERON

    {"project":"LitCovid-PD-FMA-UBERON","denotations":[{"id":"T3","span":{"begin":4539,"end":4543},"obj":"Body_part"}],"attributes":[{"id":"A3","pred":"fma_id","subj":"T3","obj":"http://purl.org/sig/ont/fma/fma9712"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-PD-UBERON

    {"project":"LitCovid-PD-UBERON","denotations":[{"id":"T3","span":{"begin":4539,"end":4543},"obj":"Body_part"}],"attributes":[{"id":"A3","pred":"uberon_id","subj":"T3","obj":"http://purl.obolibrary.org/obo/UBERON_0002398"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-PD-MONDO

    {"project":"LitCovid-PD-MONDO","denotations":[{"id":"T45","span":{"begin":443,"end":453},"obj":"Disease"},{"id":"T46","span":{"begin":503,"end":511},"obj":"Disease"},{"id":"T47","span":{"begin":529,"end":538},"obj":"Disease"},{"id":"T48","span":{"begin":561,"end":580},"obj":"Disease"},{"id":"T49","span":{"begin":1332,"end":1340},"obj":"Disease"},{"id":"T50","span":{"begin":1516,"end":1529},"obj":"Disease"},{"id":"T51","span":{"begin":1561,"end":1569},"obj":"Disease"},{"id":"T52","span":{"begin":1945,"end":1955},"obj":"Disease"},{"id":"T53","span":{"begin":3140,"end":3150},"obj":"Disease"},{"id":"T54","span":{"begin":4482,"end":4492},"obj":"Disease"},{"id":"T55","span":{"begin":4601,"end":4611},"obj":"Disease"},{"id":"T56","span":{"begin":4845,"end":4848},"obj":"Disease"},{"id":"T57","span":{"begin":5111,"end":5121},"obj":"Disease"}],"attributes":[{"id":"A45","pred":"mondo_id","subj":"T45","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A46","pred":"mondo_id","subj":"T46","obj":"http://purl.obolibrary.org/obo/MONDO_0100096"},{"id":"A47","pred":"mondo_id","subj":"T47","obj":"http://purl.obolibrary.org/obo/MONDO_0005249"},{"id":"A48","pred":"mondo_id","subj":"T48","obj":"http://purl.obolibrary.org/obo/MONDO_0005108"},{"id":"A49","pred":"mondo_id","subj":"T49","obj":"http://purl.obolibrary.org/obo/MONDO_0100096"},{"id":"A50","pred":"mondo_id","subj":"T50","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A51","pred":"mondo_id","subj":"T51","obj":"http://purl.obolibrary.org/obo/MONDO_0100096"},{"id":"A52","pred":"mondo_id","subj":"T52","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A53","pred":"mondo_id","subj":"T53","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A54","pred":"mondo_id","subj":"T54","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A55","pred":"mondo_id","subj":"T55","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"},{"id":"A56","pred":"mondo_id","subj":"T56","obj":"http://purl.obolibrary.org/obo/MONDO_0005812"},{"id":"A57","pred":"mondo_id","subj":"T57","obj":"http://purl.obolibrary.org/obo/MONDO_0005550"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-PD-CLO

    {"project":"LitCovid-PD-CLO","denotations":[{"id":"T47","span":{"begin":311,"end":316},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T48","span":{"begin":465,"end":470},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_9606"},{"id":"T49","span":{"begin":474,"end":479},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_9606"},{"id":"T50","span":{"begin":561,"end":566},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T51","span":{"begin":861,"end":862},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T52","span":{"begin":1048,"end":1053},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T53","span":{"begin":1174,"end":1177},"obj":"http://purl.obolibrary.org/obo/CLO_0053001"},{"id":"T54","span":{"begin":1196,"end":1199},"obj":"http://purl.obolibrary.org/obo/CLO_0053001"},{"id":"T55","span":{"begin":1208,"end":1211},"obj":"http://purl.obolibrary.org/obo/CLO_0009126"},{"id":"T56","span":{"begin":1228,"end":1231},"obj":"http://purl.obolibrary.org/obo/CLO_0053001"},{"id":"T57","span":{"begin":1259,"end":1260},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T58","span":{"begin":1798,"end":1803},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T59","span":{"begin":2290,"end":2295},"obj":"http://purl.obolibrary.org/obo/CLO_0050050"},{"id":"T60","span":{"begin":2670,"end":2671},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T61","span":{"begin":2957,"end":2960},"obj":"http://purl.obolibrary.org/obo/CLO_0053001"},{"id":"T62","span":{"begin":3006,"end":3007},"obj":"http://purl.obolibrary.org/obo/CLO_0001021"},{"id":"T63","span":{"begin":3353,"end":3365},"obj":"http://purl.obolibrary.org/obo/OBI_0000245"},{"id":"T64","span":{"begin":3504,"end":3505},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T65","span":{"begin":3535,"end":3536},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T66","span":{"begin":3587,"end":3588},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T67","span":{"begin":3870,"end":3882},"obj":"http://purl.obolibrary.org/obo/OBI_0000968"},{"id":"T68","span":{"begin":4158,"end":4169},"obj":"http://purl.obolibrary.org/obo/OBI_0000968"},{"id":"T69","span":{"begin":4274,"end":4286},"obj":"http://purl.obolibrary.org/obo/OBI_0000968"},{"id":"T70","span":{"begin":4373,"end":4374},"obj":"http://purl.obolibrary.org/obo/CLO_0001020"},{"id":"T71","span":{"begin":4401,"end":4406},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T72","span":{"begin":4675,"end":4680},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T73","span":{"begin":4788,"end":4793},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T74","span":{"begin":4834,"end":4844},"obj":"http://purl.obolibrary.org/obo/OBI_0000968"},{"id":"T75","span":{"begin":4903,"end":4912},"obj":"http://purl.obolibrary.org/obo/BFO_0000030"},{"id":"T76","span":{"begin":5029,"end":5034},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T77","span":{"begin":5619,"end":5624},"obj":"http://purl.obolibrary.org/obo/NCBITaxon_10239"},{"id":"T78","span":{"begin":5853,"end":5855},"obj":"http://purl.obolibrary.org/obo/CLO_0037066"},{"id":"T79","span":{"begin":6006,"end":6008},"obj":"http://purl.obolibrary.org/obo/CLO_0037066"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-PD-HP

    {"project":"LitCovid-PD-HP","denotations":[{"id":"T6","span":{"begin":529,"end":538},"obj":"Phenotype"},{"id":"T6","span":{"begin":529,"end":538},"obj":"Phenotype"}],"attributes":[{"id":"A6","pred":"hp_id","subj":"T6","obj":"http://purl.obolibrary.org/obo/HP_0002090"},{"id":"A6","pred":"hp_id","subj":"T6","obj":"http://purl.obolibrary.org/obo/HP_0002090"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-sentences

    {"project":"LitCovid-sentences","denotations":[{"id":"T77","span":{"begin":0,"end":15},"obj":"Sentence"},{"id":"T78","span":{"begin":16,"end":82},"obj":"Sentence"},{"id":"T79","span":{"begin":83,"end":171},"obj":"Sentence"},{"id":"T80","span":{"begin":172,"end":265},"obj":"Sentence"},{"id":"T81","span":{"begin":266,"end":425},"obj":"Sentence"},{"id":"T82","span":{"begin":426,"end":494},"obj":"Sentence"},{"id":"T83","span":{"begin":495,"end":754},"obj":"Sentence"},{"id":"T84","span":{"begin":755,"end":912},"obj":"Sentence"},{"id":"T85","span":{"begin":913,"end":1020},"obj":"Sentence"},{"id":"T86","span":{"begin":1021,"end":1144},"obj":"Sentence"},{"id":"T87","span":{"begin":1145,"end":1361},"obj":"Sentence"},{"id":"T88","span":{"begin":1362,"end":1543},"obj":"Sentence"},{"id":"T89","span":{"begin":1544,"end":1844},"obj":"Sentence"},{"id":"T90","span":{"begin":1845,"end":2108},"obj":"Sentence"},{"id":"T91","span":{"begin":2109,"end":2151},"obj":"Sentence"},{"id":"T92","span":{"begin":2152,"end":2366},"obj":"Sentence"},{"id":"T93","span":{"begin":2367,"end":2568},"obj":"Sentence"},{"id":"T94","span":{"begin":2569,"end":2671},"obj":"Sentence"},{"id":"T95","span":{"begin":2672,"end":3007},"obj":"Sentence"},{"id":"T96","span":{"begin":3008,"end":3104},"obj":"Sentence"},{"id":"T97","span":{"begin":3105,"end":3323},"obj":"Sentence"},{"id":"T98","span":{"begin":3324,"end":3463},"obj":"Sentence"},{"id":"T99","span":{"begin":3464,"end":3701},"obj":"Sentence"},{"id":"T100","span":{"begin":3702,"end":3855},"obj":"Sentence"},{"id":"T101","span":{"begin":3856,"end":4090},"obj":"Sentence"},{"id":"T102","span":{"begin":4091,"end":4188},"obj":"Sentence"},{"id":"T103","span":{"begin":4189,"end":4243},"obj":"Sentence"},{"id":"T104","span":{"begin":4244,"end":4525},"obj":"Sentence"},{"id":"T105","span":{"begin":4526,"end":4702},"obj":"Sentence"},{"id":"T106","span":{"begin":4703,"end":4861},"obj":"Sentence"},{"id":"T107","span":{"begin":4862,"end":4894},"obj":"Sentence"},{"id":"T108","span":{"begin":4895,"end":5122},"obj":"Sentence"},{"id":"T109","span":{"begin":5123,"end":5408},"obj":"Sentence"},{"id":"T110","span":{"begin":5409,"end":5589},"obj":"Sentence"},{"id":"T111","span":{"begin":5590,"end":5695},"obj":"Sentence"},{"id":"T112","span":{"begin":5696,"end":5735},"obj":"Sentence"},{"id":"T113","span":{"begin":5736,"end":6097},"obj":"Sentence"}],"namespaces":[{"prefix":"_base","uri":"http://pubannotation.org/ontology/tao.owl#"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}

    LitCovid-PubTator

    {"project":"LitCovid-PubTator","denotations":[{"id":"123","span":{"begin":465,"end":470},"obj":"Species"},{"id":"124","span":{"begin":474,"end":479},"obj":"Species"},{"id":"125","span":{"begin":581,"end":587},"obj":"Species"},{"id":"126","span":{"begin":443,"end":453},"obj":"Disease"},{"id":"127","span":{"begin":503,"end":511},"obj":"Disease"},{"id":"128","span":{"begin":529,"end":538},"obj":"Disease"},{"id":"129","span":{"begin":561,"end":577},"obj":"Disease"},{"id":"135","span":{"begin":1133,"end":1137},"obj":"Gene"},{"id":"136","span":{"begin":1699,"end":1706},"obj":"Species"},{"id":"137","span":{"begin":1332,"end":1340},"obj":"Disease"},{"id":"138","span":{"begin":1516,"end":1526},"obj":"Disease"},{"id":"139","span":{"begin":1561,"end":1569},"obj":"Disease"},{"id":"141","span":{"begin":2523,"end":2524},"obj":"Gene"},{"id":"143","span":{"begin":3140,"end":3150},"obj":"Disease"},{"id":"147","span":{"begin":4819,"end":4823},"obj":"Gene"},{"id":"148","span":{"begin":4493,"end":4499},"obj":"Species"},{"id":"149","span":{"begin":4577,"end":4584},"obj":"Species"},{"id":"152","span":{"begin":5667,"end":5673},"obj":"Species"},{"id":"153","span":{"begin":5111,"end":5121},"obj":"Disease"}],"attributes":[{"id":"A123","pred":"tao:has_database_id","subj":"123","obj":"Tax:9606"},{"id":"A124","pred":"tao:has_database_id","subj":"124","obj":"Tax:9606"},{"id":"A125","pred":"tao:has_database_id","subj":"125","obj":"Tax:9606"},{"id":"A126","pred":"tao:has_database_id","subj":"126","obj":"MESH:D007239"},{"id":"A127","pred":"tao:has_database_id","subj":"127","obj":"MESH:C000657245"},{"id":"A128","pred":"tao:has_database_id","subj":"128","obj":"MESH:D011014"},{"id":"A129","pred":"tao:has_database_id","subj":"129","obj":"MESH:D001102"},{"id":"A135","pred":"tao:has_database_id","subj":"135","obj":"Gene:118"},{"id":"A136","pred":"tao:has_database_id","subj":"136","obj":"Tax:9606"},{"id":"A137","pred":"tao:has_database_id","subj":"137","obj":"MESH:C000657245"},{"id":"A138","pred":"tao:has_database_id","subj":"138","obj":"MESH:D007239"},{"id":"A139","pred":"tao:has_database_id","subj":"139","obj":"MESH:C000657245"},{"id":"A141","pred":"tao:has_database_id","subj":"141","obj":"Gene:4137"},{"id":"A143","pred":"tao:has_database_id","subj":"143","obj":"MESH:D007239"},{"id":"A147","pred":"tao:has_database_id","subj":"147","obj":"Gene:118"},{"id":"A148","pred":"tao:has_database_id","subj":"148","obj":"Tax:9606"},{"id":"A149","pred":"tao:has_database_id","subj":"149","obj":"Tax:9606"},{"id":"A152","pred":"tao:has_database_id","subj":"152","obj":"Tax:9606"},{"id":"A153","pred":"tao:has_database_id","subj":"153","obj":"MESH:D007239"}],"namespaces":[{"prefix":"Tax","uri":"https://www.ncbi.nlm.nih.gov/taxonomy/"},{"prefix":"MESH","uri":"https://id.nlm.nih.gov/mesh/"},{"prefix":"Gene","uri":"https://www.ncbi.nlm.nih.gov/gene/"},{"prefix":"CVCL","uri":"https://web.expasy.org/cellosaurus/CVCL_"}],"text":"Empirical model\nOur analysis sample includes 304 prefecture-level cities in China. We exclude Wuhan, the capital city of Hubei province, from our analysis for two reasons. First, the epidemic patterns in Wuhan are significantly different from those in other cities. Some confirmed cases in Wuhan contracted the virus through direct exposure to Huanan Seafood Wholesale Market, which is the most probable origin of the virus6. In other cities, infections arise from human-to-human transmissions. Second, COVID-19 cases were still pneumonia of previously unknown virus infections in people’s perception until early January so that Wuhan’s health care system became overwhelmed as the number of new confirmed cases increased exponentially since mid-January. This may have caused severe delay and measurement errors in the number of cases reported in Wuhan, and to a lesser extent, in other cities in Hubei province. To alleviate this concern, we also conduct analyses excluding all cities in Hubei province from our sample.\nTo model the spread of the virus, we consider within-city spread and between-city transmissions simultaneously (Adda 2016). Our starting point is yct=∑s=114αwithin,syc,t−s+∑s=114αbetween,s∑r≠cdcr−1yr,t−s+∑s=114ρszt−s+xctβ+𝜖ct, where c is a city other than Wuhan, and yct is the number of new confirmed cases of COVID-19 in city c on date t. Regarding between-city transmissions, dcr is the log of the distance between cities c and r, and ∑r≠cdcr−1yrt is the inverse distance weighted sum of new infections in other cities. Considering that COVID-19 epidemic originated from one city (Wuhan) and that most of the early cases outside Wuhan can be traced to previous contacts with persons in Wuhan, we also include the number of new confirmed cases in Wuhan (zt) to model how the virus spreads to other cities from its source. We may include lagged yct, yrt, and zt up to 14 days based on the estimates of the durations of the infectious period and the incubation period in the literature7. xct includes contemporaneous weather controls, city, and day fixed effects8. 𝜖ct is the error term. Standard errors are clustered by province.\nTo make it easier to interpret the coefficients, we assume that the transmission dynamics (αwithin,s, αbetween,s, ρs) are the same within s = 1,⋯ ,7 and s = 8,⋯ ,14, respectively, but can be different across weeks. Specifically, we take averages of lagged yct, yrt, and zt by week, as y¯ctτ=17∑s=17yct−7τ−1−s, y¯rtτ=17∑s=17yrt−7τ−1−s, and z¯tτ=17∑s=17zt−7τ−1−s, in which τ denotes the preceding first or second week. Our main model is 1 yct=∑τ=12αwithin,τy¯ctτ+∑τ=12αbetween,τ∑r≠cdcr−1y¯rtτ+∑τ=12ρτz¯tτ+xctβ+𝜖ct.Model A\nWe also consider more parsimonious model specifications, such as the model that only considers within-city transmissions, 2 yct=∑τ=12αwithin,τy¯ctτ+xctβ+𝜖ct,and the model where the time lagged variables are averages over the preceding 2 weeks, yct=αwithin114∑s=114yc,t−s+αbetween114∑s=114∑r≠cdcr−1yr,t−s+ρ114∑s=114zt−s+xctβ+𝜖ct.Model B\nThere are several reasons that y¯ctτ, y¯rtτ, and z¯tτ may be correlated with the error term 𝜖ct. The unobserved determinants of new infections such as local residents’ and government’s preparedness are likely correlated over time, which causes correlations between the error term and the lagged dependent variables. As noted by the World Health Organization (2020b), most cases that were locally generated outside Hubei occurred in households or clusters. The fact that big clusters give rise to a large number of cases within a short period of time may still be compatible with a general low rate of community transmissions, especially when measures such as social distancing are implemented. Therefore, the coefficients are estimated by two-stage least squares in order to obtain consistent estimates on the transmission rates in the population.\nIn Eq. 2, the instrumental variables include averages of daily maximum temperature, total precipitation, average wind speed, and the interaction between precipitation and wind speed, for city c in the preceding third and fourth weeks. Detailed discussion of the selection of weather characteristics as instruments is in Section 3.2. The timeline of key variables are displayed in Fig. 1. The primary assumption on the instrumental variables is that weather conditions before 2 weeks do not affect the likelihood that a person susceptible to the virus contracts the disease, conditional on weather conditions and the number of infectious people within the 2-week window. On the other hand, they affect the number of other persons who have become infectious within the 2-week window, because they may have contracted the virus earlier than 2 weeks. These weather variables are exogenous to the error term and affect the spread of the virus, which have been used by Adda (2016) to instrument flu infections9.\nFig. 1 Timeline of key variables\nAnother objective of this paper is to quantify the effect of various socioeconomic factors in mediating the transmission rates of the virus, which may identify potential behavioral and socioeconomic risk factors for infections. For within-city transmissions, we consider the effects of local public health measures (see Section 5 for details) and the mediating effects of population density, level of economic development, number of doctors, and environmental factors such as temperature, wind, and precipitation. For between-city transmissions, apart from proximity measures based on geographic distance, we also consider similarity in population density and the level of economic development. To measure the spread of the virus from Wuhan, we also include the number of people traveling from Wuhan. The full empirical model is as follows:\n3 yct=∑τ=12∑k=1Kwithinαwithin,τkh¯ctkτy¯ctτ+∑τ=12∑k=1Kbetween∑r≠cαbetween,τkm¯crtkτy¯rtτ+∑τ=12∑k=1KWuhanρτkm¯c,Wuhan,tkτz¯tτ+xctβ+𝜖ct,where h¯ctkτ includes dummies for local public health measures and the mediating factors for local transmissions. m¯crtkτ and m¯c,Wuhan,tkτ are the mediating factors for between-city transmissions and imported cases from Wuhan."}