数学 統計に詳しい人が語るコロナウイルス

1 名前：１３２人目の素数さん [2020/02/29(土) 02:18:41 ID:twdO677Q.net] 東大数学科卒の元官僚はこう分析してるが、お前らはどうなると思う？ www.zakzak.co.jp/soc/news/200220/dom2002200003-n2.html 中国国外感染者の中国国内との比率をみると、 １月２０日の数字公表以降は、０・８〜２・６％で比較的安定している。 これは、新型肺炎の感染者のほとんどは中国国内、それも湖北省に集中しているからだ。 ちなみに中国国外での感染者数は、中国国内の１・１％だ（２月１６日現在）。 本コラムで紹介したが、現時点では、最終的な中国国内の感染者数は２０万人超と筆者は推計している。 となると、中国国外の感染者は数千人程度になるだろう。 中国国外のうち日本の比率は１割弱なので、日本の感染者数は数百人程度であろう。 その場合、死者も数人から１０人程度になるだろう。 こうした推計をすると、今の感染者は氷山の一角だと思われるが、今後の増加ペースはどうなるだろうか。 新型コロナウイルスの検査は簡単に行えるので、今後、日本での感染者数は増えていくだろう。 ある時点ではそれがネズミ算的に増えるかのように思える局面もあるだろうが、 筆者の推計が正しければ、現時点ではせいぜい数百人が一つのメドだ。 772 名前：１３２人目の素数さん [2020/07/03(金) 01:23:32.12 ID:DYGsod0z.net] 新規感染者が増えかかっているがどう見る？ 773 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 02:35:31.91 ID:90y63y3Z.net] >>727 西浦氏が再評価される。 774 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 10:13:27.22 ID:hNK4yhyJ.net] https://egg.5ch.net/test/read.cgi/sisou/1593487605/133-135 すまん寝＠秩序回復・財務省廃止・反財政再建 @sumannne 【新型コロナウィルス】藤井さん、そのグラフ変じゃないですか？ https://himorjp.blog.fc2.com/blog-entry-179.html こっちで改めて計算したところ、「自粛要請は滅茶苦茶効果あった」という結論になっちゃったんで、藤井さん、どうにかした方が良いと思うよ。多分専門家はもうツッコんですらくれないだろうし。 ちなみにデータソースは感染研に1か月遅れですがありましたので、怪しげなデータを使っているわけではありません。 むしろ藤井さんが使ってるデータソースが不明。 藤井さんの統計のいじり方で一番分からんのが、対数グラフを取るでもなく、普通に差分するでもなく、「対前日比」という謎の計算で「速度」「加速度」を名乗ることなんだよなぁ。 この程度のオーダーで対数グラフ取る意味が分からんし、実際藤井さんは対数取ってるわけじゃないので、そうであれば変化率も加速度も差分しないと正しく見えないと思うが。 対数グラフを取る、もしくは対前日比でモノを言うのは、日ごとに指数増大していることがほぼ明らかな状況でなら分かるが、少なくともデータ上はそうではないから適切な計算方法とはちょっと思えない。 藤井さんの記事の「加速度」、縦軸を見てもらえばわかるが上下が0.01の幅で動いてる。 計算は「対前日比」でやってるのでほとんど0ですよこれは。 数値計算で、厳密解が0のものが数値誤差でちょっと荒れてるのを拡大顕微鏡で見て騒いでるような感じがする。 (実際、指数増大しているわけでないものを「比」で取ったら本質的に比率1しか出ません。) 極めて不適切だと思う。 まぁ、今回はさすがに時系列データに相関分析かけたりはしなかったから、その点は成長してるんじゃないっすか。 「対前日比」がいかにおかしいかの例題。 2次多項式のオーダーで増大する場合でも、「対前日比」を取っちゃうと増減が全然分からない。 https://pbs.twimg.com/media/Eb7GF0JUMAEnVH-.png https://pbs.twimg.com/media/Eb7GHp4U4AEvAyt.png 775 名前：１３２人目の素数さん [2020/07/03(金) 10:13:49.21 ID:GtC8ue70.net] 西浦モデルはもう使われないんじゃないの？ とりあえず夜の飲食業は営業停止で。 776 名前：１３２人目の素数さん [2020/07/03(金) 10:16:46 ID:GtC8ue70.net] >>729 藤井聡太が迷惑しそうな名前のオヤジだなｗ もう藤井恥 777 名前：に改名したほうがよろしい。 [] [ここ壊れてます] 778 名前：１３２人目の素数さん [2020/07/03(金) 21:04:53.02 ID:I3C9PAZr.net] 数学系でかつ鉄道趣味のある人で通勤電車の中にどのくらい元気な感染者が存在しうるか なんて計算してる人っているんだろうか？　（元気な感染者＝感染力の強い時期の人） 　あるいはそんな研究あるんだろうか？　　　　　　-素朴な疑問です　他意はありません 779 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 21:52:40 ID:UvmxXx4m.net] https://forbesjapan.com/articles/detail/32245 上の記事、新型インフルエンザで申し訳ないが、感染研が2011年に電車に1人乗ってただけでも感染爆発するというシミュレーション出してる。 780 名前：１３２人目の素数さん [2020/07/03(金) 22:04:30.84 ID:GtC8ue70.net] コロナの飛沫感染に関しては、マスクをしてるだけで感染力は 大幅に下がるんでねぇの？ インフルエンザとの違いはそこだね。 781 名前：１３２人目の素数さん [2020/07/03(金) 22:05:15.79 ID:GtC8ue70.net] したがって、電車に乗る人はマスク着用を義務付けるべきなんだよな。 782 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 22:45:08.73 ID:90y63y3Z.net] >>733 今は削除されているのでWebarchiveにある。 曰く 首都圏の鉄道に新型インフルエンザを発症した人が1人乗ったと仮定して、まったく対策を取らなかった場合、重症軽症を含めた感染者数は1週間で12万人に拡大するというシミュレーションがあります。 https://web.archive.org/web/20090915203236/https://www.terumo-taion.jp/fever/pandemic/06.html 783 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 22:47:39.63 ID:90y63y3Z.net] コロナが始まって日本人はマスクを買い増したが、アメリカ人が買い増したのはライフル銃の弾。 マスクに関してはこれが面白かった。 A modelling framework to assess the likely effectiveness of facemasks in combination with ‘lock-down’ in managing the COVID-19 pandemic https://royalsocietypublishing.org/doi/10.1098/rspa.2020.0376 Identifying airborne transmission as the dominant route for the spread of COVID-19 https://www.pnas.org/content/early/2020/06/10/2009637117 https://www.pnas.org/content/pnas/early/2020/06/10/2009637117/F3.medium.gif 飛沫の可視化 https://youtu.be/UNHgQq0BGLI 784 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 22:56:23.10 ID:90y63y3Z.net] >>729 ここのK値というのも、意味が理解できなかった。 https://asahi.5ch.net/test/read.cgi/newsplus/1593743783/ 「（感染増の）波のとらえ方に、K値を採用しようと思っています。中野先生のおっしゃるK値モデルは、 感染の収束速度を計算して感染者数を想定するうえでは、かなり正確だと思う」 K値について詳しい説明はこちらで確認 https://www.dailyshincho.jp/article/2020/07020603/?all=1 https://www.dailyshincho.jp/article/2020/07020603/?all=1&page=2 785 名前：１３２人目の素数さん mailto:sage [2020/07/03(金) 23:50:48 ID:UvmxXx4m.net] >>738 K値ってのは簡単に言うと、適当に数値の比をとったもの。疫学的な意味はない。 786 名前：１３２人目の素数さん mailto:sage [2020/07/04(土) 07:13:01 ID:kFD5yPDP.net] >>739 直線でに線形回帰みたいなもの？ 787 名前：１３２人目の素数さん mailto:sage [2020/07/04(土) 07:55:47.68 ID:2sU7d1SU.net] >>740 適当に線引いたらそれっぽい数字が出たという数値 788 名前：１３２人目の素数さん mailto:sage [2020/07/04(土) 08:21:47.85 ID:kFD5yPDP.net] こういうパラメータも既知の分布によく当て嵌まるかで決定なのだろうな　(>681参照) Distributional fits to key COVID-19 distributions. https://i.imgur.com/SJ6QS3b.png https://www.medrxiv.org/content/10.1101/2020.03.03.20028423v3.full.pdf 789 名前：１３２人目の素数さん mailto:sage [2020/07/04(土) 13:22:16.78 ID:Y1vg+l8y.net] K値ってのは１週間の感染者数を総感染者数で割っただけなんだが、そうすると既に大流行が終わったところに2波目が来ても不当に小さい値になってしまう。 そこで適当なところで「リセット」出来ることになっているのだが、それが恣意的な時刻で可能なのが最大の問題。 発明者は「ほら、2波目を予測できてるでしょ」と言うのだが、いやいやそれ都合の良いとこでリセットかけてるだけやんとしか。 790 名前：１３２人目の素数さん mailto:sage [2020/07/04(土) 20:49:26.54 ID:kFD5yPDP.net] >>716 https://www.youtube.com/watch?v=xWT3g-oGFkY 15分辺りで２キットの検査の違いを解説している。 791 名前：１３２人目の素数さん [2020/07/04(土) 21:33:00.37 ID:YBorZJOe.net] >>733 ありがとうございます コロナの現状（マスクあり）でも誰か やっててもよさそうなもんですね 不安を煽るとかで規制あるのかな 792 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 08:07:43.49 ID:4xu+FlwC.net] >>743 解説ありがとうございます。 都合よくあう点を選んで接線を引いて予測どおりというみたいなものという理解でいいですか？ 793 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 08:09:26.87 ID:4xu+FlwC.net] >>745 最悪の事態がカイヒできたら叩かれるから自主規制しているのでは。 794 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 11:21:40.50 ID:rpUuAKzr.net] 10/14って何人から感染したんだろうなぁ？ どうやれば計算できるだろ？ 再生産数から計算できる？ >> 京都市は先月市内の飲食店で開かれた小規模のパーティーで新型コロナウイルスの集団感染＝クラスターが発生したと発表しました。これまでに男女12人の感染が確認され、市が濃厚接触者などの調査を進めています。 クラスターが発生したのは先月27日に京都市内の飲食店で開かれたパーティで参加した14人のうち今月2日から4日にかけて男女10人の感染が確認されたということです。 また、二次会が開かれたバーの従業員2人の感染も確認され、関係する感染者は合わせて12人に上っています。 京都市は、目安となる「同一の場所で5人以上の感染」が確認されたとして、二次会での感染も含めて新たなクラスターが発生したとしています。 市は濃厚接触者などの調査を進めています。 << https://www3.nhk.or.jp/news/html/20200705/k10012497081000.html 795 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 11:25:37.47 ID:by0UA+MF.net] >>746 何も分かってない気がする。 >>748 再生産数は平均値だから計算できないよ。（再生産数の推計に対して、実際に感染させた人数がバラついていることからｸﾗｽﾀｰ対策が本丸であることを西浦は見抜いた） 796 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 13:11:55.13 ID:rpUuAKzr.net] >672のデータ https://www.researchsquare.com/article/rs-29548/v1 From the empirical offspring distribution and fitted negative binomial distribution shown in Figure 2B, we estimated an observed reproductive number (R) of 0.58 (95% CI: 0.45 - 0.71) and dispersion parameter (k) of 0.45 (95% CI: 0.31 - 0.76). を使って乱数発生させてシミュレーションしてみた。 mu=0.58 size=0.45 (prob = size/(size+mu)) Rt=rnbinom(1e5,size=size,mu=mu) # reproductive number hist(Rt,breaks = 'scott',freq=F,ann=F) sim <- function(n=10){ infector=sample(1:n,1) # how many infectors? infectee=sum(sample(Rt,infector)) # how many infectees? if(infectee==n) return(infector) # is the number of infectees n? else return(NA) } s=replicate(1e6,sim()) spreader=s[!is.na(s)] hist(spreader,freq=F,ylab='',axes=F,breaks='scott') ; axis(1) BEST::plotPost(spreader,xlim=c(0,1 797 名前：0)) mean(spreader) HDInterval::hdi(spreader)[1:2] sum(spreader==1)/length(spreader) 期待値は > mean(spreader) [1] 8.032594 95％CIは > HDInterval::hdi(spreader)[1:2] lower upper 4 10 1人のスーパースプレッダーから広がった確率は > sum(spreader==1)/length(spreader) [1] 0.001580299 [] [ここ壊れてます] 798 名前：１３２人目の素数さん mailto:sage [2020/07/05(日) 20:02:50.90 ID:rpUuAKzr.net] >>750 バグに気づいたので撤回。 infectee==10ではなくてinfector+infectee==10だな。 デバッグ版 R=0.58 # mean of reproductive number k=0.45 # dispersion parameter (prob = k/(k+R)) # its probability Rt=rnbinom(1e5,k,mu=R) # random numbers of negative binomial distribution hist(Rt,breaks = 'scott',freq=F,ann=F) # show its histgram sim <- function(n=10){ # simulation infected=0 # initial value while(infected!=n){ # while infected is unequal to n infector=sample(n,1) # prior discrete uniform distribution of infector number infectee=sum(sample(Rt,infector)) # number of infectee infected=infectee+infector # number of infected } return(infector) # when n infected, return infector number } spreader=replicate(1e5,sim()) # simulation & calculation hist(spreader,freq=F,ylab='',axes=F,breaks='scott') ; axis(1) HDInterval::hdi(spreader)[1:2] # 95% credibility interval BEST::plotPost(spreader) # graph with 95%CI & mean summary(spreader) sum(spreader==1)/length(spreader) # the probability of single super-spreader まあ、パーティー中に感染させる人数に再生産数を流用していいかは疑問ではあるが、 実効結果は https://i.imgur.com/ApDSW1J.png 95% CI > HDInterval::hdi(spreader)[1:2] # 95% credibility interval lower upper 3 9 中央値、平均値、四分位値 > summary(spreader) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.000 6.000 7.000 6.844 8.000 10.000 1人のスーパースプレッダーの確率 > sum(spreader==1)/length(spreader) # the probability of single super-spreader [1] 0.00104 799 名前：１３２人目の素数さん mailto:sage [2020/07/08(水) 05:37:59 ID:uI62U+we.net] >>751 Rt値はそのままで良いのだろうか。まあ、東京で計算したらそれぐらいだから、今の環境的に言ってそれで良いか。 というか、Rtが変化しても、dispersion parameterは0.1のままで良いのかな？ 800 名前：１３２人目の素数さん mailto:sage [2020/07/08(水) 11:05:25 ID:XD7Ql8W/.net] >>752 >751の論文に95%CIが書いてあるのでそれを使って階層モデルを作ろうかと思って原著を読んだら95%CIはブートストラップ法で算出したと書いてあった。 Rtの方は期待値がCI境界の平均なので正規分布を期待したのだが。 dispersion parameterは非対称の分布みたい。 801 名前：１３２人目の素数さん mailto:sage [2020/07/08(水) 20:15:37 ID:I3BoIViR.net] >>752 レスありがとうございます。 パラメータの分布を勝手に設定して階層モデルもどきでシミレーションしてみた。数字のお遊びでしょうけど。 From the empirical offspring distribution and fitted negative binomial distribution shown in Figure 2B, we estimated an observed reproductive number (R) of 0.58 (95% CI: 0.45 - 0.71) and dispersion parameter (k) of 0.45 (95% CI: 0.31 - 0.76). Rは正規分布、kはガンマ分布として、期待値と95%CIが上記に適合するようなパラメータを ニュートン・ラフソン法で算出すると 正規分布の標準偏差sd＝0.0663277494 > pnorm(0.71,0.58,0.0663277494)-pnorm(0.45,0.58,0.0663277494) [1] 0.95 ガンマ分布の形状パラメータsh＝24.38466,尺度パラメータ=0.01845 > pgamma(0.76,sh,scale=sc)-pgamma(0.31,sh,scale=sc) [1] 0.95 が得られた。 各パラメータが上記の分布に従うとして R ~ normal(mean=0.58,sd=sd) k ~ gamma(shape=sh,scale=sc) Rt ~ nbinom(dispersion=k,mean=R) のモデルでシミレーション 既述のコードに RR=rnorm(1e5,R,sd) # R ~ dnorm(R,sd) kk=rgamma(1e5,sh,scale=sc) # k ~ dgamma(shape=sh,scale=sc) RRt <- function(x){ # Rt ~ dnbinom(k,mu=R) k=sample(kk,1) R=sample(RR,1) rnbinom(x,k,mu=R) # pick x random Rt numbers } を追加して while loop内のinfecteeの数を infectee=sum(RRt(infector)) # number of infectee に変更。 最初の10人に何人が感染していたかの95％推定 > HDInterval::hdi(spreader)[1:2] # 95% credibility interval lower upper 3 9 > BEST::plotPost(spreader,xlim=c(1 802 名前：,10)) # graph with 95%CI & mean https://i.imgur.com/LKu4hKD.png > summary(spreader) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.000000 6.000000 7.000000 6.894006 8.000000 10.000000 > sum(spreader==1)/length(spreader) # the probability of single super-spreader [1] 0.001023 予想通り、ほとんど、結果は変わらず。乖離した方がおかしい。 [] [ここ壊れてます] 803 名前：１３２人目の素数さん mailto:sage [2020/07/09(木) 22:44:44 ID:dYeNIQef.net] >>754 この結果からは、 再生産数はもっと高い and/or 無症状の市中感染が蔓延している と考えられる 804 名前：１３２人目の素数さん [2020/07/10(金) 02:12:10 ID:u3owlDvV.net] 東京で224人　日本全国でも350人の新規感染者が出てますが このスレ的にはどうすればいい思いますか。 要するに ・このまま放置 ・夜の街になんらかの規制かける ・緊急事態宣言 とかいろいろ手段はあると思いますが 805 名前：１３２人目の素数さん [2020/07/10(金) 03:23:48.78 ID:hTe5Oog6.net] K値予測って、実効再生産数が1未満の地域で、突発的にたくさんの感染者が生まれたときの感染者数推移だろ 東京は実効再生産数が1を超えるので、感染者が右肩上がり 東京以外も、東京から感染者が継続的に供給され続けるので、K値予測は当たらない 806 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 06:03:42 ID:uMeV9IwL.net] 全く納得できませんというpeerのコメント K値の胡散臭さだけは伝わってきた。 https://pubpeer.com/publications/2D5A3D7BC3C6E9EAF6F7FEF9563C8B 807 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 08:20:42.04 ID:++P0OH09.net] >>756 数学民としては放置して、推移がちゃんとSIR方程式に従うのかを確認したい。 808 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 10:04:28 ID:uMeV9IwL.net] >>759 ワクチンや治療薬が開発されたときのパラメータを追加すればいいんじゃないかなぁ？ SEIR MODEL dS(t)/dt = mu*(N-S) - b*S(t)*I(t)/N - nu*S(t) dE(t)/dt = b*S(t)I(t)/N - (mu+sig)*E(t) dI(t)/dt = sig*E(t) - (mu+g)*I(t) dR(t)/dt = g*I(t) - mu*R + nu*S(t) mu:自然死亡率 b:感染率(S->I) nu:ワクチン有効率(S->R) sig:発症率(E->I),g:回復率(I->R) 今の時期に国外から持ち込まれるCOVID-19はパラメータが異なるように思える。 809 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 17:37:19.00 ID:opwaioSQ.net] >>756 素晴らしい。なんと素晴らしい馬○発見機会を得られたのだろうと喜べばよい。 強く疑われる場合にのみ検査して出された感染発覚者人数と 見舞金目的（？）と思われるような集団にさえ検査を行い、カウントされた感染発覚者人数を 同レベルで比較するコメンテーターや、無調整で統計処理を施そうとする人を見つけたら、ラッキーと思えばよい。 彼らの発言は信用できないという強い証拠を得られたのだから。 810 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 17:46:44.67 ID:9VixCbDU.net] >>756 30代以下の感染者が80%以上を占めているわけだから慌てる必要は無い 現状の対策を継続維持だ 811 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 18:08:08.15 ID:uMeV9IwL.net] >>762 SARS罹患した若者の長期予後は芳しくないと耳にしたことがある。 SARS-COV-2の場合の長期予後は未知だから楽観視できないと思う。 812 名前：１３２人目の素数さん [2020/07/10(金) 20:04:48.04 ID:PNuztIXJ.net] SIRモデルにしろ、はしかモデルにしろ 感染症数理モデルはウィルス拡散防御のためにある公衆衛生.医学の道具で、 方程式だけをいじっても素人は素人。 最� 813 名前：ﾟ　医療現場無縁の関西の似非専門家 ど素人が声高く徘徊してますね。 ・都のデータの注意書きも読めない 　IPS細胞のY京大教授 ・医療システム無視でPCR検査信奉 　のH京大名誉教授 ・医療公衆衛生無知でウィルス実験と 　混同するM京大准教授 ・核物理が専門なのにグラフを 　いじってK値をとなえるN阪大教授 　不思議にお維新に近い考え方です。 西浦北大教授がまともでしょうし、 生物の数理モデルがベースでしょ。 故山口京大教授の数理モデルも考え方も 関西から絶滅したのだろう。 [] [ここ壊れてます] 814 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 20:22:08 ID:uMeV9IwL.net] 西浦氏が山中氏と対談 “市中感染が広がり始めている認識” https://www3.nhk.or.jp/news/html/20200710/k10012508491000.html 西浦先生は正しかったからな 815 名前：１３２人目の素数さん [2020/07/10(金) 20:39:56.75 ID:hTe5Oog6.net] K値予測は、実効再生産数が0.7前後くらいの地域で、一時的要因により感染者が急増して域内での感染が蔓延、 その後収束していく様子を描いてるのでは？ 日本でいえば、東京以外の大都市部 鹿児島は、基本的な実効再生産数がかなり低く、クラスター対策もしっかりできたので、 大量に出た感染者は、K値予測よりかなりはやくストンと直滑降で収束 東京は、実効再生産数が継続的に1を超えてるので、K値予測が適応不可 東京以外の大都市部は、東京から継続的に感染者が供給されてるので、 「一時的要因」ではないのでK値予測が当てはまらない 816 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 20:42:04.82 ID:hTe5Oog6.net] 東京はすでに、感染者の急増と、市中感染の増加により、クラスター対策が不可能になったので、 大規模規制でおさえるしかない状況 817 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:08:25.69 ID:9VixCbDU.net] 3月4月はまだ正体がよく分かってなかったこともあり非常事態宣言を出したけど、今は分かってるからね 撲滅は目指さない。共存戦略でいく 風邪やインフルエンザウイルスと同じようにね。既存の4種のコロナウイルスとも同様だ 急速な終息を目指すのは賢い戦略では無いからね 医療体制の崩壊を防ぐこと 院内感染を防止すること 老人施設などクラスター防止 オフィスなどのクラスター防止 60歳以上の高齢者の重症化対策に重点を置く これでいいんだよ 重症化対策がこの病気の肝 若者層は無症状がほとんどだし感染クラスターを把握していく今の体制で十分 たとえ1000人規模になっても30歳以下が80%以上で上記施設でクラスターが抑え込めていれば現状の対処でよい 818 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:09:32.25 ID:uMeV9IwL.net] >>766 K値って原子核の崩壊定数からの類推なのかな？ 819 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:14:22.02 ID:NckxUv6S.net] >>768 そのいいとこ取りの解ってどんな解？ 820 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:32:40.79 ID:9VixCbDU.net] 逆に、 院内感染急増 老人施設クラスター急増 60歳以上の感染者急増 入院患者急増 重症者急増 といった事態が生じたら速やかに非常事態宣言だな 要は感染者数だけを強調するのはナンセンスだ メディアで報道しないといけないのは 入院者数 重症者数 死亡者数 821 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:35:03.25 ID:NckxUv6S.net] >>771 イヤ、数学板なんだから数理防疫学の話してよ。 どんな解の事言ってんの？ 822 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:46:21.32 ID:9VixCbDU.net] 現状の 「感染者出ました！」 「感染者増えました！」 と大騒ぎするのは知能が足りてないよ メディアがレポすべきなのは 医療現場がいかに対処したか？ 医療現場がいかに重症患者を救ったか？ 治療レシピはなにが最適か？ 医療体制はどうあるべきか？ 感染者数で一喜一憂するなよ。バカみたい 823 名前：１３２人目の素数さん mailto:sage [2020/07/10(金) 21:50:12.89 ID:osz7X9Cr.net] >>773 他所でやって 824 名前：１３２人目の素数さん [2020/07/11(土) 01:47:29 ID:0/AYQWhD.net] >>761 人数はともかく、増大のトレンドにあることは確かだけどね。 東京以外でも、症状が出てる感染者が確実に増えてる。 825 名前：１３２人目の素数さん mailto:sage [2020/07/11(土) 09:11:06.67 ID:z+lU7R6q.net] >>734 ウイルスから守ると言う意味では、特に布マスクや、不織布マスクでも隙間があると無力らしい https://www.asahi.com/articles/ASN7430PVN72UBQU00B.html 各社が夏にぴったりとか言って、隙間だらけの布マスクを販売してるけど、 電車の中で、マスク無しの感染者が咳・くしゃみをすると相当ヤバそう。。 826 名前：１３２人目の素数さん mailto:sage [2020/07/11(土) 09:25:46 ID:TZvEYwn/.net] >>776 私のマスクがあなたを守り、あなたのマスクが私を守る。 NEJMへ投稿された動画 https://i.ytimg.com/vi/UNHgQq0BGLI/default.jpg 827 名前：１３２人目の素数さん [2020/07/11(土) 12:05:53.64 ID:0/AYQWhD.net] >>776 通常のマスクに防御の役割が低いのはとっくの昔にわかってたこと。 その記事はなにを今更って感じで論外。大西君も今頃そんなこと やってるなんて、、、。恥を知ってほしい。 >>734は感染させないためのマスクという役割を前提として指摘ですよ。 828 名前：１３２人目の素数さん mailto:sage [2020/07/11(土) 16:32:55.60 ID:SN7B5Z3H.net] コロナに感染させたマウスが入ったケージと、非感染のマウスが入ったケージを用意し、 送風機で、感染側から非感染側に風を送ったところ、非感染マウスの一部が感染した。 非感染側のケージに布マスクに相当するような障壁を置いたところ、感染割合が約2/3に、 感染側のケージに布マスクに相当するような障壁を置いたところ、感染割合が約1/3にまで減った。 という論文が出された、というテレビ放送を見たことがある。つまり、非感染者が布マスクを つけること（＝自衛目的）により、感染を1/3防ぐことができ、感染者が布マスクをつけること （＝感染拡大抑止目的）により、感染を2/3防ぐことができるということだ。 布マスクの隙間サイズが、コロナウィルスのサイズの数百倍〜数千倍あるから、 布マスクの効果は限定的だという話は以前から指摘されているが、それを覆す報告だ。 ウィルスが単独で飛翔することは希で、大概飛沫に乗って移動する。 従って飛沫の通過／進入を防ぐことができるなら、その分の効果は期待できると言うことなのだろう。 もしかすると、布マスクの最大の効用は、口や鼻を直接触る機会を減らすことにあるかもしれない。 だとすると、上のようなマウス実験例を示すまでも無く、「布マスクでも効果はある」と断言できる。 花粉症が蔓延している日本では、春先にマスクを着ける姿は違和感無く受け入れられている。 ファクターＸは、多くの要因の集合体と思われるが、一要因に「マスク」も含まれるだろう。 829 名前：１３２人目の素数さん [2020/07/12(日) 06:05:28 ID:JrtRNErT.net] >>769 >K値 物理なら標準的な基本理論に 基づき実験し測定値を実験式で 表して誤差を考察するけどK値は物理と 無関係でしょ。 感染症素人が個人的にグラフや資料を 見てピピっとお告げがあったんだろう 競馬予想や株価預言でも見られる 現象かな、まあ阪大の保証もないし 雨合羽にK値ではネタしかならないな 830 名前：１３２人目の素数さん [2020/07/12(日) 10:18:36.16 ID:JrtRNErT.net] >>780 続き >SIR微分方程式と原子核崩壊 ・SIRモデル。 　感染症数理モデルは� 831 名前：Vンプルで 　集団モデルのSIRが基本モデルだ。 　微分方程式は(ミクロな)個人間の接触　 　が均等より非線型項βSIがあり 　非線型。 ・原子核崩壊。 　放射平衡は 微分方程式d N1/dt=-λ1N1 　が基本で相対論的場の量子論. 　断面積から導かれる。 　親→娘→孫の線型方程式。 ・SIRモデルの局所線型化。 　SIRモデルで初期t≒0の時 　S(t)≒S(0)=1 で線型方程式になる 　SIRモデル(t≒0)と放射平衡は 　共に線型微分方程式として 　対応できる ・マスク効果は個人間接触の相互作用. 　断面積パラメータで扱え 　集団のSIRモデルで感染率パラメータ 　となる。 　ミクロ相互作用の接触は密度に依存 　する 以上は数理モデルの概論です 　 [] [ここ壊れてます] 832 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 10:41:01.76 ID:9O0Ntz/e.net] >>780 レスありがとうございます。 K値とは合羽値の頭文字出会ったとｗ >694のκ値とは別物。 833 名前：１３２人目の素数さん [2020/07/12(日) 10:57:43 ID:JrtRNErT.net] >>782 K値は根拠不明 834 名前：１３２人目の素数さん [2020/07/12(日) 14:27:41.83 ID:yl7MdC0U.net] 工学の人が好きそうな、いきあたりばったり的なパラメータだよね＞K値 理学系なら演繹的に捻り出すモデルを考えろよ、と言いたい。 835 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 14:54:44 ID:qAF2CiMg.net] モデルを考えている間に感染が広がったり経済的損失が拡大してしまったら役に立たない 病気の原因はわかっても患者が死んだら意味がない 学問的には意味があるだろうけど実践的には意味がない モデルを考えるのは後からでもいい 836 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 15:17:13 ID:Wf1trDX+.net] 日本に必要なのは、都道府県別実効再生産数を1未満に保つ、 仮に1を上回ったら2週間以内に1以下に下げる っていう数値目標を加えること 実効再生産数が1を超えると2週間後に感染者が目に見えて増大、4週間後には急拡大する 感染者数とかだけで見てると先手先手をうった対応が出来ず後手後手にまわるのは いまの東京を見れば明白 東京は5月後半に実効再生産数が継続的に1を超えて、 感染拡大するのは明白だったのにろくな対策しなかった 837 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 15:19:19 ID:Wf1trDX+.net] 実効再生産数の目標には外部要因も考慮する必要がある たとえば、自分ところの都道府県は実質的には実効再生産数が1以下で、 自分ところの都道府県内での感染は抑制されてるが、 他県、たとえば東京からどんどん感染者が流入してくるので数値的には1を超えた、みたいな場合ね 838 名前：１３２人目の素数さん [2020/07/12(日) 17:30:49.38 ID:PYIegNFg.net] >SIRモデルと再生産数 潜伏期間が0や住民間で均等な接触 を仮定し集団を対象にしたシンプルな 非線型微分方程式。 生物の基本的な数理モデル ・防疫戦略の基本的な道具の一つ ・インドボンベイのペスト流行でも 　流行曲線を再現した ・体内のウィルス感染モデルにも 　使用されている。 ・初期に非感染者S≒1の場合 　d I/dt=(β-γ) I =(Ro-1)γI 　基本再生産数Ro≡β/γ 　Ro>1の時　感染者 Iが増加 ・実効再生産数R 　感染原因が住民間の接触だから　 　接触を減らせばRが減少する。 　8割接触が減少すれば 　R=Ro×(1-0.8)=0.2R だからマスクをして三密はさける。 (物理公式と同じく意味もわからず 公式へ数値を入れ結果を振り回さない、 Rは指標の一つでデータの取り方に 依存する>マスコミ素人学者政治家) ・都道府県単位で平均された 　データ指標には注意が必要。 　ウィルス拡散は放射能拡散と 同様に県境がないし業種交通網 　年齢他にも依存する ・検査対象データは検査数対象 　特性に依存するし東京以外は 　戦略的検査数が少なく過ぎると 　思われる 　抗体検査では大阪>東京>宮城 　だった ・SIRモデルから潜伏期間EのSEIR モデル、進行波解がある(SIR＋拡散 項)モデル、再帰性があるインフル 　エンザモデルなどがあり個人行動 　をシュミレーションしたibmモデル 　も以前からある。 　新型インフルエンザ対策として 　政府へ提言されていた 　マスク備蓄もね 839 名前：１３２人目の素数さん [2020/07/12(日) 17:31:39.31 ID:PYIegNFg.net] >>787 検査していますか？ 840 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 17:43:16.95 ID:PLqlj++l.net] >>788 海外から流入するパラメータの違うIを加えたモデルってあるかな。 841 名前：１３２人目の素数さん [2020/07/12(日) 18:02:36 ID:PYIegNFg.net] >>790 移流項＋SIRモデルやibmモデルかな。 SIRモデルは人口一定の閉鎖地域だから 感染者Iの流入流出項を加える。 (物理化学モデルの類似なら 閉じ込められた分子が衝突して化学 反応をしている。 更に外部から分子が流入するとか) 西浦教授が海外からの水際対策を説明 するため玉を流してぶつけていたが 西浦教授に聞いてください。 842 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 18:06:06 ID:4t109+lN.net] >>787 その数が分かっているなら、Rtの分子から流入分は除くべき。 西浦先生の計算では、輸入分の患者はそう取り扱っていたと思う。 843 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 18:18:56.02 ID:PLqlj++l.net] >>791 レスありがとうございます。 SEIRに人口増加の項を加えたモデルはみたことがあるけど IBMって今はレノボではと思ったら別のIBMだったｗ こういうのがみつかりました。 https://bradduthie.github.io/blog/individual-based-models-in-r/ 理解できるかわからないけど読んでみます。 844 名前：１３２人目の素数さん mailto:sage [2020/07/12(日) 21:07:07.01 ID:PLqlj++l.net] RにIBMのパッケージがあった。 https://cran.r-project.org/web/packages/ibm/index.html 845 名前：１３２人目の素数さん [2020/07/13(月) 09:41:22 ID:dpBf9py+.net] >>785 だから、実践なんか二の次で真理を求めるのが理学。 846 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:17:11.63 ID:ux2hxf89.net] K値による予測を使うべきでない理由 katukawa.com/?p=6228 847 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:21:00.33 ID:FAV+lxAc.net] >>795 真理って言っても確率的に大雑把にしか言えないと思う 個々の人から排出される飛沫にウイルスがどの程度含まれて それかが別の人の体内にどの程度入って その体内でウイルスがどう動いて免疫がどう働いて 結果的に症状や他の人に感染させるかとか 具体的な事は複雑で判らないと思う 真理にはたどり着けない 848 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 19:33:06.81 ID:tR49RdlH.net] まぁ何をどう言い繕ってみても再生産数を1未満、最低でも1以下にしない限り社会は崩壊する。 再生産数が1より大きい状態では何をどうやっても社会を維持することはできない。 849 名前：１３２人目の素数さん [2020/07/13(月) 20:39:03.48 ID:loo9MSbm.net] >>796 >K値というものが物理学者によって 提唱されて、大阪などでは利用されて いるようです 阪大核物理 Nが提唱し大阪府がK値を 利用していると明記してください。　 京大RIMSのMグループはIUTが数学 のabc予想から天気予報の改善まで 有効と主張しています。 関西の数理科学 850 名前：はアホくさく悲惨です がコロナ対策は命を守る問題です。 数学と物理は感染症分野は素人だから 生物数理モデルから個人的に学ぶべき であくまで素人の範囲でしょう。 生物数理モデルは故山口昌哉京大教授 がカオスを含め解析し弟子も活躍し 今では世界的に正統派です。 西浦北大教授の感染症対策の仕事や話も 山口先生の仕事の発展があると思って います。 [] [ここ壊れてます] 851 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 20:46:24.24 ID:DPFA8Q9h.net] All analyses were performed using R version 3.6.3 (R Foundation). なんて書いてあるんで読んでみたけど 割り算に使っただけみたいだ。 https://jamanetwork.com/journals/jama/fullarticle/2768351 対照群もなんにもない。小学生に夏休みの自由研究みたいだ。 852 名前：１３２人目の素数さん mailto:sage [2020/07/13(月) 20:48:12.18 ID:DPFA8Q9h.net] >>800 平均値と標準偏差の計算に使っただけだろうな。 853 名前：１３２人目の素数さん mailto:sage [2020/07/14(火) 14:27:26.92 ID:cML8ukUt.net] K値は基本役に立たない 一番役に立つのは実効再生産数 実効再生産数をモニタリング指標、および数値目標にして管理するのが最適 K値予測が当たるのは、一時的に実効再生産数が1を超えて、また1以下に戻る、みたいな場合 もちろん、実効再生産数が1を超えたら行政がすぐい対応して1以下に減らす対策をとることが前提 行政が1か月以上実効再生産数を減らす対策をとらない東京には当てはまらない 854 名前：１３２人目の素数さん [2020/07/14(火) 22:48:00.75 ID:MNivzfqp.net] SIRモデルによるCov19対策 ・法令システム.医療システム 　行政システムによる。 　SIRモデルから状況分析と予測より 感染者集団Ｙを減少させる つまりd Ｙ/dt<0にする。 SIRモデルで( S.I.R)→(X. Y. Z) とした。 ・SIRモデルは集団免疫モデルで 　対策がない場合を予測する。 　対策を考える場合も使える 　第1式　新規感染者の増加　 　第2式　感染対策の効果 　　　　=新規感染者-隔離感染者 　第3式　医療システム ・対策のための実効再生産数R 　定性的な式 　個人間は均等な接触から R→Ro×(個人間の接触回数割合) =(感染率)/(回復率)]× (個人間の均等な　　接触回数割合) SIRモデルのβは広義の(感染率)× (個人間の均等な接触回数割合) ・第2式で dＹ/dt<0のためには 　β感染率→小　接触回数→小 　γ隔離率→大　 　結果　Rは減少する γ→大の対策について ・第2式　 d Ｙ/dt=βS Y-γＹ 戦略的PCR検査対策。 隔離率には抗体獲得.ワクチンと 検査隔離がある。 抗体集団検査が約大阪015>東京0.1> 宮城0.01パーセントだから集団免疫 なしかつワクチンがない。 クラスターが予想される地域店などへ 特にPCR検査を増やし感染者を 病院.ホテルへ隔離することは隔離者が 増加し Yが減少する。 ・第3式医療システム d Z/dt=γY隔離感染者 検査隔離(γ→大)で 隔離者増大より病床数が満杯になり 結果医療システムが崩壊するから 軽症ならホテル隔離へ。 東京方針はRを小さくする。 ・大阪K値について 　データへ指数関数を当てはめた 現象論に基づくから対策に繋がらない 　戦略的PCR検査なし 本音は 経済優先でコロナはただの風邪路線 855 名前：１３２人目の素数さん [2020/07/15(水) 11:05:27.21 ID:+c3dKonF.net] >>797 単純な原理から、いかに自然現象をより高い精度で近似できるかって ことが、自然現象を理解するってことでしょ。科学における「真理」 という言葉は、そういうただし書きつきだと理解しないと。 そこらへんは数学とは事情が異なる。 自然科学において、絶対的な真理など空疎な言葉。そんなものは、 宗教屋が唱えるおまじないでしかない。 856 名前：１３２人目の素数さん mailto:sage [2020/07/15(水) 14:49:16 ID:nLcbkn+C.net] >>804 同じようにある条件下で成り立つ法則を数式などで数量も含めて表現したものと思う ニュートン力学は量子領域や宇宙領域では成り立たないだろうし 量子領域の物理で人間サイズを扱おうとするとシュレーディンガーの猫みたいになると思う ウイルスと人間の免疫の活動の結果として健康に大きな影響なければ普通の風邪などと同様に無視していいからミクロな視点も無視はできない ウイルスが体内でどのように増殖するのかも知るのも真理を知ることにはなると思う マクロはミクロの積み重ねだから 個人の行動変容が感染拡大に影響するのも確かなこと 同じウイルスと人間の活動をどんな観点で記述するかの違いの様に思う 同じ現象だから互いに整合するはずだと思う 857 名前：１３２人目の素数さん mailto:sage [2020/07/15(水) 17:41:50.43 ID:Y5eCgaZb.net] K値予測 katukawa.com/wp-content/uploads/fig6-1024x574.png 858 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 05:27:26 ID:f7C3/cjh.net] PCR検査、アメリカでは200人だかの検体をまとめて検査してるそうだ。 つまり検査結果が陰性なら全員がセーフ、陽性なら誰か陽性の人が1人以上いるぞ！ってこと 複数の検査を 検査1,2,3,4,5,6,7,8,9,10 検査a,b,c,d,e,f,g,h,i,j のように受けることでシロの人の洗い出しが素早くできたりはしないでしょうか またはその逆にクロのひとの洗い出し 859 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 07:13:17.08 ID:C2pOqxoK.net] >>807 サンプルを混ぜる手法はバランスパズル・偽コイン問題やiPS細胞の実験手法で有名よね www.jscf.org/jscf/katudou/Walk_Again_2008/08_10_5-1.pdf 860 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 08:22:41.73 ID:muripj9V.net] >>807 >>633の式を当てはめると、アメリカは陽性率の事前確率を25/100万ぐらいに見積もってるのだな。 861 名前：１３２人目の素数さん [2020/07/16(木) 09:19:06 ID:dzMbDhwI.net] >>805 マクロなシステムの振る舞いをミクロな構造から説き起こそうとする 演繹的な方法論には限界があるが、それでもそれしかないからね。 モデル化ってのまさにそういうことでしょ。 マクロなシステムの予測に役立つk値のようなパラメータを経験的に 見出すような、帰納的なやり方もありなのは確かにそうなんだろう。 だけど、熱力学が統計力学で置き換えられたように、そこでとどまっ ていては現象を理解した気にはなれんのよね。 やはり普遍的な基本法則から（たとえ数値シミュレーションの結果 であっても）演繹的に導かれるような説明ができないと理解した気 になれない。 862 名前：１３２人目の素数さん [2020/07/16(木) 09:22:20 ID:dzMbDhwI.net] >>807 へー、実際にやってるんだ。 863 名前：１３２人目の素数さん [2020/07/16(木) 14:37:18 ID:tJZRomem.net] >>803 訂正 ばつ.大阪015>東京0.1> 宮城0.01 まる.大阪0.17>東京0.1>宮城0.03 864 名前：１３２人目の素数さん [2020/07/16(木) 18:04:59.93 ID:d2bliECv.net] K値モデルによる予想。 >8府県(大阪兵庫京都愛知神奈川 埼玉千葉福岡) 7月9日ごろにピークアウト →大阪府認定のk値だが7月15日 　大阪61人では、大風呂敷を広げて 　根拠に欠ける曲線を入れても予測 　能力に全く欠ける。 　現象論的にも過去に検証がなく 　思いつきのレベルだが大阪府が　 　K値を全国へばらまいた 　第1波の厚労省集団抗体検査で 　大阪府が約0.17パーセントから 　推測する。 　府民882万人より15000人の潜在 　感染者がいて隔離感染者2126 人　　 (7.14)だから大阪府に潜在感染者が　　 　12.8千人もいた。 　医療システムが脆弱では戦略的PCR 　検査が不可能でk値やメディア出演 　のパフォーマンスよりその場を 　しのいできた。 　大阪夜の街ホストクラブで感染は 　徳島の指摘から発覚した。 　大阪府は関空の水際の規制緩和を 　要望など経済優先に前のめり 　よく知らんが都構想の都合がある 　のだろうか 865 名前：１３２人目の素数さん [2020/07/16(木) 18:11:15.37 ID:d2bliECv.net] >>813 つづき 移流項があるSIRモデル。 ・移流項は個人が一様性等方の運動を 　して個人間の感染,相互作用が仮定 される、 　ブラウン運動も一様等方の運動。 　移流項d∇^2(S、 I、R)より 　∂S/∂ t=-βSI＋d(∇^2)S 　 ∂ I/∂t=βSI-γY＋d(∇^2 )I 　∂R/∂t=γY ＋d(∇^2)R まず移流項なし閉鎖地域のSIR方程式を 解くこと。 解き方は力学系.近似解.数値計算などが あり力学系ベクトル場から閾値の存在が すぐわかる。 Murray数理生物学は感染症の章に ペスト狂犬病などがあり 空間.時空間の発展モデルがある。 更にロジスティック曲線もある。 念のためゴンペルツ曲線は全くない。 866 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 22:14:21 ID:lhXKpfMQ.net] >>807 検体をまとめて調査する手法は、検査能力が足りないときに検査数拡大したいときによく使われる手法 中国では、武漢の集団PCR検査で使ってた 韓国も、3月ころの検査数足りないときにやってた しかし200人ってのはちょっと信憑性かけるのでは？ アメリカの陽性率だと、200人もまぜたらほぼ全部陽性になって検査にならない 867 名前：１３２人目の素数さん mailto:sage [2020/07/16(木) 22:19:18.21 ID:lhXKpfMQ.net] 日本は学校なんかの全生徒検査とかがあれば、 検体混ぜるのを使えばいいとおもうよ 感染者と同じクラスの児童や担任、クラス外で親密な人→検体まとめずに個別検査、 それ以外→検体混ぜて検査 みたいな感じで ただし、紙とFAXで検査してる日本じゃ、検体ミックスなんてやったら保健所がパンクしたりミスがたくさん出そう ちゃんとコンピュータ化してれば、検体ミックスが有効に使える 868 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 05:10:32.66 ID:l8XG6+rW.net] 水中でも感染力を失わないならやばそう。 https://www.medrxiv.org/content/10.1101/2020.06.17.20133504v1 869 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 13:35:13 ID:nAA2dWKD.net] K値モデル完全失敗、西浦予想に近い形になる K値モデルは、感染者が増えたらあわてて政府・行政が対策するなら、K値モデルに近くなる 現実は、感染者増えても自治体も政府も無対策どころかさらに感染者ふやす対策をとる 870 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 15:50:29.96 ID:3kM9UPoB.net] >>815 >>816 毒ワイン問題の帰結から類推すると、陽性と出るか陰性と出るか1/2の確率になる� 871 名前：ｭらいの 人数分を混ぜて検査するのが、最適な方法ではないだろうか。 [] [ここ壊れてます] 872 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:10:22.26 ID:l/H5WeF7.net] >>819 >>419 873 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:30:29.34 ID:PvhkH04x.net] 誰かk値理論の論文読んで解説してくれん？ なんかアホらしくて探す気にも読む気にもならん。 1) 感染者数の推移についてどんなモデル曲線を採用してるのか、どんなパラメータを持ってるのか 2) そのパラメータをどんなモニタリング指標を使ってどう統計処理して推定するのか？ もちろん2)はk値とかいうのを使うんやろうけど 874 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 16:58:10.33 ID:l/H5WeF7.net] ツイッターかnoteかで解読してるヤツおらんかった？ 875 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:04:19.18 ID:PvhkH04x.net] なんかどうせデタラメやろうから調べる気にすらならん。 検索して読む時間全部無駄になる予感しかせん 876 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:08:29.39 ID:l/H5WeF7.net] 実際に大阪とかはK値モデル採用して意思決定してるんだから、その予測が出来るから無駄ではないと思うけど。 877 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:17:55.08 ID:PvhkH04x.net] つまりは大阪がどんなデタラメ理論に基づいて動いてるかわかる事に意味があると？ イヤ、ちょっとご勘弁w 878 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 17:48:59 ID:l/H5WeF7.net] >>825 発言権を得て大阪以外も追随する可能性があるし、トンデモ理論の成り立ちとかわかるから面白いと思うけどね。（大阪専門家委員にはウイルス学の宮坂いるんだけどこの人もちょっとアホだからK値に感心している。あと山中伸弥教授もいたようだから、K値と接触して宣伝する可能性があって危険。） 879 名前：１３２人目の素数さん [2020/07/17(金) 18:13:37 ID:WeNCWXbE.net] >>802 >行政が1か月以上実効再生産数を 減らす対策をとらない東京には 当てはまらない デマですね。 あなたはなぞのK値でも信奉してください 東京都は戦略的にPCR検査を増やし 隔離感染者を増大していますね。 これは基本的なSIRモデルで有効な 対策で実効再生産数を小さくします 880 名前：１３２人目の素数さん [2020/07/17(金) 19:43:00.71 ID:Ew4KkMJN.net] >>821 >>796にK値のポイントが整理されて いますね。 「論文」も検索できるし 881 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 21:58:21.84 ID:4GCjkvh6.net] >> 882 名前：828 なんかアホらしすぎて読む気にならんのよねぇ。 ただ例え便所の落書きとはいえ、k値理論に口出しするなら、やっぱ読まんといかんかなぁと言う気もするし。 なんでこんな怪しさ満載の理論採用されてるんかねぇ？ そもそも実績としてもことごとく予測値外してるのに [] [ここ壊れてます] 883 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 22:53:16.49 ID:l8XG6+rW.net] >>821 >796にでてくるゴンペルツ曲線 https://en.wikipedia.org/wiki/Gompertz_function ここに利用法として Examining disease spread と書いてあるな。 884 名前：１３２人目の素数さん mailto:sage [2020/07/17(金) 23:19:19.44 ID:CMxFmiKu.net] >>830 thx 未知パラメータ３つか。 コレをK値で推定すると。 しかしコレ増加しっぱなし、減少しっぱなしやん。 それでみんなどこから起点日に取るかの恣意性があるって言ってんのか。 それとも起点日もパラメータにしてちゃんと推計値出すのかな？ 読んだ？と思われる人のレスだと無さそうだなww 885 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 00:58:41 ID:nqMg34Lo.net] >>827 東京はクラスター対策がいい加減すぎて、PCR検査で実効再生産数を十分低くすることは難しい 中国並みに厳格なクラスター対策しない限り、東京のようなザルクラスター対策では 大して実効再生産数下げられない 886 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 01:04:31 ID:nqMg34Lo.net] ちなみに北京は、6月中旬からはじまった感染クラスターで、 普段の感染対策と、大量のPCR検査を組み合わせて抑え込みに成功した 高リスクグループから� 887 名前：�に毎日数十万件とかいう検査をおこなった結果抑え込めたわけで、 いまは低リスクグループの大量検査を実施中で、すでに検査数が1000万件を超えてる [] [ここ壊れてます] 888 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 07:39:06 ID:kzNofEnA.net] >>831 ゴンペルツ曲線って新規発生数じゃなくて累積患者数を表しているんじゃないかな？ ゴンペルツ分布というのがあるらしいが 新規発生数=pdf 累積患者数=cdfとすると パラメータによってpdfは単調にも一峰性にもなるみたい。 https://en.m.wikipedia.org/wiki/Gompertz_distribution 889 名前：１３２人目の素数さん [2020/07/18(土) 07:49:49 ID:b9ODbCJm.net] >>832 中国は共産党独裁国家だし 大阪は検査もろくにしなかった。 法令システムでは 措置法は規制がなくザルで厚労省は 食中毒並みのPCR検査義務と休業命令 も通達しない。 国、政府、国会の法改正責任が重大なのにGO TO キャンペーンへ、 まあコロナはただの風邪路線なのだろう 890 名前：１３２人目の素数さん [2020/07/18(土) 07:55:01 ID:b9ODbCJm.net] >>829 論文を読むとインチキ手口がわかるよw 891 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 08:49:38.59 ID:bK/VOAbE.net] >>834 そうかもね。よんだ？ >>836 自信が確信になってもなぁw 892 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 09:17:12 ID:bK/VOAbE.net] >>830 しかしなぁ sir方程式の解はピーク時点は左右対称になんかならんからな。 もちろん現実世界では再生産数は社会の対策で変動してるからsirの解通りには変動しないけど。 なのでどこでピークアウトするかは理論的に予測できるハズはほとんどないんだけどな。 893 名前：１３２人目の素数さん [2020/07/18(土) 10:06:46 ID:b9ODbCJm.net] >>837 K値論文の仮定が根拠不明だ。 SIRは数値計算を行なっている。 実際の測定値は検査などによる隔離感染者数だな 感染者数Iは推定で対策よりγが 変わればSIR予測も変わる。 894 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 11:49:40.36 ID:ILbvGgBu.net] >>839 まぁまだ論文読んでないけど(そして読む気も消え失せたけどw)大体インチキなのは確定したな。 モデルの曲線に何とるにせよ、そのモデルの曲線はなんかの意味もって選定しないと話にならん。 sirのモデルではsir方程式 d/dt I = β SI - γ I (感染者の増加分は感受性宿主数と感染者数に比例し、減少分(死亡、回復する数)は感染者数に比例する をとる。 「そやろな」と納得いく式だし、実績もある。 この方程式に従えば感受性宿主数がγ/βに到達するまでピークアウトしない。 にもかかわらず実際の第一波でははるかに早い段階でピークアウトしてる。 それは感染率βが定数ではないからだ。 行動制限などによって接触率が下がれば感染率βも下がる。 第一波が早期にピークアウトしたのは緊急事態宣言下で市民が接触率を下げた効果。 つまりどこでピークアウトするかは市民の生活様式によるものでとても現在の新規感染者数から予測できるようなものじゃない。 もうそれを現在の新規感染者すうから計算できるK値なるもので予測できると言ってる時点でインチキ確定と言っていいやろ。 もう全く読む気が失せたww 895 名前：１３２人目の素数さん mailto:sage [2020/07/18(土) 11:51:05.18 ID:ILbvGgBu.net] あ、達するまでというか減るまでね。 896 名前：１３２人目の素数さん [2020/07/18(土) 14:36:31 ID:+nayeGtY.net] K値論文 Indicator K to analyze the transition of COVID-19 spread In this study, we introduce a new indicator called the K value defined by K(d) = 1−N(d − 7)/N(d), where d is the number of days from the reference d 897 名前：ate, and N (d) and N (d − 7) are the total number of infected people on days d and (d − 7), respectively. Since N (d) is greater than N (d − 7) during the period from the initiation of spread to convergence, K takes a value between 0 and 1. ( infected people 感染者) k = 1 + 2.88K′ (Fig. S1). この公式を図へあてはめたのか、 図から帰納的に公式を導いたのか、 特に2.88の由来はなんだろう？？ [] [ここ壊れてます] 898 名前：１３２人目の素数さん [2020/07/18(土) 23:28:48.51 ID:b9ODbCJm.net] K値について調べてみました https://tatsuharug.com/kvalue 899 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 01:35:49 ID:m6+VaLKl.net] >>842 Indicator K の K はなんの K ? 900 名前：１３２人目の素数さん [2020/07/19(日) 06:37:20 ID:GNq9W41i.net] KはあまカッパのKからともいわれている が、よくしらない。 阪大か大阪府へきいてみたら 901 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 07:29:35 ID:w6BaQv4H.net] >>842 R0 = 2.5であとは適当にパラメータ設定したSEIRモデルにK値を重ねてグラフにしてみた。 https://i.imgur.com/5AsiQxW.png どうやって予想に使うんだろう？ 902 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 08:00:38 ID:w6BaQv4H.net] >>846 変数名が別のデータと被っていたので rm(list=ls())でメモリクリアしてから 作図し直した。 https://i.imgur.com/KpwhZ0G.png 903 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 08:24:34 ID:w6BaQv4H.net] >>847 K値の基準となる日を１週間ずつずらしてK値の変動をグラフにしたみた。 https://i.imgur.com/fkN9JLg.png 904 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 21:02:46 ID:iNIKliXH.net] K値モデルが正しく、なにもしなくて収束するなら、 アメリカやブラジルはとっくに収束してるわ K値モデルがダメな点は、ただの一学者の主張じゃなくて、 大阪府知事が信じて取り入れてしまった点だな やっぱ実効再生産数がいちばん大事 905 名前：１３２人目の素数さん mailto:sage [2020/07/19(日) 21:20:04.76 ID:iNIKliXH.net] なぜいまだに政府も東京都も大阪府も、 実効再生産数をモニタリング指標および政策目標につかわないのかまったく不明 東京と大阪以外は実効再生産数0.7、東京と大阪は実効再生産数0.8 これくらいを政策目標にすればいいのに 一時的に実効再生産数が0.5以下になるまで規制して感染者を大きく減らし、 その後は実効再生産数0.8目標でいいとおもう 906 名前：１３２人目の素数さん [2020/07/19(日) 21:56:54.03 ID:GNq9W41i.net] >実効再生産数 山中京大教授は実効再生産数のエクセル用公式を見つけて数値いれて威張って いますが、 実効再生産数は感染症数理モデルに基づきますのでマネすることはおやめください 907 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:09:27.09 ID:Ze2tCmmK.net] 疫学者が頑張って計算する正確な実効再生産数じゃなくても、 日々の確認陽性者数をエクセルに入れて処理する簡単な実効再生産数でも十分役に立つんだよ ネットで公表されてる簡易モデルは、実際の確認感染者数の増減とだいたい一致してる 908 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:13:48.19 ID:Ze2tCmmK.net] 日本の都道府県別実効再生産数 ttps://rt-live-japan.com/ ttps://covid.gutas.net/ アメリカの州ごとの実効再生産数 ttps://rt.live/ 世界の実効再生産数 ttps://covid.gutas.net/world 疫学者ががんばって計算する実効再生産数じゃなくても、 このレベルの簡易モデルで、十分通用するのが今回の新型コロナだよ 日本の場合、感染者がかなり少ない都道府県はまともに算出されないが、 それ以外はだいたい妥当な値になってる 909 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 00:30:32.49 ID:Ze2tCmmK.net] ドイツやイタリアは実効再生産数ベースに対策を考えてるのでは？ 実効再生産数が上がらないように対策してる感じがする 日本は、日本全国および東京都の実効再生産数が継続的に1を超えてるのに、 感染対策をやらずに自粛解除だけ進めていったのがこの事態を引き起こした 本来、自粛解除しても実効再生産数が１を超えないように感染対策をした上で 自粛解除をする必要があった 910 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 01:56:42.64 ID:L7jTb9sR.net] そう、なのに日本だとわけのわからんモニタリング指標に従って対策立ててる アホじゃないかと 911 名前：１３２人目の素数さん [2020/07/20(月) 08:24:49.04 ID:IvKiJrmj.net] 感染者数　成長曲線の3パターン ・指数関数型 マルサス人口論　原子核崩壊　 　官邸　アベノミクスカンフル剤 　コロナ対策より 　株価が指数関数型の上昇希望 ・ロジスティック曲線 　dX/dt=N(1-X/N)X SIモデル=SIRモデルの第1式.第2式 (隔離率γ=0)の感染者増加曲線。 PCR検査より感染者の隔離対策は 過去に実証された有効対策であり 　SIRモデルに基づく ・ゴンペルツ曲線 　高齢者の人口減少 　腫瘍ガン細胞の体積増加 　大阪府地方で流行っているK値 　=感染者数の新興占いの式。 　神奈川も流行り始めた模様 912 名前：１３２人目の素数さん [2020/07/20(月) 08:44:21 ID:IvKiJrmj.net] 東京都のコロナ対策はSIRモデルが 重要な役割をしている。 913 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 08:55:29 ID:mgzJmLOO.net] SIRって言ってもRは当分の間ほぼ0のままだから、事実上SIモデルだよな。 914 名前：１３２人目の素数さん [2020/07/20(月) 09:31:31 ID:IvKiJrmj.net] 積極的にPCR検査しない道県は SIモデル。 府はガン細胞増殖のK値モデル 915 名前：１３２人目の素数さん mailto:sage [2020/07/20(月) 10:24:10.81 ID:jF/Tht4x.net] γは回復率やろ 平均回復日数の逆数 916 名前：１３２人目の素数さん [2020/07/20(月) 10:47:58 ID:AYBv/aE8.net] γは隔離率とも回復率ともいう。 抗体ができる、死亡、隔離より感染しない者を表わす係数 917 名前：１３２人目の素数さん [2020/07/20(月) 11:45:14 ID:AYBv/aE8.net] ・SIRモデル 閉鎖区域の総人口Nが一定。 S＋ I＋R= N 一定。 感染可能者S→潜期→ウィルス排出 →発病I→一部回復(抗体).隔離.死亡 などR、 ウィルス排出→発病の期間0とする。 第1式　dX/dt=-βXY 第2式dY/dt=βXY-γY 第3式d Z/dt=γY. ・発展SIRSモデル 　一部R→Sがある場合 918 名前：１３２人目の素数さん mailto:sage [2020/07/23(木) 02:43:39 ID:iTEunCvL.net] 大阪モデル完全崩壊だな K値唱える人がいるのはいいが、それを政治が採用するとかやばすぎる 実効再生産数を見るのが最適 ただし、実効再生産数はクラスター対策が不可能になったら突然急騰するので、 クラスター対策できる上限には注意が必要 919 名前：１３２人目の素数さん [2020/07/23(木) 18:41:45.80 ID:bfbrr4ZK.net] 実効再生産数は数理モデルによる。 数学として SIモデルからロジスティック方程式 となり感染率βを小さくする対策 新しい日常不要不急不急の外出禁止 が必要だ。 更にオイラー差分を使い離散ロジスティック方程式が得られる。 シャルコフスキーの定理、リ.ヨーク の定理のカオスの意味で β Nよりカオスになる場合があるが、 感染症数理モデルとしては感染症 の基礎知識が必要だろう。 920 名前：１３２人目の素数さん [2020/07/24(金) 18:33:22.00 ID:g3DPzG/z.net] K値って、大阪だけじゃなく神奈川も使ってるのか 知事がアホだとこういったオカルトに走り出す 921 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 21:41:46.70 ID:E42oSbW+.net] 世間は無駄に大騒ぎだが以下の事実がある 感染者増と全く比例してない おまえら全く割出せなかっただろ？ つまり今後の予測もギブアップだよな 死者数 7月*1日　*1人 7月*2日　*1人 7月*3日　*1人 7月*4日　*0人 7月*5日　*0人 7月*6日　*1人 7月*7日　*2人 7月*8日　*1人 7月*9日　*1人 7月10日　*0人 7月11日　*0人 7月12日　*0人 7月13日　*0人 7月14日　*2人 7月15日　*1人 7月16日　*0人 7月17日　*0人 7月18日　*0人 7月19日　*0人 7月20日　*3人 7月21日　*1人 7月22日　*1人 922 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 22:44:07.67 ID:leZfEHQl.net] おいおまえら また新しいの出てきたぞ K値　改め　Rw値 wwwww 923 名前： [] [ここ壊れてます] 924 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 23:11:05 ID:kQJgVFx5.net] 医療提供体勢が整っているから死者数が少なくなるんじゃね？ 925 名前：１３２人目の素数さん [2020/07/24(金) 23:11:27 ID:2uKiwAiC.net] 死亡者は重症感染者数による。 第1波の重症者は回復か死亡より 激減したが、今後第2波が GO TOキャンペーンで拡大し 若年軽症から感染した重症者 (特に高齢者)が増加するだろう。 東洋経済の実効再生産数公式 西浦教授監修とのこと 7月23日 全国　1.4 ・東京1.20神奈川1.19埼玉1.19千葉1.1 ・大阪2.46兵庫2.23京都1.5 ・愛知9.51 岐阜3.88 ・北海道1.29 福岡 K値予想の大阪 新規感染者149人 陽性率13パーセント。 大阪と官邸はコロナはただの風邪路線「GO TOキャンペーン」 全国へコロナウィルスばらまき 926 名前：１３２人目の素数さん mailto:sage [2020/07/24(金) 23:13:52 ID:kQJgVFx5.net] 薬やその使い方も少しわかってきただろうし 血栓に注意が必要ということもわかってきただろう サイトカインストームにステロイドを使うとか 回復した人の血漿を治療に使うとか 対処法が変われば結果も変わるだろう 927 名前：１３２人目の素数さん [2020/07/24(金) 23:14:10 ID:2uKiwAiC.net] >>869 追加 福岡3.63 928 名前：１３２人目の素数さん [2020/07/25(土) 01:04:16.54 ID:3rBi8tDq.net] >>870 確かに。 対症療法とはいえ、それなりに対策が立てられるようになってんのかもね。 アメリカもUKも新規死者数の減少幅のほうが、感染者の減少幅より大きい （致死率が半分くらいに下がってる？） 929 名前：１３２人目の素数さん [2020/07/25(土) 22:57:57.28 ID:M+D5uJMO.net] 発展形も含め基本のSIRモデルが コロナ対策の基盤で ・第3式から隔離者増は医療システム内。 医療崩壊がおきるとコロナに限らず、 死亡者増加を招く。 野党は日本でも独並みのPCR検査7万件をと喚いているが、現実の医療システム 無視で机上の妄想。 与論島などベッド数が少ないし 医療崩壊危機がすぐおきる。 また個人の任意PCR検査は高額だから 政府か共産党、 まあ国会議員が寄付すれば良い 930 名前：１３２人目の素数さん [2020/07/25(土) 23:04:39.05 ID:M+D5uJMO.net] ・移流項つきSIRモデルは感染経路 対策。 ・GOTOキャンペーン推進の、官邸や 大阪関西愛知などは東京問題 東京由来などと非科学的な発言を している。 東京発新幹線は名古屋関西方面の 感染爆発地域も東北長野新幹線方面では 落ち着いていて実効再生産数もこれを 示している。 第1波後厚労省の集団抗体検査は 宮城0.03<東京0.1<大阪0.17パーセント。 大阪関西愛知などPCR検査も過少で 第1波が地域で燻っていても見逃して いただけだ ・海外からの感染者入国管理は 政府の責任だが成田空港も毎日感染者 が発見されている。 水際対策は行動アプリ義務もなく ズブズブの上更に緩めようとしている。 西浦教授の試算で毎日10人の感染者が 入国すれば3ヶ月で感染爆発だ。 ・野党国会議員は緊急事態宣言中に 新宿歌舞伎町で歓楽するバカ丸出し だが、横田基地の兵隊さんも歓楽し ている？ まず防衛省は米軍が軍人従業員の感染 を管理しているか情報を開示すべきだ。 野党は沖縄の米軍基地だけに興味が あるらしい、、 931 名前：１３２人目の素数さん [2020/07/25(土) 23:07:21.66 ID:M+D5uJMO.net] 国際空港は関空も中部国際空港も あるだろ 932 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 08:31:58.69 ID:g9nf//ah.net] 君ら肝心なことを認識してないけど、感染自体を敵視するのは間違いだからね コロナウイルスは4種類あり、今回が5種類目 コロナウイルスは風邪ウイルスであり、新型コロナウイルスもまた 933 名前：風邪ウイルス 風邪にはワクチンも治療薬も無い だが新型コロナウイルスにはすでに治療薬があるしワクチンも出てくる予定 既存のコロナウイルス同様、共存していく時代に入った 感染自体で大騒ぎするのは間違いだ 重症化に陥らないよう対処するためのノウハウを周知していくほうが大事 [] [ここ壊れてます] 934 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 08:46:16.62 ID:4EvbZFoF.net] >>876 アホか 防疫に治療法もへったくれもあるか むしろ防疫学ではrecover(回復)するか死亡(removed)するかなんか区別せずどっちもR、まぁ再感染しない場合だけど 実際治療法がいくらできても現状の倍化時間が10日前後≒再生産数2ちょっとの状態なら１ヶ月もしたら医療崩壊するわ 935 名前：１３２人目の素数さん [2020/07/26(日) 11:36:41.02 ID:SV5vsEiJ.net] 現状で比較すると感染者数でドイツレベル、死者数だとまだ ドイツの数分の一なんだよね。 経済活動を戻しても現状レベルで推移できればいいんだけど、 ちょっと怪しいよねぇ。マスク着用の義務化とかで、新しい 日常のガイドラインを守るような強制力を行使しないと無理 かもしれん。飛沫を人に向けてばら撒かなきゃいいだけなん だけどねぇ。 936 名前：１３２人目の素数さん [2020/07/26(日) 11:38:01.20 ID:SV5vsEiJ.net] >>878 累積じゃなくて、新規の数での比較ね。 オーストラリアの新規感染者の推移も日本と似てるけど、 向こうは死者数がかなり増えてる。 937 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 11:43:38.72 ID:p4sxlIGA.net] >>878 どんな対策が有効か数学畑にはわからんし、数学板にはそんな話でイキって意見言う奴おらんやろ ホンマに全員が公衆衛生上正しい行動したら再生産数が1切るかどうかなんか数学でわかるハズがない 数学でわかるのは現状の再生産数が変わらなければ１ヶ月後どんな事態になるのかとか 938 名前：１３２人目の素数さん [2020/07/26(日) 12:09:13.87 ID:SV5vsEiJ.net] >>880 同じ話が4月の段階でもあったわけだけど、そのときと何か違うか？ 実効再生産数が現状のままなら、ってのは、単に現状の指数関数的推移を 外挿するとどうなるかっていう超単純な議論と大して変わらんのでは？ 公衆衛生対策が数理モデルにどう反映するかを考えなきゃ面白くもなん ともない。 939 名前：１３２人目の素数さん [2020/07/26(日) 13:27:15.33 ID:fbUsmkh7.net] 大阪府のK値は全く感染症の根拠に欠ける 感染症の予測に使うことがトンデモ だろう 940 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 15:34:13.59 ID:p4sxlIGA.net] >>881 4月でも5月でもコレから先もずっと数学でできるのは感染者数がどのように推移するのか予測することしかできん 現在の再生産数がどれくらいなのかをモニタリング指標(新規感染者数)を用いて推定値を出し、予測値を出す、コレが全て。 941 名前：１３２人目の素数さん [2020/07/26(日) 16:35:58.67 ID:SV5vsEiJ.net] >>883 で、４月の感染ピークの予測は当たったんだっけ？ 西浦先生の予測値はどうだったっけ？ 942 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 17:14:37.06 ID:hvMjKhGS.net] >>883 指標値をこの値にしたらこうなるという予測も立てられる それを実現するのは個人や組織の行動による あと行政が行動を上手く誘導できるかどうか 税制とか給付金とか罰則とか会見などでの演説とか 943 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 17:46:33.47 ID:UGJIzTTp.net] >>884 当たったるよ ほぼ再生産数は2.5前後で間違いないよ >>885 そう、それが西浦先生がやった感染率を8割削減したらどのくらい影響が出るのか、7割ならどうかという予測 その数値を見ながら政治がどう判断するかは数学の預かり知るところではない 944 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 18:37:06.63 ID:hO4a28mo.net] >>886 え、2.5だった？ 945 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 19:28:46.58 ID:FDP28PFB.net] >>887 色々それらしい解(新規感染者数の増減と抗体検査の結果に付合する値)は見つけたけど2.5はないと到底抗体検査と符合しない 多分もっと大きい 946 名前：１３２人目の素数さん [2020/07/26(日) 20:33:01.96 ID:TwUfJ85a.net] 　基本のSIRモデルで第1式から 　SIモデルはd I/dt=βSIでβが集団S 　と集団Iの接触相互作用の大きさを 　表す。 ・高校化学で平衡状態の化学反応 　A＋B→Cから反応速度は d[C]/dt=a[A][B]となる。 接触率β⇔ 反応速度定数a 　[⚫]は濃度を表す。 　実際に感染症数理モデルには 　質量作用の羅患率の用語がある。 ・更に簡単な化学反応の分子衝突 　モデルは一様等方な分子運動の 　弾性衝突(均等な衝突) だ。 　高校物理で気体分子運動論があり 　PV=zRTの式があるが、 　浸透圧でもなりたつ。 ・SIRからの感染症の対策。 　数理モデルの感染⇔化学反応の 　対応から考えると、三密を避ける 　(濃度を小さくする)マスクをする 　(反応断面積σを小さくする) ・高校数学化学物理のポイントを 　身につけている人たちには、 　接触8割削減が感染を防ぎ結果実効 　再生産数を下げると納得するだろう 　まあロックダウンとは、インフル 　エンザ対策でも学級封鎖をします 947 名前：１３２人目の素数さん [2020/07/26(日) 20:35:08.92 ID:TwUfJ85a.net] 大阪府のK値のように思いつきの 　ゴンペルツ曲線からガラクタ 　のシステムへ数値を入れても結論は 　ガラクタ対策にすぎない。 　アインシュタインの1905年ブラウン 　運動.原子の実在も読んだこともない だろうがPV=zRTからはじまる。 　K値は感染症の熱素？ 　これではマッハに失礼すぎるなw 　 948 名前：１３２人目の素数さん [2020/07/26(日) 21:04:52.86 ID:SV5vsEiJ.net] >>886 1.5くらいじゃねーの？2.5もあるか？ 西浦予測の影響はなかったってのが結論じゃなかったっけ？ それ以前に収束が始まってたんじゃなかったっけ？ 949 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 21:18:45.37 ID:TKDcUAL0.net] >>891 無理 東京の場合 仮に3月上旬くらいから拡大が始まって緊急事態宣言から再生産数を1未満に落として5月末で一旦収束、そして総罹患者数が14000になる コレラの条件を満たす解は3月中の拡大期の再生産数が2.5未満なんかでは到底届かない 950 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 21:23:40.24 ID:hvMjKhGS.net] 一部で新型コロナウイルスに暴露しても 自然免疫で対応できる間は発症しないし抗体も獲得しないという説がある 自然免疫だけでウイルスを除去出来る人は抗体もできないから抗体検査でも陽性にならない 951 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 21:34:50.81 ID:TKDcUAL0.net] >>893 いわゆる山中先生が言うところのfactorXが自然免疫によるものとする説だな その説が正しいかどうか数学的にわかるハズないのでそれをここで議論してもスレ荒れるだけだから無駄だろ そもそもその説が正しいなら再生産数はさらに大きいことになるしw 952 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 22:18:02.44 ID:fB//8lKq.net] >>892 自分の計算では、最も高い頃で2ぐらいだけどなあ。 2が一か月続いたら、SerialIntervalを5日として6世代、64倍になる。4月上旬に5人、それから1か月後もまだ200人に達していないから、少なくともその期間平均は2を超えていないはず。 953 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 22:26:14.19 ID:TKDcUAL0.net] >>895 オレは西浦先生の評価に異論はないけど、異論があるならココに書いてみればいい ココは数学板 キミの出した答えが合ってるか検証できる人間はゴロゴロいる 3/10くらいのときの感染者xxx人、 そこからxx日目まで再生産数はxx、 そこからxx日目まで再生産数はxx、 そこからxx日目まで再生産数はxx、 ‥‥ と書けば感染者数がどう推移するかこの板の住人なら誰でも計算できるし、その解が新規感染者数の認知件数や抗体検査による総罹患者数の値と矛盾しないかチェックしてもらえる それで西浦先生の推定がおかしいという結論がでるなら数学板のお墨付きで他の板にも書けるよ 954 名前：１３２人目の素数さん [2020/07/26(日) 23:09:26.16 ID:PbHdIAnK.net] >>895 実行再生産数と基本再生産数を混同してるのでは？ 955 名前：１３２人目の素数さん mailto:sage [2020/07/26(日) 23:34:46 ID:mF0auFFv.net] 君らに良いこと教えてあげるけど 元々ビタミンDはコロナウイルスに効果があるエビデンスがあった 新型コロナウイルスにも効果があることが世界各国の臨床試験で判明してる 罹患するとビタミン・ミネラルを猛烈に消費し欠乏状態になる なかでも抗ウイルス作用があるビタミンDと亜鉛が顕著で、亜鉛が欠乏すると味覚障害になる 新型コロナウイルスで味覚障害が出るのは亜鉛が欠乏するから 味覚障害からの回復が遅い人は普段の食生活に問題があることを示す（亜鉛摂取不足） 普段からビタミンDや亜鉛の血清濃度が高い人は感染しずらいし感染しても無症状で済むし後遺症も出ない 956 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 00:06:00.51 ID:sQ0hXnoT.net] 重症患者は ビタミンD欠乏 亜鉛欠乏 カリウム欠乏 低アルブミン グルタチオン欠乏 が起こってることが分かっている なかでもビタミンDの臨床試験は多くの国で行われ論文が上がってる ちなみに日本からは全く上がってこない（これ、日本がかなり深刻な問題を抱えてることを示してる） 分かる人は分かったと思うが、ビタミンDと亜鉛は新型コロナウイルスの予防に使えるんだよ 957 名前：１３２人目の素数さん [2020/07/27(月) 10:12:17 ID:zGP/Kefd.net] NHKスペシャル。 国立感染研　鈴木基センター長 7/17時点で7/7の東京のRtは1.2。 昨日今日が300なら400から500まで 行ってもおかしくないが倍々に増えることはないと推測。 心配する報道が増えて行動が抑制されてきている。 一方、関西は東京より高い1.7から2 で流行拡大リスクが高い 東洋経済　実効再生産数Rt 7月25日 全国　1.33 東京1.12 埼玉 1.37 千葉0.96 神奈川0.97 大阪2.02兵庫1.56 福岡3.16 愛知3.48 岐阜2.82 Rtから見ると確かに、東日本都市が 埼玉以外が全国1.33未満。 西日本都市が感染拡大で、第1波の残存 ウィルスと国際空港からの輸入が 疑われる。 政府と東京都以外の道府県はGo To キャンペーンの移動推進政策だが 現況に基づかず西日本都市を中心に 全国へウィルス拡散となっている 958 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 11:11:34 ID:JaT+AB6z.net] >>894 感染者数と重傷者数の割合が変わるから検証するときに考慮しないといけなくなるし 自然免疫でウイルスがなくなるなら基本再生産数の値も変わるのでは？ 959 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 11:16:10 ID:JaT+AB6z.net] >>900 クラスターか感染経路不明の市中感染かでも対策は変わると思う Rtにはその区別はないと思うけど 960 名前：１３２人目の素数さん [2020/07/27(月) 11:44:07.43 ID:mjBkePcb.net] >>902 第1波後大阪府の集団抗体検査 では0.17パーセントあり ざっくりと府内で15000人が 感染していた。 961 名前：１３２人目の素数さん [2020/07/27(月) 12:03:45.11 ID:mjBkePcb.net] >>903 詳しくいえば 厚労省による大阪府の集団抗体検査 962 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 12:21:02.25 ID:0LgObjZ6.net] >>901 感染したけど抗体できなくで治るなら感染者数は抗体獲得者数より大きくなるので 963 名前：抗体検査の予測値三県で29000人より大きくなります 再生産数2.5で29000人に届かないのだから29000より大きくなればさらに再生産表式は大きいということになる [] [ここ壊れてます] 964 名前：１３２人目の素数さん [2020/07/27(月) 14:37:26.30 ID:s+Y7T6ku.net] >>892 1.5くらいじゃね？って言ってるのは現在の値だよ。3月、4月の話じゃなくて。 965 名前：１３２人目の素数さん [2020/07/27(月) 14:40:58.64 ID:s+Y7T6ku.net] >>904 抗体検査は特異度が低いからなぁ。 それは上限だろ。いや、抗体ができなかったり、消えちゃったり するから上限とも言えんか。まあ、そのくらいいても不思議は ないけど、はるかに少ない可能性もはるかに多い可能性もある。 で、実質意味無しｗ 966 名前：１３２人目の素数さん [2020/07/27(月) 14:49:11.54 ID:wv4NRHrQ.net] >>907 定量検査 967 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 14:52:09.80 ID:3aL20tW8.net] >>906 現在値どっかで月毎6倍らしいからそんなもんかもね 968 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 17:27:30.04 ID:JaT+AB6z.net] >>905 再生産数は一人が何人に感染させるかという値だと思う 自然免疫で治癒するなら 他の人に感染させなかったり ウイルスを体内に取り込んでも治癒して他の人に感染させない可能性があるという事 自然免疫で治癒する人は再生算数の計算から除外するなら影響ないけど 自然免疫だけで治癒するという仮説が正しいなら 自然免疫で治癒可能な集団とそうでない集団に分けて考える必要もあるかなと思う 969 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 17:39:40.00 ID:ZGna9uAs.net] >>910 違う 自然免疫で撃退できる分がどうこうとか再生産数なんかに関係しません 例の高橋とかいうバカの仮説があつてようが間違ってようが3月当時の再生産数が2.5未満の解なんて存在しません あるなら見つけて論文発表なりなんなりして下さい 970 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 17:56:23 ID:JaT+AB6z.net] >>911 自然免疫で治癒できる人が100%だったら実効再生算数は0になるのでは？ 自然免疫と実効再生算数は関係あると思うけど 971 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 18:22:27 ID:bRN6Z2LH.net] >>912 なるわけないやん？ 意味わかって書いてるの？ 972 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 18:26:56 ID:bRN6Z2LH.net] >>912 ああ、もしかして自然免疫で一人も写らなかったら再生産数0って？ そりゃそうだwwww そういう意味ではない ともかく再生産数 Sβ(t)/γ が2.5より小さい解見つけたら論文にでもなんでもすればよろしい 973 名前：１３２人目の素数さん [2020/07/27(月) 18:55:03 ID:s+Y7T6ku.net] >>914 おっしゃってることが支離滅裂ですよ。 974 名前：１３２人目の素数さん [2020/07/27(月) 19:11:44 ID:IOEE9iOc.net] >自然免疫で治癒する 自然免疫で治癒=ワクチン投与者 とみなす。 SIRモデルで初期条件は ワクチンを投与された者を 初期の隔離感染者roとする場合、 全人口N=so＋lo＋ro =一定 対策も考え基本再生産数 Ro≡βN/γとすると形式的に 変わらない。 一方、初期条件がN'=N-ro=so＋lo ＋0=一定とすれば基本再生産数が 小さくなりRtも小さくなる。 どちらも方程式の解は同じになりそう 975 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 20:09:15 ID:bRN6Z2LH.net] >>915 わからんけど、その自然免疫で撃退できる場合には再生産数の計算式が変わるという理論を論文にでもして出せばいい sirにおける再生産数は β(t)S(t)/γ ＝ 単位時間あたりに一人の感染者が生み出す二次感染者の数×平均罹患日数 別の定義なり方程式なり新たに作りたいならどうぞご勝手に 976 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 20:55:01 ID:ck4ARvK3.net] 問題 PCRプーリングしてスクリーニング検査する際、検体をいくつまとめて検査するのが検査数を減らせるのか、 それぞれ以下の陽性率について算出せよ 陽性率0.01%、0.1%、1%、10% ただし、1回目の検査はプーリング法を使い、2回目の検査はプーリング使わずに全員個別検査をおこない、 検体採取数は最大2回とする 977 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 22:06:48 ] [ここ壊れてます] 978 名前：.80 ID:I+okm8ek.net mailto: >>918 >>633 productlogはランベルトのW関数な。 [] [ここ壊れてます] 979 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 22:11:54.18 ID:+OVelfiF.net] n検体を同時に検査するとする。 n検体全てが陰性である確率は　(1-p)^nで　この場合、1回の検査で検査終了 n検体の中に、陽性が含まれている確率は　1-(1-p)^nで　この場合、n+1回の検査が必要。 n検体を同時に検査する場合に、必要な検査回数の期待値は、1×(1-p)^n + (n+1)×{1-(1-p)^n}　であり、 1検体当たりに必要な検査回数の期待値は、これに　(1/n)　を書ければ良い。 つまり、f(x)=(1/x){(1-p)^x + (x+1)×{1-(1-p)^x}} = 1 + 1/x -(1-p)^x を最小にするxを求める問題に帰着する p=0.1 の時　x=3.75458... ;　f(3;p=0.1)=0.604333、f(4;p=0.1)=0.5939　なので、4検体が最適 p=0.01 の時　x=10.5162... ;　f(10;p=0.01)=0.195618、f(11;p=0.01)=0.195571　なので、11検体が最適 p=0.001の時　x=32.1271... ;　f(32;p=0.001)=0.0627589、f(33;p=0.001)=0.0627804　なので、32検体が最適 p=0.0001の時　x=100.501...;　f(100;p=0.0001)=0.0199506613、f(101;p=0.0001)=0.019950656　なので、101検体が最適 980 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 22:16:01.91 ID:JaT+AB6z.net] 自然免疫で治癒する人が多ければ 集団免疫を獲得したのと同様になるだろうから 値を求める数式は同じでもその値には影響すると思うな 有効な免疫を持っていない人だけの集合 有効な免疫を持っている人だけの集合 それぞれ50%が混在した集合 で実効再生算数を考えると値は違うものになるのでは？ 劇場とかの空間に有効免疫を持っている人が多ければ その人達が壁になって飛沫感染は防げる可能性が上がる ドアノブなどを介した接触感染には影響ないかもしれないけど 981 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 22:24:01.56 ID:JaT+AB6z.net] 抗体を獲得した人は自然免疫では治癒しなかったのだから 抗体獲得率の測定値や再生算数などと比較することで 自然免疫で治癒する人がどの程度の割合で存在するか推定出来そうな気がする 982 名前：１３２人目の素数さん mailto:sage [2020/07/27(月) 22:40:42.66 ID:BVT2YP9e.net] >>921 ああ、つまり絶対罹患しない人が幾ばくかいる感染症なら総罹患者数はもっと落ちる可能性があると。 それはそうだな。 自然免疫で治癒する人がいて、その人が何回ウイルスの被曝を浴びても感染しないならその人はRに入るな。 しかし結局実効再生産数が β(t)S(t)/γ である事に変わりはないし、その場合のI(t)が満たす方程式 d/dt I(t) = β(t)S(t)I(t)-γI(t) も変わらない 仮に絶対罹患しない人の割合かpだけあるならその場合の実効再生産数は(1-p)β(t)S(t)/γ (有効率×接種率=pのワクチンを打った状態)になるけど、結局I(t)<<S(t)におけるその解の振る舞いは感染率が(1-p)β(t)の感染症に対する方程式と同じになる なのでもしそんな人が割合pでいるなら3月中の再生産数が2.5という表現は再生産数×(1-絶対罹患しない人の割合)が2.5というべきかもね。 まぁそれなら西浦先生の総罹患者数の予想にも影響は出てくるな まぁアホ高橋の仮説を入れればだけどなww 983 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 00:56:10 ID:p32FeE4v.net] あまり検体混ぜすぎると感度が落ちるので、 10検体くらいが現実的？ 血液検査なら、20検体ミックスとか可能だと思う 984 名前：１３２人目の素数さん [2020/07/28(火) 09:47:33.57 ID:/xPjTU2o.net] >>921 自然免疫で感染が防げて、そういう人が多いのなら、たとえば ミニシアターなんかの集団感染での感染率からある程度わかり そうな気もする。報道による情報が少なすぎてよくわからんけど、 濃厚接触した人の何割が感染したんだろうねぇ？ 985 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 10:30:12 ID:WVF7YRda.net] >>925 まぁ調べようないね そもそもモデル曲線のパラメータをボコボコ増やしていけばそりゃいくらでも表現力 986 名前：は高まるんだからドンドン実測値に合わせていける でも神ならぬ人間にはそんな実測不能のパラメータだらけの曲線を持ってきてもしょうがない 結局モデル曲線の母数の推定に使えるモニタリング指標なんて限りがあるし、その少ない情報の中で可能な限りせめて全体の傾向だけでも捉える事ができないかと考えられたのがsirなんだから 「こういう調査すれば××がわかってより正確な予測ができる」なんてものはいくらでもあるだろうけどないものねだりしたってしょうがない まぁsirは病理的パラメータとして回復率しか見ないけど実際には病理的に決定できるパラメータがもっとたくさんあるけとな しかし日本人の実際の生活において自然免疫だけで撃退できる比率はムズイんだろうと思う もちろん研究してる人もいるだろうけど しかしどれくらい被曝してどれくらい発症する人がいるのかなんて人体実験でもしないと出ないだろうし [] [ここ壊れてます] 987 名前：１３２人目の素数さん [2020/07/28(火) 12:03:34 ID:/xPjTU2o.net] >>926 調べようはあるんじゃないの？ 濃厚接触者の罹患率は簡単に調べられるでしょ。 俺たち部外者には調べようがないだけのこと。 988 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 12:17:39 ID:WVF7YRda.net] >>927 まぁ当局の人は外野が知らない指標はいっぱい持ってるやろな しかしそれでsirでは感染性宿主に繰り入れられてる人数のうち、実際にはワクチン接触してるレベルで自然免疫だけで撃退できる割合なんて母数を推定するのはむりな気はする 家族とか職場での接触の度合いなんて個人差ありすぎやろ？ 989 名前：１３２人目の素数さん [2020/07/28(火) 14:55:30 ID:/xPjTU2o.net] >>928 >個人差ありすぎやろ？ それでも意味があるのが統計調査。 990 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 16:00:15.94 ID:3tKaz31X.net] >>928 意味ない それに本当に意味があるのか、誤差と言っていい程度の誤差なのかは数学では出ない そこから先は得られた結論が実際にどれくらいズレてるのかフィールドワークで疫学の専門家たちが実地に調べてきたのが現在の防疫学 素人が紙の上でヒョコヒョコ計算して「こんなんどう？」と言ってる時点では仮説にしか過ぎない それを分からずに大阪に持ち込んで大顰蹙買ってるのが阪大のアホ物理学者だよ 991 名前：１３２人目の素数さん [2020/07/28(火) 16:20:28 ID:/xPjTU2o.net] そういえば、googleがスマホの位置情報から、人の移動量を調べることが できるって言ってるらしいけど、西浦説の裏付けには使ったんだろうか？ 人の移動を8割削減しないと感染爆発とか言ってたやつね。 992 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 18:18:01.96 ID:AUVo4ZRB.net] >>931 日本の移動量は欧米のロックダウン相当の減少と聞いたな。 993 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 19:15:00.14 ID:koi0QwxX.net] >>931 移動じゃない。接触量だ。 移動が半分になれば接触量は1／4になるから、ほぼ達成。 994 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 19:26:56.02 ID:Tev5pFvA.net] >>931 テレビではヤフーから提供されたデータを使って 緊急事態宣言前の日の人出がどの程度になっていたかを数値やグラフで表示してた 渋谷とかの繁華街だけだけど 995 名前：１３２人目の素数さん mailto:sage [2020/07/28(火) 21:12:33.53 ID:MZcVolQH.net] >>933 当初は西浦さんも移動量（人出）8割減つってたよね。 接触量減が二乗に比例するんじゃないかって話について当人はのらりくらり。その後のRの下がり方みると、二乗に比例すると考え外出5割減→接触25%≒8割減でよかったらしい。 普通にわからんならわからんつってほしかったわ。 996 名前：１３２人目の素数さん mailto:sage [2020/07/29(水) 05:04:52 ID:kvGiW0JL.net] >>933 接触量っていっても、条件によって重みづけしないと 換気のいい所で双方マスクした状態での接触と 997 名前：、 換気の悪い室内でノーマスクで10人でカラオケや飲み会したのでは、 感染確率がまったく違うからな 接触時間、接触距離の他に、 濃厚接触係数？みたいな数値が必要 マスクの有無、換気の有無、声出しの有無等で係数を変えて重みづけしないと [] [ここ壊れてます] 998 名前：１３２人目の素数さん [2020/07/30(木) 00:14:52 ID:dea9G3z1.net] 東京はピークアウトしそうじゃないか？ 地方の大都市部はまだ増加中みたいだけど。 999 名前：１３２人目の素数さん [2020/07/30(木) 00:48:17.91 ID:ustZkzWN.net] >>916 式も書いておく ・実効再生産数 対策なしでβとγは定数とする。 σ≡β/γとして S(t)=so exp[-σ R(t)]≦soより Rt=σ S(t)≦σso= Ro Ro→小なら Rt→小 ・自然治癒=ワクチン投与の場合 　t=0の時　 　感染者が極めて少ない 　0<lo<<1 とする。 基本再生産数Ro≡σso≒ σN またN'=N-ro≒so 　R'o/ Ro=σN'/σN=(N-ro)/N =1-ro/N≡1-p 流行しないためには 　つねにR't<R'oよりR'o<1 　(1-p)Ro<1 1-1/ Ro<p ・Cov 19のRoをSARSの3程度 　　と仮定すれば 　　2/3<p 初期にワクチン投与者(自然治癒者) 　が全人口の67パーセント以上なら 　流行しない 1000 名前：１３２人目の素数さん [2020/07/30(木) 01:06:31.43 ID:ustZkzWN.net] 大阪府はK値のおつげより7月9日ごろに ピークアウトしたはずだ。 政府大阪愛知共々GO TO キャンペーン で全国へバカンス 新型肺炎ウィルス共々に 1001 名前：１３２人目の素数さん mailto:sage [2020/07/30(木) 04:58:08.17 ID:mCTKNYbJ.net] >>937 今日の数値を見てからだな 土日のに出来なかった分の検査結果が出るのが木曜で一気に増えるから平均を押し上げる 1002 名前：１３２人目の素数さん mailto:sage [2020/07/30(木) 05:07:36.95 ID:bjYtrpOr.net] 今日は検査数も多かったから多いらしい。 1003 名前：１３２人目の素数さん [2020/07/30(木) 07:26:28.01 ID:+O9O8n76.net] 大阪はPCR検査が5日待ちと聞いたが、 つい最近まで検査数自体が異様に 少ない、医者が死亡した、K値、 闇につつまれている 1004 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 11:14:04 ID:sNNGYlEF.net] 【コロナ抗体検査】 宇都宮市の感染者数は推定６３７８人、市が把握していた陽性者の２７７倍 [影のたけし軍団★] https://asahi.5ch.net/test/read.cgi/newsplus/1596158248/ >> 調査は、地域の人口ベースでの感染状況を把握し、感染拡大防止策に役立てるのが狙い。 市の協力を得て、市内１千世帯（２２９０人）を無作為で抽出。 ６〜７月に検査したところ、参加者７４２人のうち感染歴を示す抗体の陽性者は３人だった。 この時点の感染率は０・40％だが、統計学的に精度の高い数字に補正すると１・23％になるとしている。 << >統計学的に精度の高い数字に補正 ってなんだろうね？ 各種信頼区間を求めると > binom.ci(3,742) method x n mean lower upper 1 agresti-coull 3 742 0.004043126685 0.0007875167910 0.012407591505 2 asymptotic 3 742 0.004043126685 -0.0005227594808 0.008609012850 3 bayes 3 742 0.004710632571 0.0006700883674 0.009654411476 4 cloglog 3 742 0.004043126685 0.00115452 1005 名前：03237 0.011224715356 5 exact 3 742 0.004043126685 0.0008345672857 0.011769890197 6 logit 3 742 0.004043126685 0.0013045836463 0.012458623458 7 probit 3 742 0.004043126685 0.0012220936604 0.011703448731 8 profile 3 742 0.004043126685 0.0010103960338 0.010450783974 9 lrt 3 742 0.004043126685 0.0010330128348 0.010444714948 10 prop.test 3 742 0.004043126685 0.0010447711937 0.012807577663 11 wilson 3 742 0.004043126685 0.0013759603328 0.011819147963 > [] [ここ壊れてます] 1006 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 11:25:53.62 ID:uMo8LokM.net] 参加者人数742人は2290人からp=742/2290の二項分布に従った数値と設定するのかな？ 1007 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 15:18:45.72 ID:ImynyFWQ.net] >>944 やってみたけど、全く無意味だったｗ https://i.imgur.com/2b3l4zP.png 1008 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 15:19:06.30 ID:ImynyFWQ.net] この統計処理云々には疑問があるな。 偽陽性が多い検査らしいから 検査感度の分布を平均0.9標準偏差0.1、特異度を平均0.8　標準偏差0.1に設定して 742人中3人陽性でパラメータの分布をMCMCでだすと こんな感じになるのでこの結論とは全く逆だな。 https://i.imgur.com/KVwSwwa.png 1009 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 16:36:10 ID:I5zVtTcu.net] >>943 無作為抽出したけど、協力してくれた人の年齢・職業分布に偏りがあって、それを補正したんじゃない？ 世論調査でよくやるやつ。 1010 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 21:52:50.14 ID:xsqMU8gU.net] PCR検査の発明者は生前、「これは感染症検査には使うな」と言っていた そして昨年亡くなってしまった 彼が健在なら動きは違ったものになっただろう メディアが報道しないPCR検査の秘密 https://pbs.twimg.com/media/EeLojPZVoAEjCVA.jpg 1011 名前：１３２人目の素数さん mailto:sage [2020/07/31(金) 23:22:42.85 ID:ElpU0iaW.net] >>948 アホはツイッターに帰れ 1012 名前：１３２人目の素数さん [2020/08/01(土) 00:41:16.20 ID:DbvvNTta.net] >>943 これまでの抗体検査と同じで、意味無し推定だな。 1013 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 03:31:20 ID:ksq2JqO+.net] 宇都宮は抗体検査の陽性者すくなすぎて、もっと検査人数を拡大しないといまいちわからんな また抗体検査は、独立して開発・製造された2社の抗体検査キット両方使って判定しないと、 いまいち信頼できない A者の抗体検査キットでもB社の抗体検査キットでも陽性 A者の抗体検査キットでもB社の抗体検査キットでも陰性 A者の抗体検査キットとB社の抗体検査キットで結果が違う これらを考慮しないと 1014 名前：１３２人目の素数さん [2020/08/01(土) 10:29:25 ID:MOVlgcIX.net] >>943 なぜ二項分布ですか？ 1015 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 10:47:24 ID:85irbUct.net] >>952 別人だけど非復元抽出で陽性か陰性の2通りだからだと思う 1016 名前：１３２人目の素数さん [2020/08/01(土) 10:49:24 ID:kSBhb5Su.net] どうみても弱毒化しているように見えますので、 俺は弱毒化（致死率の低下）についてモデルを 考えてみます。 つか、コロナウィルスってどんどん変化するらしいので、 すぐに性質変わるんですよね？ ウィルスはド素人なので全くわかりませんけど・・・。 1017 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 11:07:59.20 ID:85irbUct.net] >>954 患者の年齢層とか検査対象範囲の拡大とかウイルスそのものの要因以外に起因する変化という見解もあるみたいですよ 1018 名前：１３２人目の素数さん [2020/08/01(土) 11:21:14 ID:DbvvNTta.net] 年齢層が特にね。 ４月の感染拡大と比べると明らかに高齢者の割合が減ってるから、 そのせいじゃね？＞致死率の低下 たとえば、７０代の感染者数、死者数の推移で比べないと弱毒化 とは判断できん。 1019 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 11:21:39 ID:6WqfoPjj.net] >>952 1の目の出る確率が1/6のサイコロを100回振って1の目のでる回数の分布は二項分布じゃないかな？ ポアソン分布の方がいい？？ 1020 名前：１３２人目の素数さん [2020/08/01(土) 11:28:19 ID:kSBhb5Su.net] >>955 >>956 レスありがとうございます。 同じ属性（年齢）で比較できるか データを探してみます。 1021 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 12:00:13 ID:HQxy6OrC.net] >>954 確かに、第2波の方が感染者数が多くても死者が少ない国が多いんだよなあ。日本もそんな感じ。 ところがオーストラリアでは感染者数に比例して、あるいはそれ以上に死んでる。 https://www.worldometers.info/coronavirus/country/australia/ これ、感染力に気温の影響は無くても、重症化率には効いてくるという話なんじゃね？ 1022 名前：１３２人目の素数さん [2020/08/01(土) 14:25:09 ID:z4AwLSVw.net] >>957 珍しい事象はポアソン分布だね。 数セミ5月号で最尤度推定は2項分布ですが。 宇都宮の件は山中教授同様に エクセル公式へぶっこんだのか？ 統計推定方法がわからない。 1023 名前：１３２人目の素数さん [2020/08/01(土) 14:27:09 ID:XiLTKLm0.net] >>960 あと抗体検査の系統誤差は 計算しているのかな？ 1024 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 15:34:53.02 ID:akw9FJQf.net] >>960 上限がないのがポアソンであるのが二項分布では？ 1025 名前：１３２人目の素数さん [2020/08/01(土) 15:47:14.76 ID:TPZxQKbP.net] >>962 二項分布の根拠は？ 1026 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 16:00:50.79 ID:85irbUct.net] ポアソン分布は珍しい事がある時間内に何回発生するかを確率変数とした分布だと思う 確か馬に蹴られて死ぬ兵士の人数を考えていて発見した分布 問題にしている抗体検査の場合の確率変数は調査時点以前に感染したことがある人数 非復元抽出だと厳密には超幾何分布かもしれないけど 全体の個数が多ければ復元しなくても復元しないことの影響が小さくなって二項分布で近似できるんじゃないかな 1027 名前：１３２人目の素数さん mailto:sage [2020/08/01(土) 16:14:33.60 ID:85irbUct.net] あと パラメータが n と p = λ/nである二項分布において、λ を一定に保ったまま n を無限大に近づけると、その分布は平均 λ のポアソン分布に近づく。 時間の場合は時間間隔は無限に小さくなるからnが無限大になるかも 空間の場合も量子レベルとかのほんの少しの違いを区別するなら無限大とみなして良いかもしれない 人数の場合はnが無限大と言えるくらい大きければポアソン分布になるかも 1028 名前：１３２人目の素数さん [2020/08/01(土) 18:33:08.31 ID:voACzb9W.net] まず感染現象と類推するにしても 古典的な場合で、全く量子や量子統計 は関係ない。 中心極限定理からは近似が ポアソン分布→二項分布→正規分布 が考えられる。 正規分布は(一様等方な自由粒子の 衝突モデルなど)均等な接触による 感染を表し感染者は均等に感染させる。 βが集団平均のパラメーターだから 正規分布の感染は蔓延した市中感染も 表す。 一方、クラスターで1人が多数に感染させたり感染者が感染させないなど不均等 場合、分布を無理やり平均場で近似させるだけでは不十分だ。 感染症へ素人だが、両方の場合を分けて分析するな。 K値関係は論外として 1029 名前：１３２人目の素数さん mailto:sage [2020/08/02(日) 06:26:41 ID:5m0cS9Rk.net] PCR検査を調べた人間はわかるが インフルにも反応する クラミジアにも反応する もともとコロナウイルス専用の検査ではないんだよ これが感染爆発の秘密だ https://i.imgur.com/6J0qvev.jpg 1030 名前：１３２人目の素数さん mailto:sage [2020/08/02(日) 07:03:58.68 ID:koevBDz6.net] >>967 プライマーに共通部分がなければ誤� 1031 名前：沛oは理論的にありえん。 [] [ここ壊れてます] 1032 名前：１３２人目の素数さん mailto:sage [2020/08/02(日) 10:55:07 ID:ZtgCF6NZ.net] >>967 PCRは遺伝子指定して使う検査方法 だからいろんなウイルスの検出に使用できるんよ それを勘違いした人が一定数いるらしい 新型コロナをターゲットにしたら新型コロナしか検出せんよ 1033 名前：１３２人目の素数さん mailto:sage [2020/08/02(日) 16:34:11.82 ID:jZXxVIQr.net] >>969 抗体検査とごっちゃにしてる人居るよね 1034 名前：１３２人目の素数さん mailto:sage [2020/08/02(日) 17:57:32 ID:/N/GqpmS.net] インフルとクラミジアウイルス(?)の遺伝子配列と相同なプライマー使ってたら検出されるだろうけど、ちゃんと配列は選んでるでしょ。 ミスマッチ起こったら仕方ないけど 1035 名前：１３２人目の素数さん [2020/08/03(月) 09:09:21.10 ID:nP+Sm0D1.net] 東京はそろそろピークアウトしそうな感じ？ 連日ワイドショーで騒いでくれるおかげで、実効再生産数が 漸減してほぼ１になってるみたい。 「新しい日常」に馴染めず旧態依然の生活を続けようとする 馬鹿者たちが感染を拡げてるが、ようやくそれじゃだめだと 気づくころ？ 1036 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 10:19:10 ID:7Vi8/EtE.net] >>972 まだ無理やろ 再生産数計算してるサイト見たらまだ1.3ほどあるみたいやし https://rt-live-japan.com/ 1037 名前：１３２人目の素数さん [2020/08/03(月) 11:05:13.14 ID:eJM6VgsU.net] >>973 そのサイトは計算の根拠がわかりません。 東洋経済では　8月1日時点 全国　1.43 東京都0.99 大阪1.64 愛知1.93 福岡　1.3 沖縄　4.1 東京都の実効再生産数は全都の指数 です、 更にクラスターの地域や職種などの 情報も必要でしょう 東洋経済 「実効再生産数= 直近7日間の新規陽性者数／その前 7日間の新規陽性者数）^（平均世代 時間／報告間隔）」 平均世代時間は5日、報告間隔は7日 https://toyokeizai.net/sp/visual/tko/covid19/ 1038 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 11:11:53 ID:C8O7DWVU.net] そもそも「今の」再生産数が出てる段階でおかしいわな。どちらも。 1039 名前：１３２人目の素数さん [2020/08/03(月) 11:28:42.92 ID:eJM6VgsU.net] 速報値だからな 1040 名前：１３２人目の素数さん [2020/08/03(月) 11:59:22.14 ID:fIfeS9d2.net] 学術の巨大掲示板群 - アルファ・ラボ ttp://x0000.net 数学 物理学 化学 生物学 天文学 地理地学 IT 電子 工学 言語学 国語 方言 など 1041 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 12:27:13 ID:xgcgj0Hi.net] >>974 一番下に参照したサイトとか書いてあるよ 1042 名前：１３２人目の素数さん [2020/08/03(月) 12:29:56.16 ID:eJM6VgsU.net] >>978 説明して 1043 名前：１３２人目の素数さん [2020/08/03(月) 13:34:21.90 ID:nP+Sm0D1.net] >>975 速報値だから別に構わん。 感染から陽性判明までの時間差があるから過去の値だし、 事例ごとに時間差を補正しないと正確ではないことも みなさん承知でしょ。 1044 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 13:49:14.05 ID:7Vi8/EtE.net] >>974 実効再生産数の推定方法は使うモデルでもちがうし一概になんとも言えない 例えばsirにしても実効再生産数は感染率と回復率の両方で日々の拡大の比率が決まる事を利用するので、回復率の概算値の取り方でもズレてくる その東洋経済の計算式もかなりデータ処理する人によって恣意的に変わりうるパラメータを含んでいるのでなんとも言えないけど、監修してるのが西浦先生なら割とどちらにもバイアスのかかってない式だと言えるかもしれない でもどのみち再生産数は「投げたボールが自然に落ちてくる」ような法則で決まるものではなく、市民の集団心理のような数学だけではでは決まらないものなのでハッキリとした事が言えるような性格のものではないとしか言えないと思う 少なくとも今0.99になったのがなんらかの対策の結果によって全体の傾向として下がり始めているのか、たまたま、ココ数日の数値データを見て市民の意識が警戒する方にふれているのかはまだなんとも言えないと思う 前者ならこのまま1未満が持続されるだろうけど、後者で不急ではあるが不要でない用事を後回しにしてるだけなら、また増えてくるかもしれない こんなの多分数週間単位で数字見ないとわかりっこないと思う 1045 名前：１３２人目の素数さん [2020/08/03(月) 16:57:55.99 ID:nP+Sm0D1.net] >>981 だんだん小さくなっていくトレンドははっきり見て取れるんじゃない？ で、実際、感染者の増加傾向の倍加時間も長くなってるように見える。 ま、連日あれだけ騒がれてりゃ、都民も飲み会自粛したいりして気を つけるだろうからね。 逆に言えば、７月に入ってバーっと感染拡大したのは、規制解除で浮か れて油断したからでしょ。 1046 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 17:14:26.84 ID:7Vi8/EtE.net] >>982 だからそのトレンドが方程式上ないんだよ 今日少ないから明日も少ないとは方程式の形からは言えない あくまでβ(t)が一定値を取る傾向にあると言うのは数学上の法則ではなく経験則に過ぎない まぁそう考えても過去の経験上、数週間〜数ヶ月単位のトレンドを見る限りにおいては大きくはハズレないと言うだけ β(t)がこのところ下がってるから明日も下がるかなんてとても言えない 実際今日東京かなり多いみたいだし 再生産数のトレンドはやはり数週間単位でないとわからないと思う しかも今は7月の4連休の影響と梅雨明けの影響をダブルで受けてるから8/10くらいまではかなり値が揺れると思う 1047 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 17:40:14.51 ID:C8O7DWVU.net] >>980 いや、だから事例ごとに補正した値を出すべきなんじゃね？遅れも考慮して。手法は確立してるし。 もちろんある程度信頼できる値は1、２週間前のものしか求まらないから、１週間前のRtは1.3でしたと発表するしかないが。 1048 名前：１３２人目の素数さん [2020/08/03(月) 19:15:01.96 ID:qIXOsgp1.net] >>984 モデルから説明してくれ 1049 名前：１３２人目の素数さん mailto:sage [2020/08/03(月) 21:24:43 ID:rU5XKbyf.net] >>985 色々あるけど西浦のやり方で良いよ。 1050 名前：１３２人目の素数さん [2020/08/04(火) 00:10:14.54 ID:QOYlLXf+.net] >>983 >今日少ないから明日も少ないとは方程式の形からは言えない んなこたわかってるよ。予測ができりゃ誰も苦労しないわけで。 単純に長期トレンドを外挿するのが関の山だが、それで大概は あたりだったりする。 1051 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 00:21:42.31 ID:buqp/rtg.net] >>967 そんなデマ信じるなよ 現在厚労省が認可してるPCR検査キットは、新型コロナウイルス以外に反応した例は確認されてない 1052 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 00:33:47.31 ID:OPfgl3ou.net] >>987 あてもんしたいなら好きにすれば良い ココの板の住人はそんなあてもん興味ない K値中野とか宮沢は好きそうだけど 1053 名前：１３２人目の素数さん [2020/08/04(火) 00:46:28.14 ID:QOYlLXf+.net] 俺も「ココの板の住人」だが、あてもんも必要だよ。 先のことはわからない、で済ませちゃなにも面白くない。 過去を分析するのは、未来を予想するため。 1054 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 00:57:17.62 ID:OPfgl3ou.net] >>990 なんの根拠もない数字が当たってなんでそれで未来が予測できるんだよ？ それこそまさにK値理論で中野がやらかしてる事やろ？ 1055 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 09:15:15 ID:kB52B1VU.net] >>715 メキシコのデング熱ワクチンの話が業界では有名。 1056 名前：１３２人目の素数さん [2020/08/04(火) 10:46:52.87 ] [ここ壊れてます] 1057 名前： ID:QOYlLXf+.net mailto: >>991 大きな変動要素が見当たらなければ、単なる外挿でだいたいあうはず。 夏休み、お盆休みの影響が出るまではこの調子じゃね？ [] [ここ壊れてます] 1058 名前：１３２人目の素数さん [2020/08/04(火) 10:57:46.29 ID:NqlKcFnC.net] 次スレ 数学 統計に詳しい人が語るコロナウイルス ☆2 https://rio2016.5ch.net/test/read.cgi/math/1596506253/ 1059 名前：１３２人目の素数さん [2020/08/04(火) 11:44:49.84 ID:EnxLnyf8.net] >>993 株価の予言でがんばれ 1060 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 13:18:34.30 ID:UTSYjGjo.net] >>993 だから合わせても意味ないっての 中野はそれがわかってないんだって あるモデルを作ってそれが現実の値にあってるかどうかなんて数学的にはわからない それがちゃんと予言性を持つモデルか、設定されてる仮説が妥当かどうかはモデルから出された預言値と実測値が合うかどうかを何度も何度も検証して実証される sirとかはいろんな民族のいろんな社会のいろんな疫病で検証されて「β(t)が数週間くらいのオーダーでは一定」という仮説による予言性がわりと信頼ができると検証されてきた β(t)が変化するようなモデルで信頼にたる予言性のあるモデルなんかないやろ？ 1061 名前：１３２人目の素数さん [2020/08/04(火) 13:55:10 ID:QOYlLXf+.net] 意味はないだろうね。単なる外挿だから、天気の長期予報みたいなもの。 どっちみちβは予測できないんだから、モデルを細かく云々しても 予想に限れば意味ないのは同じ。 株価のほうは市場にもろにフィードバックされるから、たとえ正しく 予想できても公言できない罠ｗ 1062 名前：１３２人目の素数さん [2020/08/04(火) 13:58:57 ID:QOYlLXf+.net] >>996 >β(t)が数週間くらいのオーダーでは一定 β(t)のトレンドが数週間くらいのオーダーでは一定としてもええやんか。 1063 名前：１３２人目の素数さん mailto:sage [2020/08/04(火) 14:22:44.10 ID:UTSYjGjo.net] >>398 それが一定としていい実証がないっての もちろん一定として予測だしたけりゃ出せばいい 誰も支持しないだけ 1064 名前：１３２人目の素数さん [2020/08/04(火) 14:59:06 ID:NqlKcFnC.net] 次スレ 数学 統計に詳しい人が語るコロナウイルス ☆2 https://rio2016.5ch.net/test/read.cgi/math/1596506253/ 1065 名前：1001 [Over 1000 Thread .net] このスレッドは１０００を超えました。 新しいスレッドを立ててください。 life time: 157日 12時間 40分 26秒 1066 名前：過去ログ ★ [[過去ログ]] ■ このスレッドは過去ログ倉庫に格納されています

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef