高校数学の質問スレ(医者・東大卒専用) Part438 (898レス)
高校数学の質問スレ(医者・東大卒専用) Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前次
1-
新
通常表示
512バイト分割
レス栞
抽出解除
必死チェッカー(本家)
(べ)
自ID
レス栞
あぼーん
863: 132人目の素数さん [sage] 2025/06/03(火) 00:31:29.12 ID:CfA5PBxZ x = c(-0.86, -0.3, -0.05, 0.73) n = c(5, 5, 5, 5) y = c(0, 1, 3, 5) (data = cbind(x, n, y)) (response = cbind(y, n - y) ) results = glm(response ~ x, family = binomial(link="logit") ) summary(results) b=MASS::mvrnorm(1e5,results$coef,S=vcov(results)) # log(p/(1-p)) = b1+b2*x LD = function(b1,b2,p=0.5) (log(p/(1-p)) - b1)/b2 LD50=mapply(LD,b[,1],b[,2]) hist(LD50) quantile(LD50,c(0.025,0.5,0.975)) HDInterval::hdi(LD50) http://rio2016.5ch.net/test/read.cgi/math/1723152147/863
864: 132人目の素数さん [sage] 2025/06/03(火) 23:40:25.34 ID:CfA5PBxZ k=10 m=400 Nmax=1000 sim = \(){ N=sample(m:Nmax,1) if(max(sample(1:N,k))==m) return(N) } N=NULL while(length(N)<1000) N=c(N,sim()) summary(N) quantile(N,c(0.5,0.95)) http://rio2016.5ch.net/test/read.cgi/math/1723152147/864
メモ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
Google検索
Wikipedia
ぬこの手
ぬこTOP
0.029s