- 【ウハも】 開業医達の集い 28診 【粒も】
707 ：卵の名無しさん[sage]：2020/05/18(月) 12:07:13.22 ID:t88GB6TH - >>705
裏口容疑者って不親切なのか？
|
- 【ウハも】 開業医達の集い 28診 【粒も】
708 ：卵の名無しさん[sage]：2020/05/18(月) 12:07:39.21 ID:t88GB6TH - >>706
インバウンド頼みじゃないからだろうね。
|
- 【ウハも】 開業医達の集い 28診 【粒も】
709 ：卵の名無しさん[sage]：2020/05/18(月) 12:12:50.21 ID:t88GB6TH - >裏口容疑者って不親切なのか？
お金が大好きだから親切だけど、馬鹿だから>699の答は出せない。
自分で出せなれば子供に計算してもらえばいいのに。
そうか、息子も裏口容疑者か（爆笑ｗｗｗｗ）
|
- 【ウハも】 開業医達の集い 28診 【粒も】
712 ：卵の名無しさん[sage]：2020/05/18(月) 13:02:22.87 ID:t88GB6TH - >>699
# 裏口容疑者発症のdelay日後にキャバ嬢が発症
# キャバ嬢→裏口容疑者感染の確率
をstanを使ってMCMCしてグラフにしてみた。
https://i.imgur.com/LuMtiVY.png
|
- 臨床統計もおもしろいですよ、その3
168 ：卵の名無しさん[sage]：2020/05/18(月) 13:03:08.05 ID:t88GB6TH - stancode=
"
data{
real onset_delay;
real ln_par1;
real ln_par2;
}
parameters{
real <lower=0> d_incubation;
real <lower=0> c_incubation;
}
transformed parameters{
real infection_delay = onset_delay + d_incubation - c_incubation;
}
model{
d_incubation ~ lognormal(ln_par1,ln_par2);
c_incubation ~ lognormal(ln_par1,ln_par2);
}
"
model=stan_model(model_code = stancode)
fn.stan <- function(delay){
dataList=list(onset_delay=delay,ln_par1=ln.par1,ln_par2=ln.par2)
fit=sampling(model,data=dataList)
ms=rstan::extract(fit)
mean(ms$infection_delay < 0)
}
fn.stan(2)
onset_delays=0:20
y=sapply(onset_delays,fn.stan)
plot(onset_delays,y, ylab='Pr[ Infected Later ])',axes=F) ; axis(1)
|
- 【ウハも】 開業医達の集い 28診 【粒も】
713 ：卵の名無しさん[sage]：2020/05/18(月) 13:05:06.02 ID:t88GB6TH - >>712
stan のコードも埋め込んだRのスクリプトはここに置いた。改良歓迎。
https://egg.5ch.net/test/read.cgi/hosp/1583407025/168
|
- 【ウハも】 開業医達の集い 28診 【粒も】
715 ：卵の名無しさん[sage]：2020/05/18(月) 13:53:12.38 ID:t88GB6TH - >>712
正確には、キャバ嬢の方が裏口容疑者より先に感染していた確率だな。
感染させたかどうかは別の話になるから。
|
- 【ウハも】 開業医達の集い 28診 【粒も】
716 ：卵の名無しさん[sage]：2020/05/18(月) 13:56:39.78 ID:t88GB6TH - >>714
>699の答をサクッとかけば、裏口容疑者でもRが使えるんだなと少しは見直したんだが。
|
- 底辺私立医大を卒業した医者って頭悪いよね？ Part16
464 ：卵の名無しさん[sage]：2020/05/18(月) 19:34:28.64 ID:t88GB6TH - 学生時代には女子医大の学生には息子が大変お世話になりましたなぁ。
|
- 【ウハも】 開業医達の集い 28診 【粒も】
721 ：卵の名無しさん[sage]：2020/05/18(月) 19:37:20.72 ID:t88GB6TH - >>718
>699の答をサクッとかけばいいのに？
アホなの？子供に計算法を教えてもらったら？
子供も裏口容疑者なのか？
|
- 【ウハも】 開業医達の集い 28診 【粒も】
722 ：卵の名無しさん[sage]：2020/05/18(月) 19:39:10.22 ID:t88GB6TH - >>717
答が出せないのは恥ずかしいぞ。
こういうアホが職場で最初に発症した職員を感染源だと断定するんじゃないかなぁ？
|
- 底辺私立医大を卒業した医者って頭悪いよね？ Part16
465 ：卵の名無しさん[sage]：2020/05/18(月) 19:42:37.78 ID:t88GB6TH - 今日は内視鏡４件で終了。１時間程度の勤務で４万と楽だったが１０件はないといつものペースが乱れるなぁ
|
- 底辺私立医大を卒業した医者って頭悪いよね？ Part16
466 ：卵の名無しさん[sage]：2020/05/18(月) 19:46:12.14 ID:t88GB6TH - 裏口容疑者って統計処理の基本的な処理もできないアホなんだなぁ。
お金が好物の開業医スレに合わせた問題にしたんだが、誰も答が出せない。
必ずしも感染順に発症するわけではないというのを定量的に議論できる教養を欠いているね。
COVID19の潜伏期間の論文
https://www.nejm.org/doi/full/10.1056/NEJMoa2001316
結論は
#--- incubation period ---
# from Li et al NEJM 2020
# lognormal mean = 5.2
ln.par1 = 1.434065
ln.par2 = 0.6612
ある裏口容疑者が新型コロナ肺炎に罹患したとする。
行動調査によって発症前にキャバクラに行っており接客したキャバ嬢が裏口容疑者発症の２日後に発症していたことがわかった。
キャバ嬢は裏口容疑者から移されたと主張して１億円の賠償を求めている。
潜伏期間には幅がありキャバ嬢から移された可能性もあると主張してその確率を計算して賠償金を値切りたい。
いくら値切れるか計算せよ。
|
- 底辺私立医大を卒業した医者って頭悪いよね？ Part16
467 ：卵の名無しさん[sage]：2020/05/18(月) 20:47:32.80 ID:t88GB6TH - Temporal dynamics in viral shedding and transmissibility of COVID-19
https://www.nature.com/articles/s41591-020-0869-5
のRのコード
https://github.com/ehylau/COVID-19/blob/master/Fig1c_Rscript.R
と
西浦モデルのコード
https://nbviewer.jupyter.org/github/contactmodel/COVID19-Japan-Reff/blob/master/scripts/C.%20Calculating%20the%20Rt%20in%20Stan.ipynb
から発症間時間(serial interval)の分布を重ねてみた。
https://i.imgur.com/vrnra5F.png
> c(Mode,mean,var) # Xi He
[1] 3.630563 5.982241 20.728221
> c(Mode,mean,var) # Nishiura
[1] 4.259623 4.830129 4.940682
|
- 【ウハも】 開業医達の集い 28診 【粒も】
725 ：卵の名無しさん[sage]：2020/05/18(月) 20:49:28.43 ID:t88GB6TH - >>724
読んで自分で計算できる頭脳はないだろ。
|
- 底辺私立医大を卒業した医者って頭悪いよね？ Part16
472 ：卵の名無しさん[sage]：2020/05/18(月) 23:00:10.17 ID:t88GB6TH - >>471
学生時代には女子医大の学生には息子が大変お世話になりましたなぁ。
俺の後輩がコロナをもらわないようにと心温まる配慮だなｗ
|