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

---

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

---

563: 01/14(火)18:27 ID:KQCqu9ZF(1) AAS
"
Fくんが２日前（４８時間とする）に喫食した食材でノロウイルスによる食中毒が発生したとする。
画像ﾘﾝｸ

のデータから現在無症状の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]')

---

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

---

ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ

---

ぬこの手 ぬこTOP 0.039s