[過去ログ] 数学 統計に詳しい人が語るコロナウイルス (1002レス)
上下前次1-新
抽出解除 レス栞
このスレッドは過去ログ倉庫に格納されています。
次スレ検索 歴削→次スレ 栞削→次スレ 過去ログメニュー
22(1): 2020/03/11(水)13:23:09.57 ID:hVKkfTiV(5/26) AAS
>>19
病人を100人集めたら10人が左利きであった。
左利き試験の感度は10%、右利き試験の感度は90%
ただ、これだけのこと。
88: 2020/03/22(日)08:39:00.57 ID:liILqu/N(2/6) AAS
>>85
プログラムを組んで計算させてみた(有病率の事前分布は一様分布を仮定)
感度0.7 特異度0.9として
100人に1人陽性だった場合の有病率の期待値、中央値、最頻値、95%信頼区間
> sim(100,1)
mean median mode lower upper
0.0196 0.0166 0.0101 0.0004 0.0464
省12
157: 2020/03/28(土)15:42:06.57 ID:BJlezchp(3/5) AAS
>>156
n(=10)人の中にi人の感染者がいるとき無作為にm(=2)人を選ぶ。
選ばれた2人の中に少なくとも一人の感染者がいる確率をP[x]として、
n個からr個選ぶ組み合わせの数をChoose(n,r)で表すと
P[xi]=1- choose(10-x,2)/choose(10,2)
xを0から10まで変化させて、
Σx*P[x]/(ΣP[x])で
省1
165: 2020/03/29(日)10:47:28.57 ID:WogCQeQk(3/9) AAS
>>163
そうみたいですね。
> data.frame(有名人=1:10,期待値=sapply(1:10,function(x) fn(100,x)$mean))
有名人 期待値
1 1 67.00000
2 2 62.75000
3 3 60.20000
省19
182(1): 2020/03/31(火)14:50:29.57 ID:ncBHjUEo(1/2) AAS
>>180
感染率の程度、感度・特異度の値の精度の言及無しに出された結論に、ほとんど説得力は無い。
230: 2020/04/09(木)11:57:57.57 ID:cpUiyQfM(1/2) AAS
これこそ、国を守るってことだな。
キューバ政府は2020年3月24日(火)以降、外国人観光客の入国を認めないと発表しました。
キューバに入国できるのは居住者のみで、これには出張で短期または長期滞在する外国人も含まれます。
加えて、キューバに入国する人は全員14日間隔離されます。
商業活動は維持されますが、商船や航空機の乗務員の移動は制限されます。キューバの領空は封鎖されません。
今回の制限措置は30日間有効で、状況に応じて延長される可能性もあります。
318(1): 2020/04/16(木)00:30:18.57 ID:oinCpTGH(2/5) AAS
>>304
> あんたは感染してもいないのに感染の危険にさらされてもいいのかね?
> 運が悪かったと諦めろと?
それは検査後の扱い、処遇の方法論の問題であって、
検査それ自体がもたらすリスクではない。
例えば、偽陰性の可能性があることを被験者に伝えれば、
検査期間が被験者に陰性のお墨付きを与えていることにもならない。
省4
346(1): 2020/04/17(金)02:11:00.57 ID:7ANIVQwh(1) AAS
>>345
穿った質問に素直に時間かける阿呆が居ると思ってる奴
399: 2020/04/23(木)17:36:48.57 ID:3vcirHk0(6/6) AAS
>>397
これであってるかな?
> # 期待値
> integrate(function(x) x*pdf(x),0,100)$value
[1] 37.13
> # 50人以上の確率
> integrate(pdf,50,100)$value
省4
415(2): 2020/04/24(金)22:01:25.57 ID:v55OWzbu(9/11) AAS
>>413
事前分布をJeffereyにしたら、
> js=PCRj4(67,4,SEN=0.5,SD1=0.2,SPC=0.9,SD2=0.1,N.ITER=1e6)$js
mean lower upper
1.4972e-01 6.3120e-14 8.5733e-01
> summary(prev)
Min. 1st Qu. Median Mean 3rd Qu. Max.
省3
493: 2020/05/11(月)11:02:15.57 ID:+2HA8Yn9(1) AAS
>>492
レスありがとうございます。
確かにSEIRモデルやSIRIモデルだと
「感染してから発表の数字に載るまでの日数の分布」って考慮されてないね。
497(1): 2020/05/14(木)08:00:23.57 ID:kec+XbRE(1/9) AAS
>>494
これにある、
JapaneseDataCOVID19
ってどこかからダウンロードできるんだろうか?
570: 2020/05/21(木)11:39:47.57 ID:hAZkHjNF(1/21) AAS
西浦モデルの再生算数の事前分布を変化させてグラフにしてみた。
西浦モデルでのデフォルト
画像リンク[png]:i.imgur.com
事前分布を一様分布(0,10)近似の正規分布で近似させた場合
画像リンク[png]:i.imgur.com
再生算数の平均0、標準偏差1の場合
画像リンク[png]:i.imgur.com
省2
583: 2020/05/21(木)17:26:10.57 ID:BjgAHZWs(1) AAS
まぁこのスレは用語がめちゃくちゃだからなぁ。
602(1): 2020/05/22(金)08:35:38.57 ID:DQesskhT(1/9) AAS
>>599
あれ3/3になったのか。
それだと、95%CrIは0.5559329 1
事前分布にはJeffereysで計算
661: 2020/05/30(土)10:10:43.57 ID:oTxR9Gdg(1) AAS
>>659(自己解決)
mu = 1.434065
sg = 0.6612
pdf1 <- function(x) dlnorm(x,mu,sg)
pdf2 <- function(y) dlnorm(y,mu,sg)
f <- function(x,y) pdf1(x+y)*pdf2(y)
vf=Vectorize(f,vectorize.args = 'y')
省3
692(2): 2020/06/16(火)10:25:31.57 ID:IaE5USJ0(1) AAS
厚労省がやってた抗体検査も結果出た。
外部リンク[pdf]:www.mhlw.go.jp
925(1): 2020/07/28(火)09:47:33.57 ID:/xPjTU2o(1/4) AAS
>>921
自然免疫で感染が防げて、そういう人が多いのなら、たとえば
ミニシアターなんかの集団感染での感染率からある程度わかり
そうな気もする。報道による情報が少なすぎてよくわからんけど、
濃厚接触した人の何割が感染したんだろうねぇ?
上下前次1-新書関写板覧索設栞歴
スレ情報 赤レス抽出 画像レス抽出 歴の未読スレ AAサムネイル
ぬこの手 ぬこTOP 0.032s