14

おそらく単純な質問ですが、Field の「Discovering Statistics Using R」の 532 ページに示されているようなエラーグラフを作成しようとしました。

コードはhttp://www.sagepub.com/dsur/study/DSUR%20R%20Script%20Files/Chapter%2012%20DSUR%20GLM3.Rにあります。

line <- ggplot(gogglesData, aes(alcohol, attractiveness, colour = gender))
line + stat_summary(fun.y = mean, geom = "point") + 
stat_summary(fun.y = mean, geom = "line", aes(group= gender)) + 
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = 0.2) + 
labs(x = "Alcohol Consumption", y = "Mean Attractiveness of Date (%)", colour = "Gender")

同じグラフを作成しました。私のy軸変数には4点しかありません（離散スケール、1-4です）が、y軸には点1.5、2、2.5があり、線が変化します。

問題は、これらの点とグラフが何を表しているかということです。重要な部分はstat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = 0.2)、そのグループとそのレベル（x軸）の観測数ですか？それらは周波数ですか？それともプロポーションですか？

私はこれを見つけましたhttp://docs.ggplot2.org/0.9.3/stat_summary.htmlしかし、それは私を助けませんでした

ありがとうございました

4

2 に答える 2

15

これは、 83ページのggplot2の本に書かれていることですmean_cl_boot()

Function          Hmisc original        Middle Range
mean_cl_boot() smean.cl.boot() Mean Standard error from bootstrap

smean.cl.boot()Hmisc パッケージからのものだと思いますがmean.cl.boot()、ggplot2 のように名前が変更されています。

Hmiscパッケージの元の関数の定義は次のとおりです。

smean.cl.bootは、正規性を仮定せずに母平均の信頼限界を取得するための基本的なノンパラメトリック ブートストラップの非常に高速な実装です。

于 2013-07-01T23:09:13.990 に答える