R で zeligverse パッケージを使用して二変量ロジット モデルを実行しています。P(Y1=1)、P(Y2=1)、P(Y1=1、Y2=0) に対する独立変数の影響を計算したいと考えています。 P(Y1=1,Y2=1)、P(Y1=0,Y2=1)、P(Y1=0,Y2=0)、P(Y1=1|Y2=0)、その他すべての条件付き確率 ( Y1 と Y2 は従属変数で、どちらも 0 または 1 です)。また、独立変数ごとに、これらの確率に関連するすべての限界効果も必要です。
このパッケージ (または、より適切に機能する場合は別のパッケージ) でそれらを見つける方法を知っていますか?