さまざまな因子レベル (サンプル) と加重平均 (重み = カバレッジ) の変数 (対立遺伝子固有の式) を表示しようとしています。
私はいくつかのサンプルデータを作成しました:
set.seed(2)
x <- sample(c("A","B","C"), 100, replace=T)
y <- rnorm(100)
w <- ceiling(rnorm(100,200,200))
df <- data.frame(x, y, w)
library(ggplot2)
ggplot(df, aes(x=factor(x), y=y, weight=w)) +
geom_point(aes(size=w)) +
stat_summary(fun.y=mean, colour="red", geom="point", size=5)
（そして、プロットも投稿しようとしましたが、まだ十分なポイントがありません）。
これは問題なく動作しますが、重み付けされていない平均を示しています...
library(plyr)
means <- ddply(df, "x", function(x) data.frame(wm=weighted.mean(x$y, x$w),
m=mean(x$y)))
means
x wm m
1 A 0.00878432 0.11027454
2 B -0.07283770 -0.13605530
3 C -0.14233389 0.08116117
そのため、代わりに「wm」値を赤い点として表示しようとしています-ggplot2を使用しています。「weight = ..」を正しく使用しているに違いないと思いますが、今はあきらめています...
誰かが助けてくれることを本当に願っています。