底辺私立医大を卒業した医者って頭悪いよね？ Part10

底辺私立医大を卒業した医者って頭悪いよね？ Part10at HOSP

底辺私立医大を卒業した医者って頭悪いよね？ Part10 - 暇つぶし2ch456:j・最頻値なしで信頼区間でけ求めるのはカッコ悪いなぁ。ド底辺シリツ裏口調査団が100人を順次調査した。裏口判明人数をそのまま公表はヤバすぎる結果であったため、連続して裏口がみつかった最大数は4人であったとだけ公表した。公表結果が正しいとして裏口入学人数の期待値、最頻値、及び95%信頼区間を述べよ。スクリプトに計算追加して期待値・最頻値も表示させることにした。 seq2pCI <- function(N,K,alpha=0.05){ vp=Vectorize(function(p)seqNp(N,K,p)-seqNp(N,K+1,p)) curve(vp(x),lwd=2,bty='l') ; abline(h=0.05,lty=3) peak=optimize(vp,c(0,1),maximum=TRUE)$maximum mean=integrate(function(x)x*vp(x),0,1)$value/integrate(function(x)vp(x),0,1)$value lwr=uniroot(function(x,u0=alpha) vp(x)-u0,c(0.01,peak))$root upr=uniroot(function(x,u0=alpha) vp(x)-u0,c(peak,0.99))$root c(lower=lwr,mean=mean,mode=peak,upper=upr) } seq2pCI(100,4)

次ページ

続きを表示

1を表示