20/08/04 10:57:33 NqlKcFnC.net
前スレ
数学 統計に詳しい人が語るコロナウイルス
ｽﾚﾘﾝｸ(math板)
2:１３２人目の素数さん
20/08/04 16:31:14.51 Tyz6QTpz.net
>>1
乙
3:１３２人目の素数さん
20/08/04 21:57:56.55 QOYlLXf+.net
過去数週間β(t)が減少し続けているという実態があるのに、
「β(t)は今後数週間一定と想定していい」けど、
「β(t)は今後数週間減少しつづけると想定しちゃ駄目」
という理屈がわからん。
4:１３２人目の素数さん
20/08/04 22:58:16.22 esYlvxwz.net
SIRモデルの解析解。
・カーマック.マッケンドリックの
KM近似解.
βとγは正定数と仮定する。
R(t)をσ=β/γによる テーラー展開を使い
dR(t)/dt=asech^2(bt-c)
a、b、cは正定数
R(t)を死亡数として
1905-6年インドポンペイのペスト
流行を表すことができる。
a.b.cのパラメーターは外挿による。
・ケンドリックの解.
β(R(t)) 、γ正定数と仮定する。
テーラー展開を使わずに
dR(t)/dt=a'sech^2(b't-c')
を求めた。
・KM近似解は対称でテーラー展開の
範囲で成立する.
ケンドリック解は非対称だが
適用範囲が広い
今後も発展形を含めSIRモデルは
Cov19対策の基盤モデルだろう
5:１３２人目の素数さん
20/08/04 23:01:39.43 esYlvxwz.net
>>4
続き
・SIRモデルのβ(t)を感染症から
化学反応で類推すれば、
dI/dt=βSI⇔化学反応⇔自由分子の
気体衝突.
β(t)⇔σ<u>=反応断面積×平均速度
σの1次近似が弾性衝突で一定。
しかし平均速度が遅い.静止でも
唾を飛ばし近接した会話などは
反応断面積が極めて大きいし
マスク着用は小さくなる。
・天気予報の基礎方程式.
ナビエストークス方程式は
低気圧の発生発達乱流まで予測する
が、大域的なSIRモデルは
局所的に発生するクラスターを
予測できない。
クラスターの予測には
具体的な感染例を分析し類推とも
比較しシンプルなβ(t)の個人行動
感染モデルを構築する必要がある
と思う。
パーコレーションも参考になる?
6:１３２人目の素数さん
20/08/04 23:04:47.55 esYlvxwz.net
>β(t)は今後数週間減少しつづけると想定
大阪府ではハッタリＫ値に続き
イソジンでうがいが流行ってる話ですか、
口内のウィルスを一時的に減少させ
唾液をPCR検査すれば偽陰性
になるかもしれませんが、
体内のウィルスは変化しません。
(飲酒の検査で誤魔化す手口と同様
です。)
偽陰性の陽性者が市中を徘徊すれば
β(t)は増大します。
7:１３２人目の素数さん
20/08/04 23:29:56.51 esYlvxwz.net
>>4
×ケンドリック ⚪ケンドール
8:１３２人目の素数さん
20/08/05 00:57:59.93 v2n/Q9hi.net
>>6
引用文となんの関係があるんだよ。
あんた頭おかしいわ。
9:１３２人目の素数さん
20/08/05 01:01:03.12 xlvgOX8A.net
>>8
引用文が意味不明で
あしからず
10:１３２人目の素数さん
20/08/05 01:04:39.26 xlvgOX8A.net
偽陰性が増加なら形式上β(t)は減少し続けますね
実際は増加しますが
11:１３２人目の素数さん
20/08/05 01:17:50.70 ZposI61J.net
>>6
朝に検査をしてるが検査前にうがいはしていない
一晩経っての結果だから案ずる必要は無い
口腔内のウイルス細菌が身体に与える影響は大で口腔内から減少すると増殖率が大幅ダウンするよ
体内のウイルスが変わらないというのは根拠が無い
12:１３２人目の素数さん
20/08/05 01:36:03.27 CePUW9QK.net
歯茎から細菌とかが入って血管を通じて体内に広がるという話は聞いたことあるな
新型コロナウイルスでどうなのかは未確認だけど
13:１３２人目の素数さん
20/08/05 01:58:55.42 BYree24Q.net
さっさと政府や自治体は、簡易計算でいいから実効再生産数を毎日発表しろよ
・簡易モデルで毎日発表
全国の実効再生産数
各都道府県の実効再生産数
・詳細モデル(疫学者が計算)で週1回発表
全国の実効再生産数
人口や感染者数が多い自治体の実効再生産数
こんな感じでいいから
14:１３２人目の素数さん
20/08/05 06:52:47 JITqXXPB.net
>>11
>検査前にうがいはしていない
根拠がない
15:１３２人目の素数さん
20/08/05 09:58:40.21 qaIyIyoy.net
>>14
朝に検査する前にうがいをしてない事は発表してる
イソジンうがいは感染者が飛沫感染で他の人に移すのを防ぐのに効果がありそうに思うが、予防効果がある訳ではないと発表をした
あくまで重症化率を下げる可能性があるという発表だった
口腔内の消毒で体内のウィルスが少なくなるわけがないというのが一般的な常識だと思われるが
夜に口腔内のウィルスが唾液と共に肺に入る可能性があるので軽症者から重症化率を下げるのではないかという発表だった
他の医師の見解によると誤嚥性肺炎は口腔内ケアをすると発生率を抑えられるので、そのようなメカニズムが働いている可能性も否定はできないとの事
ただ検査対象数も少ないし、二重盲検も出来てないだろうから信憑性は劣る
16:１３２人目の素数さん
20/08/05 10:55:03 yjhNI4md.net
>>15
あくまでイソジンによるうがいですね。
>夜に口腔内のウィルスが唾液と共に肺に入る可能性があるので軽症者から重症化率を下げるのではないか
歯磨きや水他によるうがいと比較したの
可能性なら
陰性証明が問題になっているこの時期は
イソジン悪用の可能性も大きい
17:１３２人目の素数さん
20/08/05 21:08:17 SucOmwcd.net
>>9
あんたの頭が悪し。
18:１３２人目の素数さん
20/08/05 21:13:52 SucOmwcd.net
>>11
確かに。
特に舌に多いと言われるACE2受容体にとりついて増殖したウィルスを
減らせば、回復が早まっても不思議はないね。新型コロナにのみ効く
やりかたかもしれん。研究の価値はありそ。
19:１３２人目の素数さん
20/08/05 21:19:51 SucOmwcd.net
>>12
ふむふむ。口腔内で増殖したウィルスが血流に入り込み、血管壁にとりついて
血栓ができるという悪い流れにも予防効果があるかもね。
理屈の上ではメリットがあっても不思議はない。
20:１３２人目の素数さん
20/08/06 00:51:53 F98ygAwQ.net
> 研究の価値はありそ
イソジン信者さんですか
41人が対象と少ない上にうがいに
イソジンを使用と歯磨き.水や他の
うがい薬の使用と比較していないな
まあ 大阪府は夜の街でイソジンの使用
を勧めていたからPCR検査で
偽陰性が増えそ
表面上陽性率が低下しても
実際はβの増大として現れる
21:１３２人目の素数さん
20/08/06 01:51:05 7UeKJ9tJ.net
悲報
イソジンうがい薬の効果とは
大阪府のあのK値を絶賛する
京都大学ウイルス再生医科学研究所
宮沢孝幸准教授にまで「あまり期待
しない方がいい」とまでいわれる始末
URLﾘﾝｸ(news.yahoo.co.jp)
22:１３２人目の素数さん
20/08/06 13:12:17.94 K6lZ2KP4.net
あまり期待してないけど、効き目がないとも言い切れないからね。
とはいえ、重症化率が低い現況では、統計的に意味のかる臨床実験
するのも難しそう。アメリカでやってくれりゃいいのになｗ
23:１３２人目の素数さん
20/08/06 14:14:16.38 GF7YbrPw.net
>被験者の数が少ないということも
あるし、実験手順が不完全だ。
ヨードはタンパクとくっつくので、PCRの酵素にもくっつく可能性が
ある。
つまり口の中にヨードが残っている
場合、PCR検査の感度が下がってしまう可能性もある。
その影響が調べられていないし、
そのコントロールをして実験をしなければならないので、科学論文として
成立しない。
>確かに喉にいるウイルスを殺すことは可能だろうが、肺の奥や鼻腔でウイルスが増えていれば、喉のウイルスを減らしても何にもならない
ここまであのk値を絶賛するM京大准教授にすら言われているのに
いまだアホイソジン信者とは
お維新か、似非科学思考か、
PCR検査すり抜け目的か、
株価操作まで可能性がある
24:１３２人目の素数さん
20/08/06 19:20:37.08 +74uujMK.net
会見の直後に放送された「辛坊治郎ズームそこまで言うか！」(2020.08.04)の解説が秀逸
「そりゃー殺菌力の強いヨードで口内を洗浄すれば、舌の上で増殖しているウィルスが
死ぬ（注意）のは当たり前で、それ以外の何か新しい知見があったのだろうか？」（要約）
25:１３２人目の素数さん
20/08/06 20:20:49 4ondg9iD.net
>>23
>>>確かに喉にいるウイルスを殺すことは可能だろうが、
影響が大きいのは喉じゃなくて、舌だろうね。
しかし、誰も「確かな事実だ」なんて言ってないのに
必死で叩きたがるあんたの気がしれんわ。
なんだよ、信者ってｗ
あほはお前だ。
26:１３２人目の素数さん
20/08/06 20:25:10 4ondg9iD.net
ちゃんと治験して調べてみりゃ白黒つく話。
全否定することもないだろ。全否定しなきゃ実害が
出るっていうのなら話は別だけどね。
そういえば、マスク有害説ってのもあったな。WHOが
言いふらしてたわけだが、結果的に大間違いだった。
27:１３２人目の素数さん
20/08/06 21:38:07.55 loeZ7Mf+.net
>>25
>確かに喉にいるウイルスを殺すことは
可能だろうが、
>しかし、誰も「確かな事実だ」なんて言ってないのに
上はK値信者のの宮澤京大准教授の 発言だよ
宮澤京大准教授とあらそってくれw
京大でもイソジンの対照実験を
したのね
イソジンでアレルギー反応をする
人も多いなあ、
イソジン信者にはわからんだろうが
28:１３２人目の素数さん
20/08/06 21:46:15.86 loeZ7Mf+.net
βのはなしの前に
関西はK値だの、イソジンだの、
神様がいっぱいいて内紛がすごいな
まあ冷ややかに見ているんだが
29:１３２人目の素数さん
20/08/06 22:09:02 w8yiRRv5.net
>>28
発想することは否定しないけど信じるアホが多いのが関西かなぁ？
30:１３２人目の素数さん
20/08/07 00:53:55.05 uU9kwSRk.net
× 発想 ⚪妄想
31:１３２人目の素数さん
20/08/07 09:30:00 RxNjPnvX.net
>>27
＞上はK値信者のの宮澤京大准教授の 発言だよ
おまえ、とんでもない低能だな。
宮澤は「ヨードうがいでコロナが軽症化する」ことを
確かな事実だなんて言ってないだろ。その逆で、懐疑的
だろ。単に、「確かな」っていう文言が一致してるだけじゃん。
おまえ、どんだけ間抜けなんだよｗってか、キチガイだなｗ
32:１３２人目の素数さん
20/08/07 09:31:21 RxNjPnvX.net
>>28
関西でひとくくりにする馬鹿はコロナにかかって死んだほうが
いいんじゃない？馬鹿は死ななきゃ治らないっていうからね。
33:１３２人目の素数さん
20/08/07 11:58:41.08 W7LEcOm3.net
>>32
安倍ちゃん死んじゃうの？
34:１３２人目の素数さん
20/08/07 15:01:27.38 PanBrIc0.net
しかしこの話題数学板でやっても答えなんか出るハズもなし
35:１３２人目の素数さん
20/08/07 16:26:46.21 uU9kwSRk.net
>>31
なんだ
おまえさんは野蛮なお維ソ新の可能性
が大だなあー
トンデモK値の宮澤京大准教授にまで
「科学論文として成立しない」と
言われる始末
36:１３２人目の素数さん
20/08/07 16:29:15.73 uU9kwSRk.net
>>34
元々偽陰性とβの関連ですが
37:１３２人目の素数さん
20/08/08 00:07:35.52 HZRN3c99.net
>>35
維新に投票したことなどないわ、ヴォケ。
誰が論文として成立するかどうかなんて話してんだよ。
ピンボケもいいとこだよ。
おまえは、単に維新嫌い、関西嫌いで、ただただ
貶したいだけのゲス野郎にすぎん。
38:１３２人目の素数さん
20/08/08 03:03:43.41 B9yfy30n.net
関西が好かれる要素なんてそれこそどこにあるんだ？。
じゃぶじゃぶ広告屋に陰に陽にゼニ払い込んでる自信でもあるのか？。
39:１３２人目の素数さん
20/08/08 08:24:52 KrlNqCZd.net
>論文として成立するかどうかなんて話してんだよ。
prerpintでなくしっかりした査読を
通過した論文に基づいた話ではないね。
PCR検査で陽性判断が真の陽性でなく
偽陰性の増大によると、統計上陽性率
が減少しても真の感染者は対策なしで
闇に潜りβが増加する。
40:１３２人目の素数さん
20/08/08 13:33:00 TyH5eM+e.net
>>32
トンカスだろどうせｗ
41:１３２人目の素数さん
20/08/08 13:34:44 TyH5eM+e.net
コロナでバタバタ人死んでるのに
「岩手とか一人コロナになっただけでどんだけ騒ぐんだか
やっぱ田舎はこえーよな」
とか言っちゃうのがトンカスクオリティｗｗ
42:１３２人目の素数さん
20/08/08 14:06:12.08 HZRN3c99.net
>>39
別に査読論文に投稿する前に発表したっていいだろうし、
そうでないとまったく信用できないっていうような、
STAP細胞なみに荒唐無稽な話でもないからな。
まあ、PCR検査逃れに使われると困るけど、そんな手合は
ごく少数だろう。バレたら袋叩きに合うし。
43:１３２人目の素数さん
20/08/08 15:08:29.04 Ij1ACVOU.net
>PCR検査逃れに使われると困るけど、
そんな手合は ごく少数だろう。
イソジンうがいは府知事の推薦による
44:１３２人目の素数さん
20/08/08 20:22:27 VEtuiSLO.net
>>37
トンデモK値を絶賛したが、一方で
イソジンうがい「効果」ではごく普通に
評価した宮澤准教授は関西の京大だ。
45:１３２人目の素数さん
20/08/08 20:47:18 UDl/zNu1.net
>>44
その人って臨床医？
46:１３２人目の素数さん
20/08/08 21:26:56.28 dUEkPzLp.net
>>45
獣医だったハズ
47:１３２人目の素数さん
20/08/08 21:35:07.44 VEtuiSLO.net
>>45
吉村府知事は素人
48:１３２人目の素数さん
20/08/09 00:42:24.57 R5lJWQH4.net
>>44
だからなんなんだよ。
お前は宮澤もディスってるだろうが、関西だからってことで。
49:１３２人目の素数さん
20/08/09 04:41:20 H5ze4A8T.net
専門がイソジン生産なんだ
府知事語録
あとは必要があれば(イソジン薬の)
生産ラインを増強しなきゃならない
と思いますから
50:１３２人目の素数さん
20/08/09 04:55:40 H5ze4A8T.net
>>48
カッパ値さん乙
51:１３２人目の素数さん
20/08/09 07:42:46.68 R/8hHENx.net
経済重視にするには、クラスター対策がしっかりできる範囲内に感染者数を抑えるのが必須
クラスター対策がっかり出来れば、かなり実効再生産数を下げられるので、
最小限の規制・自粛で済む
いまの日本は、田舎を除いてクラスター対策が不可能になってるので、
経済的ダメージが大きくなる
52:１３２人目の素数さん
20/08/09 12:35:55.96 R5lJWQH4.net
>>51
もうクラスター潰しではおっつかないほど拡がっちゃってるからね。
ロックダウンも事実上不可能だし、経済を回しながら感染を抑える
ためには、屋内でのマスク着用義務化くらいしかないだろうな。
飲食のときもマスク着用を義務化して、口に食べ物、飲み物を放り
込む時だけマスクを外してすぐに戻すようにさせるべき。
これが広まれば、感染は間違いなく抑え込める。
（犬、猫間での感染もあるから、完全収束は難しいだろうけど）
53:１３２人目の素数さん
20/08/09 12:42:32.65 Ym6X0bb6.net
>>52
イヤ、もう行動制限以外止まらんよ
54:１３２人目の素数さん
20/08/09 13:56:25.39 R5lJWQH4.net
なんで？マスクしてりゃ、飛沫は飛ばんのだが？
55:１３２人目の素数さん
20/08/09 14:34:33.47 qVuZZJDA.net
水のうがいでも飛沫感染を低減できる効果があると検証できたら外食でのクラスター化はだいぶ減らせそうではある
客全員がまずトイレに行ってうがいはトイレスペース的に無理だろうから
フィンガーボウルみたいに食事前に水とうがいした水を吐き出す蓋つきのボウルを並べる
食事の配膳の時にそのボウルを店員が回収って流れなら新しい食事マナーとしてできなくはない
そもそも大元のうがいにどれ程の効果があるかによるが
56:１３２人目の素数さん
20/08/09 20:11:53.81 x003anze.net
マスクを過信しすぎ。
布マスクと医療用マスクで性能が
異なるし熱中症予防で外す場合もある。
少人数でも長々と会話しないこと
結局
戦略的にクラスターをいかに潰すかだ
沖縄の松山地区も一部の店に集中してる
せめて食中毒並みの法規制による
休業要請、空港など水際や米軍基地からの輸入感染阻止、
一方GO TOトラベルで大量の観光客が
満ち溢れているから本土へウィルス
拡散が心配
何やっているんだろ
57:１３２人目の素数さん
20/08/09 21:07:57.48 uWnPYTDr.net
>>54
飛沫感染しかしないならそうかもしれんが、おそらく飛沫感染だけではない
もちろん直接、関節接触感染もあるし、噂されてる空気感染もあるだろ
残念ながら全感染経路のどれくらいをそれぞれの感染経路で起こってるかなんて相当病理学的データが揃わないとわからないから数学的に推定などできようもないが
58:１３２人目の素数さん
20/08/09 21:09:07.81 R5lJWQH4.net
>>56
＞マスクを過信しすぎ。
軽んじすぎるよりはいいだろ。
マスクがツバの飛散を防ぐ効果は絶大だぞ。布マスクでも9割
がたは減らせるだろ。不織布（医療用マスクじゃないよ）だと
99%くらいいくんじゃね？（根拠薄弱ｗ）
とにかく口からツバを飛ばさないという目的に特化してマスク
をさせるのはそんなに難しくはないはず。それで感染はかなり
防げるはず。
クラスターなんてもう増えすぎちゃって、いまさらつぶすった
ってどうやんのよ？無理、無理。
59:１３２人目の素数さん
20/08/09 21:11:04.01 R5lJWQH4.net
>>57
空気感染で拡がるなら満員電車でとっくに蔓延してるはずだよ。
手洗い励行はあたりまえの話で、これはもうすでに実行されてる。
マスクをすれば手洗いしなくていいって話ではないよ。
60:１３２人目の素数さん
20/08/09 22:01:16.50 uWnPYTDr.net
>>59
それはよくある誤解
別の板ならともかくココでさすがにそのレスは笑われる
61:１３２人目の素数さん
20/08/10 00:22:16 2F2W2dQW.net
空気感染（空気感染とは言っていない）
62:１３２人目の素数さん
20/08/10 01:22:24.04 u28LsAPn.net
>>60
何言ってんのあんた。笑われてんのはあんただよ。
>関節接触感染
言い忘れたが、これもないだろうなｗ
63:１３２人目の素数さん
20/08/10 01:41:05.03 K8IUSDCp.net
>>62
バカだねぇ
どうやって現時点でコロナの感染経路についてそれぞれがどのくらいであるか定量的に評価できんの？
そんな事できるはずないでしょ？
なんでマスクのみで再生産数が一未満に抑えられるなどという結論が導けるんだよ
君がそういう仮説を信じたいなら信じとけばいい
しかしココは数学板、証明可能でないものを断定する者などいない
それから50台のガン患者がコロナに感染してると知っててフィリピンパブに行ってホステス感染させたのが関節接触感染だよ
64:１３２人目の素数さん
20/08/10 04:02:40.99 ssilZkGN.net
感染を100％防がなくても、感染をある程度防いで、実効再生産数を1未満にするだけで勝手に収束するんだよ
まあ1だといつまでも減らないし、突発的なクラスターや、外部からの持ち込みの影響を防げないので、
現実的には、実効再生産数が0.7以下になるように規制強化がいい
居酒屋・バー・カラオケ・クラブ・キャバ・ホスト・ダンス・パブ・スナック屋内空間では飲食時以外マスク義務化、
マスク外してるときは会話してはいけないってやれば、
勝手に収束するよ
65:１３２人目の素数さん
20/08/10 07:20:03.86 QEEJaVWL.net
>>64
だから収束するには仮に今実効再生産数を1.5だとすれば、再生産数を33%強減らせは均衡状態になり、40%減らせば実効再生産数は実効再生産数は0.9となり収束へ向かう
なので本当に「マスクをすれば収束する」というには「マスクをしていれば防げる感染が33%より大きい」と証明しないとダメでしょ？
そんなのどうやって証明するの？
できるわけないやん？
自分の願望を語りたければいくらでも語ればいいけど、それは願望に過ぎない
科学だというなら証明しなよ
現時点で我々がコロナについてわかってる事なんかほんの少ししかないやろ？
66:１３２人目の素数さん
20/08/10 09:29:47.57 u28LsAPn.net
>>63
疫学的な証拠がなきゃなんの証明もできないのは当たり前だよ、馬鹿。
数学以前の問題。
否定する証拠も肯定する証拠もないってだけの話。
空気感染の可能性が低く、接触感染はすでに手洗いで抑えられると
すれば、飛沫を飛ばさないマスク装着の有無が最重要ファクターに
なりうるのは馬鹿でもわかるだろ。信じる信じないってレベルの
話じゃなくて、蓋然性の問題。
しのごのいわずに社会実験してみりゃ間違いかどうかはすぐに答え
がでる。それだけのこと。
67:１３２人目の素数さん
20/08/10 09:43:44.50 m4tOcnal.net
>>66
オレが否定してるのはマスクや手洗いなどを徹底すれば再生産数を必ず1以下にできるとなんの根拠もなく断定している点
その可能性は否定するつもりはないが、必ず1以下にできるなどと断定できるハズがない
もうそのレベルでわからんならこの板で書き込む資格ないよ
68:１３２人目の素数さん
20/08/10 11:09:29.38 rn6MEtT3.net
松江市で91人感染
市によると、感染が確認された生徒
86人のうち、寮生が80人で、6人が
自宅から通っていたという。
生徒は全員がサッカー部員で、
感染者が発生した寮はサッカー部の
寮だった。
教員も同部関係者。部員約30人に
発熱や頭痛などの症状があるが、
いずれも軽症という。
サッカー部は7月23～25日に大阪府、
8月3、4日に鳥取県、4～6日に
香川県に練習試合などで遠征していたという。
69:１３２人目の素数さん
20/08/10 11:10:28.49 eWhEvjwV.net
統計学は統計学
数学ではない
応用計算の一種
70:１３２人目の素数さん
20/08/10 11:26:02.35 u28LsAPn.net
>>67
>なんの根拠もなく断定している点
断定できるわけないのは自明なんだから、レトリカルな
表現に噛み付いてもしょうがないだろ。子供の喧嘩じゃ
あるまいし。
>>69
激しく同意するわ。
71:１３２人目の素数さん
20/08/10 11:26:33.65 rn6MEtT3.net
情報幾何も数理統計も数理モデルも
数学です
例 双対接続 確率 反応拡散方程式
72:１３２人目の素数さん
20/08/10 11:30:06.64 u28LsAPn.net
>>68
結構な感染力だな。教員も感染してるってことは、
生徒との接し方を検証すべきだろうね。お互い、
飛沫を受ける状況にあったのか、なかったのか。
常時マスクして会話してたのなら、マスク有効説も
疑わざるを得なくなるが、まあ、してないだろうな。
73:１３２人目の素数さん
20/08/10 11:45:31.14 lKfbDbiO.net
>>70
なんで断定できるのが自明なんだよ
今誰がコロナの感染経路のうちマスク、手洗いで防げる経路の割合をpとしてpの値がわかるんだよ？
アホか
74:１３２人目の素数さん
20/08/10 13:01:15 Cq2ctg1u.net
>>72
レスに触らないで
75:１３２人目の素数さん
20/08/10 13:13:21 6R7zFEVs.net
なんか前スレと比べて質落ちてない？
語り尽くしちゃった？
76:１３２人目の素数さん
20/08/10 13:19:26 Cq2ctg1u.net
これではねぇ
↓
>統計学は統計学
数学ではない
応用計算の一種
>激しく同意するわ
77:１３２人目の素数さん
20/08/10 13:34:40.55 tu2mFW9U.net
>>75
数式も何も出さないレスしか出てこないんだもん
78:１３２人目の素数さん
20/08/10 13:39:19.07 Cq2ctg1u.net
数式の根拠も大切さ
79:１３２人目の素数さん
20/08/11 07:17:10 nTcwhTDK.net
関節接触感染w
軟骨接触感染は確認されまちたか?w
冷感接触感染は？w
80:１３２人目の素数さん
20/08/11 07:17:44 nTcwhTDK.net
直感接触感染とかもいいそうw
81:１３２人目の素数さん
20/08/11 07:45:11.32 LrKhZAdv.net
島根90人クラスターあっても、
感染経路不明じゃなくクラスター対策で一網打尽にできる感染は、疫学的には大して影響ない
感染経路不明だったり、クラスター対策が追い付かない感染が疫学的に危険
82:１３２人目の素数さん
20/08/11 14:42:17.72 Z41QhYaR.net
今はPCR検査の偽判断率ってどんなもんなの？
83:１３２人目の素数さん
20/08/11 18:07:40.49 FOYNfiZF.net
>>77
なんか、1人のおかしな人が乱入しただけでスレが崩壊する例が見られたな。まさに感染病。
84:１３２人目の素数さん
20/08/11 23:17:43.47 JqZPCil+.net
>>83
スルースキルのない瞬間湯沸かし器みたいな人がいるとそうなる。
85:１３２人目の素数さん
20/08/12 00:11:22.46 0pruDery.net
>数学統計に詳しい人が語るコロナ
ウィルス2
看板に偽りありだった
学習すら分かっていない
86:１３２人目の素数さん
20/08/12 01:02:38 Fqsv6mTi.net
>>73
＞なんで断定できるのが自明なんだよ
「断定できるわけないのは自明」って書いたんだけど？
アホなのはあんただろｗ
87:１３２人目の素数さん
20/08/12 01:06:01 Fqsv6mTi.net
>>85
じゃあ、もうこんなスレは放置で、あらたに、
「数学統計に詳しい人が語る株価予測」
でスレ立て願います。
88:１３２人目の素数さん
20/08/12 01:29:23 1v69O89i.net
バカばっか
89:１３２人目の素数さん
20/08/12 04:57:42 0/wbQTbU.net
>>87
イソジンと株価予測なら
インサイダー疑惑が指摘されている。
90:１３２人目の素数さん
20/08/12 08:24:52.23 9h7qwJ2K.net
>>71
追加
特異点 漸近展開
91:１３２人目の素数さん
20/08/12 10:22:56 Fqsv6mTi.net
>>89
そういうこと書いて名誉毀損で訴えられても知らんよ。
92:１３２人目の素数さん
20/08/12 11:24:53 9S0Sad4r.net
>>91
テレビ番組のコメンターの発言を
報道したからね、
調査するかは知らんけど
93:１３２人目の素数さん
20/08/12 11:31:35 9S0Sad4r.net
URLﾘﾝｸ(news.yahoo.co.jp)
94:１３２人目の素数さん
20/08/12 11:42:22 Fqsv6mTi.net
>>92
うっかり引用すれば、コメンテータと同罪で有罪になるよ。
95:１３２人目の素数さん
20/08/12 16:19:19 KrQ981jo.net
インサイダーの呼称はイソサイダーならいいの？
96:１３２人目の素数さん
20/08/12 17:54:30 qSeYF8jt.net
今月の数学セミナー、西浦がK値を痛烈に批判してるな。いやK値とは一言も書いてないが、読めばまあ分かる。
97:１３２人目の素数さん
20/08/12 21:15:58.21 Fqsv6mTi.net
数学セミナーってまだ生きてるんだ。
書店で見かけないようなきがして、廃刊になってんのかと思った。
98:１３２人目の素数さん
20/08/12 21:31:16.04 FTjTXnav.net
>>96
役に立つ情報ありがとうございます。
早速、アマゾンで注文しました。
高校数学程度の知識で理解できるかな？
Rとstanくらいは使えます。
99:１３２人目の素数さん
20/08/12 21:46:14.83 d4hU/MI9.net
＊感染症数理モデル元年に機構と外挿の狭間に立つ……西浦 博 8
＊感染症数理モデル入門／SIRコンパートメントモデルと
基本再生産数……小林鉄郎・西浦 博 14
＊集団免疫論を超えて……稲葉 寿 19
＊国内の流行初期のデータによる予測とその評価……國谷紀良 26
＊大規模流行の発生確率にまつわる数理……安齋麻美・西浦 博 30
＊新型コロナウイルスの生体内感染動態の定量化とその応用／
数理科学で挑む新型コロナウイルス感染症治療の確立
……岩波翔也・キム・クァンス・藤田泰久・江島啓介・岩見真吾 33
＊感染症流行の収束にまつわる数理……ナタリー・リントン・西浦 博 40
明日本屋で見よっと
100:１３２人目の素数さん
20/08/12 21:57:49.52 qSeYF8jt.net
>>98
いや、K値の数式を出してここがダメだ、みたいな記事じゃないから、期待するものとは違うかも。
喧嘩を売ってるわけでもないが、まあK値批判だろうなと分かる感じ。
101:１３２人目の素数さん
20/08/14 10:52:08 qijqdH+j
102:.net
103:１３２人目の素数さん
20/08/14 11:11:34.63 qijqdH+j.net
>>101
本文
「言い換えると2次の時間項を含む指数
関数がAIDSの流行曲線であり、この式
はAIDSの流行曲線が正規分布に従う
と想定していることにほかならない」
104:１３２人目の素数さん
20/08/14 11:16:31 6FR6kR83.net
改行アホ
105:１３２人目の素数さん
20/08/14 12:35:47.98 KqAf97KY.net
>>1 アラビア人の操った嘘だよ
学術の巨大掲示板群 - アルファ・ラボ URLﾘﾝｸ(x0000.net)
数学 物理学 化学 生物学 天文学 地理地学
IT 電子 工学 言語学 国語 方言 など
106:１３２人目の素数さん
20/08/14 12:42:09.82 9CprJwpw.net
>>102
本文 追加
「特に、リアルタイムで流行が進展
する中で、幼弱な政治に巻き込まれる
形で、不完全な数理モデルがあたかも
流行を確実に記述した法則のように
流布され、放置すると人の命に
かかわってしまう場合には科学的に
妥当な議論をして介入し、より適切な
未来へ導かないといけない」.
107:１３２人目の素数さん
20/08/14 20:31:04.66 Y1vrMBqV.net
>>105
それって、西浦さんが書いた文章なの？
思い上がりもはなはだしいな。駄目だこりゃ。
単純な数理モデルでたいして個体数も多くない複雑系を扱って、
正確な予測なんかできるわけもなし。科学的に妥当な議論を
つきつめれば、なんにも決定できないよ。
108:１３２人目の素数さん
20/08/14 22:37:53.92 sCLy2vas.net
>>106
数セミ9月号を読んでください
以後レスは不要です
109:１３２人目の素数さん
20/08/15 00:34:18.66 8U0p4jnM.net
無意味なレスにレスのつけようも無し。
110:１３２人目の素数さん
20/08/15 07:14:04.95 NzuwPV4R.net
かにはこうらににせてあなをほる
111:１３２人目の素数さん
20/08/15 09:37:59.78 8U0p4jnM.net
バカが墓穴を掘ってるだけだよ
112:１３２人目の素数さん
20/08/15 11:37:05.14 QXFGkfCx.net
数セミにはシンプルなモデルSIRと発展.
SEIRを詳しく述べている。
SEIRモデルより夏のピークが予測
された。
新規感染者報告数のデータと数値計算から、緊急事態宣言中に接触率が86パーセント削減され、幸運なことに予測が
外れた
113:１３２人目の素数さん
20/08/15 14:28:51 8U0p4jnM.net
>緊急事態宣言中に接触率が86パーセント削減され
結果論だよね。
問題は、接触率の削減とされているものが、外出自粛によるのか、マスク着用の効果なのか、手洗い励行の効果なのか、それぞれが
どのくらいの割合で効いてるのかわからないこと。
まあ、外出自粛が一番効くのは当たり前なんだけど、にしてもね。
114:１３２人目の素数さん
20/08/15 15:31:57 AjOk7ftg.net
↑ バカ(笑い)
115:１３２人目の素数さん
20/08/15 20:43:43.80 Vld/Bi5L.net
前スレよりコピペ（備忘録)
"
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)
"
SEIRモデルのパラメータ
SEIR2 <- function(
# Parameters
contact_rate = 10, # number of contacts per day
transmission_probability = 0.01, # transmission
116:probability beta = contact_rate * transmission_probability, # tranmission rate infectious_period = 20, # infectious period gamma = 1 / infectious_period, # Prob[infected -> recovered] latent_period = 5, # latent perior sigma = 1/latent_period, # The rate at which an exposed person becomes infective mu = 0, # The natural mortality rate nu = 0 , # vaccination moves people from susceptible to resistant directly, without becoming exposed or infected. Ro = beta/gamma, # Ro - Reproductive number. # Initial values for sub-populations. s = 99, # susceptible hosts e = 0, # exposed hosts i = 1, # infectious hosts r = 0, # recovered hosts # Compute total population. N = s + i + r + e, # Output timepoints. timepoints = seq (0, 365, by=0.5), ... )
117:１３２人目の素数さん
20/08/15 21:55:42.28 Ed9URwl6.net
数学セミナーで西浦は
「専門内外から色々な提案をいただいてありがたく思う」
って書いてあるけど、本文読み進めれば全くそんなこと思ってない、
自分のモデル以外はただの雑音って思ってることがはっきりわかる
「人文科学を専門とする人からはモデルの無理解からくる
適切でない批判もあった」と書いてあるが、間違いなく「人文科学」だけじゃないだろう
118:１３２人目の素数さん
20/08/15 23:30:51 pRrhCaOX.net
数セミ9月号
「感染症は、細菌やウィルスなど病原体への暴露がなければ絶対に羅患する
ことがないのが特徴である。
新型コロナウィルス感染症のように
人から人へ伝播するものについては
1人で部屋にこもっていれば絶対に
感染することがない。」
SIモデル(ロジスティック曲線).
SIRモデル.SEIRモデルは接触に
よる感染率βがある。
一方
K値は感染の接触概念がない。
・K値=1-(1週間前の累積感染者数 )/
(当日の累積感染者数
119:１３２人目の素数さん
20/08/16 01:15:57 cru0Ap6a.net
ところが、このウィルスはペットの犬猫にも感染する。
よって逆もあって不思議はない。エビデンスがー、って
言われればないみたいだけど、実験してみりゃわかるはずｗ
一人で部屋にこもっていても、外界と接触せずに生きること
は不可能だしね。
現実はそんなに単純ではない。一方で、大雑把な振る舞いが
極めて単純な数式で表せたりするから面白い。
120:１３２人目の素数さん
20/08/16 03:18:28 nZfBDZvz.net
ひきこもり
121:１３２人目の素数さん
20/08/16 10:39:39.74 cru0Ap6a.net
物理学だと経験式にあとから理論が追いついて
くるということもあるね。
ウィーンの公式とかケプラーの法則とか。
122:１３２人目の素数さん
20/08/16 17:37:22 9QGL0lLZ.net
数学セミナーの記事読んでて、理論易学者の認識なんてこんなもんなのかと失望したわ
現状の数理モデルに根本的に足りていないものを認識できていない
123:１３２人目の素数さん
20/08/17 13:04:41.86 aY5B6Gwq.net
こんなスレあったんだな
再感染も考慮すると健常Sと回復Rのとこに免疫ロス的な項が必要そう
ウィルスが変異して別物になるとSがリセットされて最初からになるしややこしいな
124:１３２人目の素数さん
20/08/17 16:56:32 4OWeWZ6i.net
易学なら合羽うらない
125:１３２人目の素数さん
20/08/17 18:01:35 h5k/PqLF.net
>>121
前スレよりコピペ
SEIRモデルに再感染(免疫を失ってR からSへの変遷）があるとしてグラフ化してみた。
1万人に1人の感染者と9人の保菌者がいるとして、
再感染率０％のとき
URLﾘﾝｸ(i.imgur.com)
再感染率０．１％（１０００人に１人の割合で免疫を失う）場合
URLﾘﾝｸ(i.imgur.com)
再感染率０．１％で外出自粛等で接触が８割減になった場合
URLﾘﾝｸ(i.imgur.com)
126:１３２人目の素数さん
20/08/18 18:33:18 vcGXmCXC.net
だいたい3ヶ月で免疫無くなるって話聞いたけど時間経過でR→Sになる式とか難しそう
127:１３２人目の素数さん
20/08/18 18:39:25 QilPUchd.net
>>124
R2とか作って2段階にするだけで良い。3、4段階目も要るかも知れんが。
128:１３２人目の素数さん
20/08/19 10:11:07 DU+CzCt1.net
>>124
>123が準拠したモデルはこれ
SEIRS MODEL
dS(t)/dt = mu*(N-S) - b*S(t)*I(t)/N - nu*S(t) + rho*R(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(t) + nu*S(t) - rho*R(t)
mu:自然死亡率 b:感染率(S->I)
nu:ワクチン有効率(S->R) sig:発症率(E->I),g:回復率(I->R)
rho:再感染率(R->S)
数値計算は計算機に任せた。
129:１３２人目の素数さん
20/08/21 09:24:03.36 t1p7TC6P.net
基本のSIRモデルの定理。
t≧0で初期条件が
S(0)≧0 I(0)≧0 R(0)≧0
→S(t)≧0 I(t)≧0 R(t)≧0
理由は初期条件が(S(0).0 .0)の方程式の解がI(t)=0 R(t)=0のS(t)軸に閉じ込められるや解の一意性などから。
SEIRやSEIRSモデルも同様だ
130:１３２人目の素数さん
20/08/21 21:45:00 a3WQGnu8.net
数学統計に詳しくないんだが、日本は少ない検査ながらもせっせと聞き取り調査をやって濃厚接触者を特定し陽性者を発掘するという丁寧緻密な事をやっているっていうのに、その一日千数百人の陽性者の発掘がどれだけ再生産数を押し下げているのか、試算すら無いのは残念だ。
131:１３２人目の素数さん
20/08/21 23:32:23 5qiPpY9M.net
もう困難むだ
132:１３２人目の素数さん
20/08/24 06:44:05.41 CGsqsz2y.net
>>129
すでに市中に潜行蔓延していそうだな。
133:１３２人目の素数さん
20/08/24 12:18:43.82 JCsBPlbo.net
>>128
統計学的には既に日本の感染は収束期に入ってる。
以下のグラフからも明らか。
URLﾘﾝｸ(imgur.com)
理由はいくつか考えられるが、
感染可能な人の中では
すでにある程度感染してしまっている、
ということになると思う。
この「感染可能な人」というのは、抗体、抵抗力、
社会的行動などいろいろ考えられるということ。
134:１３２人目の素数さん
20/08/24 12:29:09.42 6DjeFZ7m.net
なんで数式使うつもりないのに数学板に書き込もうとするんだろう
やめてほしい
135:１３２人目の素数さん
20/08/24 13:09:03.64 Qp4hyvjZ.net
数式さえあれば内容がある十分条件になるわけないだろカス
136:１３２人目の素数さん
20/08/24 13:16:06.58 6DjeFZ7m.net
このスレもうダメだ
いいスレだっだのに
137:１３２人目の素数さん
20/08/24 13:47:42.03 Qp4hyvjZ.net
受験偏差値操作に詳しくても疫学的エビデンスにならんしな
138:１３２人目の素数さん
20/08/24 14:04:36.46 Ppw9W3Cg.net
家族という集団を最小のクラスターとして、職場、学校、地域と
だんだん大きなクラスターを考え、クラスターサイズ、クラスター
間の感染確率をモデル化するっていう方法論はすでにあるのかな？
要するに、個体間の実効再生算数はクラスターの規模ごとに違っ
てて、さらにクラスター間の実効再生産数的なものが別にあると
いう階層性のあるモデルね。
139:１３２人目の素数さん
20/08/24 17:04:30.50 PTaqhFtT.net
前スレはめっちゃ勉強になったのに
140:１３２人目の素数さん
20/08/24 20:24:19 Ppw9W3Cg.net
もうネタは出尽くしてんだよ。
そもそもお前の勉強のためのスレじゃないし。
141:１３２人目の素数さん
20/08/24 22:27:13.24 IhCpglkj.net
前スレまだ見てないから気になるな
まあ時間空いたら�
142:ｩるか 緊急事態宣言で感染率が時間変動してるのがめんどいな あといまくそ暑いから外出減ってるのも感染者数減少の原因 気候変動みたいな効果も感染率に入れると面白い 必要あるかと言われればなんとも言えんが あと自粛継続者と、累積的に鬱憤たまって爆発して活発化した人のモデルも入れると面白いかも 累積の表現は積分に確率密度関数埋め込んで遊んでいこう
143:１３２人目の素数さん
20/08/24 22:33:30.62 PTaqhFtT.net
少し前に東京と大阪で重症者の定義が違うって話あったけど、それをパラメータにしてる人いたらたまったもんじゃないのだろうか
それとも修正する手段があるのかな
144:１３２人目の素数さん
20/08/24 22:34:31.83 IhCpglkj.net
>>136
エージェントモデルでミクロ的にやってみるのもいいんじゃないかな
詳しいことはしらんが
将棋盤みたいなマスのなかで人をごちゃごちゃ動かして、接触したら感染or変化なしでやるとか
145:１３２人目の素数さん
20/08/24 23:31:15.96 hv97BsCm.net
第3回大阪府新型コロナウィルス対策
本部専門家会議
オブザーバー 中野貴志
大阪大学核物理研究センター教授
第二波の解析結果に基づいた新しい
波の探知について
投げ上げ運動と感染者数推移を対比させている。
投げ上げ運動の原因は物体と地球間の重力相互作用であり、感染は感染者可能者.集団と感染者.集団間のウィルス相互作用だが、K値は感染概念も感染率もない。
URLﾘﾝｸ(www.pref.osaka.lg.jp)
146:１３２人目の素数さん
20/08/25 00:43:20.82 xbb2QmaK.net
>>142
×感染者可能者 ⚪感染可能者
147:１３２人目の素数さん
20/08/25 05:55:32 xxyjr/6K.net
K値は凄いモデルだよ。感染者数のピークの日にちさえ分かれば、いつ感染者数が下落に転じるのかを100%予測できるから。
148:１３２人目の素数さん
20/08/25 06:11:43.11 VuCPdy9T.net
>>144
単なる外挿じゃないの？
149:１３２人目の素数さん
20/08/25 07:06:23.76 CVCNerWF.net
覚えたての言葉無理して使わない
150:１３２人目の素数さん
20/08/25 07:10:45.16 xxyjr/6K.net
>>145
まさかマジレスが来るとは思わなんだ。
151:１３２人目の素数さん
20/08/25 07:31:38.79 StLkbsjA.net
>>146
やはり中挿の方がエロい響き
152:１３２人目の素数さん
20/08/25 10:00:13.13 TXc1/wYK.net
俺は内挿を極めた奥挿が得意技です。
153:１３２人目の素数さん
20/08/25 10:25:31.04 kHfIUPk/.net
外挿・内挿ほど、誰でも使うことは出来るが適用可能かどういう問題があるかなど基本的なことを誰も知らないテクニックはないよな
154:１３２人目の素数さん
20/08/25 10:49:16.93 TXc1/wYK.net
そうそう。相手次第なんだよね。
クリ派かどうかとか好き好きがあるから。一律に同じ技使ってもだめ。
155:１３２人目の素数さん
20/08/25 18:11:36 Nb8QQuEX.net
URLﾘﾝｸ(mainichi.jp)
156:１３２人目の素数さん
20/08/25 23:17:28.50 TXc1/wYK.net
>>152
いやらしい記事だな。
アンチ維新がみえみえで、知事の足を引っ張る気まんまんｗ
157:１３２人目の素数さん
20/08/25 23:57:31.45 4W3iOin5.net
変に知事を持ち上げている記事のほうがいやらしい記事に見えるよ。
支持する先によるってことでしょ。
158:１３２人目の素数さん
20/08/26 09:59:13.48 8ae+cQFx.net
そうだな
俺もそう思う
159:１３２人目の素数さん
20/08/26 10:08:53.08 4v9zM8x8.net
>>154
だから、それがミエミエなのがダメだってことだよ。
偏向報道が多すぎ。バカかと思うわ。バカだけどｗ
是々非々で批判なり称賛なりすればいいのに、ウガイ
薬を紹介したことにそこまでケチつけなくてもいいだろ。
ただの風邪でもやることなんだから、covidでも効きめ
がある可能性がないとはいえない。
160:１３２人目の素数さん
20/08/27 17:09:00.88 l14UPOa/.net
K値は失敗も含めて数周遅れのトップ
ランナーだな。
行列や漸化式(差分)と万有引力の法則は
高校で学んだはずだが
161:１３２人目の素数さん
20/08/29 10:43:16 WibkBkgR.net
400円で250円買ったらお釣りはいくらのスレニュー速に立ってるけど
記事って記者が理解できる、理解したことが書かれるんだな
162:１３２人目の素数さん
20/08/29 19:08:13.59 /bR8Ix+Q.net
400円札で払えばお釣りは150円だな。
163:１３２人目の素数さん
20/09/01 11:23:08.35 ygl7IezF.net
超過死亡8月計算分の報告。Rのコードも載ってて、結構楽しめそう。
URLﾘﾝｸ(www.niid.go.jp)
解説は7月分が詳しい。
URLﾘﾝｸ(www.mhlw.go.jp)
164:１３２人目の素数さん
20/09/01 15:14:06 aDwZhxer.net
公式のコロナ死者数をさっぴいたら、結局は誤差の範囲なんだろ？
165:１３２人目の素数さん
20/09/01 19:12:43.53 2qjbTlF5.net
1245
学コン・宿題ボイコット実行委員会@gakkon_boycott 9月1日
#拡散希望
#みんなで学コン・宿題をボイコットしよう
雑誌「大学への数学」の誌上で毎月開催されている学力コンテスト(学コン)と宿題は、添削が雑で採点ミスが多く、訂正をお願いしても応じてもらえない悪質なコンテストです。(私も7月号の宿題でその被害に遭いました。)このようなコンテストに参加するのは時間と努力の無駄であり、参加する価値はありません。そこで私は、これ以上の被害者を出さないようにするため、また、出版社に反省と改善を促すために、学コン・宿題のボイコットを呼び掛けることにしました。少しでも多くの方がこの活動にご賛同頂き、このツイートを拡散して頂ければ幸いです。
URLﾘﾝｸ(twitter.com)
(deleted an unsolicited ad)
166:１３２人目の素数さん
20/09/02 23:18:51.01 62Rt+r6G.net
>>160
有用な情報ありがとうございます。
Rのコードには自分の使ったこともないlibraryが記載されていてハードルが高そう。
167:１３２人目の素数さん
20/09/07 17:49:47.59 UgWSgnVg.net
数セミ9月号
小林西浦
「新型コロナウィルスにおいては夜の街クラスターとその他で二次感染
の頻度が異なることが分かっており
どのクラスターをどのように行動制限
すればどこまで総感染者数を減らす
ことができるかが計算できる」
ワクチン対策も平均でなく夜の街
クラスター集団など陽性率感染率の
データが必要だ
168:１３２人目の素数さん
20/11/16 16:04:02.99 qsoJOq7H.net
高校入試にも類題が出題されたが、
ある地域で中学生と高校生について
インターネットの利用を調査した。
調査人数は中学生100人と高校生150人である。
インターネットを利用している
人数は中学生30人 高校生40人
であった。
中学生と高校生ではどちらが
インターネットを利用しているか？
マスコミ.ワイドショー.評論家.
ある有名大学教員の解答例
40人は30人より多い。
よって高校生
169:１３２人目の素数さん
20/11/17 12:31:49.25 UJMPy762.net
米Ｊ＆Ｊ、コロナワクチン治験中断 参加者が原因不明の病気
>>
同社は９月、新型コロナワクチンの治験の最終段階に当たる第３相試験に着手したと発表。最大６万人を対象に安全性や有効性を確認するとしていた。
<<
URLﾘﾝｸ(www.jiji.com)
一方、
>>
ロシアのワクチンは旧ソ連の人工衛星にちなみ「スプートニクＶ」と名付けられた。
ランセットに掲載された論文によると、７６人を対象に第１相および第２相の治験を行ったところ、全員が新型コロナウイルスに対する抗体を獲得したという。
...
<<
URLﾘﾝｸ(www.cnn.co.jp)
76人に中ワクチン投与による致死的副作用は０らしい。
J&Jの治験が再開されて6万例でこの奇病1例だけと仮定して発生割合を1/60000とする。
ロシアのワクチンではこの奇病の報告はないから発生割合は0/76
ここで問題 ：
奇病発生率1/60000と0/76では、どちらが奇病発生が少ないといえるか？
定量的に結論を示せ。
170:１３２人目の素数さん
20/11/17 12:33:53.20 UJMPy762.net
>>165
χ二乗検定では有意差でないな。
> prop.test(c(30,40),c(100,150))
2-sample test for equality of proportions with continuity correction
data: c(30, 40) out of c(100, 150)
X-squared = 0.18601, df = 1, p-value = 0.6663
alternative hypothesis: two.sided
95 percent confidence interval:
-0.08934679 0.15601345
sample estimates:
prop 1 prop 2
0.3000000 0.2666667
171:１３２人目の素数さん
20/11/17 12:40:11.82 UJMPy762.net
>>165
有意差はでないが、リスク比（ネットジャンキーになるリスクｗ）やリスク差をグラフにすると
リスク１が中学生。７割り程度の信憑性で中学生の方が利用が多いと言える。
URLﾘﾝｸ(i.imgur.com)
172:１３２人目の素数さん
20/11/18 15:27:50.16 c7c/Z7+t.net
イーロン・マスクが検査受けて2回陰性、2回陽性だったそうだ。ベイズ的に言って陽性の確率は・・・
URLﾘﾝｸ(srad.jp)
173:１３２人目の素数さん
20/11/20 23:12:44.57 4+ut799/.net
偽陰性が出る確率のほうがはるかに高いはずだから
たぶん陽性だろうな。
174:１３２人目の素数さん
20/11/21 21:06:09.97 nPl35AWs.net
ＳＩＲみたいな微分方程式モデルは何の役にも立たない
感染状況がひろがれば政府も個人のレベルでも対策を強化するから
あてはまらないのが自明なモデルで机上の空論を振りかざすのは愚か
175:１３２人目の素数さん
20/11/21 23:43:14.33 2LkwSkD5.net
今日の新コロナ感染者数が2586と記録更新とのこと。コロナは只の風邪かを死者数から計算してみた。
URLﾘﾝｸ(www3.nhk.or.jp)
176:-all/ によると、 日本のコロナ死亡者数/感染者数=1976/130941=0.01509077で約1.5% 全世界だと1372182/57564083=0.02383747で約2.4% 数が多いのでこれくらいの差でもカイ二乗検定で有意差がでる。 X-squared = 414.56, df = 1, p-value < 2.2e-16 日本での死亡率が1976/130941として、本日の感染者2586の39人が死亡する計算になる 95%信頼区間を求めてみた,、どの計算式でも３７～４１程度という結果になった。 > data.frame(method=method,lower=lower*x,upper=upper*x) method lower upper 1 agresti-coull 37.35312 40.76990 2 asymptotic 37.31710 40.73235 3 bayes 37.33211 40.74723 4 cloglog 37.34536 40.76098 5 exact 37.33503 40.77027 6 logit 37.35338 40.76965 7 probit 37.34889 40.76477 8 profile 37.34155 40.75695 9 lrt 37.34141 40.76518 10 prop.test 37.34385 40.77959 11 wilson 37.35351 40.76950 2019年の交通事故による死亡者の数は前年比317人減の3215人で、3年連続で統計が残る1948年以来の最少を更新した。 https://www.nippon.com/ja/japan-data/h00625/ ということなので、このペースが続けば 年間の死者は14244.(95%CI 13627 14881）なので交通事故の死亡者の4.43倍という計算になるな。
177:１３２人目の素数さん
20/11/22 01:41:06.21 PV/vygCU.net
大ざっぱに、日本の人口1億2400万。近年は、毎年約90万人生まれ、約130万人子房しているようだ。
平均的な寿命は100年いかないのだから、1% 強が毎年旅立つのは自明。
つまり、現在生きている日本人をランダムに1000人選べば、一年たつと、
何らかの原因で、1%に当たる10人程は子房している。
(ランダムでは無く、高齢者に偏った1000人を選べば、何倍にも増える。）
「感染者」13万人の中で、約2000人がコロナ氏とされているようだが、
ランダムに非感染者13万人を選んだとしても、1% 強に当たる1400人程は、
一年経てば子房している。8ヶ月では、1000人程だ。
コロナ氏とされる約2000人のうち、約1000人は、非コロナ的子房と考えるのが妥当。
（「感染者」の中の高齢者の割合で調整は必要。）
真のコロナ氏者数は、雑煮をのどに詰まらせて旅立つ数と同等かもしれない。
s://yomidr.yomiuri.co.jp/article/20191217-OYTET50009/
178:１３２人目の素数さん
20/11/22 21:36:01.18 qa0hlJK8.net
>>166
0/0⇒1/2 ぽぃから
1/60000⇒2/60002 ∴約30000分の1
0/76⇒1/78 ∴78分の1
結論、モピロン
0/76 ＜ 1/60000
定量的でなくヤマ感でも直ぐに解るが
179:１３２人目の素数さん
20/11/23 07:07:46.83 snAB1ETt.net
>>173
感染してから死亡まではせいぜい数週間だろうから、
>コロナ氏とされる約2000人のうち、約1000人は、非コロナ的子房と考えるのが妥当。
の1000人は多過ぎではないかな？
180:１３２人目の素数さん
20/11/23 08:27:00.28 snAB1ETt.net
>>166
ロシアワクチンの奇病発生率をＲＳ
米Ｊ＆Ｊワクチンの奇病発生率ＪＪ
としてリスク比ＲＳ/ＪＪ分布を図示（乱数発生させて計算）すると
URLﾘﾝｸ(i.imgur.com)
99.7％の確信でＪＪの方が奇病発生が少ないと言える。
181:１３２人目の素数さん
20/11/23 09:13:12.92 vbilpfMV.net
>>176
Rで1000マンコのβ分布乱数を使って計算
> mean(rbeta(1e7,1,77)/rbeta(1e7,2,60000)>1)
[1] 0.9974359
182:１３２人目の素数さん
20/11/23 09:37:35.51 vbilpfMV.net
今週日曜日で東京の陽性数は391という。今週の陽性数の平均値の予想をしてみる。
東京コロナ 日 月 火 水 木 金 土 計
06/28の週 060 058 054 067 107 124 131 0601
07/05の週 111 102 106 075 224 243 206 1067
07/12の週 206 119 143 165 286 293 290 1502
07/19の週 188 168 237 238 366 260 295 1752
07/26の週 239 131 266 250 367 462 472 2187
08/02の週 292 258 309 263 360 461 429 2372
08/09の週 331 197 188 222 206 389 385 1918
08/16の週 260 161 207 186 339 258 256 1667
08/23の週 212 095 182 236 250 226 247 1448
08/30の週 148 100 170 141 211 136 181 1087
09/06の週 116 077 170 149 276 187 226 1201
09/13の週 146 080 191 163 171 220 218 1189
09/20の週 162 098 088 059 195 195 270 1067
09/27の週 144 078 212 194 235 196 207 1266
10/04の週 108 066 177 142 248 203 249 1193
10/11の週 146 078 166 177 284 184 235 1270
10/18の週 132 078 139 150 185 186 203 1073
10/25の週 124 102 158 171 221 204 215 1195
11/01の週 116 087 209 122 269 242 294 1339
11/08の週 189 157 293 317 393 374 352 2075
11/15の週 255 180 298 493 534 522 539 2821
11/22の週 391 *** *** *** *** *** *** 0391
前週の土曜日より減っているから今週は減少トレンドだという主張と日曜日で過去最高だから上昇トレンドという主張があったのでこんな処理をしてみた。
このデータの、日曜日の陽性数を横軸にその週の平均値を縦軸にとって線形回帰直線を描く。
URLﾘﾝｸ(i.imgur.com)
相関係数はMultiple R-squared: 0.6577, Adjusted R-squared: 0.6397 とあまり高くはないが、
日曜日の陽性数が391であったときの平均値の
予測値と95%信頼区間は
fit lwr upr
404.671 335.0596 474.2823
他に何かいい方法あるかな？
183:１３２人目の素数さん
20/11/23 10:10:53.07 e5v2JgoN.net
>>178
たとえ上昇傾向でも、２階微分がマイナスならピークに近づきつつ
あるってことで、減少トレンド（？）にあると言えなくもない。
一週間の移動平均の対数グラフの傾きをみるとよくわかる。
現状では、まだピークに近づいてはいなさそう（傾き一定のまま）？
URLﾘﾝｸ(web.sapmed.ac.jp)
184:１３２人目の素数さん
20/11/23 11:05:15.87 4qqIyMUZ.net
0/76と1/60000の比較だと
>定量的でなくヤマ感でも直ぐに解る
と仰せの通りなので問題を変えてみた。
米Ｊ＆Ｊ、コロナワクチン治験中断 参加者が原因不明の病気
同社は９月、新型コロナワクチンの治験の最終段階に当たる第３相試験に着手したと発表。最大６万人を対象に安全性や有効性を確認するとしていた。
URLﾘﾝｸ(www.jiji.com)
一方、
ロシアのワクチンは旧ソ連の人工衛星にちなみ「スプートニクＶ」と名付けられた。
ランセットに掲載された論文によると、７６人を対象に第１相および第２相の治験を行ったところ、全員が新型コロナウイルスに対する抗体を獲得したという。
...
URLﾘﾝｸ(www.cnn.co.jp)
76人に中ワクチン投与による致死的副作用は０らしい。
ここで問題
J&Jの奇病は100例集計時点での報告と仮定すると発生割合は1/100
ロシアのワクチンではこの奇病の報告はないから発生割合は0/76
このままでは、
J&Jの方が奇病発生確率が高いと判断されてしまう。
あと何例の奇病発生なしの治験者を集めればロシアのワクチンと奇病発生の確率は同等と主張できるか？
185:１３２人目の素数さん
20/11/23 11:06:13.09 d18tkwNk.net
平均値は一つの指標にすぎない。
東京都は郊外の多摩地区もあれば
23区もあるし、データもある。
186:１３２人目の素数さん
20/11/23 11:42:34.12 4qqIyMUZ.net
>>179
レスありがとうございます。
２階微分
187:、つまり、加速度で判定するということと理解しました。 １週間の移動平均とその対数をグラフにしてみました。 https://i.imgur.com/11jbMIq.png https://i.imgur.com/ixjywer.png
188:１３２人目の素数さん
20/11/23 13:09:33.69 86FAyFV7.net
k値まであと一歩
189:１３２人目の素数さん
20/11/23 13:19:16.19 mVuKIVea.net
>>180
ベータ分布での期待値が等しくなるには総治験者数が154になればいいのだが、
リスク比が１以上になる確率が0.5になるのはJ&Jの治験者総数が185のときという結果になった。
リスク比の確率密度分布を図示すると
URLﾘﾝｸ(i.imgur.com)
190:１３２人目の素数さん
20/11/23 16:07:05.50 e5v2JgoN.net
>>182
そのグラフだと移動平均になってないと思いますよ。
あと>>179で言ってる２階微分は対数の２階微分です。すみません。
真数値では下に凸（２階微分も正）でも、対数値だと上に凸みたいな。
191:１３２人目の素数さん
20/11/23 17:06:29.80 4qqIyMUZ.net
>>185
レスありがとうございます。
根本的な理解が誤っているかもしれないのでご教授お願いしたいのですが、
1週間の移動平均はi日目の陽性者数をa[i]とすると
a[1]+a[2]+a[3]+a[4]+a[5]+a[6]+a[7])/7 ,(a[2]+a[3]+a[4]+a[5]+a[6]+a[7]+a[8])/7, a[3]+a[4]+a[5]+a[6]+a[7]+a[8]+a[9])/7, ...
で計算したのですが、間違っていますか？
192:１３２人目の素数さん
20/11/23 17:18:00.46 4qqIyMUZ.net
>>185
プログラムを見直したら初歩的なミスがありました。
ご指摘ありがとうございます。
移動平均を描画しなおしてみました。
URLﾘﾝｸ(i.imgur.com)
193:１３２人目の素数さん
20/11/23 18:04:41.50 4qqIyMUZ.net
1週間移動平均の自然対数をグラフにすると、
URLﾘﾝｸ(i.imgur.com)
移動平均の一日毎の差を計算して、その値の一日毎の差を計算してグラフにしてみた。
URLﾘﾝｸ(i.imgur.com)
どう解釈すればいいのかわからないなぁ。
194:１３２人目の素数さん
20/11/23 21:39:18.72 e5v2JgoN.net
平滑化してみたら？
195:１３２人目の素数さん
20/11/23 23:58:03.06 mVuKIVea.net
>>189
レスありがとうございます。ご助言に沿って平滑化してみました。
URLﾘﾝｸ(i.imgur.com)
これでみると、いまは横ばいという判断みたいです。
196:１３２人目の素数さん
20/11/24 11:50:15.23 tZUZOwaP.net
>>188
day＞125頃から直線的に増加してる
対数とって、直線的増加のようだし、
無限回微分して0でもやはり、
指数関数的に増加してるような
気がする。でもそろそろ減るかな
先週は割と温かい。からその
影響が今週位に現れるとみる。
197:１３２人目の素数さん
20/11/24 17:24:54.89 GE4SvM8c.net
>>190
2階微分がゼロってことだとしたら、横ばいではなく、１次関数的変化。
対数が時間の１次関数ってことは、倍加時間が一定の指数関数的増大中って
ことでしょ。
198:１３２人目の素数さん
20/11/25 11:05:30.60 VLlaV3CX.net
>>192
なるほど今は自然の成り行きってことですな。
199:１３２人目の素数さん
20/11/25 15:36:18.40 Wori00xb.net
自然の成り行きというか、
対数二次微分が０ってことは、介入が影響してないってことだな。
まあ、マスク装着とGOTOという介入の影響が相殺しているという解釈ができなくもないが。
200:１３２人目の素数さん
20/11/25 16:31:57.71 HZLk1M5C.net
あいまいな言葉ではあるけど、 勢いが横ばい といわれるとなんとなく納得。
201:１３２人目の素数さん
20/11/26 14:32:08.43 d8k7sQX5.net
助けてください
次の２つの命題を満たす互いに異なる５つの実数は存在しないことを示せ
(1)どの実数も残りの4つの和より小さい
(2)任
202:意に２数a,bを選んだときa<=bなら2a<=bである
203:１３２人目の素数さん
20/11/26 15:47:45.56 bwV9Oh45.net
>>196 えっ？
シミュレーションしようと思ったけど
それ以前に問題文の意味が多様で不明
とはいえ、反例を示すと
a=0.5
b=0.6
c=0.7
d=0.8
e=0.9
という5つの異なる実数は存在する。
モチロン小数も100％実数だと思われる
by 読解力の練習できた。ありがとう
204:１３２人目の素数さん
20/11/26 20:37:56.51 nTBSTLPq.net
>>197
a=0.5
b=0.6
で
a<=bなら2a<=bである
が成立しないから反例になっていないのでは？
205:１３２人目の素数さん
20/11/26 20:48:22.83 nTBSTLPq.net
次の２つの命題を満たす互いに異なる５つの実数は存在しないことを示せ
(1)どの実数も残りの4つの和より小さい
(2)任意に２数a,bを選んだときa<=bなら2a<=bである
"
rm(list=ls())
f <- function(x) length(unique(x))==5 # 異なる実数
g <- function(x){ # (1)どの実数も残りの4つの和より小さい
flg=NA
for(i in 1:5){
if(x[i] <- sum(x[-i])){
flg=TRUE
}else{
flg=FALSE
break
}
}
return(flg)
}
# (2)任意に２数a,bを選んだときa<=bなら2a<=bである
'%=>%' = function(P,Q) !(P & !Q) # PならばQ
pm=gtools::permutations(5,2)
h <- function(x){
sub <- function(i){
a=x[i[1]]
b=x[i[2]]
(a <= b) %=>% (2*a <= b)
}
all(apply(pm,1,sub))
}
sim <- function(x) f(x) & g(x) & h(x)
sim(c(0.5,0.6,0.7,0.8,0.9))
206:１３２人目の素数さん
20/11/26 20:50:48.54 nTBSTLPq.net
> sim(c(0.5,0.6,0.7,0.8,0.9))
[1] FALSE
になったので
乱数発生させて反例を探してみた。
flg=FALSE
while(!flg){
x=runif(5)
flg=sim(x)
}
> x
[1] 0.021129725 0.073322952 0.001661962 0.181783657 0.540831530
> sim(x)
[1] TRUE
桁を丸めても成立するか確認
> (y=round(x,2))
[1] 0.02 0.07 0.00 0.18 0.54
> sim(y)
[1] TRUE
207:１３２人目の素数さん
20/11/26 20:55:37.08 nTBSTLPq.net
バグ発見
if(x[i] <- sum(x[-i])){
は
if(x[i] < sum(x[-i])){
だな
反例は撤回。
208:１３２人目の素数さん
20/11/26 20:56:42.51 nTBSTLPq.net
>>200
flg=FALSE
while(!flg){
x=runif(5)
flg=sim(x)
}
が無限ループに入ったので、どうやら成立するみたいだな。
209:１３２人目の素数さん
20/11/27 08:15:22.17 mPM/C8GR.net
>ID:nTBSTLPq
誤って書き込まれたスレ違いの書き込みに延々と何やってんの？
他スレでとっくに証明されてるよ。
3連休の検査数減少のせいかもしれんけど、そろそろピークアウト
しそうな気配があるな。あと４，５日様子を見ないと判断できないが。
210:１３２人目の素数さん
20/11/27 08:43:28.41 anGa5WFp.net
>>203
プログラム組んで体感するのが楽しいんだね。
収束することが証明されていてもグラフを書いて収束するのがみたい。
**分布に従うとか言われてもヒストグラムに分布曲線を重ねたくなるのと同じ。
中に何があるかわかっていてもタンクトップをずり下げたくなるようなものｗ
211:１３２人目の素数さん
20/11/28 23:16:05.02 71BcBuYQ.net
やはりそろそろピークアウトしてきた感じだな。
212:１３２人目の素数さん
20/12/02 22:30:09.49 LbYHBA85.net
>>205
近隣の職場でクラスター発生したりして、ピークアウトした実感はないな。
むしろ、ますます、近づいてきたという印象。
213:１３２人目の素数さん
20/12/02 22:47:59.65 LbYHBA85.net
>>178
URLﾘﾝｸ(stopcovid19.metro.tokyo.lg.jp)
に、データが出たので線形回帰での予測と比較してみる
11/22の週 319 314 186 401 481 570 561 2832
> 2832/7
[1] 404.5714
予測値が
予測値 95%CI下限 95%CI上限
404.671 335.0596 474.2823
うまくあたっているなぁ。
214:１３２人目の素数さん
20/12/03 00:17:20.95 yHlja/PL.net
>>206
あんたの近傍だけの話をされてもなぁｗ
現状はほぼ横ばいだろ。ピークに達した感がある。
215:１３２人目の素数さん
20/12/03 10:09:02.55 wLLqGYqv.net
GOTOトラベルだと
40,000,000 / 172 = 232,558人に一人感染
日本の人口12600万/感染者15万=840人に一人感染
GOTOトラベル参加は新コロナ感染リスクを何分の１に減らすことができるか？
95%信頼区間付きで計算せよ。
216:１３２人目の素数さん
20/12/03 10:09:35.85 wLLqGYqv.net
>>208
全国各地でクラスター発生が報道されてきたようにおもうが。
217:１３２人目の素数さん
20/12/03 10:33:15.16 yHlja/PL.net
>>218
数値で示さないと無意味。印象論だったらなんとでも言える。
218:１３２人目の素数さん
20/12/03 10:36:30.92 yHlja/PL.net
すでに全国的に感染が蔓延している状態なので
GOTOトラベルが感染地域の拡大に寄与してる部分
はあまりないんじゃないかな。
ただ、グループで旅行するとそのグループ内で
感染するということはあるだろうな。
219:１３２人目の素数さん
20/12/03 14:51:41.68 dEm2SN/9.net
>>211
>現状はほぼ横ばいだろ。ピークに達した感がある。
も印象論では？
220:１３２人目の素数さん
20/12/03 15:18:30.08 dEm2SN/9.net
>>210
警戒レベルを引き上げる自治体が増えてきているよなぁ。
221:１３２人目の素数さん
20/12/03 18:49:09.99 yHlja/PL.net
>>213
グラフみろよ。現時点では間違いなく横ばい。
実行再生産数もほぼ１。
ただ、これが踊り場ですぐまた上昇に転じる可能性もあるので、
ピークアウトかどうかは単なる印象。
222:１３２人目の素数さん
20/12/07 22:17:19.54 CU+4ckWm.net
>>215
同業者の職場でも新コロナが発生してきて徐々に身近に迫ってきたなぁと実感している。
223:１３２人目の素数さん
20/12/08 00:12:15.30 hzCu4Q+1.net
およそ科学的ではないな。
収束しない限り、時間が立てば身の回りで見聞きする
確率があがるのは当たり前。
224:１３２人目の素数さん
20/12/08 05:39:33.58 bQW+eNK4.net
西浦 博 他4名
新型コロナウイルスと闘うために数学にできること
がkindleからでてるな。アマゾンのリンクが貼れない・・・
225:１３２人目の素数さん
20/12/08 07:38:28.46 4aN8y6Kd.net
URLﾘﾝｸ(www.ama)
zon.co.jp/dp/B08P9Z3N7K
226:１３２人目の素数さん
20/12/09 17:48:25.06 nv/1M1Kw.net
こんな本がでるようだ
理論疫学者・西浦博の挑戦-新型コロナからいのちを守れ! (単行本) (日本語) 単行本 ? 2020/12/8
数学の本ではなさそうな
抜粋とてこんな記述が
・・・僕には脅迫状が届き、生まれて初めて殺害予告を受けました。一番緊迫した頃には、厚労省と新橋のビジネスホテルの間を歩くだけなのに警察の方に護衛してもらったことすらありました。
227:１３２人目の素数さん
20/12/09 19:58:27.56 YBZ1yNkx.net
>新型コロナからいのちを守れ!
自分の命がやばい、ってか？ｗ
228:１３２人目の素数さん
20/12/09 20:03:51.35 TYzqqYYm.net
今年のインフルエンザの確率はほぼ
zeroだから、モピロン、
ウイルス性のカゼの症状即ち、
倦怠感かつのどがヘンな感じなら、
99.9％に限りなく近い可能性で
コロナ感染症だと思われる。
多分、知らないうちに
地球人の半分位感染してて
その内に、偽陰性になってるかもね
229:１３２人目の素数さん
20/12/09 21:16:46.12 YBZ1yNkx.net
そうならいいんだけどね。
230:１３２人目の素数さん
20/12/10 04:54:10.87 LcOGRiI3.net
西浦
231:さんによるK値の批判はあるけど 逆にK値側から、SIRモデルなど理論疫学で標準的なアプローチへの批判ってあるの？
232:１３２人目の素数さん
20/12/10 10:21:05.07 oBWH18jr.net
SIRモデルも、K値も、どっちもハズレ
やっぱり、回帰分析だろ。
R言語等でシミュレーションやって
カラフルな📈にすると
分かりやすくと思う。というか
SIRモデルは、多分完璧にわかったぁ
一次化学反応論の
k1 k2
A→B→C ソックリでIがBに相当
でその積分形は、A0はAの初期値なら
B＝A0[k1/(k1-k2)][e^(-k2t)-e^ (-k1t)]
K値モデルの微分形は簡単ぽぃけど
積分形がわからなかった。
でも初期の指数関数的増加を
正確に予測出来る気がする。
現実はカオス的に複雑で
A:健康 B:感染 C:あの世 とおくと
k1 k3
A↔B→C
k2
しかもk1 k2 k3は定数でなく確率変数
233:１３２人目の素数さん
20/12/10 11:00:01.89 oBWH18jr.net
続編
んまっ、
ゴンペルツ曲線にしろ、
ロジスティクス曲線にしろ、
I＝S0[k1/(k1-k2)][e^(-k2t)-e^ (-k1t)]
のような曲線にしろ、
無限日後に感染者数の時間微分は、
モチロン、Zeroに限りなく近づく
実無限日後⇒Zeroになるハズ
で安心するだろう
それにしても、このウイルス
感染&免疫回避性能 半端なくヤヴァィ
234:１３２人目の素数さん
20/12/10 11:20:08.16 cxBeQaI7.net
普通の風邪のコロナウィルスのやっかいさと同じなのでは？
235:１３２人目の素数さん
20/12/10 11:47:34.80 oBWH18jr.net
消息不明の怪情報だと、
コウモリに感染するが人に感染しない
コロナウイルスだったのが、
数年前に、SARS遺伝子が
そのコロナウイルスに組み込れ、
あのコロナウイルスは人の気道に
感染力を持つようになっとか
現時点で消息不明の情報なハズ
236:１３２人目の素数さん
20/12/10 12:22:37.08 cxBeQaI7.net
有り得る話だけどな。
武漢のウィルス研究所で作られた可能性はそんなに低くはない。
故意に拡散させた可能性はないとは思うが、事故であれなんで
あれ漏出させた責任を中国政府とるべきだろうな。
237:１３２人目の素数さん
20/12/10 12:27:53.61 oBWH18jr.net
【新型コロナウイルスは割と安全】
SARS-CoV2は普通ぢゃないコロナだ
SARS-CoV2には亜種がいくつもある
今までのとは
違うから、そもそも新型だ。
新型のは普通のとは異なり、
サイトカインストリームという
免疫暴走の重篤肺炎が
短期間に発生することが珠にある。
表面S-タンパクが、怪しい
これを回避できれば、ま、軽いタダの
インフルエンザ程度だろう
多分。
新型以外コロナウイルス以下に示す
危ないコロナウイルス
SARS-CoV
MERS-CoV
ショボぃ風邪コロナウイルスHCoV-229E
HCoV-OC43
HCoV-NL63
HCoV-HKU1
238:１３２人目の素数さん
20/12/10 19:24:14.59 cxBeQaI7.net
ストリームじゃなくてストームだろｗ
武漢ウィルスでいいよ。ややこしいから。
239:１３２人目の素数さん
20/12/10 19:32:32.84 oBWH18jr.net
ご指導ありがとう。
HuhanV.のイチバン怖いのは、
HuhanV.自体というより、
HuhanV.と免疫細胞の相互的な
事象で、IL6が宿主が発射され
宿主により、宿主を破壊することに
ある。肺炎により死亡する
ここがヤバい
240:１３２人目の素数さん
20/12/10 20:48:28.27 Whv7aYu3.net
一昔前はセカンドアタックセオリーとかいわれていたんだが
URLﾘﾝｸ(www.jaam.jp)
241:7.html
242:１３２人目の素数さん
20/12/11 16:21:59.05 N4X2zCt3.net
URLﾘﾝｸ(youtu.be)
243:１３２人目の素数さん
20/12/12 21:15:11.20 xixiDIl1.net
ウィルス干渉説が本当なら
インフル流行させればコロナ殲滅できる
はあり？
244:１３２人目の素数さん
20/12/13 19:41:04.73 ZoF1CxDr.net
新コロに比べてインフルは感染力低すぎて不可能じゃね？
245:１３２人目の素数さん
20/12/15 21:55:06.09 SiGqyW5Y.net
コロナが流行してるのは、
インフルエンザが流行してないから
とのトンデモが現れる予感がしてきた
246:１３２人目の素数さん
20/12/15 23:48:34.33 GfZkwHJf.net
新型コロナよりインフルエンザの方が感染力高いことはすでにわかっている
247:１３２人目の素数さん
20/12/16 00:09:59.66 FKLRpKYB.net
>>238
Timeline: How The U.S. Reached 200,000 Deaths From COVID-19
トランプは２月の段階で
It's also more deadly than your, you know, even your strenuous flus.
と発言して正しい認識をしたたなぁ。
URLﾘﾝｸ(www.youtube.com)
１分１５秒辺りから
248:１３２人目の素数さん
20/12/22 00:05:01.18 kYmv3i0r.net
やはりカッパ値はかんさいの風土やまい
249:１３２人目の素数さん
20/12/24 16:38:01.60 nYQkp2Hz.net
東京 日 月 火 水 木 金 土
06/28 060 058 054 067 107 124 131 計0601
07/05 111 102 106 075 224 243 206 計1067
07/12 206 119 143 165 286 293 290 計1502
07/19 188 168 237 238 366 260 295 計1752
07/26 239 131 266 250 367 462 472 計2187
08/02 292 258 309 263 360 461 429 計2372
08/09 331 197 188 222 206 389 385 計1918
08/16 260 161 207 186 339 258 256 計1667
08/23 212 095 182 236 250 226 247 計1448
08/30 148 100 170 141 211 136 181 計1087
09/06 116 077 170 149 276 187 226 計1201
09/13 146 080 191 163 171 220 218 計1189
09/20 162 098 088 059 195 195 270 計1067
09/27 144 078 212 194 235 196 207 計1266
10/04 108 066 177 142 248 203 249 計1193
10/11 146 078 166 177 284 184 235 計1270
10/18 132 078 139 150 185 186 203 計1073
10/25 124 102 158 171 221 204 215 計1195
11/01 116 087 209 122 269 242 294 計1339
11/08 189 157 293 317 393 374 352 計2075
11/15 255 180 298 493 533 522 539 計2820
11/22 391 314 186 401 481 570 561 計2904
11/29 418 311 372 500 533 449 584 計3167
12/06 327 299 352 572 602 595 621 計3368
12/13 480 305 460 678 821 664 736 計4144
12/20 556 392 563 748 888 *** *** 計3147
1日の感染者数にベンフォードの法則*)が成り立つかを試したみた。
先頭の数字の頻度は
URLﾘﾝｸ(i.imgur.com)
x
1 2 3 4 5 6 7 8 9
59 54 20 11 15 7 7 5 2
成立しているようにみえる。
*)
URLﾘﾝｸ(ja.wikipedia.org)
ベンフォードの法則（ベンフォードのほうそく、Benford's law）とは、自然界に出てくる多くの（全てのではない）数値の最初の桁の分布が、一様ではなく、ある特定の分布になっている、という法則である。
250:１３２人目の素数さん
20/12/25 16:14:51.97 Mak8FdjO.net
>>241
しかしχ^2検定でp=0.000356。
まあ、時系列値だから仕方ないか。
251:１３２人目の素数さん
20/12/25 16:48:59.06 mzF7SS2F.net
>>241
ペンフォードの法則が成立するという事は
人為的な操作がされていなそうという解釈でいいんだよな？
252:１３２人目の素数さん
20/12/25 18:01:34.98 SYK9UM4r.net
ついに数秘術か？
253:１３２人目の素数さん
20/12/25 23:54:04.45 fnHPIfYT.net
URLﾘﾝｸ(ja.m.wikipedia.org)ベンフォードの法則
不正検出目的では2桁以降も使うって書いてある
254:１３２人目の素数さん
20/12/27 18:28:07.55 wgykouqU.net
>>241
横軸に日曜日の陽性者数、縦軸に1週間の総数をグラフにしてみ�
255:ｽ。 https://i.imgur.com/onZOdkc.png 今日（１２/２７日曜日）の感染者数708人で線形回帰すると 今月の総数の予測は > predict(lm,new,interval = 'confidence') fit lwr upr 1 5779.214 5196.727 6361.701 各曜日の予測を計算してみると 日 fit lwr upr 1 708 708 708 月 fit lwr upr 1 507.6391 458.0591 557.219 火 fit lwr upr 1 593.0795 489.3003 696.8586 水 fit lwr upr 1 904.2954 763.8492 1044.742 木 fit lwr upr 1 996.9032 842.5825 1151.224 金 fit lwr upr 1 1000.532 873.6893 1127.375 土 fit lwr upr 1 1068.765 948.5815 1188.948 1000人を超えるのは金曜日と出た。 来週になったら、どのくらいあたったか検討してみよう。
256:１３２人目の素数さん
20/12/27 18:29:28.12 wgykouqU.net
>>246（追記）
Rのコードはここに置いた。
ｽﾚﾘﾝｸ(hosp板:395番)
257:１３２人目の素数さん
20/12/27 21:09:27.59 wgykouqU.net
>>245
一連の数値nで始まる数に遭遇する確率＝log10(1+1/n)を使って2桁目をグラフにしてみた。
棒グラフが実測値、黒線がベンフォードの法則での値
URLﾘﾝｸ(i.imgur.com)
2桁めで一様分布に似てきた。
258:１３２人目の素数さん
20/12/27 23:05:55.14 wgykouqU.net
>>241(今月→今週に訂正)
横軸に日曜日の陽性者数、縦軸に1週間の総数をグラフにしてみた。
URLﾘﾝｸ(i.imgur.com)
今日（１２/２７日曜日）の感染者数708人で線形回帰すると
今週の総数の予測は
> predict(lm,new,interval = 'confidence')
fit lwr upr
1 5779.214 5196.727 6361.701
各曜日の予測を計算してみると
日 fit lwr upr
1 708 708 708
月 fit lwr upr
1 507.6391 458.0591 557.219
火 fit lwr upr
1 593.0795 489.3003 696.8586
水 fit lwr upr
1 904.2954 763.8492 1044.742
木 fit lwr upr
1 996.9032 842.5825 1151.224
金 fit lwr upr
1 1000.532 873.6893 1127.375
土 fit lwr upr
1 1068.765 948.5815 1188.948
1000人を超えるのは金曜日と出た。
来週になったら、どのくらいあたったか検討してみよう。
259:１３２人目の素数さん
20/12/28 22:12:27.05 nslNf4hr.net
>>241
日曜日と各曜日の回帰直線をグラフ化。
縦の破線は日曜日=708
URLﾘﾝｸ(i.imgur.com)
260:１３２人目の素数さん
20/12/29 00:17:44.79 9f4Zs7Qm.net
今日は481人か
月 fit lwr upr
1 507.6391 458.0591 557.219
一応、95%信頼区間458-557には入っているな。
明日の予想
火 fit lwr upr
1 593.0795 489.3003 696.8586
261:１３２人目の素数さん
20/12/29 12:42:25.83 BqqMc2cN.net
目的論的解釈からウイルスは生存と拡散を目的に弱毒化する方に変異すると言われるんだが、
すでに無症状の患者がウイルスをばらまくように変異しているわけだから、これ以上弱毒化する方向に変異するのは考えにくいと思う。
つぎの変異の進化は感染拡大するように宿主の行動を制御することだな。
カタツムリの寄生虫ロイコクロディウムみたいに。
即ち、感染すると宿主の活動性が増して外出しやすくなるようにする。
この変異の進化はすでに起きているかもしれんな。
GOTOキャンペーン利用者ってまさにそれかも。
262:１３２人目の素数さん
20/12/29 16:50:37.14 BqqMc2cN.net
>>251
東京都 新型コロナ 856人感染確認 火曜日の発表として最多更新. 2020年12月29日 16時08分 新型コロナ 国内感染者数
線形回帰での予想数の95%信頼区間を越える増加だなぁ。
263:１３２人目の素数さん
20/12/29 18:14:50.45 9f4Zs7Qm.net
>>253
予想を外したので曜日別の相関係数を計算してみた。
> for(i in 2:7) cat(day[i],' : ',summary(LM[[i]])$r.squared,'\n')
月 : 0.9042094
火 : 0.7020253
水 : 0.7968849
木 : 0.76971
金 : 0.8385985
土 : 0.8681291
264:１３２人目の素数さん
20/12/29 23:30:25.19 PIHDW0il.net
ジジイここにも湧いてたのか
URLﾘﾝｸ(hissi.org)
265:１３２人目の素数さん
20/12/29 23:31:39.04 PIHDW0il.net
医者どころか無職だろお前。
266:１３２人目の素数さん
20/12/30 00:32:41.76 uEoEZYjH.net
>>254
感染モデルもいろいろあるが、どれもたいがい指数関数的変化
が卓越するので、線形予想じゃだめなんじゃね？
267:１３２人目の素数さん
20/12/30 00:54:44.32 u3OsTh8
268:S.net
269:１３２人目の素数さん
20/12/30 07:58:53.92 7I9M4Lp6.net
>>256
無職で生活できたらいいなぁ。
最近は基幹病院ですら救急搬送を断るのでうちみたいな個人病院への救急搬送が増えた。
今月は救急受け入れの歩合給だけで10万円を越えたのでクリスマスにT-boneステーキを買って家族で食べた。
まあ、一見さんの発熱は断れと院長から指示が出ているから何でもかんでも救急を受けなくていい。
270:１３２人目の素数さん
20/12/30 08:03:50.41 7I9M4Lp6.net
>>257
レスありがとうございます。
データを対数化して線形回帰した方がいいということでしょうか？
素のデータでやった>207でよくあたっていたから自分でもびっくりした。
271:１３２人目の素数さん
20/12/30 09:50:24.34 u3OsTh8S.net
>>259
URLﾘﾝｸ(hissi.org)
お前毎日5chしかやってねーだろ。
272:１３２人目の素数さん
20/12/30 09:51:06.60 u3OsTh8S.net
今日は7時起きかい？
273:１３２人目の素数さん
20/12/30 14:12:31.01 uEoEZYjH.net
>>259
みんなが貰えた給付金10万円に追加で医療従事者の給付金ってのもあったよね。
274:１３２人目の素数さん
20/12/30 14:13:11.21 uEoEZYjH.net
>>260
そゆこと。
275:１３２人目の素数さん
20/12/30 15:41:16.79 yK/NkJWh.net
>>264
データを対数化して>249相当の値をだしてみた。
日 fit lwr upr
1 708 708 708
月 fit lwr upr
1 472 378 590
火 fit lwr upr
1 524 376 730
水 fit lwr upr
1 836 573 1222
木 fit lwr upr
1 831 620 1115
金 fit lwr upr
1 885 686 1142
土 fit lwr upr
1 922 737 1153
276:１３２人目の素数さん
20/12/30 15:47:18.46 yK/NkJWh.net
>>263
うちは、コロナ非対応病院だから、給付金５万円。
夏場に防護服着ての内視鏡検査は体力を消耗した。
277:１３２人目の素数さん
20/12/30 16:38:25.36 uEoEZYjH.net
そうなんだ。
俺の友人はコロナ対応病院なので10万もらってたが、
本人は無関係な部局でいつもと変わらぬ仕事だった
ので、気が引けると言ってた。
278:１３２人目の素数さん
20/12/30 17:23:54.05 u3OsTh8S.net
>>266
対応病院じゃなくても20万貰えた人いましたけどね。
やっぱりエセ医者かな？
279:１３２人目の素数さん
20/12/30 17:25:06.68 uEoEZYjH.net
曜日ごとのばらつきが大きいから一週間平均で予測したほうがいい。
ここ3週間は4週間で倍増ペースの比較的ゆるやかな指数関数的増加
をしている。明日は1000人超えそう。
280:１３２人目の素数さん
20/12/30 20:07:41.05 yK/NkJWh.net
>>268
コロナ病院にバイトに行っていたら２０万でたはず。
281:１３２人目の素数さん
20/12/30 20:12:03.34 yK/NkJWh.net
>>265
対数化しないほうが当たるような印象。
日 fit lwr upr
1 708 708 708
月 fit lwr upr
1 508 458 557
火 fit lwr upr
1 593 489 697
水 fit lwr upr
1 904 764 1045
木 fit lwr upr
1 997 843 1151
金 fit lwr upr
1 1001 874 1127
土 fit lwr upr
1 1069 949 1189
明日は997人 95%CI [843-1151]の予想。
282:１３２人目の素数さん
20/12/30 23:52:01.76 uEoEZYjH.net
>>271
12/3からの4週間分の木曜データを使って、エクセルで指数近似させると
明日の予測値は1100と出た。線形近似だと1033。ってか、日毎の予想
なんかしても意味ないよ。曜日ごとの変動が大きすぎる。
283:１３２人目の素数さん
20/12/31 00:06:28.13 qiffllOG.net
一週間の移動平均値をデータとして12/28までの19日分を使うと、
線形近似で1000人を超えるのは来年の1/14で、指数近似だとそれ
が1/9に早まる。
284:１３２人目の素数さん
20/12/31 00:54:32.88 ZNO2TkSk.net
そもそも大本営発表みたいなデータが信用できないから予測もへったくれもないと思いませんか？
285:１３２人目の素数さん
20/12/31 11:09:12.91 qiffllOG.net
絶対値ではなく、変動量を見る目安として使うことに問題ない。
286:１３２人目の素数さん
20/12/31 11:47:18.12 ZNO2TkSk.net
鼻くそほじりながら思いついた数字でも問題ない？
287:１３２人目の素数さん
20/12/31 12:52:20.55 uj9CqJOS.net
>>250のグラフ見てみたので、
感想文を書いてみたぁぁ。
新規、新規の、新規感染者数が
時間経過で、ほぼ線形的増加だ。
ぢゃ、入院患者も線形的増加となるか
自問自答してみた。
入院期間は定数 2week とおくと
d入院数/dt = k*t ∵kは定数∵線形増加
d退院数/dt = k*(t-2)
∴d入院患者/dt = 2*k = K （Kも定数)
入院患者 = K*t^2 (積分定数は無視）
どうやら、二次関数的に増加するハズ
もし、新規感染者が指数関数的なら
さらに不味いけどどうにかなるだろう
288:１３２人目の素数さん
20/12/31 13:01:43.05 uj9CqJOS.net
>>277の訂正
あれ、間違えた 改訂
入院患者数は、
✕ 二次関数的に増加する
○ 線形的増加に増加する
でも、その内に無限大になるな
でも、誰かどうにかするだろう
289:１３２人目の素数さん
20/12/31 16:55:30.49 qiffllOG.net
>>276
馬鹿は鼻くそでも食ってろ。
290:１３２人目の素数さん
21/01/01 01:15:15.76 REQb9+Ix.net
>>273
直近一週間ぐらいで増加率が大きくなってるな。
このままいくと、1/4には7日移動平均が1000人
を超える。20日で倍増のペース。
291:１３２人目の素数さん
21/01/01 13:03:12.48 5WRDmcLC.net
数学科の脳足りん共
お前らに予測は無理無理かたつむり
292:１３２人目の素数さん
21/01/02 20:43:48.17 1Rj7ie4D.net
>>246
今月の総数の予測は
> predict(lm,new,interval = 'confidence')
fit lwr upr
1 5779.214 5196.727 6361.701
の結果をみてみると
今週の総数は
> thisweek=c(708, 418, 856, 944, 1337, 783, 814)
> sum(thisweek)
[1] 5860
日曜日の陽性数でわりとその週の総数が当たるみたいだな。
今週のデータを加えて、日曜日の総数と週の総数の相関係数を計算してみると
Residual standard error: 349.7 on 25 degrees of freedom
Multiple R-squared: 0.9296, Adjusted R-squared: 0.9267
よい相関関係にあるな。
293:１３２人目の素数さん
21/01/02 20:48:23.49 1Rj7ie4D.net
グラフにして相関を実感してみる。
URLﾘﾝｸ(i.imgur.com)
294:１３２人目の素数さん
21/01/02 20:49:18.46 1Rj7ie4D.net
>>282
今月の総数→今週の総数の タイプミス
295:１３２人目の素数さん
21/01/02 22:41:03.64 xs0uSabE.net
>>284
相変わらずタイプミスだらけだな爺さんよ。
さすが非医だね。