一昨日(こちら)、たくさんの受験者の結果を全部使うと、使わない場合より、個々人の真の正答力を、より正確に推定できるという話を書いた より正確に推定する方法をJames-Stein Estimatorと言った 「ふーん」と思いやすい例はないかと考えた こんな風に考え…

こちらから 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){ # 変なデータとして、正規分布と指数分…