高校数学の質問スレ（医者・東大卒専用） Part438 (899ﾚｽ)
高校数学の質問スレ（医者・東大卒専用） Part438 http://rio2016.5ch.net/test/read.cgi/math/1723152147/
上
下
前次
1-
新
通常表示
512ﾊﾞｲﾄ分割
ﾚｽ栞
抽出解除
必死ﾁｪｯｶｰ(本家)
(べ)
自ID
ﾚｽ栞
あぼーん
861: １３２人目の素数さん [sage] 2025/06/02(月) 11:14:31.22 ID:GMuHFUYr 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) #summary(results) -results$coef[1]/results$coef[2] library(MASS) # mvrnorm を使うため # 推定された係数と共分散行列 beta_hat = coef(results) (vcov_matrix = vcov(results)) # 多変量正規乱数を生成（β0, β1） set.seed(123) # 再現性のため samples = mvrnorm(n = 10000, mu = beta_hat, Sigma = vcov_matrix) # 各サンプルから LD50 を計算 LD50_samples = -samples[,1] / samples[,2] # 信頼区間（95%） CI = quantile(LD50_samples, probs = c(0.025, 0.975)) # 結果表示 cat("シミュレーションによるLD50の95%信頼区間:\n") print(CI) http://rio2016.5ch.net/test/read.cgi/math/1723152147/861
862: １３２人目の素数さん [sage] 2025/06/02(月) 12:04:32.64 ID:GMuHFUYr # データ x <- c(-0.86, -0.3, -0.05, 0.73) n <- c(5, 5, 5, 5) y <- c(0, 1, 3, 5) # JAGSモデル model_string <- " model { for (i in 1:N) { y[i] ~ dbin(p[i], n[i]) logit(p[i]) <- beta0 + beta1 * x[i] } # 事前分布（非情報的） beta0 ~ dnorm(0.0, 0.001) beta1 ~ dnorm(0.0, 0.001) # LD50の定義 LD50 <- -beta0 / beta1 } " # JAGSに渡すデータ data_jags <- list( x = x, n = n, y = y, N = length(y) ) # 初期値 inits <- function() { list(beta0 = rnorm(1, 0, 1), beta1 = rnorm(1, 0, 1)) } # モデル作成と実行 model <- jags.model(textConnection(model_string), data = data_jags, inits = inits, n.chains = 3) update(model, 1000) # バーンイン # サンプリング samples <- coda.samples(model, variable.names = c("beta0", "beta1", "LD50"), n.iter = 10000) # 結果表示（LD50の95%信用区間） summary(samples) LD50_samples <- as.matrix(samples)[, "LD50"] quantile(LD50_samples, probs = c(0.025, 0.975)) plot(samples) http://rio2016.5ch.net/test/read.cgi/math/1723152147/862
ﾒﾓ帳
(0/65535文字)
上
下
前次
1-
新
書
関
写
板
覧
索
設
栞
歴
ｽﾚ情報
赤ﾚｽ抽出
画像ﾚｽ抽出
歴の未読ｽﾚ
AAｻﾑﾈｲﾙ
Google検索
Wikipedia
ぬこの手
ぬこTOP
1.529s*