場面設定
当サイトは女性の訪問者が少ないので、女性の訪問者を増やしたいと思っている。
サイト作りの参考にしようと女性の新規訪問の割合が多いランディングページをリストアップしたい。
そこで新規セッション率で降順にソートをかけると、セッション数10〜20くらいのページばかりが上の方にきた。
分母の数が小さいときに割り算値を額面通りに受け取ってよいものか。
そういうときには区間推定をして、信頼下限で判定してやろう。
区間推定では小さく見積もればこれくらい、大きく見積もればこれくらい、という値を分母の数に応じて出すことができる。
実践
信頼度95%の信頼区間を棒で、新規率を点で表したプロットが以下である。
女性の新規率の信頼下限が 0.8 以上のページのみを抜き出している。
点の大きさでセッション数も合わせて表現した。
セッション数が多くなるとき、信頼区間の幅は狭くなることがわかる。
#パッケージをまとめて読み込みlibrary(RGA)library(dplyr)library(tidyr)library(cowplot)#パッケージ読み込み終わり#データ取得 authorize() prof <-list_profiles() dat_ga <-get_ga(profileId = prof$id[1], start.date ="2016-10-01", end.date ="2016-10-31", dimensions ="ga:landingPagePath,ga:userGender", metrics ="ga:sessions,ga:newUsers")#データ取得終わり#信頼区間を計算 CI <-t(mapply(function(x,n)binom.test(x=x,n=n)$conf.int, dat_ga$newUsers,dat_ga$sessions)) LPdf <-data.frame(dat_ga,CI)%>% rename(lower=X1,upper=X2)%>% mutate(percent_new=newUsers/sessions)#女性の新規率の信頼下限が0.8以上のページを抜き出す LPdf085 <-filter(LPdf,landingPagePath %in% landingPagePath[userGender=="female"& lower>=0.8])#以下プロット theme_set(theme_cowplot(font_family ="HiraKakuProN-W3")) ggplot(LPdf085,aes(x=landingPagePath,y=percent_new,colour=userGender))+ geom_pointrange(aes(ymin=lower,ymax=upper),position=position_dodge(width =0.3))+ geom_point(aes(size=sessions),position=position_dodge(width =0.3))+ scale_y_continuous(label=scales::percent)+ coord_flip()+ylab("新規率")+xlab("")+ geom_hline(yintercept=0.9,linetype=2)
新規率以外の比率、直帰率やコンバージョン率でも同じプロットが使える。
補足：シミュレーション
信頼度95%の信頼区間は、100回推定したら内95回くらいはそのなかに真のパラメータを含むだろう。
シミュレーションしてみる。
ここではセッション数がポアソン分布に従ってランダムに変動するとして、二項分布の成功確率パラメータの信頼区間を繰り返し構成した。
N <-rpois(10000,100) x <- rbinom(10000,N,0.8) CIsim <-t(mapply(function(x,n)binom.test(x=x,n=n)$conf.int, x,N)) cat("coverage probability is",round(mean(CIsim[,1]<=0.8& CIsim[,2]>=0.8),2))
Coverage probability（信頼区間がその範囲内に真値を含む割合）は0.96と95%に近い値になった。
- 作者:永田靖
- 出版社/メーカー:日本経済新聞社
- 発売日: 2006/01
- メディア:新書
- 購入: 2人 クリック: 7回
- この商品を含むブログ (5件) を見る