20/07/28 14:55:30 /xPjTU2o.net
>>928
>個人差ありすぎやろ？
それでも意味があるのが統計調査。
990:１３２人目の素数さん
20/07/28 16:00:15.94 3tKaz31X.net
>>928
意味ない
それに本当に意味があるのか、誤差と言っていい程度の誤差なのかは数学では出ない
そこから先は得られた結論が実際にどれくらいズレてるのかフィールドワークで疫学の専門家たちが実地に調べてきたのが現在の防疫学
素人が紙の上でヒョコヒョコ計算して「こんなんどう？」と言ってる時点では仮説にしか過ぎない
それを分からずに大阪に持ち込んで大顰蹙買ってるのが阪大のアホ物理学者だよ
991:１３２人目の素数さん
20/07/28 16:20:28 /xPjTU2o.net
そういえば、googleがスマホの位置情報から、人の移動量を調べることが
できるって言ってるらしいけど、西浦説の裏付けには使ったんだろうか？
人の移動を8割削減しないと感染爆発とか言ってたやつね。
992:１３２人目の素数さん
20/07/28 18:18:01.96 AUVo4ZRB.net
>>931
日本の移動量は欧米のロックダウン相当の減少と聞いたな。
993:１３２人目の素数さん
20/07/28 19:15:00.14 koi0QwxX.net
>>931
移動じゃない。接触量だ。
移動が半分になれば接触量は1／4になるから、ほぼ達成。
994:１３２人目の素数さん
20/07/28 19:26:56.02 Tev5pFvA.net
>>931
テレビではヤフーから提供されたデータを使って
緊急事態宣言前の日の人出がどの程度になっていたかを数値やグラフで表示してた
渋谷とかの繁華街だけだけど
995:１３２人目の素数さん
20/07/28 21:12:33.53 MZcVolQH.net
>>933
当初は西浦さんも移動量（人出）8割減つってたよね。
接触量減が二乗に比例するんじゃないかって話について当人はのらりくらり。その後のRの下がり方みると、二乗に比例すると考え外出5割減→接触25%≒8割減でよかったらしい。
普通にわからんならわからんつってほしかったわ。
996:１３２人目の素数さん
20/07/29 05:04:52 kvGiW0JL.net
>>933
接触量っていっても、条件によって重みづけしないと
換気のいい所で双方マスクした状態での接触と
997:、 換気の悪い室内でノーマスクで10人でカラオケや飲み会したのでは、 感染確率がまったく違うからな 接触時間、接触距離の他に、 濃厚接触係数？みたいな数値が必要 マスクの有無、換気の有無、声出しの有無等で係数を変えて重みづけしないと
998:１３２人目の素数さん
20/07/30 00:14:52 dea9G3z1.net
東京はピークアウトしそうじゃないか？
地方の大都市部はまだ増加中みたいだけど。
999:１３２人目の素数さん
20/07/30 00:48:17.91 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:１３２人目の素数さん
20/07/30 01:06:31.43 ustZkzWN.net
大阪府はK値のおつげより7月9日ごろに
ピークアウトしたはずだ。
政府大阪愛知共々GO TO キャンペーン
で全国へバカンス
新型肺炎ウィルス共々に
1001:１３２人目の素数さん
20/07/30 04:58:08.17 mCTKNYbJ.net
>>937
今日の数値を見てからだな
土日のに出来なかった分の検査結果が出るのが木曜で一気に増えるから平均を押し上げる
1002:１３２人目の素数さん
20/07/30 05:07:36.95 bjYtrpOr.net
今日は検査数も多かったから多いらしい。
1003:１３２人目の素数さん
20/07/30 07:26:28.01 +O9O8n76.net
大阪はPCR検査が5日待ちと聞いたが、
つい最近まで検査数自体が異様に
少ない、医者が死亡した、K値、
闇につつまれている
1004:１３２人目の素数さん
20/07/31 11:14:04 sNNGYlEF.net
【コロナ抗体検査】 宇都宮市の感染者数は推定６３７８人、市が把握していた陽性者の２７７倍 [影のたけし軍団★]
ｽﾚﾘﾝｸ(newsplus板)
>>
調査は、地域の人口ベースでの感染状況を把握し、感染拡大防止策に役立てるのが狙い。
市の協力を得て、市内１千世帯（２２９０人）を無作為で抽出。
６～７月に検査したところ、参加者７４２人のうち感染歴を示す抗体の陽性者は３人だった。
この時点の感染率は０・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:１３２人目の素数さん
20/07/31 11:25:53.62 uMo8LokM.net
参加者人数742人は2290人からp=742/2290の二項分布に従った数値と設定するのかな？
1007:１３２人目の素数さん
20/07/31 15:18:45.72 ImynyFWQ.net
>>944
やってみたけど、全く無意味だったｗ
URLﾘﾝｸ(i.imgur.com)
1008:１３２人目の素数さん
20/07/31 15:19:06.30 ImynyFWQ.net
この統計処理云々には疑問があるな。
偽陽性が多い検査らしいから
検査感度の分布を平均0.9標準偏差0.1、特異度を平均0.8 標準偏差0.1に設定して
742人中3人陽性でパラメータの分布をMCMCでだすと
こんな感じになるのでこの結論とは全く逆だな。
URLﾘﾝｸ(i.imgur.com)
1009:１３２人目の素数さん
20/07/31 16:36:10 I5zVtTcu.net
>>943
無作為抽出したけど、協力してくれた人の年齢・職業分布に偏りがあって、それを補正したんじゃない？
世論調査でよくやるやつ。
1010:１３２人目の素数さん
20/07/31 21:52:50.14 xsqMU8gU.net
PCR検査の発明者は生前、「これは感染症検査には使うな」と言っていた
そして昨年亡くなってしまった
彼が健在なら動きは違ったものになっただろう
メディアが報道しないPCR検査の秘密
URLﾘﾝｸ(pbs.twimg.com)
1011:１３２人目の素数さん
20/07/31 23:22:42.85 ElpU0iaW.net
>>948
アホはツイッターに帰れ
1012:１３２人目の素数さん
20/08/01 00:41:16.20 DbvvNTta.net
>>943
これまでの抗体検査と同じで、意味無し推定だな。
1013:１３２人目の素数さん
20/08/01 03:31:20 ksq2JqO+.net
宇都宮は抗体検査の陽性者すくなすぎて、もっと検査人数を拡大しないといまいちわからんな
また抗体検査は、独立して開発・製造された2社の抗体検査キット両方使って判定しないと、
いまいち信頼できない
A者の抗体検査キットでもB社の抗体検査キットでも陽性
A者の抗体検査キットでもB社の抗体検査キットでも陰性
A者の抗体検査キットとB社の抗体検査キットで結果が違う
これらを考慮しないと
1014:１３２人目の素数さん
20/08/01 10:29:25 MOVlgcIX.net
>>943
なぜ二項分布ですか？
1015:１３２人目の素数さん
20/08/01 10:47:24 85irbUct.net
>>952
別人だけど非復元抽出で陽性か陰性の2通りだからだと思う
1016:１３２人目の素数さん
20/08/01 10:49:24 kSBhb5Su.net
どうみても弱毒化しているように見えますので、
俺は弱毒化（致死率の低下）についてモデルを
考えてみます。
つか、コロナウィルスってどんどん変化するらしいので、
すぐに性質変わるんですよね？
ウィルスはド素人なので全くわかりませんけど・・・。
1017:１３２人目の素数さん
20/08/01 11:07:59.20 85irbUct.net
>>954
患者の年齢層とか検査対象範囲の拡大とかウイルスそのものの要因以外に起因する変化という見解もあるみたいですよ
1018:１３２人目の素数さん
20/08/01 11:21:14 DbvvNTta.net
年齢層が特にね。
４月の感染拡大と比べると明らかに高齢者の割合が減ってるから、
そのせいじゃね？＞致死率の低下
たとえば、７０代の感染者数、死者数の推移で比べないと弱毒化
とは判断できん。
1019:１３２人目の素数さん
20/08/01 11:21:39 6WqfoPjj.net
>>952
1の目の出る確率が1/6のサイコロを100回振って1の目のでる回数の分布は二項分布じゃないかな？
ポアソン分布の方がいい？？
1020:１３２人目の素数さん
20/08/01 11:28:19 kSBhb5Su.net
>>955
>>956
レスありがとうございます。
同じ属性（年齢）で比較できるか
データを探してみます。
1021:１３２人目の素数さん
20/08/01 12:00:13 HQxy6OrC.net
>>954
確かに、第2波の方が感染者数が多くても死者が少ない国が多いんだよなあ。日本もそんな感じ。
ところがオーストラリアでは感染者数に比例して、あるいはそれ以上に死んでる。
URLﾘﾝｸ(www.worldometers.info)
これ、感染力に気温の影響は無くても、重症化率には効いてくるという話なんじゃね？
1022:１３２人目の素数さん
20/08/01 14:25:09 z4AwLSVw.net
>>957
珍しい事象はポアソン分布だね。
数セミ5月号で最尤度推定は2項分布ですが。
宇都宮の件は山中教授同様に
エクセル公式へぶっこんだのか？
統計推定方法がわからない。
1023:１３２人目の素数さん
20/08/01 14:27:09 XiLTKLm0.net
>>960
あと抗体検査の系統誤差は
計算しているのかな？
1024:１３２人目の素数さん
20/08/01 15:34:53.02 akw9FJQf.net
>>960
上限がないのがポアソンであるのが二項分布では？
1025:１３２人目の素数さん
20/08/01 15:47:14.76 TPZxQKbP.net
>>962
二項分布の根拠は？
1026:１３２人目の素数さん
20/08/01 16:00:50.79 85irbUct.net
ポアソン分布は珍しい事がある時間内に何回発生するかを確率変数とした分布だと思う
確か馬に蹴られて死ぬ兵士の人数を考えていて発見した分布
問題にしている抗体検査の場合の確率変数は調査時点以前に感染したことがある人数
非復元抽出だと厳密には超幾何分布かもしれないけど
全体の個数が多ければ復元しなくても復元しないことの影響が小さくなって二項分布で近似できるんじゃないかな
1027:１３２人目の素数さん
20/08/01 16:14:33.60 85irbUct.net
あと
パラメータが n と p = λ/nである二項分布において、λ を一定に保ったまま n を無限大に近づけると、その分布は平均 λ のポアソン分布に近づく。
時間の場合は時間間隔は無限に小さくなるからnが無限大になるかも
空間の場合も量子レベルとかのほんの少しの違いを区別するなら無限大とみなして良いかもしれない
人数の場合はnが無限大と言えるくらい大きければポアソン分布になるかも
1028:１３２人目の素数さん
20/08/01 18:33:08.31 voACzb9W.net
まず感染現象と類推するにしても
古典的な場合で、全く量子や量子統計
は関係ない。
中心極限定理からは近似が
ポアソン分布→二項分布→正規分布
が考えられる。
正規分布は(一様等方な自由粒子の
衝突モデルなど)均等な接触による
感染を表し感染者は均等に感染させる。
βが集団平均のパラメーターだから
正規分布の感染は蔓延した市中感染も
表す。
一方、クラスターで1人が多数に感染させたり感染者が感染させないなど不均等
場合、分布を無理やり平均場で近似させるだけでは不十分だ。
感染症へ素人だが、両方の場合を分けて分析するな。
K値関係は論外として
1029:１３２人目の素数さん
20/08/02 06:26:41 5m0cS9Rk.net
PCR検査を調べた人間はわかるが
インフルにも反応する
クラミジアにも反応する
もともとコロナウイルス専用の検査ではないんだよ
これが感染爆発の秘密だ
URLﾘﾝｸ(i.imgur.com)
1030:１３２人目の素数さん
20/08/02 07:03:58.68 koevBDz6.net
>>967
プライマーに共通部分がなければ誤�
1031:沛oは理論的にありえん。
1032:１３２人目の素数さん
20/08/02 10:55:07 ZtgCF6NZ.net
>>967
PCRは遺伝子指定して使う検査方法
だからいろんなウイルスの検出に使用できるんよ
それを勘違いした人が一定数いるらしい
新型コロナをターゲットにしたら新型コロナしか検出せんよ
1033:１３２人目の素数さん
20/08/02 16:34:11.82 jZXxVIQr.net
>>969
抗体検査とごっちゃにしてる人居るよね
1034:１３２人目の素数さん
20/08/02 17:57:32 /N/GqpmS.net
インフルとクラミジアウイルス(?)の遺伝子配列と相同なプライマー使ってたら検出されるだろうけど、ちゃんと配列は選んでるでしょ。
ミスマッチ起こったら仕方ないけど
1035:１３２人目の素数さん
20/08/03 09:09:21.10 nP+Sm0D1.net
東京はそろそろピークアウトしそうな感じ？
連日ワイドショーで騒いでくれるおかげで、実効再生産数が
漸減してほぼ１になってるみたい。
「新しい日常」に馴染めず旧態依然の生活を続けようとする
馬鹿者たちが感染を拡げてるが、ようやくそれじゃだめだと
気づくころ？
1036:１３２人目の素数さん
20/08/03 10:19:10 7Vi8/EtE.net
>>972
まだ無理やろ
再生産数計算してるサイト見たらまだ1.3ほどあるみたいやし
URLﾘﾝｸ(rt-live-japan.com)
1037:１３２人目の素数さん
20/08/03 11:05:13.14 eJM6VgsU.net
>>973
そのサイトは計算の根拠がわかりません。
東洋経済では 8月1日時点
全国 1.43
東京都0.99
大阪1.64
愛知1.93
福岡 1.3
沖縄 4.1
東京都の実効再生産数は全都の指数
です、
更にクラスターの地域や職種などの
情報も必要でしょう
東洋経済
「実効再生産数=
直近7日間の新規陽性者数／その前
7日間の新規陽性者数）^（平均世代
時間／報告間隔）」
平均世代時間は5日、報告間隔は7日
URLﾘﾝｸ(toyokeizai.net)
1038:１３２人目の素数さん
20/08/03 11:11:53 C8O7DWVU.net
そもそも「今の」再生産数が出てる段階でおかしいわな。どちらも。
1039:１３２人目の素数さん
20/08/03 11:28:42.92 eJM6VgsU.net
速報値だからな
1040:１３２人目の素数さん
20/08/03 11:59:22.14 fIfeS9d2.net
学術の巨大掲示板群 - アルファ・ラボ URLﾘﾝｸ(x0000.net)
数学 物理学 化学 生物学 天文学 地理地学
IT 電子 工学 言語学 国語 方言 など
1041:１３２人目の素数さん
20/08/03 12:27:13 xgcgj0Hi.net
>>974
一番下に参照したサイトとか書いてあるよ
1042:１３２人目の素数さん
20/08/03 12:29:56.16 eJM6VgsU.net
>>978
説明して
1043:１３２人目の素数さん
20/08/03 13:34:21.90 nP+Sm0D1.net
>>975
速報値だから別に構わん。
感染から陽性判明までの時間差があるから過去の値だし、
事例ごとに時間差を補正しないと正確ではないことも
みなさん承知でしょ。
1044:１３２人目の素数さん
20/08/03 13:49:14.05 7Vi8/EtE.net
>>974
実効再生産数の推定方法は使うモデルでもちがうし一概になんとも言えない
例えばsirにしても実効再生産数は感染率と回復率の両方で日々の拡大の比率が決まる事を利用するので、回復率の概算値の取り方でもズレてくる
その東洋経済の計算式もかなりデータ処理する人によって恣意的に変わりうるパラメータを含んでいるのでなんとも言えないけど、監修してるのが西浦先生なら割とどちらにもバイアスのかかってない式だと言えるかもしれない
でもどのみち再生産数は「投げたボールが自然に落ちてくる」ような法則で決まるものではなく、市民の集団心理のような数学だけではでは決まらないものなのでハッキリとした事が言えるような性格のものではないとしか言えないと思う
少なくとも今0.99になったのがなんらかの対策の結果によって全体の傾向として下がり始めているのか、たまたま、ココ数日の数値データを見て市民の意識が警戒する方にふれているのかはまだなんとも言えないと思う
前者ならこのまま1未満が持続されるだろうけど、後者で不急ではあるが不要でない用事を後回しにしてるだけなら、また増えてくるかもしれない
こんなの多分数週間単位で数字見ないとわかりっこないと思う
1045:１３２人目の素数さん
20/08/03 16:57:55.99 nP+Sm0D1.net
>>981
だんだん小さくなっていくトレンドははっきり見て取れるんじゃない？
で、実際、感染者の増加傾向の倍加時間も長くなってるように見える。
ま、連日あれだけ騒がれてりゃ、都民も飲み会自粛したいりして気を
つけるだろうからね。
逆に言えば、７月に入ってバーっと感染拡大したのは、規制解除で浮か
れて油断したからでしょ。
1046:１３２人目の素数さん
20/08/03 17:14:26.84 7Vi8/EtE.net
>>982
だからそのトレンドが方程式上ないんだよ
今日少ないから明日も少ないとは方程式の形からは言えない
あくまでβ(t)が一定値を取る傾向にあると言うのは数学上の法則ではなく経験則に過ぎない
まぁそう考えても過去の経験上、数週間～数ヶ月単位のトレンドを見る限りにおいては大きくはハズレないと言うだけ
β(t)がこのところ下がってるから明日も下がるかなんてとても言えない
実際今日東京かなり多いみたいだし
再生産数のトレンドはやはり数週間単位でないとわからないと思う
しかも今は7月の4連休の影響と梅雨明けの影響をダブルで受けてるから8/10くらいまではかなり値が揺れると思う
1047:１３２人目の素数さん
20/08/03 17:40:14.51 C8O7DWVU.net
>>980
いや、だから事例ごとに補正した値を出すべきなんじゃね？遅れも考慮して。手法は確立してるし。
もちろんある程度信頼できる値は1、２週間前のものしか求まらないから、１週間前のRtは1.3でしたと発表するしかないが。
1048:１３２人目の素数さん
20/08/03 19:15:01.96 qIXOsgp1.net
>>984
モデルから説明してくれ
1049:１３２人目の素数さん
20/08/03 21:24:43 rU5XKbyf.net
>>985
色々あるけど西浦のやり方で良いよ。
1050:１３２人目の素数さん
20/08/04 00:10:14.54 QOYlLXf+.net
>>983
>今日少ないから明日も少ないとは方程式の形からは言えない
んなこたわかってるよ。予測ができりゃ誰も苦労しないわけで。
単純に長期トレンドを外挿するのが関の山だが、それで大概は
あたりだったりする。
1051:１３２人目の素数さん
20/08/04 00:21:42.31 buqp/rtg.net
>>967
そんなデマ信じるなよ
現在厚労省が認可してるPCR検査キットは、新型コロナウイルス以外に反応した例は確認されてない
1052:１３２人目の素数さん
20/08/04 00:33:47.31 OPfgl3ou.net
>>987
あてもんしたいなら好きにすれば良い
ココの板の住人はそんなあてもん興味ない
K値中野とか宮沢は好きそうだけど
1053:１３２人目の素数さん
20/08/04 00:46:28.14 QOYlLXf+.net
俺も「ココの板の住人」だが、あてもんも必要だよ。
先のことはわからない、で済ませちゃなにも面白くない。
過去を分析するのは、未来を予想するため。
1054:１３２人目の素数さん
20/08/04 00:57:17.62 OPfgl3ou.net
>>990
なんの根拠もない数字が当たってなんでそれで未来が予測できるんだよ？
それこそまさにK値理論で中野がやらかしてる事やろ？
1055:１３２人目の素数さん
20/08/04 09:15:15 kB52B1VU.net
>>715
メキシコのデング熱ワクチンの話が業界では有名。
1056:１３２人目の素数さん
20/08/04 10:46:52.87
1057: ID:QOYlLXf+.net
1058:１３２人目の素数さん
20/08/04 10:57:46.29 NqlKcFnC.net
次スレ
数学 統計に詳しい人が語るコロナウイルス ☆2
ｽﾚﾘﾝｸ(math板)
1059:１３２人目の素数さん
20/08/04 11:44:49.84 EnxLnyf8.net
>>993
株価の予言でがんばれ
1060:１３２人目の素数さん
20/08/04 13:18:34.30 UTSYjGjo.net
>>993
だから合わせても意味ないっての
中野はそれがわかってないんだって
あるモデルを作ってそれが現実の値にあってるかどうかなんて数学的にはわからない
それがちゃんと予言性を持つモデルか、設定されてる仮説が妥当かどうかはモデルから出された預言値と実測値が合うかどうかを何度も何度も検証して実証される
sirとかはいろんな民族のいろんな社会のいろんな疫病で検証されて「β(t)が数週間くらいのオーダーでは一定」という仮説による予言性がわりと信頼ができると検証されてきた
β(t)が変化するようなモデルで信頼にたる予言性のあるモデルなんかないやろ？
1061:１３２人目の素数さん
20/08/04 13:55:10 QOYlLXf+.net
意味はないだろうね。単なる外挿だから、天気の長期予報みたいなもの。
どっちみちβは予測できないんだから、モデルを細かく云々しても
予想に限れば意味ないのは同じ。
株価のほうは市場にもろにフィードバックされるから、たとえ正しく
予想できても公言できない罠ｗ
1062:１３２人目の素数さん
20/08/04 13:58:57 QOYlLXf+.net
>>996
>β(t)が数週間くらいのオーダーでは一定
β(t)のトレンドが数週間くらいのオーダーでは一定としてもええやんか。
1063:１３２人目の素数さん
20/08/04 14:22:44.10 UTSYjGjo.net
>>398
それが一定としていい実証がないっての
もちろん一定として予測だしたけりゃ出せばいい
誰も支持しないだけ
1064:１３２人目の素数さん
20/08/04 14:59:06 NqlKcFnC.net
次スレ
数学 統計に詳しい人が語るコロナウイルス ☆2
ｽﾚﾘﾝｸ(math板)
1065:1001
Over 1000 Thread .net
このスレッドは１０００を超えました。
新しいスレッドを立ててください。
life time: 157日 12時間 40分 26秒
1066:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています