myscratchpad

再現率はpsychパッケージのp.rep関数を使う
この関数で計算される効果量 (dprime) やp.repはhttp://psychology.clas.asu.edu/sqabのものとなんか違う。そもそも計算式が違う。
というか、t分布で近似するときと正規分布で近似するときがあるみたい。ようわからんな…。p値の再現性は正規分布で、効果量の再現性はt分布なのか？

                SS num Df Error SS den Df         F       Pr(>F)
(Intercept) 3444.5      1       77      7 313.13636 4.537639e-07
afact        217.5      3       71     21  21.44366 1.347309e-06

p.rep(p = 0.05, n=NULL,twotailed = FALSE)
p.rep.f(F,df2,twotailed=FALSE)
p.rep.r(r,n,twotailed=TRUE)
p.rep.t(t,df,df2=NULL,twotailed=TRUE)

引数と出力だけメモ

引数
p F、tなどの統計量の有意確率。普通のp値
F F値
df t検定の自由度。サンプルサイズが違うときは第1群の自由度
df2 F値の分母の自由度。もしくはサンプルサイズの違うt検定の第2群の自由度
r 相関係数
n rの算出に使った全体のサンプルサイズ
t t検定のときにはt統計量、もしくは回帰係数の有意性検定の統計量
twotailed 片側か両側か


p.rep  再現率 (Probability of replication)
dprime  効果量 (Cohen`s d) 。一定以上のpが得られる場合if more than just p is specified
prob  F, t, rの確率。片側か両側
r.equivalent t検定でのrとtの等質性

参考文献
Cummings, Geoff (2005) Understanding the average probability of replication: comment on Killeen 2005). Psychological Science, 16, 12, 1002-1004).
Killeen, Peter H. (2005) An alternative to Null-Hypothesis Significance Tests. Psychological Science, 16, 345-353.
Rosenthal, R. and Rubin, Donald B.(2003), r-sub(equivalent): A Simple Effect Size Indicator. Psychological Methods, 8, 492-496.
Rosnow, Ralph L., Rosenthal, Robert and Rubin, Donald B. (2000) Contrasts and correlations in effect-size estimation, Psychological Science, 11. 446-453.

Society for Quantitative Analysis of Behavior
http://psychology.clas.asu.edu/sqab