In reply to Ksharp.

Hi KSharp. Thanks for writing. My code for the funnel plot assumes that you want to compare each region's rate to the average rate over all regions. In the China data, the Wuhan province experienced many more cases and deaths than other provinces. That causes the overall rate to be very high, which is why the non-Wuhan rates are displayed as much lower than the average.

As I say in the article, "you can create a similar graph by using some other rate, such as ... target rate." Let's say you want to compare the China data to the North Carolina rate of 0.02. In the IML step that computes the funnel curves, replace the line
theta = sum(Events) / sum(Trials);
with
theta = 0.02;
Also specify
%let AvgProp = 0.02;
Now when you create the funnel plot, many provinces are inside the funnel, and Wuhan is an extreme outlier.