>>529
Stanでの西浦モデルではinfecterとinfecteeが発症するまでの期間serial intervalの分布に
## Serial interval [Nishiura et al 2020 - only certain cases]
param1_SI = 2.305,
param2_SI = 5.452,
// serial interval
vector[K] gt = pweibull(param1_SI, param2_SI, K);
として使われているので、平均値などを出してみた。
乱数発生と理論値
> x=rweibull(1e5,param1_SI,param2_SI)
> mean(x) ; param2_SI*gamma(1+1/param1_SI)
[1] 4.829273
[1] 4.830129
> var(x) ; param2_SI^2*(gamma(1+2/param1_SI)-(gamma(1+1/param1_SI))^2)
[1] 4.907755
[1] 4.940682
> median(x) ; param2_SI*(log(2)^(1/param1_SI))
[1] 4.655777
[1] 4.6505
> density(x)$x[which.max(density(x)$y)] ; param2_SI*(1-1/param1_SI)^(1/param1_SI)
[1] 4.116837
[1] 4.259624
> optimise(function(x) dweibull(x,param1_SI,param2_SI),c(0,10),maximum = T)$max
[1] 4.259623
グラフにしてみた。
https://i.imgur.com/9vvCJuZ.png
正規分布で近似してもよさそうな感じだな。
数学 統計に詳しい人が語るコロナウイルス
■ このスレッドは過去ログ倉庫に格納されています
532132人目の素数さん
2020/05/16(土) 22:57:43.46ID:BxGcLzV+■ このスレッドは過去ログ倉庫に格納されています
ニュース
- ランドセルにくぎ刺される「国に帰れ」など言われ、転校を余儀なくされた海外からの転校生 仙台市教育委員会が「いじめ重大事態」認定★2 [煮卵★]
- 「いいの?前科ついちゃうよ」万引きした女子大学生を脅す 元コンビニ店長の男(54)逮捕 ★2 [煮卵★]
- 【テレ東音楽祭】長渕剛、テレ東初出演!『とんぼ』『乾杯』熱唱へ 酒井法子は『碧いうさぎ』歌唱 [Ailuropoda melanoleuca★]
- 「撮り鉄」同士のトラブルか 車のフロントガラス破壊の男逮捕 寝台特急カシオペア撮影のため東京から宮城・白石に遠征 [少考さん★]
- 【芸能】立川志らく「ヤクザの事務所が隣近所にあれば、強盗だとか闇バイトが来ない」 暴力団の“利点”指摘に違和感 [冬月記者★]
- 東京駅で切符紛失→「3倍払って」と言われ→拒否すると「警察呼ぶ」と言い始め警備5人が包囲… BD選手のトラブル報告にネット紛糾★3 [冬月記者★]
- 【地上波/DAZNほか】 FIFAワールドカップ2026 総合スレ★107【メキシコ/カナダ/アメリカ】
- 【MLB】ドジャース vs オリオールズ ★2
- 【地上波/DAZNほか】 FIFAワールドカップ2026 総合スレ★106【メキシコ/カナダ/アメリカ】
- 【MLB】ドジャース vs オリオールズ
- はません
- 巨専】
- 右翼「中国と戦争するぞおお!!」 普通の人「いややべえでしょ…貿易してるのに」 右翼「中国を気遣うのか!!」👈これ [268718286]
- 【悲報】著名精神科医・和田秀樹さん、高市早苗に「反社会性パーソナリティ障害」の診断を下してしまう [126042664]
- 木原「高市首相がG7でボッチだったり無視された事実はない!デマを流すなー😡」 [931948549]
- 【安倍悲報】株厨「資産数千万になったんだがwwwwwwwwwww」ぼく「あのさぁ…」 [279951338]
- 🏡ガイジ製造工場
- 【悲報】「フェルメール展」のチケットを取れなかったジャップ民がSNSで大暴れwwwwwwwwwwwwwwwwwwwwwwww [551743856]