高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞
ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡
100: 2024/08/19(月)23:40:01.05 ID:TuWh8Fdq(1) AAS
他全滅
177: 2024/08/25(日)16:40:08.05 ID:1ZiEG71e(6/6) AAS
>>176
それが理解できないアホだからここで発狂してんだよ
躾ければ決まったところでトイレができる犬以下
220(1): 2024/09/11(水)23:03:13.05 ID:cP+fZTFu(1) AAS
この数学教師マジか…
元小樽桜陽高校数学科教師堤伸弘の数学教師としての能力
外部ﾘﾝｸ:note.com
252: 2024/10/27(日)09:42:05.05 ID:HIz5euo6(1/3) AAS
# 一峰性非対称分布する離散量t1の95% highest probability intervalを返す
t2=sort(table(t1),decreasing=TRUE)
t3=as.numeric(names(t2))
t4=t2/sum(t2)
t5=sum(cumsum(t4)<0.95)
range(t3[1:t5])
sum(t2[1:t5]/sum(t2))
range(t3[1:(t5+1)])
sum(t2[1:(t5+1)]/sum(t2))
451: 2024/12/22(日)12:28:31.05 ID:PT9UhnTA(1) AAS
>>447
なに自分で自演してんの？馬鹿なの？
それとも自分がFラン卒だって自白したかっただけ？
840: 05/20(火)23:31:13.05 ID:gwaBTE4C(1) AAS
library(R2jags)
# データ
data <- list(
nA1 = 100, rA1 = 80, # Study1: 治療A
nB1 = 100, rB1 = 40, # Study1: 治療B
nA2 = 100, rA2 = 10, # Study2: 治療A
nC2 = 100, rC2 = 5 # Study2: 治療C
)
# JAGSモデル（textConnection使用）
model_code <- "
model {
# 尤度関数
rA1 ~ dbin(pA1, nA1)
rB1 ~ dbin(pB1, nB1)
rA2 ~ dbin(pA2, nA2)
rC2 ~ dbin(pC2, nC2)
# 治療Aの階層モデル
mu_A ~ dbeta(1, 1)
tau_A ~ dgamma(0.001, 0.001)
pA1 ~ dbeta(mu_A * tau_A, (1 - mu_A) * tau_A)
pA2 ~ dbeta(mu_A * tau_A, (1 - mu_A) * tau_A)
sigma_A <- 1 / sqrt(tau_A) # SDに変換
# 治療BとCも階層化（平均リスクを別々に推定）
mu_B ~ dbeta(1, 1)
mu_C ~ dbeta(1, 1)
pB1 ~ dbeta(mu_B * 100, (1 - mu_B) * 100) # 高い精度を仮定
pC2 ~ dbeta(mu_C * 100, (1 - mu_C) * 100)
# リスク差
RD_A1_B1 <- pA1 - pB1
RD_A2_C2 <- pA2 - pC2
RD_B1_C2 <- pB1 - pC2 # B vs Cの直接比較
}
"
# JAGS実行
jags_model <-
(textConnection(model_code),
data = data, n.chains = 3, quiet=TRUE)
update(jags_model, 3000) #, progress.bar="none")
jags_samples <- coda.samples(jags_model,
c("mu_A", "sigma_A", "RD_A1_B1", "RD_A2_C2", "RD_B1_C2", "pA1", "pA2", "pB1", "pC2"),
n.iter=10000) # , progress.bar="none")
gelman.plot(jags_samples)
plot(jags_samples)
summary(jags_samples)
jags_samples |> as.matrix() |> as.data.frame() -> js
names(js)
925: 06/22(日)13:39:18.05 ID:DRKE+jFQ(2/5) AAS
>>912
じゃあアンタが解いてみろよ
解けないアンタはfランだからさっさと消えろ
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ
ぬこの手 ぬこTOP 0.039s