Я хотел бы нарисовать доверительный эллипс (иногда также называемый эллипсом концентрации) для двух непрерывных переменных с использованием Stata.
Для этого вида графиков существует команда сообщества , названная ellip, подробное описание которой предоставлено ее автором Андерсом Александерссоном в 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?