マルチレベル解析におけるモデリング効果のプロットR

Question

私のデータセットには3つのグループがあり、グループとyの相互作用をプロットしたいと思います。

id <- c(1,1,1,2,2,2,3,3,3) 
group <- c(0,0,0,1,1,1,2,2,2) 
x <- c(20,50,30,50,65,80,20,50,60) 
y <- c(120,130,150,200,210,180,160,170,120) 

私はすでに "interaction.plot"で試してみましたが、動作しませんでした。

interaction.plot(x,group,y) 

誰かがこの相互作用をプロットするための良いR構文を持っていますか？

Answer 1

グループのclassを希望するものとは関係なくfactorに変更する必要があります。実際にはinteraction.plot()は要素の双方向の組み合わせであり、xは要因ではありません。しかし、1つが連続している場合は、interaction.plot()が助けを与えます。あなたの場合、出力は「そのようなデータとのやりとりを考えるのは愚かな考えです」と言っています。

df <- data.frame(id = id, x = x, y = y, group = as.factor(group)) 

## Base plot 
model <- lm(y ~ x * group, data = df) 
xpara <- 20:80 

plot(y ~ x, data = df, col=c(2:4)[group], pch=19) 
for(i in 1:3) lines(xpara, predict(model, data.frame(x = xpara, group = as.factor(i-1))), col = i+1) 
legend("topleft",paste(c("group0","group1","group2")), pch=19, lty=1, col=c(2:4)) 

## ggplot2 (I plotted lines and confidence intervals to interpret) 
library(ggplot2) 
ggplot(df, aes(x = x, y = y, colour = group)) + 
    geom_point(size = 4) + 
    geom_smooth(method = "lm", se = T, fullrange = T) 

[編集]モデルのクラスは predict()でサポートされている場合

、道：あなたがしたい場合は

は、しかし、（私はあなたが線形モデルを望んでいたはず）基本的に同じです。

df2 <- data.frame(id = as.factor(id), x = x, y = y, group = as.factor(group)) 
library(nlme) 

# first; make model 
lme.mod <- lme(y ~ x * group, random = ~ 1|id, data = df2) 

# second; get predicted values 
xpara <- 20:80 # make a vector for an independent variable you use as x. 
y.g1 <- predict(lme.mod, data.frame(x = xpara, group = "0", id = "1"), type="response") 
y.g2 <- predict(lme.mod, data.frame(x = xpara, group = "1", id = "1"), type="response") 
y.g3 <- predict(lme.mod, data.frame(x = xpara, group = "2", id = "1"), type="response") 

# third; draw 
plot(y ~ x, df2, col=c(2:4)[group], pch=19) 
lines(xpara, y.g1, col=2) 
lines(xpara, y.g2, col=3) 
lines(xpara, y.g3, col=4) 

## Simplificated version 
lev <- levels(df$group) 

plot(y ~ x, data = df2, col=c(2:4)[group], pch=19, ylab="y (id = "1")") 
for(i in seq.int(length(lev))) 
    lines(xpara, predict(lme.mod, data.frame(x = xpara, group = lev[i], id = 1)), col = i+1) 
legend("topleft",paste(c("group0","group1","group2")), pch=19, lty=1, col=c(2:4))