2

以前は相互作用プロットを作成していましたがggplot2、コードを以下に示します。ggvisここで、以下に示すように、出力と同じではない同じプロットを再現したいと思いggplto2ます。で同じプロットを取得するにはどうすればよいggvisですか?

library(ggplot2)
p <- qplot(as.factor(dose), len, data=ToothGrowth, geom = "boxplot", color = supp) + theme_bw()
p <- p + labs(x="Dose", y="Response")
p <- p + stat_summary(fun.y = mean, geom = "point", color = "blue", aes(group=supp))
p <- p + stat_summary(fun.y = mean, geom = "line", aes(group = supp))
p <- p  + theme(axis.title.x = element_text(size = 12, hjust = 0.54, vjust = 0))
p <- p  + theme(axis.title.y = element_text(size = 12, angle = 90,  vjust = 0.25))
print(p)

ここに画像の説明を入力

library(ggvis)
ggvis(data=ToothGrowth, x= ~as.factor(dose), y= ~len, fill= ~supp, stroke = ~supp) %>% 
  layer_points(shape=~supp) %>% 
  layer_lines(fillOpacity=0)

ここに画像の説明を入力

4

1 に答える 1

7

でこれを実装しようとするときの基本的な問題は、のようなオプションggvisがないため、異なる値の箱ひげ図を同じ座標にプロットできないことです。したがって、軸にインデックスを付けることはオプションではないようです。ただし、一意の線量値の数に等しい長さの整数インデックスを使用し、値に応じてデータの位置を手動で左または右にオフセットすることができます。position = dodgeggplot2suppxxas.factor(dose)xsupp

library(ggvis)
library(dplyr)
d <- ToothGrowth
d$xpos <- as.integer(factor(d$dose)) + ifelse(d$supp == "OJ", -.2, .2)

x = ~xposこれで、ボックスプロットを正しい位置にプロットするために使用できるようになりました。次のステップは、線で結ばれた点をプロットするために使用される平均を保持するデータを定義することです。

means <- d %>% group_by(dose, supp) %>% summarise(len = mean(len))
means$xpos <- as.integer(factor(means$dose))
means <- group_by(means, supp) # The grouping is needed for layer_paths()

グラフは次のように取得できます。

ggvis(d, x = ~xpos, y = ~len, stroke = ~supp) %>% 
    layer_boxplots() %>%
    layer_points(data = means, fill := "blue") %>%
    layer_paths(data = means)

ここxで、プロットの位置が実際の線量値ではなく 1、2、3 になるという問題があります。軸の目盛りにラベルを付け直す方法がないため、これを克服するのは簡単ではありませんadd_axis()(また、最初に 1、2、3 の代わりに実際の線量値を使用することはできませんでした。線量値 0.5 と 1 は、線量値 1 と 2 よりも互いに近い値です)。これは、個々の線量値ごとに軸を追加するという、それほど洗練されていないハックによって克服できます。この関数add_axis()は、軸のプロパティ (ラベルを含む) を変更する方法を提供しますが、プロパティは軸全体に適用されるため、軸全体に同じラベルを使用します。したがって、線量値ごとに軸を追加することで、ラベルを 1 つずつ操作できます。これは次のように見えます

ggvis(d, x = ~xpos, y = ~len, stroke = ~supp) %>% 
    layer_boxplots() %>%
    layer_points(data = means, fill := "blue") %>%
    layer_paths(data = means) %>%
    add_axis("x", title = "Dose", 
        values = c(1, 1), # For some reason values of length 1 don't work...
        properties = axis_props(labels = list(text = "0.5"))) %>%
    add_axis("x", title = "", 
        values = c(2, 2), 
        properties = axis_props(labels = list(text = "1"))) %>%
    add_axis("x", title = "", 
        values = c(3, 3), 
        properties = axis_props(labels = list(text = "2"))) %>%     
    add_axis("y", title = "Response")

または、これらにループを使用して、同じことを何度も入力する必要がないようにすることもできます

labs <- data.frame(dose = unique(d$dose))
labs$xpos <- as.integer(factor(labs$dose))

v <- ggvis(d, x = ~xpos, y = ~len, stroke = ~supp) %>% 
    layer_boxplots() %>%
    layer_points(data = means, fill := "blue") %>%
    layer_paths(data = means) %>%
    add_axis("x", title = "Dose", ticks = 0) %>%
    add_axis("y", title = "Response")

for (i in 1:nrow(labs)) {
    v <- add_axis(v, "x", title = "", values = rep(labs[i, "xpos"], 2),
        properties = axis_props(labels = list(text = as.character(labs[i, "dose"]))))
}

最終結果はこんな感じ

ここに画像の説明を入力

于 2014-12-03T11:29:12.207 に答える