R 当观测值下降时,天文学家分离变量
我使用stargazer为多个模型创建一个表。它们实际上是同一个模型,但第一个模型基于所有观测值,而另一个模型则分别基于不同的观测值。所有变量的名称都是相同的,所以让我惊讶的是,当我将表导出到Latex时,两行(一行用于虚拟变量,另一行用于交互项)是重复的 真正奇怪的是,我无法复制结果,但我将发布一个最小的工作示例。也许你可以根据我的描述来帮助我 这是我的MWE的代码:R 当观测值下降时,天文学家分离变量,r,latex,stargazer,R,Latex,Stargazer,我使用stargazer为多个模型创建一个表。它们实际上是同一个模型,但第一个模型基于所有观测值,而另一个模型则分别基于不同的观测值。所有变量的名称都是相同的,所以让我惊讶的是,当我将表导出到Latex时,两行(一行用于虚拟变量,另一行用于交互项)是重复的 真正奇怪的是,我无法复制结果,但我将发布一个最小的工作示例。也许你可以根据我的描述来帮助我 这是我的MWE的代码: library(tibble) library(stargazer) df <- as_tibble(data.fra
library(tibble)
library(stargazer)
df <- as_tibble(data.frame(first = rnorm(100, 50), second = rnorm(100, 30), third = rnorm(100, 100), fourth = c(rep(0, 50), rep(1, 50))))
model.1 <- lm(first ~ second + third + fourth + third*fourth, data = df)
model.2 <- lm(first ~ second + third + fourth + third*fourth, data = df[!rownames(df) %in% "99",])
stargazer(model.1, model.2)
现在,我将发布Latex输出,其中包括我试图用此代码段修复的错误,它似乎工作正常
当然,我想要的是这段代码所产生的代码,因为无法重现它,我感到非常愚蠢:
您可以使用系数查看模型系数的名称。确保它们是相同的,即相同的名称模型1,名称模型2,然后使用stargazer的keep语句确保得到所需的系数 在这里,以上面的例子保持选定的变量
coefficients(model.1)
#> (Intercept) second third fourth third:fourth
#> 57.27352606 0.02674072 -0.08236250 20.23596216 -0.20288137
coefficients(model.2)
#> (Intercept) second third fourth third:fourth
#> 57.06149556 0.03305134 -0.08214812 20.85087288 -0.20885718
identical(names(model.1), names(model.2))
#> [1] TRUE
我使用type=text使它对SO更友好,但我想LaTeX也是这样
stargazer(model.1, model.2, type = "text", keep=c("third","third:fourth"))
#>
#> =========================================================
#> Dependent variable:
#> -------------------------------------
#> first
#> (1) (2)
#> ---------------------------------------------------------
#> third -0.082 -0.082
#> (0.166) (0.167)
#>
#> third:fourth -0.203 -0.209
#> (0.222) (0.223)
#>
#> ---------------------------------------------------------
#> Observations 100 99
#> R2 0.043 0.044
#> Adjusted R2 0.002 0.004
#> Residual Std. Error 1.044 (df = 95) 1.047 (df = 94)
#> F Statistic 1.056 (df = 4; 95) 1.089 (df = 4; 94)
#> =========================================================
#> Note: *p<0.1; **p<0.05; ***p<0.01
但是,如果我们无法找到复制您的问题的方法,可能很难排除这是本地问题。好的,实际上解决我问题的方法是重新运行整个脚本。但我确信,如果这不起作用,你的解决方案会给我提供答案。非常感谢。