Understanding the dynamic impact of COVID-19 through competing risk modeling with bivariate varying coefficients