Stata

アドヒアランスの指標であるMPRとPDCを計算するパッケージ[Stata]

The Stata Journalを読んでいたらMPRとPDCを計算するパッケージの紹介があったので、こちらで紹介しておこう。 https://journals.sagepub.com/doi/pdf/10.1177/1536867X19893625 http://www.lindenconsulting.org/stats.html NPRとは薬剤保持率(medication …

SPSS、SAS、StataのデータをRでインポート、エクスポートをする

SPSS、SAS、StataのデータをRでインポート、エクスポートをするhavenパッケージの紹介である。 havenパッケージはtidyverseを構成するパッケージの一つである。 他の統計パッケージのデータを読み込む有名なパッケージはforeignだが、RStudioブログによると…

Stata LCA プラグインを使用した潜在クラス分析 導入偏

Stata15から標準機能で潜在クラス分析ができるようになっている。 ただ、基本的な機能しかないため、Stataの標準機能の潜在クラス分析は現在のところ実用に耐えられない。そこで、面倒な作業が必要になるが、Lanzaらが作成したStata LCA プラグインを使わざ…

潜在クラス分析の応答確率の棒グラフ[Stata]

Stataでは潜在クラス分析の応答確率の棒グラフが描ける。棒グラフというのは今まで見たことがなかった。不要な気もしないでもないが、説明するときにわかりやすいのかもしれない。 使用するのはmarginsとmarginsplotである。説明によると条件付き応答確率の…

潜在クラス分析[Stata]

Stataで潜在クラス分析ができるのでやってみた。 Stata15から標準機能で利用できるらしい。 Lanzaらのプラグインを使えばStata11から利用可能である。 www.methodology.psu.edu Stataの標準機能だと機能が足りないように思うので、Lanzaらのパッケージを結局…

ベイジアンロジスティック回帰分析・その3

HardinとHilbeのStata本にあるベイジアンロジスティック回帰分析の続き。今回は事前分布としてコーシー分布を使う方法についてである。 1回目エントリはこちら 2回目エントリはこちら Generalized Linear Models and Extensions: Fourth Edition作者: Jame…

ベイジアンロジスティック回帰分析・その2

HardinとHilbeのStata本にあるベイジアンロジスティック回帰分析の続き。前回のエントリはこちら Generalized Linear Models and Extensions: Fourth Edition作者: James W. Hardin,Joseph M. Hilbe出版社/メーカー: Stata Press発売日: 2018/06/28メディア:…

有病率の計算[Stata]

疫学調査で有病率を出す方法である。専門用語では母比率の推定という。日本語で母比率と言うとピンとこないが、英語だとpopulation rateなので、こちらの方が直感的に理解できるはずである。 Stataで有病率を推定する Stataではcii propコマンドを使用する。…

ベイジアンロジスティック回帰分析[Stata]

HardinとHilbeが書いたglmのStata本をちらちらと読んでいて少し気になったので、自分でも分析をしてみた。 Generalized Linear Models and Extensions: Fourth Edition作者: James W. Hardin,Joseph M. Hilbe出版社/メーカー: Stata Press発売日: 2018/06/28…

ベイジアンP値

重回帰分析 Stataでデータwomenwageを読み込む。女性の就労のデータのようだ。 . use http://www.stata-press.com/data/r16/womenwage.dta, clear 回帰分析をする。初めて使うデータで詳しくないし、詳しくないし、分析内容はわりと適当。有意になる変数と有…

相関係数の比較[Stata]

Stataでも相関係数の計算をしてみた。こちらのエントリ尺度水準に適した相関係数とシミュレーションをStataで行ったバージョンである。 データはこちらからダウンロード、もしくは下部にRでの作成方法を掲載しているので、そのままRで走らせると、dta形式の…

StataでLassoとリッジ回帰

Stata16でLassoが新しく使えるようになったそうだ。 https://www.lightstone.co.jp/stata/stata16_new.html#Lasso 僕はStataはあまり使わないので、新しいStataを入手していない。今手元にあるものはバージョン13で少し古く、標準機能でLassoはできない。た…