スタガのマージン/マージンをRに入れる？

Question

私はRにStataののmarginsplotを複製しようとしていますが、それでも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に複製しようとしているグラフを取得します

はあなたが私はこのグラフを複製役立つだろう任意のポインタを持っていますか？

Answer 1

# 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 

プロットを含める必要があります。