高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除 ﾚｽ栞

---

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

---

48(1): 2024/08/15(木)17:59 ID:bFfiJSUV(2/3) AAS
n=9;
j[n_] :=(
count=0;
Until[Length@Union@a==2,a=RandomChoice[Range[3],n];count++];
b=Sort@Union@a;
If[b=={1,2}, w=Count[a,2]];
If[b=={2,3}, w=Count[a,3]];
If[b=={1,3}, w=Count[a,1]] ;
{w,count}
)

For[{w,count}=j[n],w>1k=j[w];w=k[[1]];count=count+k[[2]]]
{w,count}

---

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.040s