2016-08-19 1 views
1

私はstargazerパッケージをRの回帰出力に使用しています。私は、モデルオブジェクトにはなりませんが、係数と標準誤差のベクトルのみをもたらすカスタマイズされた推定手順を持っています。 stargazerにこれらを供給し、きれいにフォーマットされた出力テーブルを得る方法はありますか?stargazer - ユーザー提供の係数とSE

例:

dep.var <- "foo" 
regressors <- c("bar", "baz", "xyz") 
vec.coeffs <- c(1.2, 2.3, 3.4) 
vec.se <- c(0.1, 0.1, 0.3) 

出力が似ているはずです:

=============================================== 
         Dependent variable:  
        --------------------------- 
           foo    
----------------------------------------------- 
bar       1.200***     
           (0.100)   

baz       2.300***   
           (0.100) 

xyz       3.400***   
           (0.300)   

----------------------------------------------- 

答えて

2

ここに1つの提案です:主なアイデアは、など、偽lmオブジェクトを作成し、次にカスタム係数を適用し、SEのことです。stargazerへ出力:

d <- as.data.frame(matrix(rnorm(10 * 4), nc = 4)) 
names(d) <- c(dep.var, regressors) 
f <- as.formula(paste(dep.var, "~ 0 +", paste(regressors, collapse = "+"))) 
p <- lm(f, d) 

stargazer(p, type = "text", 
    coef = list(vec.coeffs), 
    se = list(vec.se), 
    t = list(vec.coeffs/vec.se), 
    omit.stat = "all") 
# ================================= 
#   Dependent variable:  
#  --------------------------- 
#     foo    
# --------------------------------- 
# bar   1.200***   
#     (0.100)   

# baz   2.300***   
#     (0.100)   

# xyz   3.400***   
#     (0.300)   

# ================================= 
# ================================= 
# Note: *p<0.1; **p<0.05; ***p<0.01 
関連する問題