r - Stata边距/边距图转换为R？

我正在尝试将 Stata marginsplot 复制到 R 中，但即使在浏览 StackExchange 并试图弄清楚几周之后也无法这样做。

您是否知道如何重新创建使用 R 中的 marginsplot 创建的图？

首先，我使用以下 R 代码生成了一个可重现的数据集:

set.seed(96)
sampleData <- data.frame(
    ID = 1:200,
    outcome = sample(1:7, 200, replace = T),
    scale = sample(1:7, 200, replace = T),
    dummy1 = sample(0:1, 200, replace = T),
    dummy2 = sample(0:1, 200, replace = T))
write.dta(dataframe, "sampledata.dta")

然后，我将数据读入 Stata 并运行以下命令以获取附加图:

   reg outcome dummy1##c.scale dummy2##c.scale

   margins dummy1 dummy2, at(scale=(1(1)7))

   marginsplot, recast(line) name(outcome, replace)

本质上，当这段代码在 Stata 中运行时，我得到了我试图在 R 中复制的图形:



你有什么指示可以帮助我复制这张图吗？

最佳答案

# your data
set.seed(96)
sampleData <- data.frame(
    ID = 1:200,
    outcome = sample(1:7, 200, replace = T),
    scale = sample(1:7, 200, replace = T),
    dummy1 = sample(0:1, 200, replace = T),
    dummy2 = sample(0:1, 200, replace = T))

# all possible combinations
newData <- data.frame(scale=rep(1:7, each=4),
    dummy1=rep(c(0, 0, 1, 1), 7),
    dummy2=rep(0:1, 2*7))
newData$dummy <- paste(newData$dummy1, newData$dummy2)

a <- lm(outcome ~ dummy1*scale + dummy2*scale, data=sampleData)

b <- predict(a, newdata=newData, interval="confidence")

c <- cbind(newData, b)

require(ggplot2)

# one pane
d <- ggplot(c, aes(scale, fit, color=dummy))
d <- d + geom_line()
# d <- d + facet_grid(dummy1 ~ dummy2)
d <- d + geom_ribbon(aes(ymin=lwr, ymax=upr), alpha=0.25)
d

# using facets
e <- ggplot(c, aes(scale, fit))
e <- e + geom_line()
e <- e + facet_grid(dummy1 ~ dummy2)
e <- e + geom_ribbon(aes(ymin=lwr, ymax=upr), alpha=0.25)
e

我应该包括情节。

关于r - Stata边距/边距图转换为R？，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/38518332/