2013-08-27から1日間の記事一覧

マルチプルテスティング:James-Stein推定

こちらから my.JS.estimate <- function(x){ m <- mean(x) A1 <- (length(x)-2)/sum((x-m)^2) y <- (1-A1)*(x-m) +m return(y) } n.iter <- 1000 N <- 1000 A <-pi R.As <- R.Bs <- rep(0,n.iter) for(i in 1:n.iter){ # 変なデータとして、正規分布と指数分…