QuartoからStataを使う

仕事でStataを使う必要があったのでメモ

参考にしたページ

Stataをモダンな環境で使う https://zenn.dev/nicetak/articles/stata-modern-2025
nbstata https://hugetim.github.io/nbstata/

環境構築

pythonを使うので、事前にvenvなどで仮想環境をつくっておく
以下のコマンドで必要なパッケージとカーネルをインストールする
- 公式サイトを見るとpip install nbstataだけになっているが、自分の環境では他のも必要だった

pip install nbstata pyyaml nbformat nbclient
python -m nbstata.install

初期設定

jupyterのマジックコマンドのような形で、いくつか初期設定が可能
ここではグラフの画像形式をsvgにしておく
詳しくは
- https://hugetim.github.io/nbstata/user_guide.html#configuration-optional
- https://hugetim.github.io/nbstata/user_guide.html#set-set

%set graph_format = svg

graph_format was png, is now svg

実際に使ってみる

データの読み込み

sysuse census, clear

(1980 Census data by state)

変数一覧を表示

describe


Contains data from /Applications/Stata/ado/base/c/census.dta
 Observations:            50                  1980 Census data by state
    Variables:            13                  6 Apr 2024 15:43
-------------------------------------------------------------------------------
Variable      Storage   Display    Value
    name         type    format    label      Variable label
-------------------------------------------------------------------------------
state           str14   %-14s                 State
state2          str2    %-2s                  Two-letter state abbreviation
region          int     %-8.0g     cenreg     Census region
pop             long    %12.0gc               Population
poplt5          long    %12.0gc               Pop, < 5 year
pop5_17         long    %12.0gc               Pop, 5 to 17 years
pop18p          long    %12.0gc               Pop, 18 and older
pop65p          long    %12.0gc               Pop, 65 and older
popurban        long    %12.0gc               Urban population
medage          float   %9.2f                 Median age
death           long    %12.0gc               Number of deaths
marriage        long    %12.0gc               Number of marriages
divorce         long    %12.0gc               Number of divorces
-------------------------------------------------------------------------------
Sorted by:

チャンクオプションは*| echo: falseのような形（コードブロックが非表示になる）


Contains data from /Applications/Stata/ado/base/c/census.dta
 Observations:            50                  1980 Census data by state
    Variables:            13                  6 Apr 2024 15:43
-------------------------------------------------------------------------------
Variable      Storage   Display    Value
    name         type    format    label      Variable label
-------------------------------------------------------------------------------
state           str14   %-14s                 State
state2          str2    %-2s                  Two-letter state abbreviation
region          int     %-8.0g     cenreg     Census region
pop             long    %12.0gc               Population
poplt5          long    %12.0gc               Pop, < 5 year
pop5_17         long    %12.0gc               Pop, 5 to 17 years
pop18p          long    %12.0gc               Pop, 18 and older
pop65p          long    %12.0gc               Pop, 65 and older
popurban        long    %12.0gc               Urban population
medage          float   %9.2f                 Median age
death           long    %12.0gc               Number of deaths
marriage        long    %12.0gc               Number of marriages
divorce         long    %12.0gc               Number of divorces
-------------------------------------------------------------------------------
Sorted by:

グラフを書いてみる

都市人口率と高齢化率を計算

gen pop65_prop = pop65p / pop
gen popurban_prop = popurban / pop

散布図を表示

scatter pop65_prop popurban_prop