Collinearity appears to be the problem.
I currently have a panel data set that contains the quarterly increment of loans initiated from more than 300 cities in China over the period from 2011Q1 to 2020Q2. I want to examine the impact of COVID-19 on lending activities.
SARS-Cov-2 is assumed to affect the lending practices of all cities within China at precisely the same time. Your event dummy 'turns on' in all cities in the first and second quarters of 2020. Note, your event dummy has the same 'within-city' pattern. Your first model should look something like this:
$$
\text{Lending}_{cq} = \beta_{0} + \beta_{1}\text{Event}_{q} + \alpha_{c} + \epsilon_{cq},
$$
where you observe the lending activities in city $c$ across quarters $q$. $\text{Event}_{q}$ is $q$-subscripted. The dummy exhibits no variation across cities. Demeaning your equation (or estimating dummies for all cities) 'sweeps out' $a_c$ (i.e., the 'city' effect). This model is estimable.
However, adding time (quarter) fixed effects introduces collinearity concerns. Suppose you have quarterly data on 300 cities from 2011 onward; I assume your panel ends in the second quarter of 2020. Now say you wish to estimate a two-way fixed effects estimator using the least-squares dummy variables (LSDV) approach; this is algebraically equivalent to estimation in deviations from means. This should result in 300 separate city effects and 42 separate quarter effects. The quarter fixed effects represent the common shocks, absent the pandemic, affecting all cities. By construction, your event dummy is collinear with your last two quarter effects.
My guess is that the impact of the dummy variable "event" has already been absorbed by time fixed effects. Is that right?
Precisely. Again, $\text{Event}_{q}$ is equal to 1 in all cities in the first and second quarter of 2020. If you include a model with quarter fixed effects, then you have a unique 2020-Q1 dummy and a unique 2020-Q2 dummy specific to all cities. In other words, $\text{Event}_{q}$ is the sum of the last two quarter effects. To break the collinearity, SAS dropped your event dummy. In most software packages, the order of your variables matters inside of a regression function call. If you want an estimate for your event dummy, then I recommend dropping the last quarter. In R/Stata, for example, if you specified your event indicator before the city and quarter effects, it will default to removing your last quarter dummy (i.e., 2020-Q2). You could do this manually in SAS to obtain an estimate for your event (i.e., exposure/intervention) dummy.
In a branched setting, instead of using the dummy variable "event", I use two dummy variable, "event_level1" and "event_level2", which equals to 1 if the period is 2020Q1 and 2020Q2 respectively, to visualize the dynamic of COVID-19 impact.
As already shown, you cannot estimate the dynamics of exposure with city and quarter fixed effects. Such a model is inestimable. Your event dummies are perfectly collinear with your last two quarter effects. Thus, a two-way fixed effects estimator will fail in this setting. Your event indicators must vary over time and across units to use this estimator.
In order to estimate city and quarter fixed effects, you require either a ‘staggered’ exposure (i.e., different cities become exposed at different times) or a setting where a subset of cities were never treated. Because SARS-Cov-2 is a population-level exposure affecting all cities/provinces in precisely the same time periods, then quarter fixed effects will completely absorb your event dummy.
You can safely estimate your model using one-way city-level fixed effects or pooled ordinary least squares.
The quote below was reproduced from the comments section:
You give me some hints that I can work around multicollinearity by making the event dummy interacted with other variables, like the treated dummy or control variables, instead of including a pure event dummy ... I can adopt a two-way fixed effects model and only include interaction terms of the event dummy with other variables.
This seems appropriate.
I suppose you can call this "working around" the collinearity problem. Collinearity isn't fatal, though. As you suggested, interacting the event dummy with a covariate(s) of interest is one way to proceed. See the equation below, where I include both city (i.e., $\alpha_{c}$) and quarter (i.e., $\tau_q$) fixed effects:
$$
\text{Lending}_{cq} = \beta_{0} + \beta_{1} (\text{Event}_{q} \times x_{cq}) + \beta_{2} x_{cq} + \alpha_{c} + \tau_q + \epsilon_{cq}.
$$
Here, I introduced a covariate $x_{cq}$, which should be of some substantive interest. Note how the individual event dummy is not important, hence why it doesn't appear in the foregoing equation. What matters with respect to this specification is the effect of $x_{cq}$ on lending practices and whether it is different before and after the pandemic. This is a different question entirely. Again, we don't care about how the pandemic (i.e., the "event") affects lending activities, but rather how our covariate affects lending activities and how this is different pre- versus post-pandemic. The event dummy will be dropped, for reasons already outlined above, but this variable, by itself, isn't important anyway.