高校数学の質問スレ（医者・東大卒専用） Part438

高校数学の質問スレ（医者・東大卒専用） Part438 (979ﾚｽ)
上下前次1-新
通常表示 512ﾊﾞｲﾄ分割ﾚｽ栞
抽出解除必死ﾁｪｯｶｰ(本家) (べ) 自ID ﾚｽ栞あぼーん

ﾘﾛｰﾄﾞ規制です｡10分ほどで解除するので､他のﾌﾞﾗｳｻﾞへ避難してください｡

563: １３２人目の素数さん [sage] 2025/01/14(火) 18:27:07.83 ID:KQCqu9ZF

"
Fくんが２日前（４８時間とする）に喫食した食材でノロウイルスによる食中毒が発生したとする。
https://i.imgur.com/Aaz7iOi.png
のデータから現在無症状のFくんが今後、発症する確率を算出しなさい。
算出に必要な条件は適宜設定してよい。
"

rm(list=ls())
library(fitdistrplus)
solve=\(hours,distr=c("lnorm","gamma","weibull"),verbose=FALSE){
noro=c(
runif(2,21,24),
runif(2,25,28),
runif(16,29,32),
runif(16,33,36),
runif(14,37,40),
runif(2,41,44),
runif(5,45,48)
)
aic=NULL
for(i in 1:length(distr)) aic=c(aic,gofstat(fitdist(noro,distr[i]))$aic)
best_fit=distr[which.min(aic)]
fit=fitdist(noro,best_fit)
par=fit$estimate
if(verbose){
hist(noro,freq=F,main='Noro Virus incubation',col=4,axes=FALSE,xlab='hours',ylab='') ; axis(1)
curve(dgamma(x,par[1],par[2]),add=TRUE)
}
paste0('p',best_fit,'(hours,par[1],par[2],lower.tail = FALSE)') |> str2lang() |> eval()

}
solve(48,verbose=TRUE)

k=1e4
p48=replicate(k,solve(48,c("lnorm","gamma","weibull")))
round(HDInterval::hdi(p48),3)[1:2]
median(p48)
mean(p48)
BEST::plotPost(p48,showCurve = F,col='lightgreen',xlab='p[hour>48]')

http://rio2016.5ch.net/test/read.cgi/math/1723152147/563

上下前次1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.025s