Rでマージンエフェクトテーブルを作成するには？

Question

私は順序付き従属変数の範囲が0から2までの3つの順序付き回帰モデルを持っています。私がしたいのは、3つのモデルすべての各レベル（0,1,2）で限界効果テーブル（プロットではありません）しかし

stargazer(M1,M2, M3, type = ”text”) 

を：だから、3つのレベル0、1で限界効果を示す各持つテーブル、および2

## create a random data that is similar to my dataset 
    set.seed(100) 
    mydata <- data.frame(
     x1 = sample(c(0, 1, 2), 100, replace = TRUE), 
     x2 = sample(c(0, 1, 2, 3, 4), 100, replace = TRUE), 
     x3 = sample(c(0, 1, 2, 3, 5), 100, replace = TRUE), 
     x4 = sample(c(1:100), 100, replace = TRUE), 
     x5 = sample(c(10:1000), 100, replace = TRUE), 
     Z1 = sample(c(0, 1, 2), 100, replace = TRUE) 
    ) 


    ## makeit factor 
    mydata$Z1 <- as.factor(mydata$Z1) 

    ## My models 
    require(MASS) 

    M1<- polr(Z1 ~x1+x2+x3+x4, data=mydata, Hess = TRUE, method="logistic") 

    M2<- polr(Z1 ~x2+x3+x4+x5, data=mydata, Hess = TRUE, method="logistic") 

    M3<- polr(Z1 ~x1+x2+x3+x4+x5, data=mydata, Hess = TRUE, method="logistic") 

    ## Calculate marginal effects using the erer package 
    require(erer) 

    M1ME<- ocME(M1) 

    M2ME <- ocME(M2) 

    M3ME <- ocME(M3) 

は通常、私が使用して、たとえば、適切なテーブルを作成するために、パッケージstargazerを使用します、OcME（）からの出力は、テーブルの同じタイプを生成しませんし、また私は、各レベルでテーブルを生成することができます stargazer(M1ME$out,M2ME$out, M3ME$out, type = "text")

• あなたはありますかどのようにこれらのタイプの テーブルを生成するための任意の提案？これを手動で行うには非常に時間がかかります。

EDIT： だから望ましい出力は、以下に示す3つのテーブルを作成することである。

限界効果（注番号は、単に例示正しくない） - レベル0（Z）

========================================== 
         Dependent variable:  
       ----------------------------- 
           Z1    
        (1)  (2)  (3) 
    ------------------------------------------ 
    x1    0.301    0.302 
        (0.250)    (0.250) 

    x2    0.143  0.174  0.142 
        (0.131) (0.128) (0.132) 

    x3    0.121  0.106  0.122 
        (0.117) (0.116) (0.117) 

    x4   -0.008 -0.008 -0.008 
        (0.007) (0.007) (0.007) 

    x5      -0.00004 -0.0001 
          (0.001) (0.001) 

    ------------------------------------------ 
    Observations 100  100  100 
    ========================================== 
    Note:   *p<0.1; **p<0.05; ***p<0.01 

限界効果 - レベル1（Z）

========================================== 
        Dependent variable:  
      ----------------------------- 
          Z1    
       (1)  (2)  (3) 
------------------------------------------ 
x1    0.301    0.302 
       (0.250)    (0.250) 

x2    0.143  0.174  0.142 
       (0.131) (0.128) (0.132) 

x3    0.121  0.106  0.122 
       (0.117) (0.116) (0.117) 

x4   -0.008 -0.008 -0.008 
       (0.007) (0.007) (0.007) 

x5      -0.00004 -0.0001 
         (0.001) (0.001) 

------------------------------------------ 
Observations 100  100  100 
========================================== 
Note:   *p<0.1; **p<0.05; ***p<0.01 

限界効果--level 3（Z）

========================================== 
        Dependent variable:  
      ----------------------------- 
          Z1    
       (1)  (2)  (3) 
------------------------------------------ 
x1    0.301    0.302 
       (0.250)    (0.250) 

x2    0.143  0.174  0.142 
       (0.131) (0.128) (0.132) 

x3    0.121  0.106  0.122 
       (0.117) (0.116) (0.117) 

x4   -0.008 -0.008 -0.008 
       (0.007) (0.007) (0.007) 

x5      -0.00004 -0.0001 
         (0.001) (0.001) 

------------------------------------------ 
Observations 100  100  100 
========================================== 
Note:   *p<0.1; **p<0.05; ***p<0.01 

Answer 1

のGithub上の新しい発達パッケージには、この時にビットを取得します。 Thomas Leeperによると、かなりここに投稿しています。

https://github.com/leeper/margins

お楽しみください！