34

因子変数に従ってファセットを導入することにより、単純なフォレストプロットの例を変更しようとしています。

この構造のデータを想定すると：

test <- structure(list(characteristic = structure(c(1L, 2L, 3L, 1L, 2L
), .Label = c("Factor1", "Factor2", "Factor3"), class = "factor"), 
    es = c(1.2, 1.4, 1.6, 1.3, 1.5), ci_low = c(1.1, 1.3, 1.5, 
    1.2, 1.4), ci_upp = c(1.3, 1.5, 1.7, 1.4, 1.6), label = structure(c(1L, 
    3L, 5L, 2L, 4L), .Label = c("1.2 (1.1, 1.3)", "1.3 (1.2, 1.4)", 
    "1.4 (1.3, 1.5)", "1.5 (1.4, 1.6)", "1.6 (1.5, 1.7)"), class = "factor"), 
    set = structure(c(1L, 1L, 1L, 2L, 2L), .Label = c("H", "S"
    ), class = "factor")), .Names = c("characteristic", "es", 
"ci_low", "ci_upp", "label", "set"), class = "data.frame", row.names = c(NA, 
-5L))

そしてコードを実行します：

p <- ggplot(test, aes(x=characteristic, y=es, ymin=ci_low, ymax=ci_upp)) + geom_pointrange() +
  coord_flip() + geom_hline(aes(x=0), lty=2) + 
  facet_wrap(~ set, ncol = 1) +
  theme_bw() + 
  opts(strip.text.x = theme_text())

次のような出力を生成します。

ここに画像の説明を入力してください

これまでのところすべて良い。ただし、下のパネルから空のFactor3レベルを削除したいのですが、その方法が見つかりません。それを行う方法はありますか？

手伝ってくれてありがとう。

4

2 に答える 2

33

編集をggplot20.9.3に更新

別の解決策があります。とを使用facet_gridしspace = "free"ます; geom_point()また、とを使用geom_errorbarh()するため、の必要はありませんcoord.flip()。また、x軸の目盛りラベルは下部パネルにのみ表示されます。以下のコードでは、themeコマンドは必須ではありません。ストリップテキストを回転させて水平に表示するために使用されます。上記のデータフレームを使用するtestと、次のコードで必要なものが生成されます。

library(ggplot2)

p <- ggplot(test, aes(y = characteristic, x = es, xmin = ci_low, xmax = ci_upp)) +
   geom_point() +
   geom_errorbarh(height = 0) +
   facet_grid(set ~ ., scales = "free", space = "free") +
   theme_bw() +
   theme(strip.text.y = element_text(angle = 0))

p

このソリューションは、Wickhamのggplot2ブックの124ページの例に基づいています。

于 2012-04-18T23:38:43.843 に答える
27

次のように使用scales = "free"します。

p <- ggplot(test, aes(x=characteristic, y=es, ymin=ci_low, ymax=ci_upp)) + geom_pointrange() +
  coord_flip() + geom_hline(aes(x=0), lty=2) + 
  facet_wrap(~ set, ncol = 1, scales="free") +
  theme_bw() + 
  opts(strip.text.x = theme_text())

p

生成するもの：

ここに画像の説明を入力してください

編集：私は実際、次のようdrop = TRUEにこのソリューションの議論の方が好きだと思います：

p <- ggplot(test, aes(x=characteristic, y=es, ymin=ci_low, ymax=ci_upp)) + 
  geom_pointrange() +
  coord_flip() + geom_hline(aes(x=0), lty=2) + 
  facet_wrap(~ set, ncol = 1,  drop=TRUE) +
  theme_bw() + 
  opts(strip.text.x = theme_text())

p
于 2012-04-17T00:07:29.927 に答える