Cohenのdの算出用いるにpooled varianceを計算するスクリプト。

pooledVariance <- function(...) {
args = list(...)
n.args=length(args)
ss2=0
df=0
for(i in 1:n.args){
ss2 = ss2 + var(args[[i]])*(length(args[[i]])-1)
df = df + (length(args[[i]])-1)
}
ss2/df
}

effectsize <- function(y1,y2){
diff=mean(y1)-mean(y2)
var=(var(x1)*(length(x1)-1)+ var(x2)*(length(x2)-1))/(length(c(y1,y2))-2)
sd=sqrt(var)
diff/sd
}

パッケージを探したらやっぱりあった。
library(effsize)
?cohen.d
パッケージはエラー処理などが整ってて便利だが
coreの部分の計算は自分でもやってみないと
一次方程式も立てられいド底辺卒のようなクズになるからね。