2013-01-23から1日間の記事一覧

二次的現象なのか

特定の要因がなくても発生する事象があって、ある要因があるとその事象が引き起こされるような印象がある。その印象はどんなときにどのくらいもたらされる? 印象を強める要素は何だろう? N<-20 # 単位期間にDという有害事象が起きる確率 pD<-0.1 D<-runif(…

大小不同の瞬発力

貧血における赤血球の大小不同や悪性細胞の大小不同の定量の瞬発力 癌の細胞診(参考) # 一視野あたりの細胞数 n.cell <- 10 # 半径の平均はすべてのスライドでそろえよう r.mean <- 1 # 結構ばらばらな2x2分割表を作る library(MCMCpack) library(shape) …

放っておけば相乗的

my.2x2.dev <- function(p,q,r){ #tmp <- r*sqrt(p*q*(1-p)*(1-q)) x <- c(p*q,p*(1-q),(1-p)*q,(1-p)*(1-q)) tmp <- r*min(x) matrix(x+c(tmp,-tmp,-tmp,tmp),byrow=TRUE,2,2) } my.2x2x2.dev <- function(p,q,r){ tmp <- r*sqrt(prod(p)*prod(q)*r) tmp1 <…

あとだしじゃんけんを見破る

グリコを赤・黒の2人でやっている 赤が後だししているかもしれない ただし、赤は後だしをするゲームとしないゲームを分けている しかも、赤はいつも後だしをするわけではなく、何度かに1度だけ後だしをしているらしい 10回、グリコを繰り返して記録をつけ…

密度分布と累積密度分布を結びつける

n <- list() x <- list() k <- 3 library(MCMCpack) for(i in 1:k){ n[[i]] <- rdirichlet(1,rep(0.1,3))*1000 x[[i]] <- runif(n[[i]][1]) x[[i]] <- c(x[[i]],rnorm(n[[i]][2],0.5)) x[[i]] <- c(x[[i]],rexp(n[[i]][3])) } par(mfcol=c(2,k)) s <- sample…

2x2表:ぱっと見て感じる

2x2分割表のクイズ 実行すると、2x2表に対応する棒グラフがいくつか並べて提示される 図を先送りすると、そのフィッシャーの正確確率検定p値が表示され、最小p値の棒グラフがハイライトされることで、答え合わせができる Q『次の3つの棒グラフに対応…

瞬発力クイズ