如何用权重创建置信椭圆？

Question

我想使用Stata绘制两个连续变量的置信椭圆（有时也称为浓度椭圆）。

对于此类图，有一个称为ellip的社区贡献命令，其作者Anders Alexandersson在Stata Journal上提供了详细说明。但是，我想为椭圆应用权重，这是使用此命令无法实现的。

下面是一个可重现的示例，在该示例中，我尝试使用人口规模作为权重，将每个1000居民的婚姻数量与美国各州城市人口的百分比作图：

clear all
sysuse census

gen marpop = marriage / pop * 1000
gen urbpop = popurban / pop * 100

ellip marpop urbpop if state2!="NV", plot(scatter marpop urbpop if state2!="NV") ///
    name(ellip_noaw)
ellip marpop urbpop if state2!="NV", plot(scatter marpop urbpop [aw=pop] if state2!="NV") ///
    name(ellip_aw)

graph combine ellip_noaw ellip_aw

结果如下：

如您所见，第二个椭圆（散点图的点由总体大小加权/膨胀）可能应该移到右下角。

如何在Stata中创建加权置信椭圆？