0

我想使用 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 中创建加权置信椭圆？

4

0 回答 0