抑制R中的固定效应系数

抑制R中的固定效应系数,r,summary,lm,R,Summary,Lm,在使用summary()函数时,是否有方法抑制线性模型中固定效应的系数(例如stata中吸收()函数的等效值)。例如,我希望summary函数只输出截距和x,而不是系数和因子的标准误差: frame <- data.frame(x = rnorm(100), y = rnorm(100), z = rep(c("A", "B", "C", "D"),25)) summary(lm(y~x + as.factor(z), data = frame)) Call: lm(formula =

在使用summary()函数时,是否有方法抑制线性模型中固定效应的系数(例如stata中吸收()函数的等效值)。例如,我希望summary函数只输出截距和x,而不是系数和因子的标准误差:

frame <- data.frame(x = rnorm(100), y = rnorm(100), z = rep(c("A", "B", "C", "D"),25))
summary(lm(y~x + as.factor(z), data = frame)) 

Call:
lm(formula = y ~ x + as.factor(z), data = frame)

Residuals:
Min      1Q  Median      3Q     Max 
-2.2417 -0.6407  0.1783  0.5789  2.4347 

Coefficients:
          Estimate Std. Error t value Pr(>|t|)  
(Intercept)   -0.25829    0.19196  -1.346   0.1816  
 x              0.09983    0.09788   1.020   0.3104  
frame | t |)
(截距)-0.25829 0.19196-1.346 0.1816
x 0.09983 0.09788 1.020 0.3104


谢谢。

您可以执行以下任一操作:

mod <- lm(y~x + as.factor(z), data = frame)
coef(mod)[c("(Intercept)", "x")]
# (Intercept)           x 
#  0.12357491 -0.06430765 
coef(mod)[grepl("Intercept|x", names(coef(mod)))]
# (Intercept)           x 
#  0.12357491 -0.06430765 
coef(mod)[1:2]
# (Intercept)           x 
#  0.12357491 -0.06430765 
mod$coefficients[1:2]
# (Intercept)           x 
#  0.12357491 -0.06430765 

mod我强烈推荐
stargazer
套装:

frame <- data.frame(x = rnorm(100), y = rnorm(100), z = rep(c("A", "B", "C", "D"),25))
summary(lm(y~x + as.factor(z), data = frame)) 
model <- lm(y~x + as.factor(z), data = frame)
library(stargazer)
stargazer(model, type="text")
stargazer(model, omit = "z", type="text")

frame是否有显示整个输出的内容?罗伯特凡桑德斯说:“有什么比仅仅获取系数更普遍的东西吗?”我回答了眼前的问题。OP想要拦截和x。如果你有问题,你可以向社区提出。谢谢你的否决票。我会改的,但他的问题语无伦次。首先,他说他只想抑制固定效应,但有总结,然后他说他只想有截距和系数。我认为他的问题是我要找的第一件事。我正在寻找如何做第一件事,但这些问题的许多答案都没有做到,直到我找到了“lfe”包。似乎最好在问题下发表评论,然后等待答复,以消除这种困惑。OP检查了这个答案,所以它解决了他眼前的问题。我不太清楚你要我怎么做。我的错,我会把它改回来的。(做一个小编辑,这样我就可以)我只是否决了没有做第一部分的答案。令人惊讶的是,很难找到。。。!