2013-03-11から1日間の記事一覧

Sidak法からBonferroni法へ

N回独立な検定を行ったところ、p値の最小値としてpminを得たとする。これよりも「珍しい」事象は「p値の最小値として、pminよりも小さい値を得た」である。これは「全てのp値がpminよりも大きかった」事象の余事象だから、その確率はで与えられる。多重検定…

ベクトル化されたif

以下の関数は、ベクトルに対しては期待した動作をしない。 mysqrt <- function(x) { if (x > 0) sqrt(x) else 0 } そこで、ifelse を使って mysqrt2 <- function(x) { ifelse(x > 0, sqrt(x), 0) } のようにすれば動く。ただ、sqrt(x) が x のすべての要素に…