抑制R中的固定效应系数
在使用summary()函数时,是否有方法抑制线性模型中固定效应的系数(例如stata中吸收()函数的等效值)。例如,我希望summary函数只输出截距和x,而不是系数和因子的标准误差:抑制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 =
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检查了这个答案,所以它解决了他眼前的问题。我不太清楚你要我怎么做。我的错,我会把它改回来的。(做一个小编辑,这样我就可以)我只是否决了没有做第一部分的答案。令人惊讶的是,很难找到。。。!