Gene-set enrichment analysisの一環としてこちらでコルモゴロフ=スミルノフ検定(KS検定)の話がちらっと出た KS検定は１標本と２標本の区別がある 観察された標本の元の分布が、ある特定の分布ではない、と考えた方がよいかを判断する(１標本のKS検定) もし…

こちらにもあるように。 ２標本が同一母集団から取られているとすると、それらを併せて順序をつけると、第１群、第２群としたときに1,2,2,1,1,1,1,1,2,2,1,2,2,みたいな1,2の数値列ができるけれど、この1,2がばらばらなのか、という検定をすることでも検定可…

累積分布の比較には特徴がある スタートは0でゴールは1。必ず そのため、二つの累積分布の差を見るとき、スタートが0でゴールも0になる 経験分布２つ(もしくは経験分布と理論分布との２つ)の差を取ると、「0で始まって0で終わる」 横軸は「0で始まって1で終…

前の記事では、正規分布と指数分布との混合の場合をやった x.norm(n) x.exp(n) の混成だった。ソートしてプロットすれば、完全になめらかとは言わないが、一峰性 そこで、二峰性にしてみよう x.norm(n,-10) x.exp(n) の混成に。図にすれば、こう いわゆる密…

まず、Rのks.test()関数の中身をぱくって、２標本の場合のks.test()統計量の算出のところだけを抜き出した関数を作っておこう my.ks.stat <- function(x,y,type=c("less","greater","two.sided")){ n.x <- length(x) n.y <- length(y) w <- c(x, y) z <- cum…