計算された統計を使用してggplot2 boxplotでx軸を注文する

Question

ggplot2を使用してボックスプロットとして表示するデータがあります。基本的には2つの変数で層別化されています。ここでは（実際にそこに多くのですが、構造は同じである）データの例です：各タグの

TAG Count Condition 
A  5   1 
A  6   1 
A  6   1 
A  6   2 
A  7   2 
A  7   2 
B  1   1 
B  2   1 
B  2   1 
B 12   2 
B  8   2 
B 10   2 
C 10   1 
C 12   1 
C 13   1 
C  7   2 
C  6   2 
C 10   2 

、条件1の観測値の固定された数があり、条件2（ここではそれはです3、実際のデータではそれ以上です）。私は次のようなボックスプロット（「s」は、上記のように配置されたデータフレームである）欲しい：

ggplot(s, aes(x=TAG, y=Count, fill=factor(Condition))) + geom_boxplot() 

これは結構ですが、私はによりx軸を注文できるようにしたいです各タグのWilcoxon検定からのp値。例えば、上記のデータを用いて、値が（それぞれ、タグAのために、B、およびC）のようになります。

x軸上に順序付けA、C、Bを誘導する

> wilcox.test(c(5,6,6),c(6,7,7))$p.value 
[1] 0.1572992 
> wilcox.test(c(1,2,2),c(12,8,10))$p.value 
[1] 0.0765225 
> wilcox.test(c(10,12,13),c(7,6,10))$p.value 
[1] 0.1211833 

（最大に最小）。しかし、私はこの情報を自分のデータに追加する方法を知らない（具体的には、余分な列を追加するのではなく、単にタグレベルでp値を付ける）か、x軸注文。どんな助けでも大歓迎です。

Answer 1

ここは方法です。第1のステップは、各TAGのp値を計算することである。これを行うには、ddplyを使用してデータをTAGで分割し、数式インタフェースを使用してwilcox.testを使用してp値を計算します。 plot文は、そのp値に基づいてTAGを並べ替えます。

library(ggplot2); library(plyr) 
dfr2 <- ddply(dfr, .(TAG), transform, 
    pval = wilcox.test(Count ~ Condition)$p.value) 

qplot(reorder(TAG, pval), Count, fill = factor(Condition), geom = 'boxplot', 
    data = dfr2)