2表の共通オッズ比の推定で、独立な積二項モデルの仮定のもとではCMLEもMH推定量も一致性があるが、層内で相関がある(ベータ二項モデルなど)とMH推定量しか一致性を持たない。2020年2月14日と言ったしたものの確認したことがなかったのでシミュレーションしてみました。マルコフ連鎖よりもベータ二項分布の方が簡単そうだったので、積ベータ二項分布モデルについて検討しています(原典と微妙な違いがあるかもしれません)。100$のときを検討します。mh:Mantel–Haenszel推定量、cmle:条件付き最尤推定量、mle:無条件の最尤推定量です。Mantel–Haenszel推定量は$K$が大きければ$\phi=5$に近づいていきます。条件付き最尤推定量は$K$が大きくなっても$\phi=5$から離れたままです。無条件の最尤推定量はもともと一致性がありません。おまけ:無相関の積ベータ二項分布モデル相関がない積ベータ二項分布モデルは積二項分布モデルに等しく、Mantel–Haenszel推定量も条件付き最尤推定量も一致性があります。Procedure14.2に続いて因果推論系のプロシジャが追加されていますね。Procedure)。分解した推定ができるようです。PHモデルなど実務で使いそうなものがサポートされています。differenceの推定値、信頼区間が出力されるようになった。testのexact版が追加された。Procedure)。Statement)。Statement)。Procedure)。新情報が結構まとまっています。Procedures)。Procedure)。Procedure)。Enhancements)。C-statistic)、経時AUC・各時点のROC曲線の図示。(STAT14.1までは変量効果がある場合にNTHREADSを指定しても無視されていたらしいです)。いうパッケージ内に含まれた独自のものに変更され、他のパッケージからの拡張が容易になったようです。公式の解説があるので、これを熟読すると良いでしょう。されたので、かなり簡単に拡張が書けるようです。継承し、前処理を少し加えたものを作成しました。いう引数を追加しています。いう関数を内部で利用していて、この関数の出力を使う場合の例外処理の有無を指定できるようにしたかったので入れました。GeomStepribbon実際に処理している本体の部分です。継承すれば良いと思います。追加引数名の指定です。方で追加引数を入れた場合はここで指定していないとエラーになります。データのセットアップ処理の部分です、今回はここに例外処理を入れています。内部で順番に継承している他の処理と上述の処理が呼ばれていくのですが、今回は一部で事足りてしまいました。今後は今回触れていない部分を勉強してみようと思っています。以下のようになりました。合わせて使うために作っています。変えようかと思っています。コメントを含めて読むと、仕組みがよくわかりました。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。