X
トップページ⇒数学
1002コメント349KB
数学 統計に詳しい人が語るコロナウイルス ☆2
レス数が1000を超えています。これ以上書き込みはできません。
0002１３２人目の素数さん
垢版 |
2020/08/04(火) 16:31:14.51ID:Tyz6QTpz
>>1
乙
0003１３２人目の素数さん
垢版 |
2020/08/04(火) 21:57:56.55ID:QOYlLXf+
過去数週間β(t)が減少し続けているという実態があるのに、
「β(t)は今後数週間一定と想定していい」けど、
「β(t)は今後数週間減少しつづけると想定しちゃ駄目」
という理屈がわからん。
0004１３２人目の素数さん
垢版 |
2020/08/04(火) 22:58:16.22ID:esYlvxwz
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対策の基盤モデルだろう
0005１３２人目の素数さん
垢版 |
2020/08/04(火) 23:01:39.43ID:esYlvxwz
>>4
続き
・SIRモデルのβ(t)を感染症から
　化学反応で類推すれば、
　dI/dt=βSI⇔化学反応⇔自由分子の
　気体衝突.
　β(t)⇔σ<u>=反応断面積×平均速度
　σの1次近似が弾性衝突で一定。
　しかし平均速度が遅い.静止でも
　唾を飛ばし近接した会話などは
　反応断面積が極めて大きいし
　マスク着用は小さくなる。
・天気予報の基礎方程式.
ナビエストークス方程式は
　低気圧の発生発達乱流まで予測する
　が、大域的なSIRモデルは
　局所的に発生するクラスターを　
　予測できない。
　クラスターの予測には
　具体的な感染例を分析し類推とも
　比較しシンプルなβ(t)の個人行動
　感染モデルを構築する必要がある
　と思う。
　パーコレーションも参考になる?
0006１３２人目の素数さん
垢版 |
2020/08/04(火) 23:04:47.55ID:esYlvxwz
>β(t)は今後数週間減少しつづけると想定

大阪府ではハッタリＫ値に続き
イソジンでうがいが流行ってる話ですか、
　口内のウィルスを一時的に減少させ
　唾液をPCR検査すれば偽陰性
　になるかもしれませんが、
　体内のウィルスは変化しません。
　(飲酒の検査で誤魔化す手口と同様
　です。)
偽陰性の陽性者が市中を徘徊すれば
　β(t)は増大します。
0007１３２人目の素数さん
垢版 |
2020/08/04(火) 23:29:56.51ID:esYlvxwz
>>4
×ケンドリック　⚪ケンドール
0008１３２人目の素数さん
垢版 |
2020/08/05(水) 00:57:59.93ID:v2n/Q9hi
>>6
引用文となんの関係があるんだよ。
あんた頭おかしいわ。
0009１３２人目の素数さん
垢版 |
2020/08/05(水) 01:01:03.12ID:xlvgOX8A
>>8
引用文が意味不明で
あしからず
0010１３２人目の素数さん
垢版 |
2020/08/05(水) 01:04:39.26ID:xlvgOX8A
偽陰性が増加なら形式上β(t)は減少し続けますね
実際は増加しますが
0011１３２人目の素数さん
垢版 |
2020/08/05(水) 01:17:50.70ID:ZposI61J
>>6
朝に検査をしてるが検査前にうがいはしていない
一晩経っての結果だから案ずる必要は無い
口腔内のウイルス細菌が身体に与える影響は大で口腔内から減少すると増殖率が大幅ダウンするよ
体内のウイルスが変わらないというのは根拠が無い
0012１３２人目の素数さん
垢版 |
2020/08/05(水) 01:36:03.27ID:CePUW9QK
歯茎から細菌とかが入って血管を通じて体内に広がるという話は聞いたことあるな
新型コロナウイルスでどうなのかは未確認だけど
0013１３２人目の素数さん
垢版 |
2020/08/05(水) 01:58:55.42ID:BYree24Q
さっさと政府や自治体は、簡易計算でいいから実効再生産数を毎日発表しろよ

・簡易モデルで毎日発表
全国の実効再生産数
各都道府県の実効再生産数

・詳細モデル(疫学者が計算)で週1回発表
全国の実効再生産数
人口や感染者数が多い自治体の実効再生産数

こんな感じでいいから
0014１３２人目の素数さん
垢版 |
2020/08/05(水) 06:52:47.80ID:JITqXXPB
>>11
>検査前にうがいはしていない

根拠がない
0015１３２人目の素数さん
垢版 |
2020/08/05(水) 09:58:40.21ID:qaIyIyoy
>>14
朝に検査する前にうがいをしてない事は発表してる

イソジンうがいは感染者が飛沫感染で他の人に移すのを防ぐのに効果がありそうに思うが、予防効果がある訳ではないと発表をした
あくまで重症化率を下げる可能性があるという発表だった
口腔内の消毒で体内のウィルスが少なくなるわけがないというのが一般的な常識だと思われるが
夜に口腔内のウィルスが唾液と共に肺に入る可能性があるので軽症者から重症化率を下げるのではないかという発表だった
他の医師の見解によると誤嚥性肺炎は口腔内ケアをすると発生率を抑えられるので、そのようなメカニズムが働いている可能性も否定はできないとの事
ただ検査対象数も少ないし、二重盲検も出来てないだろうから信憑性は劣る
0016１３２人目の素数さん
垢版 |
2020/08/05(水) 10:55:03.96ID:yjhNI4md
>>15
あくまでイソジンによるうがいですね。

>夜に口腔内のウィルスが唾液と共に肺に入る可能性があるので軽症者から重症化率を下げるのではないか

歯磨きや水他によるうがいと比較したの
可能性なら
陰性証明が問題になっているこの時期は
イソジン悪用の可能性も大きい
0017１３２人目の素数さん
垢版 |
2020/08/05(水) 21:08:17.46ID:SucOmwcd
>>9
あんたの頭が悪し。
0018１３２人目の素数さん
垢版 |
2020/08/05(水) 21:13:52.45ID:SucOmwcd
>>11
確かに。
特に舌に多いと言われるACE2受容体にとりついて増殖したウィルスを
減らせば、回復が早まっても不思議はないね。新型コロナにのみ効く
やりかたかもしれん。研究の価値はありそ。
0019１３２人目の素数さん
垢版 |
2020/08/05(水) 21:19:51.65ID:SucOmwcd
>>12
ふむふむ。口腔内で増殖したウィルスが血流に入り込み、血管壁にとりついて
血栓ができるという悪い流れにも予防効果があるかもね。
理屈の上ではメリットがあっても不思議はない。
0020１３２人目の素数さん
垢版 |
2020/08/06(木) 00:51:53.04ID:F98ygAwQ
> 研究の価値はありそ

イソジン信者さんですか

41人が対象と少ない上にうがいに
イソジンを使用と歯磨き.水や他の
うがい薬の使用と比較していないな

まあ　大阪府は夜の街でイソジンの使用
を勧めていたからPCR検査で
偽陰性が増えそ
表面上陽性率が低下しても
実際はβの増大として現れる
0022１３２人目の素数さん
垢版 |
2020/08/06(木) 13:12:17.94ID:K6lZ2KP4
あまり期待してないけど、効き目がないとも言い切れないからね。
とはいえ、重症化率が低い現況では、統計的に意味のかる臨床実験
するのも難しそう。アメリカでやってくれりゃいいのになｗ
0023１３２人目の素数さん
垢版 |
2020/08/06(木) 14:14:16.38ID:GF7YbrPw
>被験者の数が少ないということも
あるし、実験手順が不完全だ。
ヨードはタンパクとくっつくので、PCRの酵素にもくっつく可能性が
ある。
つまり口の中にヨードが残っている
場合、PCR検査の感度が下がってしまう可能性もある。
その影響が調べられていないし、
そのコントロールをして実験をしなければならないので、科学論文として
成立しない。
>確かに喉にいるウイルスを殺すことは可能だろうが、肺の奥や鼻腔でウイルスが増えていれば、喉のウイルスを減らしても何にもならない

ここまであのk値を絶賛するM京大准教授にすら言われているのに
いまだアホイソジン信者とは

お維新か、似非科学思考か、
PCR検査すり抜け目的か、
株価操作まで可能性がある
0024１３２人目の素数さん
垢版 |
2020/08/06(木) 19:20:37.08ID:+74uujMK
会見の直後に放送された「辛坊治郎ズームそこまで言うか！」(2020.08.04)の解説が秀逸

「そりゃー殺菌力の強いヨードで口内を洗浄すれば、舌の上で増殖しているウィルスが
死ぬ（注意）のは当たり前で、それ以外の何か新しい知見があったのだろうか？」（要約）
0025１３２人目の素数さん
垢版 |
2020/08/06(木) 20:20:49.65ID:4ondg9iD
>>23
>>>確かに喉にいるウイルスを殺すことは可能だろうが、

影響が大きいのは喉じゃなくて、舌だろうね。

しかし、誰も「確かな事実だ」なんて言ってないのに
必死で叩きたがるあんたの気がしれんわ。
なんだよ、信者ってｗ
あほはお前だ。
0026１３２人目の素数さん
垢版 |
2020/08/06(木) 20:25:10.73ID:4ondg9iD
ちゃんと治験して調べてみりゃ白黒つく話。
全否定することもないだろ。全否定しなきゃ実害が
出るっていうのなら話は別だけどね。

そういえば、マスク有害説ってのもあったな。WHOが
言いふらしてたわけだが、結果的に大間違いだった。
0027１３２人目の素数さん
垢版 |
2020/08/06(木) 21:38:07.55ID:loeZ7Mf+
>>25
>確かに喉にいるウイルスを殺すことは
可能だろうが、
>しかし、誰も「確かな事実だ」なんて言ってないのに

上はK値信者のの宮澤京大准教授の　発言だよ
宮澤京大准教授とあらそってくれw
京大でもイソジンの対照実験を
したのね

イソジンでアレルギー反応をする
人も多いなあ、
イソジン信者にはわからんだろうが
0028１３２人目の素数さん
垢版 |
2020/08/06(木) 21:46:15.86ID:loeZ7Mf+
βのはなしの前に
関西はK値だの、イソジンだの、
神様がいっぱいいて内紛がすごいな
まあ冷ややかに見ているんだが
0030１３２人目の素数さん
垢版 |
2020/08/07(金) 00:53:55.05ID:uU9kwSRk
× 発想　⚪妄想
0031１３２人目の素数さん
垢版 |
2020/08/07(金) 09:30:00.78ID:RxNjPnvX
>>27
＞上はK値信者のの宮澤京大准教授の　発言だよ
おまえ、とんでもない低能だな。
宮澤は「ヨードうがいでコロナが軽症化する」ことを
確かな事実だなんて言ってないだろ。その逆で、懐疑的
だろ。単に、「確かな」っていう文言が一致してるだけじゃん。

おまえ、どんだけ間抜けなんだよｗってか、キチガイだなｗ
0032１３２人目の素数さん
垢版 |
2020/08/07(金) 09:31:21.49ID:RxNjPnvX
>>28
関西でひとくくりにする馬鹿はコロナにかかって死んだほうが
いいんじゃない？馬鹿は死ななきゃ治らないっていうからね。
0035１３２人目の素数さん
垢版 |
2020/08/07(金) 16:26:46.21ID:uU9kwSRk
>>31
なんだ
おまえさんは野蛮なお維ソ新の可能性
が大だなあー
トンデモK値の宮澤京大准教授にまで
「科学論文として成立しない」と
言われる始末
0036１３２人目の素数さん
垢版 |
2020/08/07(金) 16:29:15.73ID:uU9kwSRk
>>34
元々偽陰性とβの関連ですが
0037１３２人目の素数さん
垢版 |
2020/08/08(土) 00:07:35.52ID:HZRN3c99
>>35
維新に投票したことなどないわ、ヴォケ。

誰が論文として成立するかどうかなんて話してんだよ。
ピンボケもいいとこだよ。

おまえは、単に維新嫌い、関西嫌いで、ただただ
貶したいだけのゲス野郎にすぎん。
0038１３２人目の素数さん
垢版 |
2020/08/08(土) 03:03:43.41ID:B9yfy30n
関西が好かれる要素なんてそれこそどこにあるんだ？。
じゃぶじゃぶ広告屋に陰に陽にゼニ払い込んでる自信でもあるのか？。
0039１３２人目の素数さん
垢版 |
2020/08/08(土) 08:24:52.96ID:KrlNqCZd
>論文として成立するかどうかなんて話してんだよ。

prerpintでなくしっかりした査読を
通過した論文に基づいた話ではないね。
PCR検査で陽性判断が真の陽性でなく
偽陰性の増大によると、統計上陽性率
が減少しても真の感染者は対策なしで
闇に潜りβが増加する。
0040１３２人目の素数さん
垢版 |
2020/08/08(土) 13:33:00.28ID:TyH5eM+e
>>32
トンカスだろどうせｗ
0041１３２人目の素数さん
垢版 |
2020/08/08(土) 13:34:44.27ID:TyH5eM+e
コロナでバタバタ人死んでるのに
「岩手とか一人コロナになっただけでどんだけ騒ぐんだか
やっぱ田舎はこえーよな」
とか言っちゃうのがトンカスクオリティｗｗ
0042１３２人目の素数さん
垢版 |
2020/08/08(土) 14:06:12.08ID:HZRN3c99
>>39
別に査読論文に投稿する前に発表したっていいだろうし、
そうでないとまったく信用できないっていうような、
STAP細胞なみに荒唐無稽な話でもないからな。

まあ、PCR検査逃れに使われると困るけど、そんな手合は
ごく少数だろう。バレたら袋叩きに合うし。
0043１３２人目の素数さん
垢版 |
2020/08/08(土) 15:08:29.04ID:Ij1ACVOU
>PCR検査逃れに使われると困るけど、
そんな手合は ごく少数だろう。

イソジンうがいは府知事の推薦による
0044１３２人目の素数さん
垢版 |
2020/08/08(土) 20:22:27.45ID:VEtuiSLO
>>37
トンデモK値を絶賛したが、一方で
イソジンうがい「効果」ではごく普通に
評価した宮澤准教授は関西の京大だ。
0047１３２人目の素数さん
垢版 |
2020/08/08(土) 21:35:07.44ID:VEtuiSLO
>>45
吉村府知事は素人
0048１３２人目の素数さん
垢版 |
2020/08/09(日) 00:42:24.57ID:R5lJWQH4
>>44
だからなんなんだよ。
お前は宮澤もディスってるだろうが、関西だからってことで。
0049１３２人目の素数さん
垢版 |
2020/08/09(日) 04:41:20.98ID:H5ze4A8T
専門がイソジン生産なんだ

府知事語録
あとは必要があれば(イソジン薬の)
生産ラインを増強しなきゃならない
と思いますから
0050１３２人目の素数さん
垢版 |
2020/08/09(日) 04:55:40.58ID:H5ze4A8T
>>48
カッパ値さん乙
0051１３２人目の素数さん
垢版 |
2020/08/09(日) 07:42:46.68ID:R/8hHENx
経済重視にするには、クラスター対策がしっかりできる範囲内に感染者数を抑えるのが必須
クラスター対策がっかり出来れば、かなり実効再生産数を下げられるので、
最小限の規制・自粛で済む

いまの日本は、田舎を除いてクラスター対策が不可能になってるので、
経済的ダメージが大きくなる
0052１３２人目の素数さん
垢版 |
2020/08/09(日) 12:35:55.96ID:R5lJWQH4
>>51
もうクラスター潰しではおっつかないほど拡がっちゃってるからね。

ロックダウンも事実上不可能だし、経済を回しながら感染を抑える
ためには、屋内でのマスク着用義務化くらいしかないだろうな。
飲食のときもマスク着用を義務化して、口に食べ物、飲み物を放り
込む時だけマスクを外してすぐに戻すようにさせるべき。

これが広まれば、感染は間違いなく抑え込める。
（犬、猫間での感染もあるから、完全収束は難しいだろうけど）
0054１３２人目の素数さん
垢版 |
2020/08/09(日) 13:56:25.39ID:R5lJWQH4
なんで？マスクしてりゃ、飛沫は飛ばんのだが？
0055１３２人目の素数さん
垢版 |
2020/08/09(日) 14:34:33.47ID:qVuZZJDA
水のうがいでも飛沫感染を低減できる効果があると検証できたら外食でのクラスター化はだいぶ減らせそうではある
客全員がまずトイレに行ってうがいはトイレスペース的に無理だろうから
フィンガーボウルみたいに食事前に水とうがいした水を吐き出す蓋つきのボウルを並べる
食事の配膳の時にそのボウルを店員が回収って流れなら新しい食事マナーとしてできなくはない

そもそも大元のうがいにどれ程の効果があるかによるが
0056１３２人目の素数さん
垢版 |
2020/08/09(日) 20:11:53.81ID:x003anze
マスクを過信しすぎ。
布マスクと医療用マスクで性能が
異なるし熱中症予防で外す場合もある。
少人数でも長々と会話しないこと

結局
戦略的にクラスターをいかに潰すかだ
沖縄の松山地区も一部の店に集中してる
せめて食中毒並みの法規制による
休業要請、空港など水際や米軍基地からの輸入感染阻止、

一方GO TOトラベルで大量の観光客が
満ち溢れているから本土へウィルス
拡散が心配
何やっているんだろ
0057１３２人目の素数さん
垢版 |
2020/08/09(日) 21:07:57.48ID:uWnPYTDr
>>54
飛沫感染しかしないならそうかもしれんが、おそらく飛沫感染だけではない
もちろん直接、関節接触感染もあるし、噂されてる空気感染もあるだろ
残念ながら全感染経路のどれくらいをそれぞれの感染経路で起こってるかなんて相当病理学的データが揃わないとわからないから数学的に推定などできようもないが
0058１３２人目の素数さん
垢版 |
2020/08/09(日) 21:09:07.81ID:R5lJWQH4
>>56
＞マスクを過信しすぎ。
軽んじすぎるよりはいいだろ。

マスクがツバの飛散を防ぐ効果は絶大だぞ。布マスクでも9割
がたは減らせるだろ。不織布（医療用マスクじゃないよ）だと
99%くらいいくんじゃね？（根拠薄弱ｗ）

とにかく口からツバを飛ばさないという目的に特化してマスク
をさせるのはそんなに難しくはないはず。それで感染はかなり
防げるはず。

クラスターなんてもう増えすぎちゃって、いまさらつぶすった
ってどうやんのよ？無理、無理。
0059１３２人目の素数さん
垢版 |
2020/08/09(日) 21:11:04.01ID:R5lJWQH4
>>57
空気感染で拡がるなら満員電車でとっくに蔓延してるはずだよ。

手洗い励行はあたりまえの話で、これはもうすでに実行されてる。
マスクをすれば手洗いしなくていいって話ではないよ。
0062１３２人目の素数さん
垢版 |
2020/08/10(月) 01:22:24.04ID:u28LsAPn
>>60
何言ってんのあんた。笑われてんのはあんただよ。


>関節接触感染
言い忘れたが、これもないだろうなｗ
0063１３２人目の素数さん
垢版 |
2020/08/10(月) 01:41:05.03ID:K8IUSDCp
>>62
バカだねぇ
どうやって現時点でコロナの感染経路についてそれぞれがどのくらいであるか定量的に評価できんの？
そんな事できるはずないでしょ？
なんでマスクのみで再生産数が一未満に抑えられるなどという結論が導けるんだよ
君がそういう仮説を信じたいなら信じとけばいい
しかしココは数学板、証明可能でないものを断定する者などいない
それから50台のガン患者がコロナに感染してると知っててフィリピンパブに行ってホステス感染させたのが関節接触感染だよ
0064１３２人目の素数さん
垢版 |
2020/08/10(月) 04:02:40.99ID:ssilZkGN
感染を100％防がなくても、感染をある程度防いで、実効再生産数を1未満にするだけで勝手に収束するんだよ
まあ1だといつまでも減らないし、突発的なクラスターや、外部からの持ち込みの影響を防げないので、
現実的には、実効再生産数が0.7以下になるように規制強化がいい

居酒屋・バー・カラオケ・クラブ・キャバ・ホスト・ダンス・パブ・スナック屋内空間では飲食時以外マスク義務化、
マスク外してるときは会話してはいけないってやれば、
勝手に収束するよ
0065１３２人目の素数さん
垢版 |
2020/08/10(月) 07:20:03.86ID:QEEJaVWL
>>64
だから収束するには仮に今実効再生産数を1.5だとすれば、再生産数を33%強減らせは均衡状態になり、40%減らせば実効再生産数は実効再生産数は0.9となり収束へ向かう
なので本当に「マスクをすれば収束する」というには「マスクをしていれば防げる感染が33%より大きい」と証明しないとダメでしょ？
そんなのどうやって証明するの？
できるわけないやん？
自分の願望を語りたければいくらでも語ればいいけど、それは願望に過ぎない
科学だというなら証明しなよ
現時点で我々がコロナについてわかってる事なんかほんの少ししかないやろ？
0066１３２人目の素数さん
垢版 |
2020/08/10(月) 09:29:47.57ID:u28LsAPn
>>63
疫学的な証拠がなきゃなんの証明もできないのは当たり前だよ、馬鹿。
数学以前の問題。

否定する証拠も肯定する証拠もないってだけの話。
空気感染の可能性が低く、接触感染はすでに手洗いで抑えられると
すれば、飛沫を飛ばさないマスク装着の有無が最重要ファクターに
なりうるのは馬鹿でもわかるだろ。信じる信じないってレベルの
話じゃなくて、蓋然性の問題。
しのごのいわずに社会実験してみりゃ間違いかどうかはすぐに答え
がでる。それだけのこと。
0067１３２人目の素数さん
垢版 |
2020/08/10(月) 09:43:44.50ID:m4tOcnal
>>66
オレが否定してるのはマスクや手洗いなどを徹底すれば再生産数を必ず1以下にできるとなんの根拠もなく断定している点
その可能性は否定するつもりはないが、必ず1以下にできるなどと断定できるハズがない
もうそのレベルでわからんならこの板で書き込む資格ないよ
0068１３２人目の素数さん
垢版 |
2020/08/10(月) 11:09:29.38ID:rn6MEtT3
松江市で91人感染

市によると、感染が確認された生徒
86人のうち、寮生が80人で、6人が
自宅から通っていたという。
生徒は全員がサッカー部員で、
感染者が発生した寮はサッカー部の
寮だった。
教員も同部関係者。部員約30人に
発熱や頭痛などの症状があるが、
いずれも軽症という。
　サッカー部は7月23〜25日に大阪府、
8月3、4日に鳥取県、4〜6日に
香川県に練習試合などで遠征していたという。
0070１３２人目の素数さん
垢版 |
2020/08/10(月) 11:26:02.35ID:u28LsAPn
>>67
>なんの根拠もなく断定している点

断定できるわけないのは自明なんだから、レトリカルな
表現に噛み付いてもしょうがないだろ。子供の喧嘩じゃ
あるまいし。

>>69
激しく同意するわ。
0071１３２人目の素数さん
垢版 |
2020/08/10(月) 11:26:33.65ID:rn6MEtT3
情報幾何も数理統計も数理モデルも
数学です
例　双対接続　確率　反応拡散方程式
0072１３２人目の素数さん
垢版 |
2020/08/10(月) 11:30:06.64ID:u28LsAPn
>>68
結構な感染力だな。教員も感染してるってことは、
生徒との接し方を検証すべきだろうね。お互い、
飛沫を受ける状況にあったのか、なかったのか。
常時マスクして会話してたのなら、マスク有効説も
疑わざるを得なくなるが、まあ、してないだろうな。
0073１３２人目の素数さん
垢版 |
2020/08/10(月) 11:45:31.14ID:lKfbDbiO
>>70
なんで断定できるのが自明なんだよ
今誰がコロナの感染経路のうちマスク、手洗いで防げる経路の割合をpとしてpの値がわかるんだよ？
アホか
0074１３２人目の素数さん
垢版 |
2020/08/10(月) 13:01:15.17ID:Cq2ctg1u
>>72
レスに触らないで
0076１３２人目の素数さん
垢版 |
2020/08/10(月) 13:19:26.29ID:Cq2ctg1u
これではねぇ
↓
>統計学は統計学
数学ではない
応用計算の一種
>激しく同意するわ
0078１３２人目の素数さん
垢版 |
2020/08/10(月) 13:39:19.07ID:Cq2ctg1u
数式の根拠も大切さ
0079１３２人目の素数さん
垢版 |
2020/08/11(火) 07:17:10.38ID:nTcwhTDK
関節接触感染w

軟骨接触感染は確認されまちたか?w

冷感接触感染は？w
0080１３２人目の素数さん
垢版 |
2020/08/11(火) 07:17:44.46ID:nTcwhTDK
直感接触感染とかもいいそうw
0081１３２人目の素数さん
垢版 |
2020/08/11(火) 07:45:11.32ID:LrKhZAdv
島根90人クラスターあっても、
感染経路不明じゃなくクラスター対策で一網打尽にできる感染は、疫学的には大して影響ない
感染経路不明だったり、クラスター対策が追い付かない感染が疫学的に危険
0082１３２人目の素数さん
垢版 |
2020/08/11(火) 14:42:17.72ID:Z41QhYaR
今はPCR検査の偽判断率ってどんなもんなの？
0083１３２人目の素数さん
垢版 |
2020/08/11(火) 18:07:40.49ID:FOYNfiZF
>>77
なんか、1人のおかしな人が乱入しただけでスレが崩壊する例が見られたな。まさに感染病。
0085１３２人目の素数さん
垢版 |
2020/08/12(水) 00:11:22.46ID:0pruDery
>数学統計に詳しい人が語るコロナ
ウィルス2

看板に偽りありだった
学習すら分かっていない
0086１３２人目の素数さん
垢版 |
2020/08/12(水) 01:02:38.03ID:Fqsv6mTi
>>73
＞なんで断定できるのが自明なんだよ

「断定できるわけないのは自明」って書いたんだけど？
アホなのはあんただろｗ
0087１３２人目の素数さん
垢版 |
2020/08/12(水) 01:06:01.85ID:Fqsv6mTi
>>85
じゃあ、もうこんなスレは放置で、あらたに、
「数学統計に詳しい人が語る株価予測」
でスレ立て願います。
0089１３２人目の素数さん
垢版 |
2020/08/12(水) 04:57:42.63ID:0/wbQTbU
>>87
イソジンと株価予測なら
インサイダー疑惑が指摘されている。
0090１３２人目の素数さん
垢版 |
2020/08/12(水) 08:24:52.23ID:9h7qwJ2K
>>71
追加
特異点　漸近展開
0091１３２人目の素数さん
垢版 |
2020/08/12(水) 10:22:56.19ID:Fqsv6mTi
>>89
そういうこと書いて名誉毀損で訴えられても知らんよ。
0092１３２人目の素数さん
垢版 |
2020/08/12(水) 11:24:53.21ID:9S0Sad4r
>>91
テレビ番組のコメンターの発言を
報道したからね、
調査するかは知らんけど
0094１３２人目の素数さん
垢版 |
2020/08/12(水) 11:42:22.24ID:Fqsv6mTi
>>92
うっかり引用すれば、コメンテータと同罪で有罪になるよ。
0096１３２人目の素数さん
垢版 |
2020/08/12(水) 17:54:30.00ID:qSeYF8jt
今月の数学セミナー、西浦がK値を痛烈に批判してるな。いやK値とは一言も書いてないが、読めばまあ分かる。
0097１３２人目の素数さん
垢版 |
2020/08/12(水) 21:15:58.21ID:Fqsv6mTi
数学セミナーってまだ生きてるんだ。
書店で見かけないようなきがして、廃刊になってんのかと思った。
0098１３２人目の素数さん
垢版 |
2020/08/12(水) 21:31:16.04ID:FTjTXnav
>>96
役に立つ情報ありがとうございます。
早速、アマゾンで注文しました。
高校数学程度の知識で理解できるかな？
Rとstanくらいは使えます。
0099１３２人目の素数さん
垢版 |
2020/08/12(水) 21:46:14.83ID:d4hU/MI9
＊感染症数理モデル元年に機構と外挿の狭間に立つ……西浦 博　8
＊感染症数理モデル入門／SIRコンパートメントモデルと
　基本再生産数……小林鉄郎・西浦 博　14
＊集団免疫論を超えて……稲葉 寿　19
＊国内の流行初期のデータによる予測とその評価……國谷紀良　26
＊大規模流行の発生確率にまつわる数理……安齋麻美・西浦 博　30
＊新型コロナウイルスの生体内感染動態の定量化とその応用／
　数理科学で挑む新型コロナウイルス感染症治療の確立
　　……岩波翔也・キム・クァンス・藤田泰久・江島啓介・岩見真吾　33
＊感染症流行の収束にまつわる数理……ナタリー・リントン・西浦 博　40

明日本屋で見よっと
0100１３２人目の素数さん
垢版 |
2020/08/12(水) 21:57:49.52ID:qSeYF8jt
>>98
いや、K値の数式を出してここがダメだ、みたいな記事じゃないから、期待するものとは違うかも。
喧嘩を売ってるわけでもないが、まあK値批判だろうなと分かる感じ。
0101１３２人目の素数さん
垢版 |
2020/08/14(金) 10:52:08.04ID:qijqdH+j
数セミ9月号を買ってはじめの記事
から読んでいる途中、オーソドックスな視点から疑問に答えている。
説明が明確かつ合理的で行間を埋め
学ぶ必要があると思う

・機構的なモデル=現象を数理的に
　記述し数理構造を理解.と外挿的
　モデル=正しいっぽい曲線で現象論　
　的に観察データを説明.の解説は優れ
　教訓的だ。
　ロジスティック曲線は2種モデルや　
　SIモデルから導かれるし単なる
　現象論の式ではない
・人口学のFarrの法則を天然痘流行
　曲線へ応用し様々な流行曲線を
　正規分布で近似を試みたが、
　仮説が間違えていて失敗した。
しかし1980年代米国で無批判にHIV/
　AIDSの疫学予測へ応用したことが
　あるそうだ。
・AISの年間新規患者数について
　比の比が定数→差分と階差→2階
　差分方程式→2階微分方程式→
正規分布で近似→予測を間違えた
(「力学系に明るい方なら」→
　差分と階差と思う)

阪大.大阪府のK値については
明らかだろうが官邸まで扱って
いる
0102１３２人目の素数さん
垢版 |
2020/08/14(金) 11:11:34.63ID:qijqdH+j
>>101
本文
「言い換えると2次の時間項を含む指数
関数がAIDSの流行曲線であり、この式
はAIDSの流行曲線が正規分布に従う
と想定していることにほかならない」
0104１３２人目の素数さん
垢版 |
2020/08/14(金) 12:35:47.98ID:KqAf97KY
>>1 アラビア人の操った嘘だよ

学術の巨大掲示板群 - アルファ・ラボ ttp://x0000.net
数学 物理学 化学 生物学 天文学 地理地学
IT 電子 工学 言語学 国語 方言 など
0105１３２人目の素数さん
垢版 |
2020/08/14(金) 12:42:09.82ID:9CprJwpw
>>102
本文　追加
「特に、リアルタイムで流行が進展
する中で、幼弱な政治に巻き込まれる
形で、不完全な数理モデルがあたかも
流行を確実に記述した法則のように
流布され、放置すると人の命に
かかわってしまう場合には科学的に
妥当な議論をして介入し、より適切な
未来へ導かないといけない」.
0106１３２人目の素数さん
垢版 |
2020/08/14(金) 20:31:04.66ID:Y1vrMBqV
>>105
それって、西浦さんが書いた文章なの？
思い上がりもはなはだしいな。駄目だこりゃ。

単純な数理モデルでたいして個体数も多くない複雑系を扱って、
正確な予測なんかできるわけもなし。科学的に妥当な議論を
つきつめれば、なんにも決定できないよ。
0107１３２人目の素数さん
垢版 |
2020/08/14(金) 22:37:53.92ID:sCLy2vas
>>106
数セミ9月号を読んでください
以後レスは不要です
0108１３２人目の素数さん
垢版 |
2020/08/15(土) 00:34:18.66ID:8U0p4jnM
無意味なレスにレスのつけようも無し。
0109１３２人目の素数さん
垢版 |
2020/08/15(土) 07:14:04.95ID:NzuwPV4R
かにはこうらににせてあなをほる
0110１３２人目の素数さん
垢版 |
2020/08/15(土) 09:37:59.78ID:8U0p4jnM
バカが墓穴を掘ってるだけだよ
0111１３２人目の素数さん
垢版 |
2020/08/15(土) 11:37:05.14ID:QXFGkfCx
数セミにはシンプルなモデルSIRと発展.
SEIRを詳しく述べている。
SEIRモデルより夏のピークが予測
された。
新規感染者報告数のデータと数値計算から、緊急事態宣言中に接触率が86パーセント削減され、幸運なことに予測が
外れた
0112１３２人目の素数さん
垢版 |
2020/08/15(土) 14:28:51.59ID:8U0p4jnM
>緊急事態宣言中に接触率が86パーセント削減され

結果論だよね。
問題は、接触率の削減とされているものが、外出自粛によるのか、マスク着用の効果なのか、手洗い励行の効果なのか、それぞれが
どのくらいの割合で効いてるのかわからないこと。
まあ、外出自粛が一番効くのは当たり前なんだけど、にしてもね。
0113１３２人目の素数さん
垢版 |
2020/08/15(土) 15:31:57.88ID:AjOk7ftg
↑ バカ(笑い)
0114１３２人目の素数さん
垢版 |
2020/08/15(土) 20:43:43.80ID:Vld/Bi5L
前スレよりコピペ（備忘録)
"
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 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),
...
)
0115１３２人目の素数さん
垢版 |
2020/08/15(土) 21:55:42.28ID:Ed9URwl6
数学セミナーで西浦は
「専門内外から色々な提案をいただいてありがたく思う」
って書いてあるけど、本文読み進めれば全くそんなこと思ってない、
自分のモデル以外はただの雑音って思ってることがはっきりわかる

「人文科学を専門とする人からはモデルの無理解からくる
適切でない批判もあった」と書いてあるが、間違いなく「人文科学」だけじゃないだろう
0116１３２人目の素数さん
垢版 |
2020/08/15(土) 23:30:51.11ID:pRrhCaOX
数セミ9月号
「感染症は、細菌やウィルスなど病原体への暴露がなければ絶対に羅患する
ことがないのが特徴である。
新型コロナウィルス感染症のように
人から人へ伝播するものについては
1人で部屋にこもっていれば絶対に
感染することがない。」

SIモデル(ロジスティック曲線).
SIRモデル.SEIRモデルは接触に
よる感染率βがある。
一方
K値は感染の接触概念がない。
・K値=1-(1週間前の累積感染者数 )/
(当日の累積感染者数
0117１３２人目の素数さん
垢版 |
2020/08/16(日) 01:15:57.95ID:cru0Ap6a
ところが、このウィルスはペットの犬猫にも感染する。
よって逆もあって不思議はない。エビデンスがー、って
言われればないみたいだけど、実験してみりゃわかるはずｗ

一人で部屋にこもっていても、外界と接触せずに生きること
は不可能だしね。

現実はそんなに単純ではない。一方で、大雑把な振る舞いが
極めて単純な数式で表せたりするから面白い。
0118１３２人目の素数さん
垢版 |
2020/08/16(日) 03:18:28.59ID:nZfBDZvz
ひきこもり
0119１３２人目の素数さん
垢版 |
2020/08/16(日) 10:39:39.74ID:cru0Ap6a
物理学だと経験式にあとから理論が追いついて
くるということもあるね。

ウィーンの公式とかケプラーの法則とか。
0120１３２人目の素数さん
垢版 |
2020/08/16(日) 17:37:22.23ID:9QGL0lLZ
数学セミナーの記事読んでて、理論易学者の認識なんてこんなもんなのかと失望したわ
現状の数理モデルに根本的に足りていないものを認識できていない
0121１３２人目の素数さん
垢版 |
2020/08/17(月) 13:04:41.86ID:aY5B6Gwq
こんなスレあったんだな

再感染も考慮すると健常Sと回復Rのとこに免疫ロス的な項が必要そう
ウィルスが変異して別物になるとSがリセットされて最初からになるしややこしいな
0122１３２人目の素数さん
垢版 |
2020/08/17(月) 16:56:32.99ID:4OWeWZ6i
易学なら合羽うらない
0123１３２人目の素数さん
垢版 |
2020/08/17(月) 18:01:35.49ID:h5k/PqLF
>>121
前スレよりコピペ

SEIRモデルに再感染(免疫を失ってR からSへの変遷）があるとしてグラフ化してみた。

1万人に1人の感染者と9人の保菌者がいるとして、

再感染率０％のとき
https://i.imgur.com/XWzQN7I.png

再感染率０．１％（１０００人に１人の割合で免疫を失う）場合
https://i.imgur.com/EdxIhMT.png

再感染率０．１％で外出自粛等で接触が８割減になった場合
https://i.imgur.com/sDppUfY.png
0124１３２人目の素数さん
垢版 |
2020/08/18(火) 18:33:18.90ID:vcGXmCXC
だいたい3ヶ月で免疫無くなるって話聞いたけど時間経過でR→Sになる式とか難しそう
0125１３２人目の素数さん
垢版 |
2020/08/18(火) 18:39:25.59ID:QilPUchd
>>124
R2とか作って2段階にするだけで良い。3、4段階目も要るかも知れんが。
0126１３２人目の素数さん
垢版 |
2020/08/19(水) 10:11:07.54ID:DU+CzCt1
>>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)

数値計算は計算機に任せた。
0127１３２人目の素数さん
垢版 |
2020/08/21(金) 09:24:03.36ID:t1p7TC6P
基本の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モデルも同様だ
0128１３２人目の素数さん
垢版 |
2020/08/21(金) 21:45:00.18ID:a3WQGnu8
数学統計に詳しくないんだが、日本は少ない検査ながらもせっせと聞き取り調査をやって濃厚接触者を特定し陽性者を発掘するという丁寧緻密な事をやっているっていうのに、その一日千数百人の陽性者の発掘がどれだけ再生産数を押し下げているのか、試算すら無いのは残念だ。
0129１３２人目の素数さん
垢版 |
2020/08/21(金) 23:32:23.69ID:5qiPpY9M
もう困難むだ
0131１３２人目の素数さん
垢版 |
2020/08/24(月) 12:18:43.82ID:JCsBPlbo
>>128
統計学的には既に日本の感染は収束期に入ってる。
以下のグラフからも明らか。

https://imgur.com/a/yezEiuc

理由はいくつか考えられるが、
感染可能な人の中では
すでにある程度感染してしまっている、
ということになると思う。
この「感染可能な人」というのは、抗体、抵抗力、
社会的行動などいろいろ考えられるということ。
0132１３２人目の素数さん
垢版 |
2020/08/24(月) 12:29:09.42ID:6DjeFZ7m
なんで数式使うつもりないのに数学板に書き込もうとするんだろう
やめてほしい
0136１３２人目の素数さん
垢版 |
2020/08/24(月) 14:04:36.46ID:Ppw9W3Cg
家族という集団を最小のクラスターとして、職場、学校、地域と
だんだん大きなクラスターを考え、クラスターサイズ、クラスター
間の感染確率をモデル化するっていう方法論はすでにあるのかな？
要するに、個体間の実効再生算数はクラスターの規模ごとに違っ
てて、さらにクラスター間の実効再生産数的なものが別にあると
いう階層性のあるモデルね。
0138１３２人目の素数さん
垢版 |
2020/08/24(月) 20:24:19.76ID:Ppw9W3Cg
もうネタは出尽くしてんだよ。
そもそもお前の勉強のためのスレじゃないし。
0139１３２人目の素数さん
垢版 |
2020/08/24(月) 22:27:13.24ID:IhCpglkj
前スレまだ見てないから気になるな
まあ時間空いたら見るか

緊急事態宣言で感染率が時間変動してるのがめんどいな
あといまくそ暑いから外出減ってるのも感染者数減少の原因
気候変動みたいな効果も感染率に入れると面白い
必要あるかと言われればなんとも言えんが

あと自粛継続者と、累積的に鬱憤たまって爆発して活発化した人のモデルも入れると面白いかも
累積の表現は積分に確率密度関数埋め込んで遊んでいこう
0140１３２人目の素数さん
垢版 |
2020/08/24(月) 22:33:30.62ID:PTaqhFtT
少し前に東京と大阪で重症者の定義が違うって話あったけど、それをパラメータにしてる人いたらたまったもんじゃないのだろうか
それとも修正する手段があるのかな
0141１３２人目の素数さん
垢版 |
2020/08/24(月) 22:34:31.83ID:IhCpglkj
>>136
エージェントモデルでミクロ的にやってみるのもいいんじゃないかな
詳しいことはしらんが
将棋盤みたいなマスのなかで人をごちゃごちゃ動かして、接触したら感染or変化なしでやるとか
0142１３２人目の素数さん
垢版 |
2020/08/24(月) 23:31:15.96ID:hv97BsCm
第3回大阪府新型コロナウィルス対策
本部専門家会議
オブザーバー　中野貴志
大阪大学核物理研究センター教授

第二波の解析結果に基づいた新しい
波の探知について

投げ上げ運動と感染者数推移を対比させている。
投げ上げ運動の原因は物体と地球間の重力相互作用であり、感染は感染者可能者.集団と感染者.集団間のウィルス相互作用だが、K値は感染概念も感染率もない。

http://www.pref.osaka.lg.jp/attach/37375/00366623/12%20siryou6.pptx
0143１３２人目の素数さん
垢版 |
2020/08/25(火) 00:43:20.82ID:xbb2QmaK
>>142
×感染者可能者　⚪感染可能者
0144１３２人目の素数さん
垢版 |
2020/08/25(火) 05:55:32.87ID:xxyjr/6K
K値は凄いモデルだよ。感染者数のピークの日にちさえ分かれば、いつ感染者数が下落に転じるのかを100%予測できるから。
0147１３２人目の素数さん
垢版 |
2020/08/25(火) 07:10:45.16ID:xxyjr/6K
>>145
まさかマジレスが来るとは思わなんだ。
0149１３２人目の素数さん
垢版 |
2020/08/25(火) 10:00:13.13ID:TXc1/wYK
俺は内挿を極めた奥挿が得意技です。
0150１３２人目の素数さん
垢版 |
2020/08/25(火) 10:25:31.04ID:kHfIUPk/
外挿・内挿ほど、誰でも使うことは出来るが適用可能かどういう問題があるかなど基本的なことを誰も知らないテクニックはないよな
0151１３２人目の素数さん
垢版 |
2020/08/25(火) 10:49:16.93ID:TXc1/wYK
そうそう。相手次第なんだよね。
クリ派かどうかとか好き好きがあるから。一律に同じ技使ってもだめ。
0153１３２人目の素数さん
垢版 |
2020/08/25(火) 23:17:28.50ID:TXc1/wYK
>>152
いやらしい記事だな。
アンチ維新がみえみえで、知事の足を引っ張る気まんまんｗ
0154１３２人目の素数さん
垢版 |
2020/08/25(火) 23:57:31.45ID:4W3iOin5
変に知事を持ち上げている記事のほうがいやらしい記事に見えるよ。
支持する先によるってことでしょ。
0155１３２人目の素数さん
垢版 |
2020/08/26(水) 09:59:13.48ID:8ae+cQFx
そうだな
俺もそう思う
0156１３２人目の素数さん
垢版 |
2020/08/26(水) 10:08:53.08ID:4v9zM8x8
>>154
だから、それがミエミエなのがダメだってことだよ。
偏向報道が多すぎ。バカかと思うわ。バカだけどｗ

是々非々で批判なり称賛なりすればいいのに、ウガイ
薬を紹介したことにそこまでケチつけなくてもいいだろ。
ただの風邪でもやることなんだから、covidでも効きめ
がある可能性がないとはいえない。
0157１３２人目の素数さん
垢版 |
2020/08/27(木) 17:09:00.88ID:l14UPOa/
K値は失敗も含めて数周遅れのトップ
ランナーだな。
行列や漸化式(差分)と万有引力の法則は
高校で学んだはずだが
0158１３２人目の素数さん
垢版 |
2020/08/29(土) 10:43:16.90ID:WibkBkgR
400円で250円買ったらお釣りはいくらのスレニュー速に立ってるけど
記事って記者が理解できる、理解したことが書かれるんだな
0159１３２人目の素数さん
垢版 |
2020/08/29(土) 19:08:13.59ID:/bR8Ix+Q
400円札で払えばお釣りは150円だな。
0161１３２人目の素数さん
垢版 |
2020/09/01(火) 15:14:06.13ID:aDwZhxer
公式のコロナ死者数をさっぴいたら、結局は誤差の範囲なんだろ？
0162１３２人目の素数さん
垢版 |
2020/09/01(火) 19:12:43.53ID:2qjbTlF5
1245
学コン・宿題ボイコット実行委員会@gakkon_boycott 9月1日
#拡散希望
#みんなで学コン・宿題をボイコットしよう
雑誌「大学への数学」の誌上で毎月開催されている学力コンテスト(学コン)と宿題は、添削が雑で採点ミスが多く、訂正をお願いしても応じてもらえない悪質なコンテストです。(私も7月号の宿題でその被害に遭いました。)このようなコンテストに参加するのは時間と努力の無駄であり、参加する価値はありません。そこで私は、これ以上の被害者を出さないようにするため、また、出版社に反省と改善を促すために、学コン・宿題のボイコットを呼び掛けることにしました。少しでも多くの方がこの活動にご賛同頂き、このツイートを拡散して頂ければ幸いです。
https://twitter.com/gakkon_boycott/status/1300459618326388737
https://twitter.com/5chan_nel (5ch newer account)
0163１３２人目の素数さん
垢版 |
2020/09/02(水) 23:18:51.01ID:62Rt+r6G
>>160
有用な情報ありがとうございます。
Rのコードには自分の使ったこともないlibraryが記載されていてハードルが高そう。
0164１３２人目の素数さん
垢版 |
2020/09/07(月) 17:49:47.59ID:UgWSgnVg
数セミ9月号
小林西浦
「新型コロナウィルスにおいては夜の街クラスターとその他で二次感染
の頻度が異なることが分かっており
どのクラスターをどのように行動制限
すればどこまで総感染者数を減らす
ことができるかが計算できる」
ワクチン対策も平均でなく夜の街
クラスター集団など陽性率感染率の
データが必要だ
0165１３２人目の素数さん
垢版 |
2020/11/16(月) 16:04:02.99ID:qsoJOq7H
高校入試にも類題が出題されたが、

ある地域で中学生と高校生について
インターネットの利用を調査した。
調査人数は中学生100人と高校生150人である。
インターネットを利用している
人数は中学生30人 高校生40人
であった。
中学生と高校生ではどちらが
インターネットを利用しているか？

マスコミ.ワイドショー.評論家.
ある有名大学教員の解答例　
40人は30人より多い。
よって高校生
0166１３２人目の素数さん
垢版 |
2020/11/17(火) 12:31:49.25ID:UJMPy762
米Ｊ＆Ｊ、コロナワクチン治験中断　参加者が原因不明の病気
>>
同社は９月、新型コロナワクチンの治験の最終段階に当たる第３相試験に着手したと発表。最大６万人を対象に安全性や有効性を確認するとしていた。
<<
https://www.jiji.com/jc/article?k=2020101300367&;g=int

一方、

>>
ロシアのワクチンは旧ソ連の人工衛星にちなみ「スプートニクＶ」と名付けられた。
ランセットに掲載された論文によると、７６人を対象に第１相および第２相の治験を行ったところ、全員が新型コロナウイルスに対する抗体を獲得したという。
...
<<
https://www.cnn.co.jp/world/35159183.html
76人に中ワクチン投与による致死的副作用は０らしい。


J&Jの治験が再開されて6万例でこの奇病1例だけと仮定して発生割合を1/60000とする。
ロシアのワクチンではこの奇病の報告はないから発生割合は0/76

ここで問題　：

　奇病発生率1/60000と0/76では、どちらが奇病発生が少ないといえるか？
　定量的に結論を示せ。
0167１３２人目の素数さん
垢版 |
2020/11/17(火) 12:33:53.20ID:UJMPy762
>>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
0168１３２人目の素数さん
垢版 |
2020/11/17(火) 12:40:11.82ID:UJMPy762
>>165
有意差はでないが、リスク比（ネットジャンキーになるリスクｗ）やリスク差をグラフにすると
リスク１が中学生。７割り程度の信憑性で中学生の方が利用が多いと言える。

https://i.imgur.com/mgvWyd8.jpg
0170１３２人目の素数さん
垢版 |
2020/11/20(金) 23:12:44.57ID:4+ut799/
偽陰性が出る確率のほうがはるかに高いはずだから
たぶん陽性だろうな。
0171１３２人目の素数さん
垢版 |
2020/11/21(土) 21:06:09.97ID:nPl35AWs
ＳＩＲみたいな微分方程式モデルは何の役にも立たない
感染状況がひろがれば政府も個人のレベルでも対策を強化するから

あてはまらないのが自明なモデルで机上の空論を振りかざすのは愚か
0172１３２人目の素数さん
垢版 |
2020/11/21(土) 23:43:14.33ID:2LkwSkD5
今日の新コロナ感染者数が2586と記録更新とのこと。コロナは只の風邪かを死者数から計算してみた。

https://www3.nhk.or.jp/news/special/coronavirus/data-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倍という計算になるな。
0173１３２人目の素数さん
垢版 |
2020/11/22(日) 01:41:06.21ID:PV/vygCU
大ざっぱに、日本の人口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/
0174１３２人目の素数さん
垢版 |
2020/11/22(日) 21:36:01.18ID:qa0hlJK8
>>166
0/0⇒1/2　ぽぃから
1/60000⇒2/60002　∴約30000分の1
0/76⇒1/78 ∴78分の1

結論、モピロン
0/76 ＜ 1/60000

定量的でなくヤマ感でも直ぐに解るが
0175１３２人目の素数さん
垢版 |
2020/11/23(月) 07:07:46.83ID:snAB1ETt
>>173
感染してから死亡まではせいぜい数週間だろうから、
>コロナ氏とされる約2000人のうち、約1000人は、非コロナ的子房と考えるのが妥当。
の1000人は多過ぎではないかな？
0176１３２人目の素数さん
垢版 |
2020/11/23(月) 08:27:00.28ID:snAB1ETt
>>166


ロシアワクチンの奇病発生率をＲＳ
米Ｊ＆Ｊワクチンの奇病発生率ＪＪ
としてリスク比ＲＳ/ＪＪ分布を図示（乱数発生させて計算）すると
https://i.imgur.com/X92uIJv.png

99.7％の確信でＪＪの方が奇病発生が少ないと言える。
0177１３２人目の素数さん
垢版 |
2020/11/23(月) 09:13:12.92ID:vbilpfMV
>>176
Rで1000マンコのβ分布乱数を使って計算

> mean(rbeta(1e7,1,77)/rbeta(1e7,2,60000)>1)
[1] 0.9974359
0178１３２人目の素数さん
垢版 |
2020/11/23(月) 09:37:35.51ID:vbilpfMV
今週日曜日で東京の陽性数は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

前週の土曜日より減っているから今週は減少トレンドだという主張と日曜日で過去最高だから上昇トレンドという主張があったのでこんな処理をしてみた。

このデータの、日曜日の陽性数を横軸にその週の平均値を縦軸にとって線形回帰直線を描く。
https://i.imgur.com/eFouC9y.png
相関係数はMultiple R-squared: 0.6577, Adjusted R-squared: 0.6397 　とあまり高くはないが、
日曜日の陽性数が391であったときの平均値の
予測値と95%信頼区間は
fit lwr upr
404.671 335.0596 474.2823

他に何かいい方法あるかな？
0179１３２人目の素数さん
垢版 |
2020/11/23(月) 10:10:53.07ID:e5v2JgoN
>>178
たとえ上昇傾向でも、２階微分がマイナスならピークに近づきつつ
あるってことで、減少トレンド（？）にあると言えなくもない。
一週間の移動平均の対数グラフの傾きをみるとよくわかる。
現状では、まだピークに近づいてはいなさそう（傾き一定のまま）？

https://web.sapmed.ac.jp/canmol/coronavirus/japan.html?y=0
0180１３２人目の素数さん
垢版 |
2020/11/23(月) 11:05:15.87ID:4qqIyMUZ
0/76と1/60000の比較だと
>定量的でなくヤマ感でも直ぐに解る
と仰せの通りなので問題を変えてみた。


米Ｊ＆Ｊ、コロナワクチン治験中断　参加者が原因不明の病気

同社は９月、新型コロナワクチンの治験の最終段階に当たる第３相試験に着手したと発表。最大６万人を対象に安全性や有効性を確認するとしていた。
https://www.jiji.com/jc/article?k=2020101300367&;g=int

一方、

ロシアのワクチンは旧ソ連の人工衛星にちなみ「スプートニクＶ」と名付けられた。
ランセットに掲載された論文によると、７６人を対象に第１相および第２相の治験を行ったところ、全員が新型コロナウイルスに対する抗体を獲得したという。
...
https://www.cnn.co.jp/world/35159183.html
76人に中ワクチン投与による致死的副作用は０らしい。


ここで問題

J&Jの奇病は100例集計時点での報告と仮定すると発生割合は1/100
ロシアのワクチンではこの奇病の報告はないから発生割合は0/76
このままでは、
J&Jの方が奇病発生確率が高いと判断されてしまう。
あと何例の奇病発生なしの治験者を集めればロシアのワクチンと奇病発生の確率は同等と主張できるか？
0181１３２人目の素数さん
垢版 |
2020/11/23(月) 11:06:13.09ID:d18tkwNk
平均値は一つの指標にすぎない。
東京都は郊外の多摩地区もあれば
23区もあるし、データもある。
0183１３２人目の素数さん
垢版 |
2020/11/23(月) 13:09:33.69ID:86FAyFV7
k値まであと一歩
0184１３２人目の素数さん
垢版 |
2020/11/23(月) 13:19:16.19ID:mVuKIVea
>>180
ベータ分布での期待値が等しくなるには総治験者数が154になればいいのだが、
リスク比が１以上になる確率が0.5になるのはJ&Jの治験者総数が185のときという結果になった。
リスク比の確率密度分布を図示すると
https://i.imgur.com/nBGboGD.png
0185１３２人目の素数さん
垢版 |
2020/11/23(月) 16:07:05.50ID:e5v2JgoN
>>182
そのグラフだと移動平均になってないと思いますよ。
あと>>179で言ってる２階微分は対数の２階微分です。すみません。
真数値では下に凸（２階微分も正）でも、対数値だと上に凸みたいな。
0186１３２人目の素数さん
垢版 |
2020/11/23(月) 17:06:29.80ID:4qqIyMUZ
>>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, ...
で計算したのですが、間違っていますか？
0189１３２人目の素数さん
垢版 |
2020/11/23(月) 21:39:18.72ID:e5v2JgoN
平滑化してみたら？
0191１３２人目の素数さん
垢版 |
2020/11/24(火) 11:50:15.23ID:tZUZOwaP
>>188
day＞125頃から直線的に増加してる
対数とって、直線的増加のようだし、
無限回微分して0でもやはり、
指数関数的に増加してるような
気がする。でもそろそろ減るかな
先週は割と温かい。からその
影響が今週位に現れるとみる。
0192１３２人目の素数さん
垢版 |
2020/11/24(火) 17:24:54.89ID:GE4SvM8c
>>190
2階微分がゼロってことだとしたら、横ばいではなく、１次関数的変化。
対数が時間の１次関数ってことは、倍加時間が一定の指数関数的増大中って
ことでしょ。
0194１３２人目の素数さん
垢版 |
2020/11/25(水) 15:36:18.40ID:Wori00xb
自然の成り行きというか、
対数二次微分が０ってことは、介入が影響してないってことだな。
まあ、マスク装着とGOTOという介入の影響が相殺しているという解釈ができなくもないが。
0195１３２人目の素数さん
垢版 |
2020/11/25(水) 16:31:57.71ID:HZLk1M5C
あいまいな言葉ではあるけど、　勢いが横ばい　といわれるとなんとなく納得。
0196１３２人目の素数さん
垢版 |
2020/11/26(木) 14:32:08.43ID:d8k7sQX5
助けてください

次の２つの命題を満たす互いに異なる５つの実数は存在しないことを示せ
(1)どの実数も残りの4つの和より小さい
(2)任意に２数a,bを選んだときa<=bなら2a<=bである
0197１３２人目の素数さん
垢版 |
2020/11/26(木) 15:47:45.56ID:bwV9Oh45
>>196　えっ？
シミュレーションしようと思ったけど
それ以前に問題文の意味が多様で不明
とはいえ、反例を示すと

a=0.5
b=0.6
c=0.7
d=0.8
e=0.9
という5つの異なる実数は存在する。
モチロン小数も100％実数だと思われる

by 読解力の練習できた。ありがとう
0199１３２人目の素数さん
垢版 |
2020/11/26(木) 20:48:22.83ID:nTBSTLPq
次の２つの命題を満たす互いに異なる５つの実数は存在しないことを示せ
(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))
0200１３２人目の素数さん
垢版 |
2020/11/26(木) 20:50:48.54ID:nTBSTLPq
> 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
0202１３２人目の素数さん
垢版 |
2020/11/26(木) 20:56:42.51ID:nTBSTLPq
>>200
flg=FALSE
while(!flg){
x=runif(5)
flg=sim(x)
}
が無限ループに入ったので、どうやら成立するみたいだな。
0203１３２人目の素数さん
垢版 |
2020/11/27(金) 08:15:22.17ID:mPM/C8GR
>ID:nTBSTLPq

誤って書き込まれたスレ違いの書き込みに延々と何やってんの？
他スレでとっくに証明されてるよ。

3連休の検査数減少のせいかもしれんけど、そろそろピークアウト
しそうな気配があるな。あと４，５日様子を見ないと判断できないが。
0204１３２人目の素数さん
垢版 |
2020/11/27(金) 08:43:28.41ID:anGa5WFp
>>203
プログラム組んで体感するのが楽しいんだね。
収束することが証明されていてもグラフを書いて収束するのがみたい。
**分布に従うとか言われてもヒストグラムに分布曲線を重ねたくなるのと同じ。
中に何があるかわかっていてもタンクトップをずり下げたくなるようなものｗ
0205１３２人目の素数さん
垢版 |
2020/11/28(土) 23:16:05.02ID:71BcBuYQ
やはりそろそろピークアウトしてきた感じだな。
0206１３２人目の素数さん
垢版 |
2020/12/02(水) 22:30:09.49ID:LbYHBA85
>>205
近隣の職場でクラスター発生したりして、ピークアウトした実感はないな。
むしろ、ますます、近づいてきたという印象。
0208１３２人目の素数さん
垢版 |
2020/12/03(木) 00:17:20.95ID:yHlja/PL
>>206
あんたの近傍だけの話をされてもなぁｗ

現状はほぼ横ばいだろ。ピークに達した感がある。
0209１３２人目の素数さん
垢版 |
2020/12/03(木) 10:09:02.55ID:wLLqGYqv
GOTOトラベルだと
40,000,000 / 172 = 232,558人に一人感染
日本の人口12600万/感染者15万=840人に一人感染
GOTOトラベル参加は新コロナ感染リスクを何分の１に減らすことができるか？
95%信頼区間付きで計算せよ。
0211１３２人目の素数さん
垢版 |
2020/12/03(木) 10:33:15.16ID:yHlja/PL
>>218
数値で示さないと無意味。印象論だったらなんとでも言える。
0212１３２人目の素数さん
垢版 |
2020/12/03(木) 10:36:30.92ID:yHlja/PL
すでに全国的に感染が蔓延している状態なので
GOTOトラベルが感染地域の拡大に寄与してる部分
はあまりないんじゃないかな。
ただ、グループで旅行するとそのグループ内で
感染するということはあるだろうな。
0215１３２人目の素数さん
垢版 |
2020/12/03(木) 18:49:09.99ID:yHlja/PL
>>213
グラフみろよ。現時点では間違いなく横ばい。
実行再生産数もほぼ１。

ただ、これが踊り場ですぐまた上昇に転じる可能性もあるので、
ピークアウトかどうかは単なる印象。
0216１３２人目の素数さん
垢版 |
2020/12/07(月) 22:17:19.54ID:CU+4ckWm
>>215
同業者の職場でも新コロナが発生してきて徐々に身近に迫ってきたなぁと実感している。
0217１３２人目の素数さん
垢版 |
2020/12/08(火) 00:12:15.30ID:hzCu4Q+1
およそ科学的ではないな。
収束しない限り、時間が立てば身の回りで見聞きする
確率があがるのは当たり前。
0218１３２人目の素数さん
垢版 |
2020/12/08(火) 05:39:33.58ID:bQW+eNK4
西浦 博 他4名
新型コロナウイルスと闘うために数学にできること

がkindleからでてるな。アマゾンのリンクが貼れない・・・
0220１３２人目の素数さん
垢版 |
2020/12/09(水) 17:48:25.06ID:nv/1M1Kw
こんな本がでるようだ

理論疫学者・西浦博の挑戦-新型コロナからいのちを守れ! (単行本) (日本語) 単行本 ? 2020/12/8

数学の本ではなさそうな

抜粋とてこんな記述が

・・・僕には脅迫状が届き、生まれて初めて殺害予告を受けました。一番緊迫した頃には、厚労省と新橋のビジネスホテルの間を歩くだけなのに警察の方に護衛してもらったことすらありました。
0221１３２人目の素数さん
垢版 |
2020/12/09(水) 19:58:27.56ID:YBZ1yNkx
>新型コロナからいのちを守れ!

自分の命がやばい、ってか？ｗ
0222１３２人目の素数さん
垢版 |
2020/12/09(水) 20:03:51.35ID:TYzqqYYm
今年のインフルエンザの確率はほぼ
zeroだから、モピロン、
ウイルス性のカゼの症状即ち、
倦怠感かつのどがヘンな感じなら、
99.9％に限りなく近い可能性で
コロナ感染症だと思われる。
多分、知らないうちに
地球人の半分位感染してて
その内に、偽陰性になってるかもね
0223１３２人目の素数さん
垢版 |
2020/12/09(水) 21:16:46.12ID:YBZ1yNkx
そうならいいんだけどね。
0224１３２人目の素数さん
垢版 |
2020/12/10(木) 04:54:10.87ID:LcOGRiI3
西浦さんによるK値の批判はあるけど
逆にK値側から、SIRモデルなど理論疫学で標準的なアプローチへの批判ってあるの？
0225１３２人目の素数さん
垢版 |
2020/12/10(木) 10:21:05.07ID:oBWH18jr
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は定数でなく確率変数
0226１３２人目の素数さん
垢版 |
2020/12/10(木) 11:00:01.89ID:oBWH18jr
続編

んまっ、
ゴンペルツ曲線にしろ、
ロジスティクス曲線にしろ、
I＝S0[k1/(k1-k2)][e^(-k2t)-e^ (-k1t)]
のような曲線にしろ、
無限日後に感染者数の時間微分は、
モチロン、Zeroに限りなく近づく

実無限日後⇒Zeroになるハズ
で安心するだろう

それにしても、このウイルス
感染&免疫回避性能 半端なくヤヴァィ
0227１３２人目の素数さん
垢版 |
2020/12/10(木) 11:20:08.16ID:cxBeQaI7
普通の風邪のコロナウィルスのやっかいさと同じなのでは？
0228１３２人目の素数さん
垢版 |
2020/12/10(木) 11:47:34.80ID:oBWH18jr
消息不明の怪情報だと、
コウモリに感染するが人に感染しない
コロナウイルスだったのが、
数年前に、SARS遺伝子が
そのコロナウイルスに組み込れ、
あのコロナウイルスは人の気道に
感染力を持つようになっとか

現時点で消息不明の情報なハズ
0229１３２人目の素数さん
垢版 |
2020/12/10(木) 12:22:37.08ID:cxBeQaI7
有り得る話だけどな。
武漢のウィルス研究所で作られた可能性はそんなに低くはない。

故意に拡散させた可能性はないとは思うが、事故であれなんで
あれ漏出させた責任を中国政府とるべきだろうな。
0230１３２人目の素数さん
垢版 |
2020/12/10(木) 12:27:53.61ID:oBWH18jr
【新型コロナウイルスは割と安全】
SARS-CoV2は普通ぢゃないコロナだ
SARS-CoV2には亜種がいくつもある
今までのとは
違うから、そもそも新型だ。

新型のは普通のとは異なり、
サイトカインストリームという
免疫暴走の重篤肺炎が
短期間に発生することが珠にある。
表面S-タンパクが、怪しい

これを回避できれば、ま、軽いタダの
インフルエンザ程度だろう
多分。

新型以外コロナウイルス以下に示す

危ないコロナウイルス
SARS-CoV
MERS-CoV

ショボぃ風邪コロナウイルスHCoV-229E
HCoV-OC43
HCoV-NL63
HCoV-HKU1
0231１３２人目の素数さん
垢版 |
2020/12/10(木) 19:24:14.59ID:cxBeQaI7
ストリームじゃなくてストームだろｗ

武漢ウィルスでいいよ。ややこしいから。
0232１３２人目の素数さん
垢版 |
2020/12/10(木) 19:32:32.84ID:oBWH18jr
ご指導ありがとう。
HuhanV.のイチバン怖いのは、
HuhanV.自体というより、
HuhanV.と免疫細胞の相互的な
事象で、IL6が宿主が発射され
宿主により、宿主を破壊することに
ある。肺炎により死亡する
ここがヤバい
0235１３２人目の素数さん
垢版 |
2020/12/12(土) 21:15:11.20ID:xixiDIl1
ウィルス干渉説が本当なら
　インフル流行させればコロナ殲滅できる
はあり？
0237１３２人目の素数さん
垢版 |
2020/12/15(火) 21:55:06.09ID:SiGqyW5Y
コロナが流行してるのは、
インフルエンザが流行してないから
とのトンデモが現れる予感がしてきた
0238１３２人目の素数さん
垢版 |
2020/12/15(火) 23:48:34.33ID:GfZkwHJf
新型コロナよりインフルエンザの方が感染力高いことはすでにわかっている
0240１３２人目の素数さん
垢版 |
2020/12/22(火) 00:05:01.18ID:kYmv3i0r
やはりカッパ値はかんさいの風土やまい
0241１３２人目の素数さん
垢版 |
2020/12/24(木) 16:38:01.60ID:nYQkp2Hz
東京 日　 月　 火　 水　 木　 金　 土　
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日の感染者数にベンフォードの法則*)が成り立つかを試したみた。

先頭の数字の頻度は

https://i.imgur.com/DN3Fcr6.png

x
1 2 3 4 5 6 7 8 9
59 54 20 11 15 7 7 5 2

成立しているようにみえる。

*)
https://ja.wikipedia.org/wiki/%E3%83%99%E3%83%B3%E3%83%95%E3%82%A9%E3%83%BC%E3%83%89%E3%81%AE%E6%B3%95%E5%89%87
ベンフォードの法則（ベンフォードのほうそく、Benford's law）とは、自然界に出てくる多くの（全てのではない）数値の最初の桁の分布が、一様ではなく、ある特定の分布になっている、という法則である。
0242１３２人目の素数さん
垢版 |
2020/12/25(金) 16:14:51.97ID:Mak8FdjO
>>241
しかしχ^2検定でp=0.000356。
まあ、時系列値だから仕方ないか。
0243１３２人目の素数さん
垢版 |
2020/12/25(金) 16:48:59.06ID:mzF7SS2F
>>241
ペンフォードの法則が成立するという事は
人為的な操作がされていなそうという解釈でいいんだよな？
0244１３２人目の素数さん
垢版 |
2020/12/25(金) 18:01:34.98ID:SYK9UM4r
ついに数秘術か？
0246１３２人目の素数さん
垢版 |
2020/12/27(日) 18:28:07.55ID:wgykouqU
>>241
横軸に日曜日の陽性者数、縦軸に1週間の総数をグラフにしてみた。
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人を超えるのは金曜日と出た。

来週になったら、どのくらいあたったか検討してみよう。
0249１３２人目の素数さん
垢版 |
2020/12/27(日) 23:05:55.14ID:wgykouqU
>>241(今月→今週に訂正)

横軸に日曜日の陽性者数、縦軸に1週間の総数をグラフにしてみた。
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人を超えるのは金曜日と出た。

来週になったら、どのくらいあたったか検討してみよう。
0251１３２人目の素数さん
垢版 |
2020/12/29(火) 00:17:44.79ID:9f4Zs7Qm
今日は481人か

月 fit lwr upr
1 507.6391 458.0591 557.219
一応、95%信頼区間458-557には入っているな。

明日の予想
火 fit lwr upr
1 593.0795 489.3003 696.8586
0252１３２人目の素数さん
垢版 |
2020/12/29(火) 12:42:25.83ID:BqqMc2cN
目的論的解釈からウイルスは生存と拡散を目的に弱毒化する方に変異すると言われるんだが、
すでに無症状の患者がウイルスをばらまくように変異しているわけだから、これ以上弱毒化する方向に変異するのは考えにくいと思う。

つぎの変異の進化は感染拡大するように宿主の行動を制御することだな。
カタツムリの寄生虫ロイコクロディウムみたいに。

即ち、感染すると宿主の活動性が増して外出しやすくなるようにする。
この変異の進化はすでに起きているかもしれんな。
GOTOキャンペーン利用者ってまさにそれかも。
0253１３２人目の素数さん
垢版 |
2020/12/29(火) 16:50:37.14ID:BqqMc2cN
>>251
東京都 新型コロナ 856人感染確認 火曜日の発表として最多更新. 2020年12月29日 16時08分 新型コロナ 国内感染者数

線形回帰での予想数の95%信頼区間を越える増加だなぁ。
0254１３２人目の素数さん
垢版 |
2020/12/29(火) 18:14:50.45ID:9f4Zs7Qm
>>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
0257１３２人目の素数さん
垢版 |
2020/12/30(水) 00:32:41.76ID:uEoEZYjH
>>254
感染モデルもいろいろあるが、どれもたいがい指数関数的変化
が卓越するので、線形予想じゃだめなんじゃね？
0259１３２人目の素数さん
垢版 |
2020/12/30(水) 07:58:53.92ID:7I9M4Lp6
>>256
無職で生活できたらいいなぁ。

最近は基幹病院ですら救急搬送を断るのでうちみたいな個人病院への救急搬送が増えた。
今月は救急受け入れの歩合給だけで10万円を越えたのでクリスマスにT-boneステーキを買って家族で食べた。
まあ、一見さんの発熱は断れと院長から指示が出ているから何でもかんでも救急を受けなくていい。
0260１３２人目の素数さん
垢版 |
2020/12/30(水) 08:03:50.41ID:7I9M4Lp6
>>257
レスありがとうございます。
データを対数化して線形回帰した方がいいということでしょうか？

素のデータでやった>207でよくあたっていたから自分でもびっくりした。
0263１３２人目の素数さん
垢版 |
2020/12/30(水) 14:12:31.01ID:uEoEZYjH
>>259
みんなが貰えた給付金10万円に追加で医療従事者の給付金ってのもあったよね。
0264１３２人目の素数さん
垢版 |
2020/12/30(水) 14:13:11.21ID:uEoEZYjH
>>260
そゆこと。
0265１３２人目の素数さん
垢版 |
2020/12/30(水) 15:41:16.79ID:yK/NkJWh
>>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
0266１３２人目の素数さん
垢版 |
2020/12/30(水) 15:47:18.46ID:yK/NkJWh
>>263
うちは、コロナ非対応病院だから、給付金５万円。
夏場に防護服着ての内視鏡検査は体力を消耗した。
0267１３２人目の素数さん
垢版 |
2020/12/30(水) 16:38:25.36ID:uEoEZYjH
そうなんだ。
俺の友人はコロナ対応病院なので10万もらってたが、
本人は無関係な部局でいつもと変わらぬ仕事だった
ので、気が引けると言ってた。
0269１３２人目の素数さん
垢版 |
2020/12/30(水) 17:25:06.68ID:uEoEZYjH
曜日ごとのばらつきが大きいから一週間平均で予測したほうがいい。
ここ3週間は4週間で倍増ペースの比較的ゆるやかな指数関数的増加
をしている。明日は1000人超えそう。
0271１３２人目の素数さん
垢版 |
2020/12/30(水) 20:12:03.34ID:yK/NkJWh
>>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]の予想。
0272１３２人目の素数さん
垢版 |
2020/12/30(水) 23:52:01.76ID:uEoEZYjH
>>271
12/3からの4週間分の木曜データを使って、エクセルで指数近似させると
明日の予測値は1100と出た。線形近似だと1033。ってか、日毎の予想
なんかしても意味ないよ。曜日ごとの変動が大きすぎる。
0273１３２人目の素数さん
垢版 |
2020/12/31(木) 00:06:28.13ID:qiffllOG
一週間の移動平均値をデータとして12/28までの19日分を使うと、
線形近似で1000人を超えるのは来年の1/14で、指数近似だとそれ
が1/9に早まる。
0274１３２人目の素数さん
垢版 |
2020/12/31(木) 00:54:32.88ID:ZNO2TkSk
そもそも大本営発表みたいなデータが信用できないから予測もへったくれもないと思いませんか？
0275１３２人目の素数さん
垢版 |
2020/12/31(木) 11:09:12.91ID:qiffllOG
絶対値ではなく、変動量を見る目安として使うことに問題ない。
0277１３２人目の素数さん
垢版 |
2020/12/31(木) 12:52:20.55ID:uj9CqJOS
>>250のグラフ見てみたので、
感想文を書いてみたぁぁ。

新規、新規の、新規感染者数が
時間経過で、ほぼ線形的増加だ。
ぢゃ、入院患者も線形的増加となるか
自問自答してみた。

入院期間は定数　2week　とおくと
d入院数/dt = k*t ∵kは定数∵線形増加
d退院数/dt = k*(t-2)
∴d入院患者/dt = 2*k = K　（Kも定数)

入院患者 = K*t^2 (積分定数は無視）
どうやら、二次関数的に増加するハズ

もし、新規感染者が指数関数的なら
さらに不味いけどどうにかなるだろう
0278１３２人目の素数さん
垢版 |
2020/12/31(木) 13:01:43.05ID:uj9CqJOS
>>277の訂正
あれ、間違えた　改訂
入院患者数は、
✕　二次関数的に増加する
○　線形的増加に増加する
でも、その内に無限大になるな
でも、誰かどうにかするだろう
0279１３２人目の素数さん
垢版 |
2020/12/31(木) 16:55:30.49ID:qiffllOG
>>276
馬鹿は鼻くそでも食ってろ。
0280１３２人目の素数さん
垢版 |
2021/01/01(金) 01:15:15.76ID:REQb9+Ix
>>273
直近一週間ぐらいで増加率が大きくなってるな。
このままいくと、1/4には7日移動平均が1000人
を超える。20日で倍増のペース。
0282１３２人目の素数さん
垢版 |
2021/01/02(土) 20:43:48.17ID:1Rj7ie4D
>>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

よい相関関係にあるな。
0286１３２人目の素数さん
垢版 |
2021/01/02(土) 23:30:14.34ID:1Rj7ie4D
日曜日の数と週の総数の関係だから対数をとっても
相関係数はむしろ小さくなった。
明日の総数から来週の総数を予想してみるかな。
0287１３２人目の素数さん
垢版 |
2021/01/02(土) 23:53:40.12ID:+iwtte4V
短期変動予測なら二次以上の項を無視した一次関数近似で
全然問題ないと思うけど、数週間先の予測だと合わなくな
ると思うよ。
0288１３２人目の素数さん
垢版 |
2021/01/03(日) 10:50:39.85ID:8tLYm46h
>>287
レス、ありがとうございます。
一次関数でなくて、ラグランジュの補完多項式（外挿に使うのは邪道ではあるが）を使って近傍なら近似できるかと思って
最近2週間のデータを使って13次関数を作ってみた。

https://i.imgur.com/oxGa7ag.png

これで今日(day 15)の値を外挿してみると、
　627185
になって、全く役に立たないことが判明したｗ
0289１３２人目の素数さん
垢版 |
2021/01/03(日) 11:01:58.83ID:/BKPeYSh
モデル関数を適当になんの理論的根拠もなく選んでも意味はない

・ある程度の理論的根拠を持ってモデル関数を選ぶ
・実際そのモデル関数の値と実測値がどれくらい離れてるか実測する

の長い長い繰り返しでモデル関数として何をとれば実用性が高いのか積み重ねていく
その地道な繰り返し
それを続けてきた何百年もの研究結果をまるで無視しては意味ない
しかし実際に防疫学の専門家にでもなってフィールドワークやらないとそんなの無理
数学しかできない人間が机の上だけでなんかやっても意味ない
0290１３２人目の素数さん
垢版 |
2021/01/03(日) 11:21:15.75ID:8k6lQiPc
>>288
ばらつきの大きな統計量を近似するのにパラメータを増やしていく
のはいかがなものか。補間多項式を使うなんて論外でしょ。

>>289が言うように、あらかじめ理論的なモデルを想定するか、
パラメータを最小限にとどめた初等関数を使うのが筋なのでは？
0291１３２人目の素数さん
垢版 |
2021/01/03(日) 14:12:48.62ID:8tLYm46h
フィボナッチ数列15個からラグランジュの補完式を求めてグラフにしてみたら、
　https://i.imgur.com/LLS3Xnw.png
当然、解離しているのだが、16個目は割と近似していたので、新型コロナの陽性数でやったら翌日の近似ができるだろうか、計算させてみたが失敗に終わるの巻。
0292１３２人目の素数さん
垢版 |
2021/01/03(日) 16:36:07.92ID:8tLYm46h
>>286
日曜日最多の816だったので、今週の感染者総数の予想とその95%信頼区間

> predict(lm,new,interval = 'confidence')
fit lwr upr
1 6680.85 6139.249 7222.451
0293１３２人目の素数さん
垢版 |
2021/01/03(日) 16:47:49.16ID:8tLYm46h
あんまり、当たらない曜日毎の予想（まあ、競馬の予想程度に思っていただきたい）
次の1000人超えは水曜日の予想。

日 fit lwr upr
816 816 816

月 fit lwr upr
544 493 595

火 fit lwr upr
798 682 913

水 fit lwr upr
1066 935 1197

木 fit lwr upr
1301 1135 1466

金 fit lwr upr
1048 919 1177

土 fit lwr upr
1109 981 1236

相関係数は、
> for(i in 2:7) cat(day[i],' : ',summary(LM[[i]])$r.squared,'\n')
月 : 0.9086839
火 : 0.8026629
水 : 0.8646371
木 : 0.8475512
金 : 0.8404378
土 : 0.8565624
0296１３２人目の素数さん
垢版 |
2021/01/03(日) 19:17:01.41ID:Ees9cUhc
>>288
過剰適応、オーバーフィッティングになっているんだろうと思う
標本誤差=経験誤差を小さくすると汎化誤差も小さくなるけど
過剰に経験誤差を小さくすると汎化誤差は大きくなってしまう
機械学習辺りではよく知られていること
0297１３２人目の素数さん
垢版 |
2021/01/03(日) 19:19:38.23ID:Ees9cUhc
機械学習辺りでは標本データの空間から
一度より高次元空間に写して
その高次元空間で分類や回帰をして得たい結果を得る
ような事をする方法がある
カーネル法とか言われていると思う
0298１３２人目の素数さん
垢版 |
2021/01/03(日) 19:45:10.94ID:8k6lQiPc
本質的にそういうことが問題じゃないんでないの？
0299１３２人目の素数さん
垢版 |
2021/01/03(日) 21:14:43.09ID:Ees9cUhc
本質的なモデルを導き出すためにデータを基にしているんだろう
データから事情を説明するモデルを導き出す
そのときにデータに過剰に引っ張られて本質からずれたモデルになるのが過剰適応だと思う

本質的なモデルに必要な要因、変数を測定していないから欠けているかも知れないけどな
ウイルスの特性、人間の接触状況、感染者発生数、感染対策方法とかの他に必要な情報が欠けていたりするかも
0300１３２人目の素数さん
垢版 |
2021/01/04(月) 00:04:34.06ID:oaJLG6Ed
>データから事情を説明するモデルを導き出す

そういう帰納的な方法でうまくいくとは思えんな。
複雑かつ統計的バラツキが大きいシステムだからね。
仮設から演繹的にモデルを組み立てて、パラメータを
データに合わせるべきじゃね？
0301１３２人目の素数さん
垢版 |
2021/01/04(月) 00:05:07.84ID:oaJLG6Ed
仮説
0303１３２人目の素数さん
垢版 |
2021/01/04(月) 11:08:48.44ID:oaJLG6Ed
統計的ばらつきが大きいデータを使うという大前提を
忘れちゃだめだよ。
補間多項式の外挿とかまったく論外。話にならん。
オーバーフィッティングとかいうレベルの話じゃない。
0304１３２人目の素数さん
垢版 |
2021/01/04(月) 11:16:14.04ID:/qDFBZT4
せやね
そもそもが「いかにパラメータの少ないモデル曲線群を選んでフィットさせていくか？」が問題なのにフィットさせるためにパラメータを増やして行ってどうするって話
K値の「後から後から波を追加してパラメータ増やして行く」のと同じ
それでフィットしてもなんの意味もない
0305１３２人目の素数さん
垢版 |
2021/01/04(月) 15:39:14.09ID:oaJLG6Ed
東京都のここ３週間の7日間移動平均のデータをみると、
12/31が大きな外れ値になってるけど、25日で２倍になる
指数関数的増大によく一致してる。

細かく見ると、12/9-12/21の間も倍加時間25日の指数関数
でよく一致してるが、12/22-12/31の間は倍加時間が18日
と急増し、1/1-1/3は横ばいになってる。

12/31あたりがピークであってくれればいいんだが、、、
0306１３２人目の素数さん
垢版 |
2021/01/04(月) 15:53:29.56ID:oaJLG6Ed
いや、年明けの数値は年末年始の検査人数の減少による
蜃気楼かもしれんな。だとすると、明日当たりから急増
傾向になっても不思議はない。
0307１３２人目の素数さん
垢版 |
2021/01/04(月) 16:43:04.89ID:Q6XIKg7u
>>304
ベイズ統計の考え方だと新しいデータを得られたら分布を更新する

モデルの形が違っても結果的に同じ分布を得られたら問題無いのでは？
例えば、あるモデルは別のモデルの特殊な形だったり
0308１３２人目の素数さん
垢版 |
2021/01/04(月) 16:51:59.30ID:/qDFBZT4
>>307
アホか
どんどんデータ更新されて行くなら理論予言性ゼロやろ？
それでK値も批判されてるんやろがな？
上の方でもおんなじ事言われてるやん
0309１３２人目の素数さん
垢版 |
2021/01/04(月) 17:08:39.81ID:Q6XIKg7u
ニューラルネットワークはどんな関数も表現できると言われている
ある関数はニューラルネットワークのパラメータの特殊な場合とみなせる
データを元に学習したらどんな関数も表現できるけど
経験誤差を小さくすると汎化誤差が大きくなる問題がある
0310１３２人目の素数さん
垢版 |
2021/01/04(月) 17:11:25.86ID:Q6XIKg7u
>>308
データ数100で推定したパラメータと
データ数1000で推定したパラメータと
どちらの方を真の値として採用したい？

標本数は関係するんじゃないか？
あとそれまで測定していなかった変数を測定したら結果に影響することがわかったらモデルを更新するだろう
0311１３２人目の素数さん
垢版 |
2021/01/04(月) 17:22:58.41ID:/qDFBZT4
>>310
そんな事言ってるんじゃないよ
統計量は多いほどいい
しかしパラメータは少ないほどいいんだよ
そもそもパラメータとかわかってないやろ？
統計の教科書読まないで統計の話についてこれるわけないやん？
0312１３２人目の素数さん
垢版 |
2021/01/04(月) 17:44:41.95ID:Q6XIKg7u
>>311
母数だろう
正規分布だと平均と分散

分布を記述するのに必要最小限のパラメータ数はどうやって求めるの？
0313１３２人目の素数さん
垢版 |
2021/01/04(月) 18:06:24.22ID:/qDFBZT4
>>312
そんなもん理詰めだけで決まるわけないやん
たくさんの統計量でパラメータ減らせば信頼度は高まるけどモデルを単純化し過ぎてるとそれで合わなくなるし、考慮するパラメータを増やせば数値は“フィット”するようになるけど逆に表現力が増え過ぎて理論予言性が無くなる
当たり前
じゃあどのくらいがいいのか、なんて実際のフィールドワークやってみんなで考えながらいい塩梅のとこを探って行く
回復率と感染率と言う二つのパラメータのモデルもそう言う長い議論で磨かれてきた理論
しかしさりとてベストではない感染症もあるやろ
そんなもん数学の机の上の議論だけで決まるはずがない
0314１３２人目の素数さん
垢版 |
2021/01/04(月) 18:12:47.10ID:m5nY/ADt
金融トレードのテクニカル分析でもテクニカル指標は２〜３、その指標のパラメーターも少なけらば少ない程ロバストとされてるな
数が多くなるとカーブフィッティングして過去の再現性はいいけど将来の期待値は１を下回るケースが多い
0315１３２人目の素数さん
垢版 |
2021/01/04(月) 18:27:00.11ID:jIvyoIrl
AICとかが参考になるんだろうか？機械学習だとデータを学習用と評価用に分けるとかするけど
0317１３２人目の素数さん
垢版 |
2021/01/04(月) 19:37:44.90ID:Q6XIKg7u
>>315
学習用と評価用に分けることで
学習用から経験誤差を推定する
評価用から汎化誤差を推定する

経験誤差を小さくすることで汎化誤差も小さくなるけど
過剰に経験誤差を小さくすると汎化誤差が大きくなる
→オーバーフィッティング、過剰適応、過剰学習みたいに言われる
0318１３２人目の素数さん
垢版 |
2021/01/04(月) 19:42:33.23ID:/sna+/0b
情報量基準って日本人が創始者になった分野だから
その分日本では過大評価されていると思う
0319１３２人目の素数さん
垢版 |
2021/01/04(月) 20:57:46.22ID:/qDFBZT4
>>316
だから更新できる＝理論予言性がないって事やろ
お前以外全員わかってでお前だけ話通じてないのわかる？
それがわかるならアレ？なんでだろとか思わん？
0320１３２人目の素数さん
垢版 |
2021/01/04(月) 21:17:21.48ID:m0kvaJBH
深い話どね。数学に答えは見つからない
【オッカムの剃刀】　これだな。さて、

何も知り得ない地球生命体でも、
敢えて未来を統計的に想像したければ、
定数パラメータは一つKが良かろぅ
定数と一次回帰と二次回帰をそれぞれ
40％、40％、20％の重みでブレンド
で最新の4weekのデータから
4week後未来なら、ソコソコ予測出来そう

y = 0.4K + 0.4K*x + 0.2K*x^2
に残差自乗和が最小となるKを
コンピュータで無理やり求めるとかね

しかし例え、
その予測が当たったとしても人々は
その予測に基づき行動し、人々の
その予測は、ハズレる向きに行動する

これにより、変動はランダムというより
カオスな予測不可能な動きとなろう

哲学、再帰性理論に通ずる。
未来予測は数学でなく霊感が必要となろう。　悲しいことに必要不十分であるが
0321１３２人目の素数さん
垢版 |
2021/01/04(月) 21:49:53.58ID:oaJLG6Ed
現状のまま感染が広がれば、一ヶ月たらずで新規感染者が倍になる。

という程度の予測は立てないと、対策の手のうちように困るでしょ。
先回りした対策というのはそういうこと。そのために従来から使われ
ているモデルが採用されている。デタラメな思いつきではなく。
0324１３２人目の素数さん
垢版 |
2021/01/05(火) 08:04:54.26ID:+nklPzTT
先週までのデータを使ってどの曜日での総数予想があたるかをみようと
線形回帰での補正決定係数adjusted. r squaredを出してみた。

Sun Mon Tue Wed Thu Fri Sat
0.9267 0.8935 0.8856 0.9656 0.9417 0.9235 0.9394
水曜日のデータが一番、回帰の当てはまりがいいようだ。


ついでに曜日間の補正決定係数も出してみた。

Sun Mon Tue Wed Thu Fri Sat Sum
Sun 1 - - - - - - -
Mon 0.905 1 - - - - - -
Tue 0.7948 0.6906 1 - - - - -
Wed 0.8592 0.8232 0.8718 1 - - - -
Thu 0.8415 0.7614 0.9265 0.929 1 - - -
Fri 0.8341 0.8686 0.7049 0.8549 0.7829 1 - -
Sat 0.8508 0.8705 0.7269 0.8826 0.8064 0.9694 1 -
Sum 0.9267 0.8935 0.8856 0.9656 0.9417 0.9235 0.9394 1
0325１３２人目の素数さん
垢版 |
2021/01/05(火) 08:39:31.09ID:+nklPzTT
おまけ、今週の日月のデータから週の総数を線形回帰で予想してみる。


日曜日データから週の総数の予想
fit lwr upr
1 6680.85 6139.249 7222.451

月曜日データから週の総数の予想
fit lwr upr
1 10283.86 9127.747 11439.97
0326１３２人目の素数さん
垢版 |
2021/01/05(火) 09:22:22.44ID:gyHss/tq
ついでに曜日毎の予想を計算するように拡張して、今日、火曜日の予想をしてみた。

日曜日データから火曜日の予想
fit lwr upr
1 797.5329 681.8693 913.1964

月曜日データから火曜日の予想
fit lwr upr
1 1178.019 926.6623 1429.376

>324の表現を変えると

Sun Mon Tue Wed Thu Fri Sat Sum
Sun 1 0.905 0.795 0.859 0.841 0.834 0.851 0.927
Mon 0.905 1 0.691 0.823 0.761 0.869 0.871 0.893
Tue 0.795 0.691 1 0.872 0.927 0.705 0.727 0.886
Wed 0.859 0.823 0.872 1 0.929 0.855 0.883 0.966
Thu 0.841 0.761 0.927 0.929 1 0.783 0.806 0.942
Fri 0.834 0.869 0.705 0.855 0.783 1 0.969 0.923
Sat 0.851 0.871 0.727 0.883 0.806 0.969 1 0.939
Sum 0.927 0.893 0.886 0.966 0.942 0.923 0.939 1

どの曜日を使っても火曜日の値との相関はよくないみたいだな。
0327１３２人目の素数さん
垢版 |
2021/01/05(火) 09:28:07.07ID:gyHss/tq
>>324
漢字の方がみやすいな（泥タブと違って文字化けしなかった）

日 月 火 水 木 金 土 計
日 1 0.905 0.795 0.859 0.841 0.834 0.851 0.927
月 0.905 1 0.691 0.823 0.761 0.869 0.871 0.893
火 0.795 0.691 1 0.872 0.927 0.705 0.727 0.886
水 0.859 0.823 0.872 1 0.929 0.855 0.883 0.966
木 0.841 0.761 0.927 0.929 1 0.783 0.806 0.942
金 0.834 0.869 0.705 0.855 0.783 1 0.969 0.923
土 0.851 0.871 0.727 0.883 0.806 0.969 1 0.939
計 0.927 0.893 0.886 0.966 0.942 0.923 0.939 1
0328１３２人目の素数さん
垢版 |
2021/01/05(火) 12:19:03.22ID:KNfoLQVz
一連の書き込み、まったく意味ないと思うんだが。 >>327
スレが汚くなるだけなので、ほどほどに願います。
0329１３２人目の素数さん
垢版 |
2021/01/05(火) 12:20:39.66ID:KNfoLQVz
意味ないっていうのは、曜日以外にも休日かどうかという
条件が大きくデータを左右するはずなので。
0330１３２人目の素数さん
垢版 |
2021/01/05(火) 12:35:19.89ID:h+VbaBvY
緊急事態行動ωωω
言葉遊びは意味が無い
watch?v=CVOTJanrRrg
0331１３２人目の素数さん
垢版 |
2021/01/05(火) 13:24:31.42ID:+nklPzTT
>>329
＞曜日以外にも休日かどうかという
＞大きくデータを左右するはず
数値化よろしく。
0332１３２人目の素数さん
垢版 |
2021/01/05(火) 13:34:27.50ID:+nklPzTT
>>311
一概には言えないとおもう。

長方形の全周長と面積の関係を統計モデルで示すという
a bit silly but makes an important point
の解説が、統計の教科書に書いてある。

データをいくら増やしても
lurking variable(ここではアスペクト比)に気づかない統計解析が無意味だという好例。

乱数発生させて、実験したみた結果がこれ。

http://i.imgur.com/TiYDx6i.png
青線は正方形のときの値でこれが上限になっている。
0333１３２人目の素数さん
垢版 |
2021/01/05(火) 16:28:23.79ID:VzBiaGAZ
>>332
話にのポイントはそこやないって言ってんのに？
なんで他の人全員がわかってるのに自分だけが話から置いてかれてるのをおかしいと思えないんかね？
だからいつまで経ってもダメなんだよ
0334１３２人目の素数さん
垢版 |
2021/01/05(火) 17:13:39.24ID:KNfoLQVz
>>331
できりゃぁ誰も苦労しない。
本質的にできないことをやろうとしても意味がないと言ってる。
0335１３２人目の素数さん
垢版 |
2021/01/05(火) 17:16:36.87ID:+nklPzTT
>>333
何か間違っているか？

さて、明日の予想。

日曜日データから水曜日の予想
fit lwr upr
1 1065.917 935.0589 1196.775

月曜日データから水曜日の予想
fit lwr upr
1 1667.556 1408.003 1927.11

火曜日データから水曜日の予想
fit lwr upr
1 1604.5 1398.906 1810.094
0336１３２人目の素数さん
垢版 |
2021/01/05(火) 17:25:58.49ID:KNfoLQVz
そんな予想はまったく無意味だと何度言えばわかるのか >>335

過去の気温データだけから明日の気温を予想するようなもの。
馬鹿馬鹿しくて話にならん。
0337１３２人目の素数さん
垢版 |
2021/01/05(火) 17:41:25.32ID:VzBiaGAZ
>>333
もちろん間違ってる
と言うか自分以外の全ての人がオレの話理解できてて自分1人取り残されてるのが何故わからん？
0338１３２人目の素数さん
垢版 |
2021/01/05(火) 17:49:23.61ID:+nklPzTT
>>336
いや、第三波は来ないという予想が堂々と流布されていてそれを信奉していた人がいるわけだし。
体育の日は、過去のデータから晴れが多い日として選ばれたとかあるからね。
0340１３２人目の素数さん
垢版 |
2021/01/05(火) 18:45:09.59ID:KNfoLQVz
>>338
第三波が来ないなんて予測はこれっぽっちも聞いたことないが、
たとえそんな予測があったとしても、あなたの超短期予報が
まったく無意味だということとは関係ないよ。
0341１３２人目の素数さん
垢版 |
2021/01/05(火) 18:52:20.43ID:VzBiaGAZ
単に勉強不足で素頭は悪くないわけではなく、そもそも素頭が悪いんだな
ここまで引っ張ってまだ他の人が何言ってるか理解出来ないんだから
0343１３２人目の素数さん
垢版 |
2021/01/05(火) 20:40:23.19ID:VzBiaGAZ
>>342
わからないかなぁ？
君以外全員理解できて同意してるレスついてるやん？
意見もへったくれもない
数学なんだから
わかるかわからないかしかない
0345１３２人目の素数さん
垢版 |
2021/01/05(火) 21:01:25.69ID:VzBiaGAZ
>>344
さよけ
ならいいわ
ともかく数学的に考えて全然意味ないことやっててそれを他人に指摘されて全く理解できない
もうウリュウは数学から足洗った方が身のためやと思う
全く才能がない
0346１３２人目の素数さん
垢版 |
2021/01/05(火) 22:17:05.93ID:+nklPzTT
>>340
>207に書いたけど自分でも驚くほどあたったんだなぁ。
それで２匹めの泥鰌狙い。
0347１３２人目の素数さん
垢版 |
2021/01/05(火) 22:49:58.87ID:+nklPzTT
>>343
種々の 交絡因子を組み込んでいるであろうグーグルの予想もかなり外しているし、単純な線形回帰とさほど差はないようにみえるな。
0349１３２人目の素数さん
垢版 |
2021/01/06(水) 00:34:20.36ID:vePATSvQ
>>347
いまのところ、ここ一ヶ月くらい続いている倍加
時間が約25日の指数関数的増加でだいたい合ってる。
このくらいのスパンをとると1次関数近似よりは
指数関数の方がよくフィットしてることがわかるよ。

短期間なら線形近似でわりとよく合うのは、テイラー
展開の一次の項を見てると思えば当たり前の話かと。
0350１３２人目の素数さん
垢版 |
2021/01/06(水) 09:57:33.04ID:416eUeUb
過去の実績から未来を予想するのは普通するだろ。
たとえばこういう問題。

ゴルゴ13の狙撃実績は100/100で依頼するには報酬は１億とする。
ゴルゴ14の狙撃実績は10/10であるという、ゴルゴ14への報酬はいくらが適切といえるか？
狙撃成功率100%だからゴルゴ14に報酬１億を払うのが適切といえるか？
0351１３２人目の素数さん
垢版 |
2021/01/06(水) 10:28:44.34ID:vePATSvQ
>>350
ゴルゴ13は３０年ぶりで仕事をする、現在の年齢は８０歳。
一方、ゴルゴ14は現役ばりばりで、それはこの１年間の実績。

さあ、どうする？
0354１３２人目の素数さん
垢版 |
2021/01/06(水) 12:23:23.76ID:vePATSvQ
>>353
アンカー間違いか？俺は医者じゃないよ。

明らかに条件が違う事がわかってるのに、手持ちのデータだけを
比べて確率を議論しても意味がないってこと。
感染者数のミクロな変化を推計する無意味さを諭してる。
0355１３２人目の素数さん
垢版 |
2021/01/06(水) 13:57:34.49ID:CkBfrSdp
>>334
祝日を交絡因子に組み込んで回帰するだけじゃねぇの？
想像だが、Googleの予想はそれくらいは組み込んでいるだろうと思う。
それで予想数が当たる確率が増すとは思えないから手間をかける気にならん。
0356１３２人目の素数さん
垢版 |
2021/01/06(水) 13:59:32.44ID:CkBfrSdp
>>351
医師板ででた、模範解答はこれ、
　＞そんなもん、ゴルゴ14にいくら欲しんか聞いたらよろしんでんがな。
0357１３２人目の素数さん
垢版 |
2021/01/06(水) 14:06:24.64ID:CkBfrSdp
>327に依拠すると、
水曜日のデータは日月火から計算したのでは火曜日の決定係数が0.872と一番大きいので今日の予想は火曜日の回帰を採用して

火曜日データから水曜日の予想
fit lwr upr
1 1604.5 1398.906 1810.094

今日は１６００くらいと予想。
0358１３２人目の素数さん
垢版 |
2021/01/06(水) 14:53:53.98ID:vePATSvQ
>>356
言い値が適切かどうかどう判断する？
0359１３２人目の素数さん
垢版 |
2021/01/06(水) 14:54:29.84ID:vePATSvQ
競馬の予想のほうがマシだなｗ
0360１３２人目の素数さん
垢版 |
2021/01/06(水) 15:24:41.71ID:vePATSvQ
>>273,280
７日間移動平均が1000人を越えた。
結局、指数関数近似で中間の1/6だったか。
0361１３２人目の素数さん
垢版 |
2021/01/06(水) 15:25:11.92ID:Iue2otBN
とりあえず陽性者数の上限を言ってみてくださいよ
PCRのサーマルサイクラーの台数と稼働時間も有限でしょ？
そんなの専門外でめんどいからわからない？知らない？
だとすればやっぱお前らクソの役にも立たないわ
0362１３２人目の素数さん
垢版 |
2021/01/06(水) 16:08:21.09ID:CkBfrSdp
>>357
東京都 新型コロナ 1591人の感染確認 過去最多に

2021年1月6日 15時05分

予想が
fit lwr upr
1604.5 1398.906 1810.094
だったから、わりと近似したな。



明日の予想

決定係数が高い水曜日データを使って
木曜日の予想
fit lwr upr
1881.525 1712.185 2050.864
0365１３２人目の素数さん
垢版 |
2021/01/07(木) 08:55:41.52ID:+xnsY4vR
ペイズの理論はおろか普通の統計の教科書すら読んだこともない
そもそも確率論の教科書すら1ミリも知らない
0367１３２人目の素数さん
垢版 |
2021/01/07(木) 09:58:28.42ID:KvxqEhz+
>>366
医師免許の画像はネットに溢れているらしいから
医師会のネクタイピンの画像を開業医スレにあげておいたよ。
0369１３２人目の素数さん
垢版 |
2021/01/07(木) 10:26:51.46ID:TDko3Vvz
>>367
それが証明になると思ってるのが笑える。どっかから拝借してきた奴だろ？姑息な真似するね。その発想がまさしく非医の証拠。
0370１３２人目の素数さん
垢版 |
2021/01/07(木) 10:46:21.48ID:KvxqEhz+
週内変動の影響を除去するのが７日間移動平均なら
逆に週内変動を利用して曜日毎の回帰から短期予測するのもあっていいと思うね。
0371１３２人目の素数さん
垢版 |
2021/01/07(木) 15:10:48.87ID:n6NiS4+G
＜新型コロナ＞東京都で新たに2447人が感染　過去最多を大幅更新

線形回帰から随分外れたなぁ。
単純な外挿を上回る勢いでコロナ蔓延中ってことになるなぁ。
0372１３２人目の素数さん
垢版 |
2021/01/07(木) 15:11:49.03ID:n6NiS4+G
>>369
どこから拝借してくるって、盗人の発想だね。
いつもそういうことやっているから思いつくのだろうな。
0373１３２人目の素数さん
垢版 |
2021/01/07(木) 15:55:30.69ID:TDko3Vvz
>>372
ネクタイピンが証明になるとでも思ってるのか？笑
さすが非医だね。となると本人の物じゃないから、どこかから拝借したものになる。つまりお前が盗人。
ずいぶんこれ見よがしに挙げたけど本人の物という証拠はあるの？
0374１３２人目の素数さん
垢版 |
2021/01/07(木) 16:16:08.70ID:SX2cder4
>>370
>週内変動の影響を除去するのが７日間移動平均なら

そういう理由もあるけど、一週間より短いタイムスケール
で論じる意味がないからってこともあるでしょうな。
潜伏期間が一週間前後と長いし、幅もあるので、数週間単位
での変化予測ができれば十分。
0375１３２人目の素数さん
垢版 |
2021/01/07(木) 17:06:40.97ID:SX2cder4
ここで東京の検査人数を見ると、毎週土日に平日より
格段に減るのが分かる。これが日曜、月曜の陽性者数
の値を左右している。
ttps://toyokeizai.net/sp/visual/tko/covid19/

年末年始にはこの検査数減少が4日間続いた結果、1/1-1/4
の数字が低くなったと思われる。検査数が復活した1/5以降の
3日間だけをみると一週間で倍増という急激なペースで拡大を
していることになり、緊急宣言の必要性が理解できよう。

年末の酒盛りが原因かもしれんね。
0376１３２人目の素数さん
垢版 |
2021/01/07(木) 22:04:11.46ID:Oj9HGHE3
それでいつまで感染者増えるの
カーブフィッティングの予想が近いだの外れたのなんて
バカでも言えるわなあ
0377１３２人目の素数さん
垢版 |
2021/01/08(金) 00:55:43.48ID:pBLL1e65
増え方がどう変化するかだろうな。

現状では年末の飲み会で相当増えてる上に、年末年始で
検査が後回しになってた分がサージしてるかも。

それが一段落して増え方は減るだろうけど、相当強い全国
的な措置をとらないと、だらだら増えつづけるかもね。
0380１３２人目の素数さん
垢版 |
2021/01/08(金) 07:08:07.47ID:e+x8NepY
あんたは何して食ってんの？
他人の所有物を拝借していることが>369から示唆されるがｗ
0382１３２人目の素数さん
垢版 |
2021/01/08(金) 07:16:03.88ID:sT1sGJ1i
>>380
私？ただの医療関係者ですよ。
他人のものを拝借してるのはあなた。ピンバッジはどこから持ってきたの？w
0385１３２人目の素数さん
垢版 |
2021/01/16(土) 11:23:05.60ID:zWZHa8c/
2万人の被験者の半数にはワクチンを、半数には偽薬を
注射し、1ヶ月後に調べたところ、ワクチン群には100人、
偽薬群には500人の感染者が出た。
このことから、ワクチンを打たずに感染する確率と、
ワクチンを打って感染する確率との比率を信頼区間95%で
求めよ。
0386１３２人目の素数さん
垢版 |
2021/01/17(日) 11:31:20.91ID:/MrggCE/
>>385
事前確率を一様分布として乱数発生させてリスク比を計算。

> va=rbeta(k,1+100,1+9900)
> pl=rbeta(k,1+500,1+9500)
> rr=pl/va
> summary(rr)
Min. 1st Qu. Median Mean 3rd Qu. Max.
2.911 4.626 4.973 5.010 5.352 8.725
> quantile(rr,c(0.025,0.975))
2.5% 97.5%
4.041304 6.183432
>
0387１３２人目の素数さん
垢版 |
2021/01/17(日) 11:36:34.22ID:/MrggCE/
>>385
このワクチンで一人の感染者の発生を抑制する為には何人に接種する必要があるか？95%信頼区間とともに求めよ。
0388１３２人目の素数さん
垢版 |
2021/01/17(日) 11:52:02.13ID:/MrggCE/
正規分布で近似した公式を使うより
乱数発生させた方が色々な計算に対応できて(・∀・)ｲｲ!!
0389１３２人目の素数さん
垢版 |
2021/01/17(日) 11:53:50.05ID:/MrggCE/
>>385
このワクチンが感染リスクを1/10以下にする確率を求めよ、とかも
乱数発生にさせてなら簡単に計算できる。
0390１３２人目の素数さん
垢版 |
2021/01/17(日) 12:05:22.41ID:/MrggCE/
r1=100
n1=10000
r2=500
n2=10000

k=1e7
va=rbeta(k,1+r1,1+n1-r1) # vacccine
pl=rbeta(k,1+r2,1+n2-r2) # placebo
rr=pl/va
summary(rr)
quantile(rr,c(0.025,0.975))

NNT=1/(pl - va)
summary(NNT)
quantile(NNT,c(0.025,0.975))

mean(1/rr < 1/10)

> va=rbeta(k,1+r1,1+n1-r1) # vacccine
> pl=rbeta(k,1+r2,1+n2-r2) # placebo
> rr=pl/va
> summary(rr)
Min. 1st Qu. Median Mean 3rd Qu. Max.
2.932 4.626 4.973 5.009 5.353 9.562
> quantile(rr,c(0.025,0.975))
2.5% 97.5%
4.041894 6.182896
>
> NNT=1/(pl - va)
> summary(NNT)
Min. 1st Qu. Median Mean 3rd Qu. Max.
18.81 24.04 25.02 25.10 26.07 34.86
> quantile(NNT,c(0.025,0.975))
2.5% 97.5%
22.34791 28.29640
>
> mean(1/rr < 1/10)
[1] 0
>
0391１３２人目の素数さん
垢版 |
2021/01/17(日) 12:45:21.55ID:xg9gjMkS
>>386
別スレで質問したら、感染者数を正規分布で近似すると、
Fieller-Hinkley分布っていうもので表されるらしい。

ここに解析的な結果が示されている。
ttps://nshi.jp/contents/doc/nrd_ci/121030_nrd_ci.pdf

めんどくさいので、まだ計算してないｗ
0392１３２人目の素数さん
垢版 |
2021/01/17(日) 12:54:03.37ID:/MrggCE/
ワクチンで感染確率が1/5未満になる確率もこの程度
> mean(1/rr < 1/5)
[1] 0.4803701
0394１３２人目の素数さん
垢版 |
2021/01/17(日) 13:06:17.47ID:/MrggCE/
Rのfmsbというパッケージに
リスク比とその信頼区間を計算すると関数、riskratio　とかいうのがあったはず。
引数なしでriskratioと入力すればソースが出てくる。
うろ覚えだが、正規分布近似で算出していたと思う。
0395１３２人目の素数さん
垢版 |
2021/01/17(日) 13:22:13.32ID:/MrggCE/
https://rdrr.io/snippets/
で実行してみた

library(fmsb)
riskratio(500,100,10000,10000)

Disease Nondisease Total
Exposed 500 9500 10000
Nonexposed 100 9900 10000

Risk ratio estimate and its significance probability

data: 500 100 10000 10000
p-value < 2.2e-16
95 percent confidence interval:
4.041158 6.186345
sample estimates:
[1] 5

乱数発生にさせた結果とほぼ同じ。
0398１３２人目の素数さん
垢版 |
2021/01/17(日) 15:12:41.94ID:/MrggCE/
>>397
それは医師コンプのあんただろう。
医学部落ちたのか？

>385のようなデータはいずれ出てくるだろうから、自分が受けるかどうかの判断に使える。スパイク蛋白をコードするmRNAで作った抗原に対する抗体がADEを起こす可能性はないのか不安だね。
医療従事者優先は先行の人柱だから考えものだな。
0399１３２人目の素数さん
垢版 |
2021/01/17(日) 15:38:31.83ID:qi/NEKt+
>>393
どんな分布の乱数を発生させたらいいかが判っていないと結果が妥当なものにならないんじゃないか？
0400１３２人目の素数さん
垢版 |
2021/01/17(日) 16:29:32.04ID:xg9gjMkS
>ID:/MrggCE/
いろいろ計算ありがとう。

元ネタはこれ。
ttps://www.bmj.com/content/bmj/372/bmj.n124.full.pdf
陰性群の感染者が14173人中404人で、陽性群の再感染者が
6614人中44人。「修正」オッズ比が0.17（95%信頼区間が0.13
から0.24）。

実は、これでオッズ比がなんで0.17になるのかと思って（
普通に計算したら0.23）、分布入れてちゃんと期待値を
求めればそうなるのかしら、と思ったんだけど、なんか
違うみたいだね。なんらかの修正を加えた結果のようだ。

にしてもおかげさまで勉強になりました。
0401１３２人目の素数さん
垢版 |
2021/01/17(日) 18:27:37.26ID:GiUUK9x4
fmsbのoddsratioを使ってもその値にはならないなぁ。
adjustedの解説がほしいところ


Disease Nondisease Total
Exposed 44 6570 6614
Nonexposed 409 13764 14173
Total 453 20334 20787

Odds ratio estimate and its significance probability

data: a b c d
p-value < 2.2e-16
95 percent confidence interval:
0.1649139 0.3080067
sample estimates:
[1] 0.2253765
0403１３２人目の素数さん
垢版 |
2021/01/17(日) 18:59:49.55ID:GiUUK9x4
感染数は二項分布に従うとして乱数発生させてオッズ比の信頼区間を求めると
パッケージのfmsbのoddratioと似たような値になるなぁ。

> r1=409
> n1=14173
> r2=44
> n2=6614
>
> k=1e6
> neg=rbeta(k,1+r1,1+n1-r1) # negative history
> pos=rbeta(k,1+r2,1+n2-r2) # positive history
> neg.o=neg/(1-neg)
> pos.o=pos/(1-pos)
> or=pos.o/neg.o
> quantile(or, c(0.025,0.975))
2.5% 97.5%
0.1652491 0.3078512
0404１３２人目の素数さん
垢版 |
2021/01/17(日) 19:02:22.51ID:GiUUK9x4
>>402
医師免許の画像はネットに溢れているらしいから
医師会のネクタイピンの画像を開業医スレにあげておいたよ。
医師会加入したら贈られてくる。
昔は医師政治連盟の会員に自動的になったけど今は辞退できる。
0406１３２人目の素数さん
垢版 |
2021/01/17(日) 19:15:17.92ID:GiUUK9x4
有症者の方の数字もadjusted odds ratioと乖離した。

> r1=249
> n1=14173
> r2=15
> n2=6614
>
> k=1e6
> neg=rbeta(k,1+r1,1+n1-r1) # negative history
> pos=rbeta(k,1+r2,1+n2-r2) # positive history
> neg.o=neg/(1-neg)
> pos.o=pos/(1-pos)
> or=pos.o/neg.o
> quantile(or, c(0.025,0.975))
2.5% 97.5%
0.07616186 0.21299638
> summary(or)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.03219 0.11050 0.13240 0.13550 0.15710 0.37640
>
0407１３２人目の素数さん
垢版 |
2021/01/17(日) 19:17:16.87ID:GiUUK9x4
>>405
どこから拝借してくるって、盗人の発想の奴には何を示しても無理だろね。ところで医学部落ちたの？
いつもそういうことやっているから思いつくのだろうな。
0409１３２人目の素数さん
垢版 |
2021/01/17(日) 22:08:04.05ID:GiUUK9x4
>>400
ロジスティック回帰の係数で求めたのがadjusted odds ratioらしい。
https://rstudio-pubs-static.s3.amazonaws.com/138745_69f9cf4d1d8b41d781aa5f3ac4d84e35.html#contingency-table

Adjusted Odds Ratio
This is also known as model based odds ratio.
でも、既往の有無という変数１つだから単純にオッズ比を計算と同じになるのではと思った。

やってみた。

'
感染有 感染無
既往有 44 6570
既往無 409 13764
'

covid=
cbind(
replicate(44, c(1,1)),
replicate(409, c(1,0)),
replicate(6570, c(0,1)),
replicate(13764, c(0,0)))

COVID=t(covid)
colnames(COVID)=c('infection','history')
COVID19=as.data.frame(COVID)
glm=glm(infection~.,data=COVID19, family=binomial)
or=exp(glm$coef)
round(cbind(or, exp(confint(glm))),3)

やっぱり同じ。0.17にはならないな。
or 2.5 % 97.5 %
(Intercept) 0.030 0.027 0.033
history 0.225 0.163 0.304
0410１３２人目の素数さん
垢版 |
2021/01/17(日) 22:10:40.52ID:GiUUK9x4
年齢、性別、人種とかを含めてロジスティック回帰したのかも知れんが、元ネタの論文をみても明確でないな。
0411１３２人目の素数さん
垢版 |
2021/01/17(日) 22:13:24.04ID:GiUUK9x4
>>408
拝借なんて思いつくのはお前の盗人根性の反映だよ。

明日のバイト先にはautofocusの機種があるんだが0.3mm太いので使いにくい。focusがあうのに時間が微妙にかかるのでゆっくり動かさないと微小病変を見逃す。
他の医師も使いたがらないと言ってたな。内視鏡の洗浄が間に合わないときは時々使っている。
経鼻は鼻出血のリスクがあるのでコロナの時代には避けたい。
経鼻を経口でやっても苦痛は同等なので希望者には経口でやっている。ただ、拡大機能がないのとクリップが使えないので不便。
先週は月曜日が祝日だったので明日は件数が多そう。
autofocusの出番かな。
0412１３２人目の素数さん
垢版 |
2021/01/17(日) 22:19:51.40ID:GiUUK9x4
>400の元ネタのおかげで
多変量ロジスティック回帰でのadjusted odds ratioの出し方の勉強になったな。
0414１３２人目の素数さん
垢版 |
2021/01/17(日) 23:03:20.91ID:cD1/eM4D
そもそもお前のバイト()の話なんざ誰も聞いてないんだけど。医師免許はよ！
0415１３２人目の素数さん
垢版 |
2021/01/17(日) 23:22:35.89ID:cD1/eM4D
>>411
そもそも俺はお前が医者だなんて微塵も思ってないからな。自分のものではないとなるとどこからか持ってきたものに決まってるだろ。身内か何か知らんが。他人の威を借りるとはとんだ不届き者だな。
0416１３２人目の素数さん
垢版 |
2021/01/18(月) 00:08:46.83ID:O9lD5jqh
医者かどうか関心ないが
臨床に統計はいちいち持ち出して考えることは無い
個々の患者の状態から最善を尽くすだけ
統計じゃない

ぐうの音も出ない正論ですわ。
0417１３２人目の素数さん
垢版 |
2021/01/18(月) 00:53:43.58ID:00TjG9Tw
まあ、匿名掲示板で身分を云々することほど不毛なことはないからね。

決定的におかしなことでも言ってなければ、他人が自称する身分を否定
する必要はないでしょ。

逆に、何を言っても信じない人は信じないしね。
0419１３２人目の素数さん
垢版 |
2021/01/18(月) 04:26:40.65ID:TCD9nCM7
イヤ、ブロおじが馬鹿にされるのは純粋に数学的にクソだから
数学板では数学力があるかないかだけで評価される
後の事はどうでもいい
0421１３２人目の素数さん
垢版 |
2021/01/18(月) 10:21:19.52ID:00TjG9Tw
>>418
どこがおかしいのでしょうか？

>>419
数学的なことならあれこれ批判するのはむしろ建設的ですが、
身分の真偽など、まさにどうでもいいことですよね。
0423１３２人目の素数さん
垢版 |
2021/01/18(月) 17:21:42.49ID:DZTf0z3I
広島市80万とか120万とか無作為に検査したら
感染者ω10万人行くかもなωωω
0424１３２人目の素数さん
垢版 |
2021/01/18(月) 18:59:57.05ID:00TjG9Tw
相撲取りの無作為検査でも1％未満だったのに、10%もいくわけ
ないだろ。あほかいな。
0425１３２人目の素数さん
垢版 |
2021/01/18(月) 22:28:09.63ID:V4G0n17V
>>424
アメフトの件あったけど、接触するスポーツだから庶民よりも気をつけてたと思うよ。
市中はもっと高い
0426１３２人目の素数さん
垢版 |
2021/01/19(火) 00:47:34.30ID:5sxQOhOt
接触するスポーツだから感染機会も多い。

ラグビーみたいなチームスポーツだとなおさらで、
とんでもない感染率だったろ。
0428１３２人目の素数さん
垢版 |
2021/01/19(火) 07:55:47.10ID:AF/97qkc
>>416
検査のカットオフ値によって判断が変わるんだよな。
ROCくらい描けた方がいいね。
ピロリ抗体の高値陰性患者の扱いをどうするかのときに考える。
0429１３２人目の素数さん
垢版 |
2021/01/19(火) 10:26:22.39ID:5sxQOhOt
>>427
まあ、思うのは勝手だからな。

俺もおまえはIQ80くらいの馬鹿だと思ってる。
0435１３２人目の素数さん
垢版 |
2021/01/21(木) 01:21:45.49ID:wnrMDA5R
↑
なんだか頭のおかしな人が...
怖いわ。
0436１３２人目の素数さん
垢版 |
2021/01/21(木) 01:46:58.68ID:o3YD3+kW
　　　 i｀i
　　　 | .|_,.､/'i
　　　ﾉ'" ﾆヽｲ
　　r〈　　　くﾝ　　　　　　__
　　ﾄ-'r､ｨ-へ7　 _Σﾆ'ゝ=<ｰｧ
　　Y　　￣'　|'"´　　　　　 ￣｀ヽ､
　　|　　　　　|___,,.._ゝ___　　　　, 　｀ヽ､
　　|　　　　 ｲ-'/´￣/￣ﾊ｀""''ヽ､ｨ　_i
　　|　　　　 / /|　/-!､ﾚ'　|　ﾊ_ !　 ｀''〈
　　.|　　　　| / ﾚ'ー=､　　ﾚ,.=､､7　ﾊ　|
　　 |　　　 .|　| | ""　 　. 　ﾋ_ﾉｲ/ /__!ｲ　　こいつ最高にｱﾎ のAAです
　　 |､__　　|,|　 |　　ｒｰ--､ 　 "i ﾄ !ソV､
　ﾄ､,ﾄ､,___,.ｨ! ',　ﾄ､　ヽ.__,ノ 　,.ｲ　ﾊ'￣'｀ご自由にお使いください
i＼ |　　　　　 ヽ|｀「＞ｰr=i'"ﾚﾍﾚヽ|
〉　　ヽ､　　　　ｲ´　Σ>o<{ﾊヽ､ｨ'i_
ゝ､_　　｀ヽr-'"´ ､　　／　　 i　 〉,　フ
　∠　　＞ヽ. 　Σ>o<{　　_,ｲ　　ｲイ.
　　レへ,　　ﾊ　　 /　⌒　　|__,.-ﾍ.7
　　　　　 '⌒!,ﾍ／ ､　　　,　ﾄ,　　　＼
　　　　　　　　r〉　'　　 i　　 ヽ'〈´ 　　〉
　　　　　　　／　/　　　!　　ヽ　ヽ､　/
　　　　　　/　　/　　　　　　　,イ　｀く
　　　　 ／ 　 /　　　　 _rﾍo'_　　　　ヽ､
　　　rく　/　　　　　　　レiヽ_７　　ヽ　 」､
0437１３２人目の素数さん
垢版 |
2021/01/21(木) 01:57:40.20ID:wnrMDA5R
なんか、キチガイに触っちゃったみたい...
怖いわー。
0439１３２人目の素数さん
垢版 |
2021/01/28(木) 14:51:59.34ID:8i6B8AWw
コロナ終末時計
出典：英・某私立大
https://www.lboro.ac.uk/news-events/news/2020/march/coronavirus-doomsday-clock-how-close-to-midnight/

感染進行率＝（回復者数＋死者数）／感染者数
終末までの秒数（秒）＝（1−感染進行率）×60

進行率100％、残り秒数0秒で
全国民が死者もしくは感染経験者となる
感染が終息しても値は100％となるが、
ここでは考えない

回復者数のデータがない場合
潜伏期間2週間＋闘病期間2週間＝1か月
または
潜伏期間2週間＋軽傷者の闘病期間1週間＝3週間
とおいて

感染進行率
≒［総感染者数−死者数−（潜伏・闘病者数≒直近1か月の感染者数）］／総感染者数
≒［総感染者数−死者数−重傷者数−（潜伏・軽傷闘病者数≒直近3週間の感染者数）］／総感染者数
で近似できる
0441１３２人目の素数さん
垢版 |
2021/01/28(木) 15:23:51.59ID:7P00nNRm
>>428
今日、抗体価が６で旧基準だと陰性だけど新基準だと高値陰性のケースに内視鏡施行した。
犬猫を飼育していないか問診してから施行。
活動性胃炎の所見はなかったのでハイルマニも否定的だった。
0443１３２人目の素数さん
垢版 |
2021/01/29(金) 07:28:57.35ID:QcH0De8M
>>442
ピロリ抗体高値陰性で内視鏡所見は活動性胃炎のときに
鞭中ピロリ抗原を検査して陽性で除菌した経験ない？
0444１３２人目の素数さん
垢版 |
2021/01/29(金) 07:36:08.92ID:QcH0De8M
>>424
力士の検査はGOTO前だから今とは感染状況が違うと思う。
俺の居住地でも３つの病院がクラスター化した。
0445１３２人目の素数さん
垢版 |
2021/01/29(金) 08:25:51.35ID:QcH0De8M
新型コロナで死亡者数ってポアソン分布に従うだろうか？
馬に蹴られて死亡した数がポアソン分布の起源だったというのだが。
0446１３２人目の素数さん
垢版 |
2021/01/29(金) 10:33:56.47ID:Kr5BKJQ0
>>444
はあ？初場所前の検査で今年に入ってからやったはずだが。
GOTO前だと７月だぞ。何言ってるの。
0447１３２人目の素数さん
垢版 |
2021/01/29(金) 10:37:14.27ID:QcH0De8M
クラスター化した病院だと感染者が患者の場合約20%が死亡するらしい。

永寿総合病院(東京)　214人感染　43人死亡
吉田病院(北海道)　212人感染　38人死亡
（感染者には看護師も含まれているから、入院患者の死亡率は20%に達しているだろう。）
戸田中央総合病院は患者だけで148人感染、31人死亡。

問題

某院で46人の患者感染と報道されていた。
死亡する確率率を20%とすると
予想される死亡者数とその95%信頼区間はいくらになるか？

正規分布近似でなしに計算するには起こりやすい順に確率を加算していって0.95を超える値を探るしかないのかなぁ？
0449１３２人目の素数さん
垢版 |
2021/01/29(金) 10:42:03.71ID:QcH0De8M
>>446
日本相撲協会は7日、親方や力士ら全協会員に新型コロナウイルスのPCR検査を実施すると発表した。
過去に全協会員への抗体検査を実施しているが、PCR検査は初めて。1000人規模となる。

GOTO前に抗体検査しているというぞ。
0451１３２人目の素数さん
垢版 |
2021/01/29(金) 11:41:17.71ID:Kr5BKJQ0
>>449
だからどうした？
直近のPCR検査で1%未満なんだが。
0452１３２人目の素数さん
垢版 |
2021/01/29(金) 11:46:28.93ID:Kr5BKJQ0
>>447
永寿とかもそうだが、入院患者の多くは高齢者で基礎疾患あり
だから、20%の致死率はそんなもんだろうな。
0453１３２人目の素数さん
垢版 |
2021/01/29(金) 11:46:51.92ID:QcH0De8M
>447(自己レス)
数値を変えて計算できるようにRで関数化

fn <- function(
n=46, # 感染患者数
p=0.2, # 死亡確率
cl=0.95, # confidence level
print=FALSE
){
px=dbinom(0:n,n,p)
ix=order(px,decreasing = TRUE) # 確率の高い順の死亡数のindex
dx=(0:n)[ix] # 確率の高い順の死亡数
cx=cumsum(px[ix]) # 高い順の確率総和（累積和）
i=which(cx>cl)[1]
cx[c(i-1,i)]
cat('range = ',range(dx[1:i]),' conf.level =',round(cx[i],3),'\n') # dx[1:i]
cat('range = ',range(dx[1:(i-1)]),' conf.level =',round(cx[i-1],3),'\n') # dx[1:(i-1)]

nlwr=qnorm(0.05/2,n*p,sqrt(n*p*(1-p))) # 正規分布近似
nupr=qnorm(0.05/2,n*p,sqrt(n*p*(1-p)),lower=FALSE)
cat('asymptotic range =',round(nlwr,2),' ',round(nupr,2),'\n\n')
if(print){
library(HDInterval)
plot(0:n,px,type='h',bty='l',xlab='death toll',ylab='probability',col=2,lwd=3)
print(HDInterval::hdi(qbinom,size=n,prob=p))
print(HDInterval::hdi(qnorm,mean=n*p,sd=sqrt(n*p*(1-p))))}
}

> fn(46,0.2)
range = 4 14 conf.level = 0.959
range = 5 14 conf.level = 0.936
asymptotic range = 3.88 14.52

パッケージHDIntervalの値と一致するから多分、バグはないと思う。
0456１３２人目の素数さん
垢版 |
2021/01/29(金) 11:52:40.46ID:QcH0De8M
抗原検査キットはクリニックレベルまで不急しているから、ぽつぽつと陽性がでる。
インフルエンザは今シーズンは０。同僚や看護師に聞いてもインフルエンザは経験ないっていってた。
インフルエンザワクチンは必要なかったんじゃないかな。
0459１３２人目の素数さん
垢版 |
2021/01/29(金) 14:03:13.86ID:Kr5BKJQ0
>>457
>>435,437を書いた俺は識者でも指揮者でもないけど、
>>450とは別人だし、彼と意見が同じわけでもない。
0460１３２人目の素数さん
垢版 |
2021/01/29(金) 14:04:21.97ID:Kr5BKJQ0
↑
いけね、>>450じゃなくて、>>454とは別人
0462１３２人目の素数さん
垢版 |
2021/01/29(金) 15:30:39.59ID:QcH0De8M
>>458
医者が羨ましいとはね。接客業が賤業と気づいたときには時遅しだったな。
近隣の病院が３箇所クラスター化しているし、内視鏡も防護服を着ての施行。
明日は我が身。
0464１３２人目の素数さん
垢版 |
2021/01/29(金) 21:49:31.76ID:XpAxQwWj
>>463
通勤のタクシーチケットも支給される優良職場はなかなかないんだな。
去年は内視鏡休診でも全額給与が支給されたし。
昔、研修にいっていた急性期病院でもクラスター発生の報道に接した。
アクティブな病院ほどリスクが高い。
0466１３２人目の素数さん
垢版 |
2021/01/29(金) 22:19:43.81ID:jMaUU9SA
>>461
おまえが自演でないという証拠は？
0468１３２人目の素数さん
垢版 |
2021/01/30(土) 07:24:11.24ID:zGHWLTrO
>>464
積極的にコロナ診療していた民間病院のホームページから

(quote)

地域の皆さん、患者さん、ご家族の皆さん

非常事態を宣言します。
新型コロナウイルスの感染が急拡大し、極めて身近に
なってきました。院内での感染も多く認められるように
なりました。感染経路が追えなくなってきています。
救急外来も止めざるを得ません。
0469１３２人目の素数さん
垢版 |
2021/01/30(土) 08:11:29.46ID:zGHWLTrO
>>445
都道府県別の死亡者数がここにあった
https://www3.nhk.or.jp/n-data/opendata/coronavirus/nhk_news_covid19_prefectures_daily_data.csv

東京のデータだけ、抜き出すと
> TYO
[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[26] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0
[51] 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 0 1 1 1 1
[76] 7 0 0 2 5 7 0 1 4 1 4 0 2 0 5 6 3 7 5 3 6 4 0 6 6
[101] 7 0 6 2 9 3 6 15 4 5 0 5 5 11 9 0 9 7 7 9 7 11 7 4 3
[126] 3 9 7 8 9 8 4 4 3 3 2 1 0 1 0 1 2 2 0 0 0 0 0 2 1
[151] 0 0 2 1 0 3 0 0 1 2 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0
[176] 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 1 1 2 0 0
[201] 1 0 0 0 0 0 0 1 2 0 2 0 0 2 1 3 3 1 2 0 0 2 2 2 0
[226] 2 0 2 3 1 0 1 4 0 0 3 6 1 0 1 2 0 2 2 3 0 1 0 0 1
[251] 0 1 3 4 1 0 6 1 1 1 0 2 0 2 1 2 1 4 0 0 0 2 2 2 4
[276] 3 2 0 1 2 3 3 0 2 3 1 0 2 2 0 0 4 0 0 2 0 0 0 1 3
[301] 3 2 0 0 0 2 4 0 0 1 1 1 0 0 3 3 3 0 0 1 5 5 2 8 2
[326] 0 6 6 5 2 5 0 0 3 9 10 0 4 5 0 1 1 10 9 10 10 0 6 5 4
[351] 5 4 0 1 2 14 8 11 7 8 3 4 2 13 3 10 3 5 3 16 10 7 9 9 3
[376] 14 13 18 20

ヒストグラムにすると
https://i.imgur.com/3aAVZVj.png
0470１３２人目の素数さん
垢版 |
2021/01/30(土) 08:22:42.12ID:zGHWLTrO
>>469
Rのfitdistrplusを使って平均値＝実測となるポアソン分布を重ねると
https://i.imgur.com/oqmVhJ0.png

goodness of fitnessを算出させると

Chi-squared statistic: 15037.22
Degree of freedom of the Chi-squared distribution: 6
Chi-squared p-value: 0

Goodness-of-fit criteria
1-mle-pois
Akaike's Information Criterion 2184.220
Bayesian Information Criterion 2188.157

ポワソン分布には従ってないみたい。
0471１３２人目の素数さん
垢版 |
2021/01/30(土) 08:32:34.76ID:zGHWLTrO
>>469
理由はさっぱりわからんが、負の二項分布を当てはめると

https://i.imgur.com/gsWUdhr.png

なんとなく、分布が近似している。

Rのコード
library(fitdistrplus)
fTYO=fitdist(TYO,'nbinom')
summary(fTYO)
gofstat(fTYO)
plot(fTYO)

> gofstat(fTYO)
Chi-squared statistic: 8.525185
Degree of freedom of the Chi-squared distribution: 5
Chi-squared p-value: 0.1295687
ポワソン分布の時と違って、
負の二項分布に従うという帰無仮説は棄却されない。
0473１３２人目の素数さん
垢版 |
2021/01/30(土) 08:52:34.02ID:zGHWLTrO
>>471
大阪府のデータを使って負の二項分布で近似させてみると

https://i.imgur.com/Au0x2W8.png

> gofstat(fOSK)
Chi-squared statistic: 5.835746
Degree of freedom of the Chi-squared distribution: 4
Chi-squared p-value: 0.2117552

Goodness-of-fit criteria
1-mle-nbinom
Akaike's Information Criterion 1463.302
Bayesian Information Criterion 1471.177

東京よりも当てはまりが良好

ちなみに東京のAICとBICは
Goodness-of-fit criteria
1-mle-nbinom
Akaike's Information Criterion 1470.136
Bayesian Information Criterion 1478.011
0474１３２人目の素数さん
垢版 |
2021/01/30(土) 09:08:17.72ID:zGHWLTrO
>>472
医師免許証の画像ならネットにいくらでも転がっているらしいぞ。
医師会のタイピンの画像はないだろうから、開業医スレにアップしておいたよ。

良識派の感想

数学 統計に詳しい人が語るコロナウイルス ☆2
https://rio2016.5ch.net/test/read.cgi/math/1596506253/435

435 名前：１３２人目の素数さん[] 投稿日：2021/01/21(木) 01:21:45.49 ID:wnrMDA5R
↑
なんだか頭のおかしな人が...
怖いわ。
0475１３２人目の素数さん
垢版 |
2021/01/30(土) 09:26:28.98ID:zGHWLTrO
>>471
指数分布で近似しようとしたのだが、やってみたら負の二項分布の方が当てはまりがよかった。

近似した負の二項分布を図示すると
https://i.imgur.com/pchRFBr.png

１０人以上の死者がでる確率は
fTYO$estimate
Tsize=fTYO$estimate[1]
Tmu=fTYO$estimate[2]
pnbinom(10-1,size=Tsize,mu=Tmu,lower=F)
で計算できる。

こういう
https://i.imgur.com/EVaHN6x.png
トレンドを考慮しない計算であるのは当然ではあるけど。
0478１３２人目の素数さん
垢版 |
2021/01/30(土) 09:39:38.96ID:BIYwEXHb
西浦さんの最近出た「gotoがコロナ拡散に影響を与えた」とする論文読んでみたけどさ
https://www.mdpi.com/2077-0383/10/3/398


これ殆ど何も示せてなくね
因果を示せとまでとは言わないけど、意味のある相関すらないんじゃないか
旅行関連感染者が増えたって言うけど非旅行関連感染者の方が増えてるし
観光関連感染者が増えたってのもあるけど、そもそもGotoで母集団（全国民）の観光旅行客が爆増してんだからそこ考慮に入れなきゃ意味ないと思うけどそういうのもないし
そもそも時期のチョイスもちょっとおかしいとか色々あるけどそういうの置いておいてもそもそも根本的に疑似相関にすらなってないと思う


データ集めただけですよって言うならわかるけど
これって数学的にgotoと感染との関係を示唆出来てるとは思えんのよね
でもDiscussionにハッキリ、gotoが影響与えた可能性ある、って書いちゃってるし
0479１３２人目の素数さん
垢版 |
2021/01/30(土) 09:41:19.42ID:zGHWLTrO
指数分布近似でのグラフ
https://i.imgur.com/0SmMRr9.png
グラフの見た目は負の二項分布の方がよさそうだったのだが、
AICやBICは指数分布の方が小さいのでこっちの方がいいみたい。

　
Goodness-of-fit criteria
1-mle-exp
Akaike's Information Criterion 1369.557
Bayesian Information Criterion 1373.494

eTYO=fitdist(TYO,'exp')
summary(eTYO)
gofstat(eTYO)
plot(eTYO)
Trate=eTYO$estimate ; 1/Trate
curve(dexp(x,rate=Trate),0,30,bty='l')
１０人以上の死亡者がでる確率は
pexp(10,rate=Trate,lower=F)
で計算できる。

指数分布近似の方が確率が高いな。
0481１３２人目の素数さん
垢版 |
2021/01/30(土) 09:47:41.40ID:zGHWLTrO
>>478
人から人へ感染するのだから、人との接触機会が増えれば感染が増えるのは自明だと思う。
スカイダイビングでパラシュートを付けた方が生存率に統計的有意差があるという比較試験はないらしい。
0484１３２人目の素数さん
垢版 |
2021/01/30(土) 10:13:02.64ID:zGHWLTrO
>>482
二期校時代に現役で合格。
一期校は滑り止めに理一。
どちらも模試での判定がAだったから。
もし判定って過去のデータからロジスティク回帰で計算しているんだろうな。
0485１３２人目の素数さん
垢版 |
2021/01/30(土) 10:28:25.67ID:E7OVFGO9
>>484
ここだったら学歴なんか何言おうと思いのまま。嘘も言い続ければ真になるってか？
それが事実ならここまで歪みきってここで燻ってるはずないだろう。それが分からないのが非医。
0486１３２人目の素数さん
垢版 |
2021/01/30(土) 10:31:39.40ID:E7OVFGO9
可哀想だね。
それが人生のピークだったのかな？
それとも全部嘘なのかな？
0487１３２人目の素数さん
垢版 |
2021/01/30(土) 15:52:29.53ID:MIK9j1nE
>>485
ひょっとして俺より低学歴なのかｗｗｗｗ
臨床やっているから、業界ネタはいくらでも語れるぞ。
サワシリンアレルギーの患者にはPPI＋クラリス＋フラジールだと保険査定されにくい。
バイト先で2件経験したので主治医にsuggestion。
2例中１例は除菌成功。あとはグレースビットを含むレジミンだけどこっちは自費になる。
0490１３２人目の素数さん
垢版 |
2021/01/30(土) 20:59:22.36ID:xatoVPck
世論　

統計爺ってスレタイも読めない、医師免許もない、誰からも必要とされない
すこしは誇れるものはないのかね
0491１３２人目の素数さん
垢版 |
2021/01/31(日) 08:03:15.72ID:SAWCTDIo
３９００人はどうやって算出したんだろう。

1 孤高の旅人 ★ sage 2021/01/29(金) 23:30:18.25 ID:WTE/oZYZ9
広島の大規模PCR、最大3900人の陽性想定　県が計画案、検査費用10億円盛る
1/29(金) 23:18配信
https://news.yahoo.co.jp/articles/7f222da16d82a489183003c2b751077b10fb2013

　新型コロナウイルス対策で、広島県は29日、広島市内4区で取り組む任意、無料の大規模PCR検査の基本計画案を発表した。市の感染者数が減少傾向にある中、無症状者を早く見つけ出し、感染を抑え込むのが狙い。1カ月で約28万人を検査し、最大3900人の感染が判明すると想定する。緊急編成した2020年度一般会計補正予算案に事業費10億3800万円を盛り込み、2月3日開会の県議会臨時会に提案する。
0494１３２人目の素数さん
垢版 |
2021/01/31(日) 12:19:33.71ID:FdtBLtuK
>>491
不思議な予測値だね。
広島市の一週間あたりの新規感染者数の最大値が１０万人あたり４０人。
いまは減少傾向で１０人。感染後２週間はPCR陽性が続くとすれば、
２週間以内の新規感染者がひっかかるはずだからピーク値でもせいぜい
８０人。２８万人あたりになおせば２００人程度。３９００ということ
は、潜在的感染者が最大その２０倍はいると見積もってることになるね。
いくらなんでも、それはないと思うが。
0495１３２人目の素数さん
垢版 |
2021/01/31(日) 13:14:05.28ID:EUU1R5Ub
２８万人調査したときの陽性者数の95%信頼区間の上限を3900人とすると
信頼区間は正規分布近似として
> uniroot(function(p,n=280000,upper=3900) qnorm(0.975,n*p,sqrt(n*p*(1-p)))-upper,c(0,1))$root
[1] 0.0135076
陽性率1.35%に想定していることになるなぁ。
0496１３２人目の素数さん
垢版 |
2021/01/31(日) 13:31:51.13ID:EUU1R5Ub
>447のコードを流用して二項分布で計算すると
> uniroot(function(p,n=280000,upper=3900) hdi(qbinom,prob=p,size=n)[2]-upper,c(0,1))$root
[1] 0.01350607
正規分布での値とほぼ一致だな。

でも、最初から陽性率が設定できているなら、何を出すために疫学調査するのだろうかが疑問だな。
0497１３２人目の素数さん
垢版 |
2021/01/31(日) 13:37:36.29ID:EUU1R5Ub
99%信頼区間の上限が3900人としても、想定された陽性率は1.3%程度であまり変わらんな。

> uniroot(function(p,n=280000,upper=3900) qnorm(0.995,n*p,sqrt(n*p*(1-p)))-upper,c(0,1))$root
[1] 0.01338053
> uniroot(function(p,n=280000,upper=3900) HDInterval::hdi(qbinom,prob=p,size=n,cred=0.99)[2]-upper,c(0,1))$root
[1] 0.01338956
0499１３２人目の素数さん
垢版 |
2021/02/02(火) 00:25:28.86ID:u9JBgcH0
東京で14000人に「抗体」検査して、その陽性率が1.8%だからね。

広島でPCR検査して陽性率が1.4%って、まずありえないっしょ。
桁が違ってる。
0500１３２人目の素数さん
垢版 |
2021/02/02(火) 07:32:41.54ID:+Dn/x9Sw
>>498
タイピンは使わないから俺の机の引き出しに以前からずっとあるんだよ。

賤業の臨床医が羨ましいなんて医学部落ちたのかよ？
病棟をもっていると夜間休日も電話かかってきて拘束感が半端ない。感染対策委員会とか救急委員会とか会議も多くてうんざりする。

４月から内視鏡バイトをもう１日増やしてくれと頼まれた。
タクシーチケットも支給されるから引き受けた。
ナースが慣れていて仕事がしやすいというのもある。
0505１３２人目の素数さん
垢版 |
2021/02/02(火) 08:05:08.42ID:+Dn/x9Sw
>>502
今日は日曜救急当番の代休で１日休み。
昨日から抗ウイルス薬の勉強中。
抗ウイルス薬は配合薬が増えた。
TANGO試験のデータの統計処理を自分で検算したりして遊びながら論文読んでる。
算出方法が明記されていなくて95%CIが合致しないことが多いんだな。
ゾフルーザの95%CIはbootstrap法で算出と明記されていたから、これは合致しなくても納得。
0508１３２人目の素数さん
垢版 |
2021/02/03(水) 07:08:33.32ID:MNJ8JS7S
メラノーマの薬だから覚えなくていいなと思っていたら、大腸がんにまで適用拡大されていた。
おまけに同じ作用機序の新薬がでたり。
ビラフﾄﾋﾞやメクトビは阻害する遺伝子のBRAFやMEKからと分かったがトビってなんだろう？と疑問だった。
新薬のインタビューフォームにgoodを意味するヘブライ語だと記載されていた。
ベルソムラとかだとフランス語由来と簡単に想像がついたが、ヘブライ語なんぞ知らん。

インタビューフォームに薬の由来が書いてあるのを知ったのは有意義であった。
こんな感じでシグナル伝達の動画がyoutubeに沢山あるので文章を読まなくていいので楽ちん。
https://www.youtube.com/watch?v=W3ZASYwqjNc
0510１３２人目の素数さん
垢版 |
2021/02/03(水) 17:56:27.32ID:0ffS8nkj
もうちょっと下の世代もそうだった
0511１３２人目の素数さん
垢版 |
2021/02/03(水) 18:02:09.72ID:NIIMmLyz
ウリュウのジジイはそんな奴らを傍目で見て嫉妬で狂っていたと見える。
0512１３２人目の素数さん
垢版 |
2021/02/04(木) 19:58:09.56ID:/FoB3enM
日本国内で日本人初の武漢ウイルス感染が確認されたのは昨年の1月29日。
武漢からのツアー客を乗せたバスの運転手だった。
それから、ほぼ一年。きりの良い時期なので、整理してみようと思った
が、諸処事情から、2020年1月1日を起点にした一年のデータを整理。

2020年1月1日　総人口　：　12602万人 　 　 　 　 　(1)

2021年1月1日　総人口　：　12557万人
2021年1月1日　0〜4才人口　：　463万人
2021年1月1日　1才以上の人口推定値　：　12464.4万人　　 　 　 　 　(2)
2020年1月1日に生きていた人の内、2021年1月1日までに死亡した人の推定人数　：　137.6万人　 　 　 　 　(3)

2021年1月1日　国内の感染者数　：　237071人　 　 　 　 　(4)
2021年1月1日　死亡者数　：　　3540人　 　 　 　 　(5)
2021年1月1日　感染経験生存者　：　　233531人　 　 　 　 　(6)

参照元
http://www.stat.go.jp/data/jinsui/pdf/202001.pdf
http://www.stat.go.jp/data/jinsui/pdf/202101.pdf
https://www.asahi.com/special/corona/
0513１３２人目の素数さん
垢版 |
2021/02/04(木) 19:58:49.59ID:/FoB3enM
前投稿のデータを用い、2020年1月1日に生きていた12602万人を対象に
1年後の　生／死、及び感染経験の有／無　をクロス表にしてみた。

　 　 　 ＼ 　 　 生 　 　 　 　 死 　 　 　 　 　 　 　合計
感染有 　 　 　 23.35　 　 　 　 　 　 　0.354 　 　 　 23.71
感染無 　 　 　12441 　 　 　 　 　 　 　137.25 　 　 　12578.3
合 計 　 　 　12464 　 　 　 　 　 　 　137.6 　 　 　 12602

1段目のデータ(感染有)は左から、(6),(5),(4)
2段目のデータ(感染無)は左から、(2)-(6),(3)-(5),(1)-(4)を利用/計算してます。
0514１３２人目の素数さん
垢版 |
2021/02/04(木) 20:29:11.51ID:6E2Oyknc
おいジジイ
ど平日の社会人が1日中5chに張り付きとはどういうことだ？
0516１３２人目の素数さん
垢版 |
2021/02/05(金) 01:33:59.80ID:opztDdUo
>>513
だから何？感染してなくても1%は死ぬんだから、死亡率1.5%なら
たいしたことないとか？そんな比較しても意味ないでしょうに。

それよか、少子高齢化のため、長年に渡って毎年の死者数は１〜２
万人ずつ増加してるのに、2020年は減少して一昨年のレベルに戻って
んのよ。つまり、コロナで死ぬ人より、感染対策やらなにやらのおか
げで死なずにすんでる人のほうがはるかに多かったってこと。
0517１３２人目の素数さん
垢版 |
2021/02/05(金) 08:36:21.84ID:4R5y0olT
>>514
平日に毎日働かなくても生活できるからね。
新コロナ前からそういう生活している。
稼いでも税金で持っていかれるだけ。
0520１３２人目の素数さん
垢版 |
2021/02/05(金) 15:14:25.20ID:opztDdUo
>>519
12中下旬で東京0.9%
一方、感染者の累計は５万人。ってことで、やっぱりその２倍くらい
だな。無症状、軽症で終わる人の割合がそれだけ多いってことかもね。

そういえば、NHKの番組で、医療現場でマスクつけた従事者でも感染
したけど、９割近くが無症状だったという事例が紹介されてたな。
マスク着用の効果もあるかもしれん。
0521１３２人目の素数さん
垢版 |
2021/02/05(金) 15:16:05.08ID:opztDdUo
ついでにいうと、現時点で累計感染者は10万人だから、
東京の抗体陽性率は現在1.8%くらいと推定される。
0522１３２人目の素数さん
垢版 |
2021/02/05(金) 15:50:24.18ID:Z6vXAk9R
＞生保か？

こういうこというやつって
なんか生活保護受けたいのかな？
だったら発狂して生活保護受ければいいのに
もう狂ってるんだろ？
0523１３２人目の素数さん
垢版 |
2021/02/05(金) 18:31:12.13ID:fNPw/w8j
医療事務の分際で医者ぶってるって滑稽だね。
恵まれた境遇のやつが毎日ルサンチマンまみれの歪んだ書き込みするかよ。
0524１３２人目の素数さん
垢版 |
2021/02/06(土) 05:59:55.42ID:1P0wuzxR
臨床医が羨ましいのか？
リスクの高い賤業だぞ。病棟持ったらプライベートな時間がない。
まあ、バイト先には困らんが。

こういう記述は術者にしか書けんよ。
https://egg.5ch.net/test/read.cgi/hosp/1579701192/548


>>519

※ この調査では、陽性の判定をより正確に行うため、2種類の検査試薬（アボット社・ロシュ社）の両方で陽性が確認されたものを「陽性」としています。

陽性数を少なくみせたい意図を感じるな。
２種類の検査のカッパ値はどの程度なんだろ？
0525１３２人目の素数さん
垢版 |
2021/02/06(土) 07:32:11.69ID:1wyI9Vw2
>>524
何が術者だよ、笑わせんな。
お前は年中無休で5chしかしてない社会と家族のお荷物だろ。
0526１３２人目の素数さん
垢版 |
2021/02/06(土) 08:23:45.30ID:1P0wuzxR
>>525
４月から内視鏡バイトを増やしてくれと頼まれた。
今と同じくタクシーチケット出るので引き受けた。
まあ半分（ざっと所得税４割住民税１割）は税金で持っていかれる。

んで、あんたこそお荷物じゃないの？
仕事あんの？
0528１３２人目の素数さん
垢版 |
2021/02/06(土) 11:57:37.75ID:ojmboENx
>>526
あんたが医者じゃないとか無意味な難癖つけてくる異常者は
放置でいいから、コロナ関連の話をしてくれ。

不毛なやりとりはうんざりだよ。
0529１３２人目の素数さん
垢版 |
2021/02/06(土) 16:11:08.82ID:2vP6d7DR
>>513
武漢ウイルスが直接的原因で無くなった方は、報道されている数字の、27%程。
逆に言えば、新型コロナが蔓延しなくたって、コロナ死と認定された者の73%は、
別の原因での死は免れなかった。
0531１３２人目の素数さん
垢版 |
2021/02/09(火) 05:31:04.47ID:S3mmq/Em
>>528
コロナ患者の退院基準が甘すぎて
「発症日から10日間経過し、かつ、症状軽快から72時間経過した場合に退院可能」だからPCRすらなしで一般病院に押し付けて来るんだなぁ。
俺の勤務先でも受けている。
俺は主治医じゃないけど未だに微熱が続いている。
当直中に発作性心房細動が起こったので呼ばれた。
防護服来て病室に赴いて診察後は風呂に入った。
当直帯はの待機時間が多い勤務だから風呂に入ったりできるけど
コロナ対応の病院はトイレにも満足にいけないという。
0533１３２人目の素数さん
垢版 |
2021/02/09(火) 05:50:13.82ID:S3mmq/Em
HBs抗体があってもde novo肝炎みたいに再活性化することもある。
昔、勤務していた病院でリツキサンを使っていた症例に出ていたな。テノゼットもベムリディもなかった時代の話。
0534１３２人目の素数さん
垢版 |
2021/02/09(火) 06:06:57.61ID:dfx3vZX8
>>513
ありきたりの計算

> dat
[,1] [,2]
[1,] 233531 3540
[2,] 124410469 1372460
> chisq.test(dat)

Pearson's Chi-squared test with Yates' continuity correction

data: dat
X-squared = 353.87, df = 1, p-value < 2.2e-16

> fisher.test(dat)

Fisher's Exact Test for Count Data

data: dat
p-value < 2.2e-16
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.7039544 0.7525420
sample estimates:
odds ratio
0.7277541
0535１３２人目の素数さん
垢版 |
2021/02/09(火) 06:29:23.55ID:dfx3vZX8
>>513
生と死の配置をいれかえたcontingency tableにして感染した場合の死のオッズ比とリスク比を95%CI付きで出してみた。

a=3540
b=233531
c=137.6e4-3540
d=12464.4e4-233531
(dat=matrix(c(a,b,c,d),nrow=2,byrow=TRUE))
chisq.test(dat)
fisher.test(dat)
oddsratio(dat)
rateratio(a,c,a+b,c+d)

> oddsratio(dat)
Disease Nondisease Total
Exposed 3540 233531 237071
Nonexposed 1372460 124410469 125782929
Total 1376000 124644000 126020000

Odds ratio estimate and its significance probability

data: dat
p-value < 2.2e-16
95 percent confidence interval:
1.329178 1.420525
sample estimates:
[1] 1.374092

> rateratio(a,c,a+b,c+d)
Cases Person-time
Exposed 3540 237071
Unexposed 1372460 125782929
Total 1376000 126020000

Incidence rate ratio estimate and its significance probability

data: a c a + b c + d
p-value < 2.2e-16
95 percent confidence interval:
1.324104 1.414398
sample estimates:
[1] 1.368506
0536１３２人目の素数さん
垢版 |
2021/02/09(火) 07:19:36.35ID:dfx3vZX8
>437と>528のように正しい認識の方が複数あらわれたので、統計処理に専念しよう。

>535の計算は二項分布を正規分布近似して信頼区間を算出していると思う。
パッケージfmsbのコードを読むとそうなっている。

Rの
contingencyTableBF {BayesFactor} R Documentation
Function for Bayesian analysis of one- and two-sample designs
を使って乱数発生させてやってみた。MCMCでやったというとカッコいいｗ

死亡数は死亡率に従う二項分布として計算（上段）
死亡数はポアソン分布に従うとして計算（下段）

https://i.imgur.com/BE1nggk.png

どれもほぼ同じ数値。カイ二乗やフィッシャーテストの値ともほぼ同じ。
これをもって、COVID-19に感染すると死亡するリスクが３割から４割増すのですよ、と恐怖を煽るのに使える。


問題：　
(1) COVID-19感染の死亡寄与率とその９５％信頼区間を求めよ。
(2) COVID-19による死者が１人現れるには何人が感染する必要があるか？９５％信頼区間とともに求めよ。
0538１３２人目の素数さん
垢版 |
2021/02/09(火) 08:45:36.13ID:S3mmq/Em
>>536
二項分布による乱数発生で求めた分布
https://i.imgur.com/LNenVZB.png

ポワソン分布による乱数発生で求めた分布
https://i.imgur.com/ppWicYq.png

コロナ感染者が250人でるとそのうち1人は1年以内に死亡するという結果になった。

正規分布近似でのリスク差の値
> riskdifference(a,c,a+b,c+d)
equals to zero)

data: a c a + b c + d
p-value < 2.2e-16
95 percent confidence interval:
0.003532353 0.004509443
sample estimates:
[1] 0.004020898
ともほぼ一致。
0539１３２人目の素数さん
垢版 |
2021/02/09(火) 08:51:11.18ID:S3mmq/Em
>>538
>437と>528のように正しい認識の方が複数あらわれたね。

RのパッケージEpiを使うと、リスク比、オッズ比、リスク差まで95％信頼区間付きで出してくれる。
しばらく使ってなかったのでその存在を失念していた。

> Epi::twoby2(dat)
2 by 2 table analysis:
------------------------------------------------------
Outcome : Col 1
Comparing : Row 1 vs. Row 2

Col 1 Col 2 P(Col 1) 95% conf. interval
Row 1 3540 233531 0.0149 0.0145 0.0154
Row 2 1372460 124410469 0.0109 0.0109 0.0109

95% conf. interval
Relative Risk: 1.3685 1.3244 1.4140
Sample Odds Ratio: 1.3741 1.3292 1.4205
Probability difference: 0.0040 0.0035 0.0045

Asymptotic P-value: 0.0000

twoｂｙ2のコードをみるとこれは正規分布近似での信頼区間算出だな。
0540１３２人目の素数さん
垢版 |
2021/02/09(火) 11:17:18.84ID:6emMqh27
>>536
コロナの死亡率は年齢依存性が著しく高いので、
全年代を対象に分析してもあまり意味がないと
思う。
たとえば、70歳だと死亡リスクは数倍になるけど、
20歳だと0.1%しか違わないとか（数字はテキトー）。
0541１３２人目の素数さん
垢版 |
2021/02/09(火) 13:32:44.31ID:S3mmq/Em
>>540
まあ、層別に処理しないとシンプソンのパラドックスが出てくるだろうね。
新型コロナの死亡率が１〜２％くらいという話だから、
コロナ感染者が250人でるとそのうち1人は1年以内に死亡するという結果はそれほど的外れでもないと思う。
0542１３２人目の素数さん
垢版 |
2021/02/09(火) 15:54:39.61ID:6emMqh27
そもそも６０代後半になれば、コロナに感染してなくても
１年以内に死亡する確率が１％くらいあるはず。なので、
感染したら１年以内に0.4％が死ぬと言われてもどうって
ことない。
0543１３２人目の素数さん
垢版 |
2021/02/09(火) 18:57:38.71ID:S3mmq/Em
こういう記述に接すると早速グラフ化したくなる。

"
The proportion (f) of the population that must be vaccinated to end an outbreak equals (1-[1/R0])/E,
where E is the vaccine efficacy to prevent transmission.
If E = 95% (0.95) and R0 = 3, then (1-1/3)/E = 70% of the population that will need to be vaccinated.
If the vaccine actually prevents a smaller proportion of transmissible infections,
or if the virus becomes more transmissible, as the London variant is said to be,
the proportion of the population that needs to be vaccinated is greater.
If the vaccine efficacy to prevent transmission is less than (1? 1/R0),
it is impossible to eliminate an infection solely by vaccinating even the whole population.

https://www.msdmanuals.com/en-jp/professional/news/editorial/2021/01/08/14/02/covid-vaccine-update
"

https://i.imgur.com/c6vsfPA.png
0544１３２人目の素数さん
垢版 |
2021/02/09(火) 19:33:24.14ID:6emMqh27
風邪と同じで、ワクチンの効き目が弱い変異株が現れるので、
永久に感染は終息しない、に百万ペプル。
0545１３２人目の素数さん
垢版 |
2021/02/09(火) 20:11:10.00ID:Xio7PfL3
なぜ社会貢献したキャプテントムが亡くなって社会のゴミのウリュウのジジイがのうのうと生きてるのか
0546１３２人目の素数さん
垢版 |
2021/02/09(火) 20:50:34.06ID:S3mmq/Em
>>543
その頁の

The efficacy 14 to 28 days after a single dose was also roughly 90%
(2 cases in the 983 vaccine recipients and 28 cases in the 1059 placebo recipients.

というデータに日本の再生算数をモード値１．２５、分散０．２５のガンマ分布に勝手に仮定して何割にワクチン接種が必要かをグラフ化。

https://i.imgur.com/jmHK8Wh.png

３割り程度で制御できるという計算になった。
0547１３２人目の素数さん
垢版 |
2021/02/09(火) 20:54:05.66ID:S3mmq/Em
>>544
B型肝炎みたいに中和抗体とされるHBs抗体ができても免疫抑制状態で再活性化されるような振る舞いになれば駆逐されることはないだろね。
人の移動だけじゃなくて、地球温暖化で熱帯地方の風土病（例えばマラリアを媒介する蚊が繁殖するとかもあるんじゃないかなぁと思っている。
0549１３２人目の素数さん
垢版 |
2021/02/10(水) 01:44:48.18ID:RwQzsw9w
>>547
ようわからんけど、鼻風邪のワクチンをつくれば鼻風邪は
駆逐できるのか、ってことでしょ。インフルエンザがワクチン
で駆逐できないのと同様、無理っしょ。

ある程度重症化が防げれば良しとしなくてはならんのでないの？
0551１３２人目の素数さん
垢版 |
2021/02/10(水) 11:34:54.15ID:cSZfFCkj
>>549
水疱瘡のウイルスが神経節に残って免疫が低下すると帯状疱疹として発症するようなことが起こりうるのではないかという話。
de novo肝炎の予後は不良。
0552１３２人目の素数さん
垢版 |
2021/02/10(水) 13:39:06.34ID:cSZfFCkj
>>549
ポリオと痘瘡は、ほぼ駆逐じゃないかな。
日本が狂犬病の清浄国でいられるのもいつまでだろうな。
野生のアライグマに噛まれた患者に本人の希望で暴露後免疫でワクチン注射したことがある。５〜６回で合計3万5千円くらいだったかな。
0554１３２人目の素数さん
垢版 |
2021/02/10(水) 17:33:58.31ID:RwQzsw9w
>>552
ポリオは小児しか発生しなくて、発症率が極端に低いことが
大きいんじゃないかな。発症者数をワクチンで一時的に抑え
こめれば大丈夫とか。しらんけど。
天然痘はDNAウィルスだから変異が遅いんじゃね？
0555１３２人目の素数さん
垢版 |
2021/02/11(木) 06:48:55.53ID:ewUi1drK
>>554
ポリオはＲＮＡウイルスだが逆転写なしでゲノムがそのままmRNA として働くらしいので工程が少なくてすむから変異が少ないのか？
0556１３２人目の素数さん
垢版 |
2021/02/11(木) 07:56:39.19ID:ewUi1drK
>>549
経鼻インフルエンザワクチン（フルミスト）の方が理論的には注射のワクチンより予防効果が高そうなんだが、
臨床試験すると注射の方が効果が高いという。
理論と臨床が乖離する実例。

CD-28を刺激で免疫が賦活されるはず、という理論での試験がTG1412事件。
刺激じゃなくて抑制因子の抑制のヤーボイは免疫暴走するには至らなかった。
0557１３２人目の素数さん
垢版 |
2021/03/06(土) 06:18:09.70ID:Mo0cuWGb
第3波対策の緊急事態宣言は、

前半
既存株の実効再生産数を0.75くらいに抑えながらも、
変異株の実効再生産数を1.15くらいになりじわじわ増える

後半
みんな緩んできて実効再生産数が増える
既存株の実効再生産数を0.85くらいに増える
変異株の実効再生産数を1.25くらいになり急拡大

今週
一部地域で変異株が主流株に入れ替わる

来週
多くの地域で変異株が主流株に入れ替わる
0558１３２人目の素数さん
垢版 |
2021/03/06(土) 06:20:16.09ID:Mo0cuWGb
このまま実効再生産数1.2超の感染急拡大に入るのか、
それとも、東京では飲食店の営業禁止命令にまで踏み込むのか
0559１３２人目の素数さん
垢版 |
2021/03/06(土) 21:14:14.41ID:E2+4ywhN
ガン細胞みたいに免疫から攻撃されにくいウイルスとか出てくるんじゃないの？
ガン細胞はどうやって免疫から攻撃されないようにしてるか知らんけど
0560１３２人目の素数さん
垢版 |
2021/03/06(土) 21:20:21.55ID:E2+4ywhN
>>544
感染自体は無くならなくても症状が軽ければ実害は少なくなると思う
変異してワクチン効かなくて感染力も強くて毒性も強いと厄介
毒性が強過ぎると感染しにくくなるだろうからウイルス自体の数は減っていくと思う

ワクチン効かない変異ウイルスの毒性が弱ければ対処する必要は少なくなると思う
0561１３２人目の素数さん
垢版 |
2021/03/06(土) 21:21:30.56ID:E2+4ywhN
ウイルスを弱毒化するように変異させる方法があれば解決方法の１つになるかもしれない
0562１３２人目の素数さん
垢版 |
2021/03/07(日) 01:31:41.16ID:+GSW0Dow
そのかわりにワクチン打つんでしょうがね。
0563１３２人目の素数さん
垢版 |
2021/03/08(月) 06:02:38.26ID:21WHt2In
ブラジル・マナウス株P1はかなり強毒化してるもよう
英国株もすこし強毒化
0564１３２人目の素数さん
垢版 |
2021/03/08(月) 12:15:04.96ID:plZzi7vX
イスラエルはほぼ100％接種が進んでるみたいだね。
しかし、新規陽性者数は下げ止まりで実行再生産数
がほぼ1の状態が３週間くらい続いてる。
集団免疫で収束させるってのは無理なのかもね。

一方、死者数はどんどん減少し続けているので、
ワクチンの効果は出てそう。
0565１３２人目の素数さん
垢版 |
2021/03/08(月) 16:42:13.91ID:jYbKWq/s
実効再生産数を0.85と、
実効再生産数を1.25が混在してると
無限時間後、では、
実効再生産数は、1.249999…に
なると、神からのお告げを受信した
0566１３２人目の素数さん
垢版 |
2021/03/08(月) 23:09:47.71ID:E0kujJjZ
今回のワクチンは感染予防というより
症状が出ない重症化しないことを目的としたものだと思うから
陽性の人はあまり減らないかもね
0568１３２人目の素数さん
垢版 |
2021/03/17(水) 06:33:13.95ID:FENjjiOs
識者の反応

437 １３２人目の素数さん 2021/01/21(木) 01:57:40.20 ID:wnrMDA5R
なんか、キチガイに触っちゃったみたい...
怖いわー。
0570564
垢版 |
2021/03/22(月) 22:56:47.64ID:JSsZ9KFp
なんか、一旦下げ止まったからどうなることかと思ったら、
しばらくしてまたぐんぐん下がってきた。
ワクチンの効果と社会の緩みが一時的に拮抗してたのかな？
0571１３２人目の素数さん
垢版 |
2021/03/23(火) 07:28:24.60ID:E3AF3tGK
プロおじ自称医者にも関わらず期待値と応召義務が分からないことがバレてしまうww
0573１３２人目の素数さん
垢版 |
2021/05/28(金) 17:38:08.46ID:F3XksN9d
>>572
下のデータでAB型の母比率の信頼区間の上限は
95%信頼区間で11.1%
99%信頼区間で11.6%だから
上のデータは統計的に有意な差があるんじゃね？
0574１３２人目の素数さん
垢版 |
2021/05/29(土) 07:49:26.96ID:WmcfsSLD
人口10万人あたりの新規感染者数の変化のグラフを見ると
この指標が統計的に優れたものだということが
実証されているように感じる
0575１３２人目の素数さん
垢版 |
2021/05/29(土) 10:40:29.63ID:YRteybDl
「この指標」ってなんのこと？
0576１３２人目の素数さん
垢版 |
2021/05/29(土) 12:25:41.86ID:g7CBv1i/
人口10万人あたりの新規感染者数
これを統計的な指標と呼ぶのは不適切なら
その理由を教えてください
0577１３２人目の素数さん
垢版 |
2021/05/29(土) 19:03:57.69ID:YRteybDl
それは指標ではなく、基本データでしょ。

なんの指標だと言いたいの？
0578１３２人目の素数さん
垢版 |
2021/05/29(土) 19:09:04.26ID:g7CBv1i/
「特性数」の意味で「指標」を使うのは
数学では一般的
その統計的集団の特性を表す量と
単なる感染者数は分けて考える方が良い
統計学では
データと基本データを区別して使うわけでしょうか
0579１３２人目の素数さん
垢版 |
2021/05/29(土) 19:10:21.47ID:YRteybDl
なんだ、おかしなやつか。相手して損した。
0580１３２人目の素数さん
垢版 |
2021/05/29(土) 22:26:50.55ID:WmcfsSLD
>>579
指標という用語の誤用かどうかはそのお答えからは
分かりませんね
いずれにせよ
県境を越えた移動が制限されているという状況で
10万人あたりの新規感染者数というものが
多くの県でほとんど同じカーブを描いて
推移していることは
ほとんど見事と言ってよい
そのことを「指標」という言葉を使って
言ってみただけ
0581１３２人目の素数さん
垢版 |
2021/05/30(日) 00:51:10.11ID:6BrMEpm2
いや、日常語としての「指標」としてまったく意味をなさんだろ。
それが何を知るための、あるいは測るための目安になってるっていうの？
0582１３２人目の素数さん
垢版 |
2021/05/30(日) 06:49:07.89ID:zeN0ktr4
パンデミックにおいて危機的状況の度合いを測るための指標
データに直接現れる数字だから
何かを何かで割ったという
形のものではなっていないが
指標と呼んでよい重要な数字
0583１３２人目の素数さん
垢版 |
2021/05/30(日) 08:52:14.47ID:z9qtX2J+
数値がどうなったら緊急事態宣言を出して
どうなったら解除するのかが判らないと指標にはならないと思う
0584１３２人目の素数さん
垢版 |
2021/05/30(日) 09:03:36.09ID:zeN0ktr4
まさにその目安として現在使われている
0585１３２人目の素数さん
垢版 |
2021/05/30(日) 11:28:00.31ID:6BrMEpm2
人口あたりの感染者数が緊急事態宣言を判定するための指標になる、
という表現なら、それはその通りだと思うが、当たり前すぎて論評
に値しない。「統計的に優れた」という形容がどこから来るのかも
よくわからん。
0586１３２人目の素数さん
垢版 |
2021/05/30(日) 15:26:47.83ID:rRM3r1vD
そのばらつきに一定の特性が認められれば
常に何らかの利用価値がある
これは当たり前だが統計学的な考え方だろう
0587１３２人目の素数さん
垢版 |
2021/05/30(日) 15:51:03.68ID:6BrMEpm2
具体的にどういう特性で、それがなんの指標になるのか
言わないとなんとも言えん。

感染者数の数値が増大してるから感染が拡がってるって
言われても、当たり前すぎてほとんど無意味だよ。
0588１３２人目の素数さん
垢版 |
2021/05/30(日) 18:00:43.76ID:WXSKIfwS
それ自体は当然当たり前だが
「一定の特性」というものが
現れてきた時点で
当たり前ではなくなる
これに注意するかしないかは
人類の運命にかかわるかもしれない
0589１３２人目の素数さん
垢版 |
2021/05/30(日) 18:06:12.78ID:WXSKIfwS
指標らしい指標にするために
勾配をとったほうが良いことは認める
0590１３２人目の素数さん
垢版 |
2021/05/30(日) 23:56:16.72ID:6BrMEpm2
>>588
まったくナンセンス。

注意するもしないも、感染者数は感染状況を示す基本的な
データなんだから、注意するのは当たり前でしょ。
なにわけのわかんないこと言ってんの。
0591１３２人目の素数さん
垢版 |
2021/05/31(月) 08:47:13.88ID:EWJvgwFj
感染者数ではなく
10万人あたりの新規感染者数のことを言っている
0592１３２人目の素数さん
垢版 |
2021/05/31(月) 10:11:29.39ID:ZysYfXir
明確に言えない時点で指標になってない
0593１３２人目の素数さん
垢版 |
2021/05/31(月) 12:15:18.34ID:kO70dltI
>>591
ここでの感染者数というのは、人口あたりの感染者数の
話に決まってるだろ。
0594１３２人目の素数さん
垢版 |
2021/05/31(月) 12:17:44.97ID:7Qi3ka4Q
それを10万人あたりにするか
1万人あたりにするかを決定するにあたっては
統計学的な検討が必要だろう
0595１３２人目の素数さん
垢版 |
2021/05/31(月) 23:17:27.00ID:kO70dltI
>>594
10万人あたりにしようが、1万人あたりにしようが、１人あたり
にしようが、どれでも同じだよ。なんにもわかってないな。
0596１３２人目の素数さん
垢版 |
2021/06/09(水) 11:01:12.34ID:mzuxkCce
>>572
こういうのと大差ない気もするんだが。

問題
緊急事態宣言下の都道府県で
「コロナウイルスは、画数が偶数の都道府県で重症化しやすい」という説は正しいといえるか？
危険率0.05で検定せよ。

画数和 重症者数 mod2
東京 16 62 0
大阪 10 287 0
兵庫 17 78 1
京都 19 24 1
北海道 26 47 0
岡山 11 18 1
広島 15 22 1
沖縄 22 89 0
福岡 21 75 1
愛知 21 87 1

https://i.imgur.com/fuBXOv8.png
0598１３２人目の素数さん
垢版 |
2021/06/10(木) 12:24:41.45ID:rpbURcGM
感度と陽性的中率の調和平均をとった指標って医学では使われないの？
他分野ではF値って名前で使われるみたいだけど。

真陽性/(真陽性+(偽陽性＋偽陰性)/2)　ってな感じで感度、特異度のトレードオフを考慮した指標
0599１３２人目の素数さん
垢版 |
2021/06/12(土) 09:47:26.23ID:pkqO0d37
「ファイザーで139/975万の死亡例、モデルナで0/19万の死亡例らしいですが、統計的にはどちらが安全なんですか？」
と患者に質問されたら、どう答えたらいいか？
0600１３２人目の素数さん
垢版 |
2021/06/22(火) 08:08:19.10ID:aT+HIzsB
>>598
F値は
真陽性/真陽性+偽陽性　だよ
医療統計の基本中の基本


病気の判断は大抵陰性だから脳死で陰性と言っておけば99%当たるわけでそれだと意味ないからね
0601１３２人目の素数さん
垢版 |
2021/06/26(土) 14:18:07.31ID:YwxapBim
◆イングランド公衆衛生局の報告書を分析してみると

ＰＨＥは１８日（現地時間）、変異株ウイルス状況をモニタリングした報告書を公表している。この報告書は毎週出されているが、
今週の報告書の中で「接種状態に伴うデルタ株感染患者の応急治療および死亡」（表４）の部分を見ると、英国人４０８７人が接種を完了し、
１４日経過してデルタ（インド）株に感染した。このうち２６人（０．６４％）が亡くなった。反面、未接種者は３万５５２１人が感染し、このうち３４人（０．１％）が死亡した。この表だけを用いて２つを単純比較し、ワクチン接種（完了）者がデルタ株に感染したとき、未接種者よりも死亡率が６倍以上だという結論を出したのだ。

該当メディアはまた、ワクチン接種者の入院比率が未接種者に比べて高いとも主張した。同じ表４でＰＨＥは接種者のうち入院事例は４０８７人のうち８４人（２．１％）、未接種者の入院は３万５５２１人のうち５２７人（１．５％）という統計を作成した。あわせてライフサイトはワクチンの副反応「抗体依存性免疫増強（ＡＤＥ・ａｎｔｉｂｏｄｙ　ｄｅｐｅｎｄｅｎｔ　ｅｎｈａｎｃｅｍｅｎｔ）」作用を疑っている。ＡＤＥは１回目感染の時に体内で形成された抗体が２番目の感染を防御するのではなくむしろウイルスを助ける現象をいう。
https://news.yahoo.co.jp/articles/bf95706ab8e1e029e9602dcc10c0a6dc1201c021
0606１３２人目の素数さん
垢版 |
2021/07/04(日) 00:32:40.14ID:/kTQBDb5
結局pcrはどのくらい信用できるんだろう
最初の頃は有症状者のみにしてたけど、いまは五輪選手には毎日してるみたい
0607１３２人目の素数さん
垢版 |
2021/07/04(日) 08:13:46.60ID:N2VJ6Cbp
>>606
感染の診断確定はPCRで行う以外にあるのか？
症状や画像からコロナを疑ったがPCR３回めでようやく陽性というのもある。俺の勤務先の器械はCt45回までやって陽性判定。

尿瓶洗浄係は何回洗浄してんだろうな。
0608１３２人目の素数さん
垢版 |
2021/07/05(月) 18:01:31.50ID:TXloK64L
ホントにワクチン接種が原因で死亡なのか？
どなたか計算お願いします！

ワクチン接種後に死亡した高齢者が、ワクチンが原因じゃないか？
という人が多いですが、高齢者はどんどん死んでゆくので、
その数と違うかどうか検証したいです。

75歳以上の人口は、1,849万人
老衰死（ガンや肺炎、脳梗塞、心筋梗塞じゃないもの、つまり原因不明の突然死ともいえる）
は５％程度である。

ということは平均寿命が85歳とすると、1849万人のうち、10年以内に
半数ぐらい死ぬわけだから、924.5万人が10年で死ぬ。
1年間の老衰死（原因不明の突然死）は92万4500人、
1日あたり2532人である。

コロナワクチンを打つのが日に50万人の高齢者だったとすると、
1849万人のうちの50万人（ワクチン注射）と老衰（突然）死の2532人が
ちょうど同じ日にタイミングよく重なる確率を求めればいいわけである。
母集団は同じだから、1849万人のうち、任意で選んだ50万人が、
また別の復元抽出によって1849万人から選ばれた2532人と
同じ人が選ばれる確率が、
ちょうど、ワクチンが原因かどうか不明の高齢者の1日あたりの死者数、
ということになる？？？
どなたか確率を計算してください！
0609１３２人目の素数さん
垢版 |
2021/07/05(月) 19:07:23.56ID:wcfUYNNo
>>601
コロナの重症化や死亡率は年齢依存度が極端に高いので、
年齢構成を考えた比較をしないとまったく意味がない。

ワクチン接種者率は高齢者が圧倒的に高いことや、サンプルの
人数が一桁も違うことなど、いろいろ不可解。

間違った記事なんでないの？
0610１３２人目の素数さん
垢版 |
2021/07/05(月) 19:19:37.59ID:wcfUYNNo
>>608
>ちょうど同じ日にタイミングよく重なる確率を求めればいいわけである。

それは間違ってる。前日がワクチン接種日という条件つきで、
「特定の日に死ぬ確率」だけを考えればよい。

接種日を確率事象ととらえて、ある人がある特定の日をたまたま
接種日として指定してその翌日に死ぬ確率というのならそれで
それぞれの確率の積になるが、無意味でしょ。
0611１３２人目の素数さん
垢版 |
2021/07/05(月) 20:39:45.57ID:dpe6hrc6
インフルエンザワクチン接種後の死亡率と比べると有意差はあると思うな。
0612１３２人目の素数さん
垢版 |
2021/07/06(火) 08:48:41.90ID:Cki/Ij/R
>>610
>>611
ありがとうございます！
0613１３２人目の素数さん
垢版 |
2021/07/06(火) 11:23:10.06ID:17rOyoWN
>>611
何を根拠に？
接種翌日の死亡数は出てるんだから、それと接種人数あたりの
平均的な日々の死亡数を比べればいいだけじゃね？
0614１３２人目の素数さん
垢版 |
2021/07/06(火) 11:40:58.34ID:17rOyoWN
>>608
いろいろ大雑把すぎる見積もりだけど、急性心筋梗塞による突然死
だけにかぎっても60〜80歳で10万人あたり年間100人くらいになる。
https://www.mhlw.go.jp/toukei/saikin/hw/jinkou/tokusyu/sinno05/13-4.html

ある日ワクチンを接種した50万人のうち、１年以内に急死するのは
500人にもなる。つまり、ワクチンの影響がなくても、毎日1人か2人
はワクチン接種の翌日に急性心筋梗塞で死んでも不思議はない。
0615１３２人目の素数さん
垢版 |
2021/07/07(水) 14:27:29.97ID:i+OHYRE6
インフルエンザワクチン接種後の死亡が10年で20人くらいで
コロナワクチンだと200人以上じゃないの？
いずれも公式には因果関係は認められていない数だけど。
0616１３２人目の素数さん
垢版 |
2021/07/08(木) 09:55:04.30ID:1tn2baG8
コロナワクチン接種後の死者数は既に500を越えていてびっくりした。
インフルエンザワクチンの比じゃないな。
0621１３２人目の素数さん
垢版 |
2021/07/15(木) 11:08:08.07ID:gCoMz3US
>>618
ロシアでは78人の治験だった。

>>619
尿瓶ジジイとは職種の言えない医療従事者=尿瓶洗浄係のことである。
0623１３２人目の素数さん
垢版 |
2021/07/15(木) 11:14:48.65ID:gCoMz3US
>>620
母集団の人数もparameterと呼ぶなら母数と言えなくもない。
単に人数でいいと思うけど。
0624１３２人目の素数さん
垢版 |
2021/07/15(木) 11:22:00.04ID:gCoMz3US
疫学問題
1/10000で起こる致死的副作用を1%未満の間違いで
発見できる治験数はいくつか？
0626１３２人目の素数さん
垢版 |
2021/07/15(木) 21:27:58.32ID:QZAno263
>>620
騙される人の母数
だから騙される人の人数とか母比率とかなら
騙される人の分布に影響するパラメータになるかな

確率変数を騙される人数として
その分布を変える要因は何か？

統計に関する理解度がその１つって感じのことを言いたいのかな
0627１３２人目の素数さん
垢版 |
2021/07/15(木) 21:49:50.76ID:XlfWhKhm
コレの話やろ

taroyabuki.github.io/2018/01/15/misuse-of-the-word-parameter/

統計学の話する中での母数はリンク先の4の意味では使わない方がいいし避けるべき
日本語にはちゃんとこの意味もあるので「日本語としては間違ってない」とは言えるかもしれないが
0628１３２人目の素数さん
垢版 |
2021/07/16(金) 10:34:15.33ID:PDXWv2Wr
417 卵の名無しさん[sage] 2021/05/23(日) 14:16:36.03 ID:Sd8uxjIO
日当直とは別に休日発熱外来もやっているので発熱患者は診なくていいのが(・∀・)ｲｲ!!
行政から補助金がたんまり入るらしいが、勤務医に無理に押し付けると反発を買うので理事長みずから発熱外来をやっている。
銭ゲバと悪口をいう人もいるけど、こういう姿勢は交換がもてる。

>交換がもてる。

>交換がもてる。

>交換がもてる。

>交換がもてる。

>交換がもてる。

相変わらず誤字だらけでかなり耄碌してるし、化石のような絵文字使っててもう社会との関わりも何十年単位でないと思われる。

脳みそ交換してもらったら？
0629１３２人目の素数さん
垢版 |
2021/07/16(金) 15:08:10.64ID:z+sKcmDh
4週間後に東京は2000人越えとか新聞に出てたけど
8週間前に出せる予測値ではなかったのか？
0630１３２人目の素数さん
垢版 |
2021/07/16(金) 15:21:06.89ID:qOHYgDdd
ここ一ヶ月くらい綺麗な指数関数的増加で３週間で倍増くらいの
ペース。むしろ大阪のほうが増えかたは激しい
0631１３２人目の素数さん
垢版 |
2021/07/16(金) 16:14:27.69ID:ZRHF0RQK
指数関数で予測できる範囲のことなら時宜に等しいと言わざるをえない
0632１３２人目の素数さん
垢版 |
2021/07/16(金) 16:23:23.67ID:ZRHF0RQK
訂正
時宜に→児戯に
0634１３２人目の素数さん
垢版 |
2021/07/16(金) 23:51:31.41ID:by6YM2Ls
反日勢力による集団未必の故意コロナテロ発生中
札幌の反ワクチン団体の事例は氷山の一角
0635１３２人目の素数さん
垢版 |
2021/07/17(土) 17:23:46.21ID:ghfDoeWH
>>631
だれが予測の話をしてんだよ、バカ。

直近の期間できれいな指数関数的増加をしてるという素直な感想だろ。
いつまでも続かないのは当たり前。
0636１３２人目の素数さん
垢版 |
2021/07/18(日) 09:10:56.58ID:MW7lfOOz
>>635
だったら専門家ぶるのは予想
0637１３２人目の素数さん
垢版 |
2021/07/18(日) 18:35:59.98ID:buaO2lMU
おまえが勝手にそう思い込んでるだけだろ。
思い込みの強いバカって、迷惑だな。
0638１３２人目の素数さん
垢版 |
2021/07/18(日) 18:59:13.91ID:QcN72bRP
さて、明後日の数字に注目
0640１３２人目の素数さん
垢版 |
2021/07/24(土) 06:02:03.03ID:i1o8XLoL
>>639
尿瓶ジジイとは職種を言えない医療従事者のことである。どうもシリツ卒らしい。
0643１３２人目の素数さん
垢版 |
2021/07/28(水) 10:39:35.86ID:B5sfT9UY
俺は尿瓶洗浄は業務じゃないな。
BPHによる尿閉に導尿することはあっても尿瓶洗浄はしないな。
ライセンスが必要な業務じゃないから。
0645１３２人目の素数さん
垢版 |
2021/07/31(土) 21:55:54.65ID:5eL2DP8M
老害とよんで高齢者に経緯を払えない人間って親の愛情を十分に受けなかった気の毒な人間なんだろうなぁと思う。

>経緯を払えない

>経緯を払えない

>経緯を払えない

>経緯を払えない

>経緯を払えない

日本語もろくに書けない人間って親から義務教育すら十分に受けさせてもらえなかった気の毒な人間なんだろうなぁと思う。
0649１３２人目の素数さん
垢版 |
2021/08/02(月) 12:01:30.91ID:AhPOqdgh
尿瓶ジジイとは職種を言えない医療従事者のことである。どうもシリツ卒らしい。おまるの洗浄も業務みたいだな。
0651１３２人目の素数さん
垢版 |
2021/08/16(月) 11:02:44.46ID:IJhGJi3W
昨年はコロナ対策のおかげか、前年より死者数が減少するという
異例な年であったが（つまり、コロナによる死者数を相殺して
あまりある死者数減の恩恵があったようだ）、昨年後半からは、
コロナ死者数も増えて、前年より大幅に死者数が増えそう。

厚労省の最新統計をみると、前年６月から翌年５月までの１年間
の死者数はその１年前より3万5千人増えいている。この間のコロナ
による死者数は1万1千人なので、2万4千人の超過ということになる。

しかし、高齢化でほぼ年に2万人くらいの割合で死者数が増えていた
ここ５、６年の流れの中で考えるとさほどおかしな数字ではなさそう。
0654１３２人目の素数さん
垢版 |
2021/08/21(土) 12:29:34.22ID:/+lqDxvh
月毎の全死亡者数を比較したらコロナ前後で死亡者数の増加は有意なんですかね？
誰かやってみて
0655１３２人目の素数さん
垢版 |
2021/08/23(月) 01:12:07.75ID:4NOoHz58
今年は増えてんじゃないの？去年は減ってたけど。

ところで、横浜の新市長はまさに統計学に詳しい人（笑）がコロナを
語ってたら「コロナ専門家」として当選しちゃったみたいね。
医学部教授っていうからてっきり医師だと思ってたら、データサイエンス
が専門の人だとか。
0657１３２人目の素数さん
垢版 |
2021/08/23(月) 10:30:32.87ID:4NOoHz58
統計学の権威では絶対にないと思うけど、医療統計という分野では
権威かもね。しらんけど。
0658１３２人目の素数さん
垢版 |
2021/08/23(月) 10:30:56.96ID:4NOoHz58
あ、もしかして皮肉か？ >>656
0659１３２人目の素数さん
垢版 |
2021/09/22(水) 22:13:15.02ID:0ZH8Ymdz
>>61
接種日から死亡日までの統計がこれ
これでも「ワクチンのせいじゃない」と言い張るワクチン推し医とワクチン信者
でもこれ見て医者はさすがに黙ったのに
確率統計習ってない中卒のワクチン信者だけが「医者の報告が接種日から遠のくたび減る、偏ってるからだ」と
自ら『医者がワクチン死を認めず誤魔化してる』と言ってるパラドックスに気づいてない
じゃあホントはもっとワクチン死多いんじゃね？
https://pbs.twimg.com/media/E7RbEEQVoAEKNQA.jpg
0661１３２人目の素数さん
垢版 |
2021/09/24(金) 11:25:55.87ID:wJeoKQ7S
>>659
ワクチン接種した日にご飯を食べているだろうから
ご飯食べてからの日数とかでも同じになるんじゃね？
0662１３２人目の素数さん
垢版 |
2021/09/25(土) 02:04:26.94ID:yo/uOMp1
>>659
ワクチンの反応にて死ぬ人はいるだろう
ワクチンというのは、その利益が不利益より上回ってるから打つわけであって、
そのことを理解しようとしない君のような馬鹿が
世の中をひっかきまわしてるんだよ。わかる？
君はすごーく馬鹿なんだよ。

馬鹿すぎて、自分が馬鹿だってわからないだろ？
0663１３２人目の素数さん
垢版 |
2021/09/26(日) 10:03:44.99ID:tJw5/+iI
1年間の日本のコロナ死は1万人余。ワクチン打ったら死ななかったのが何人いるかね？
一方、ワクチン死は半年で1155人報告。1年間で2千数百人報告される見込み。
医者は報告義務がないから、ワクチンと関係なさそうな死亡にわざわざ報告書を書く動機はない。
関係ありそうでも死亡したことを接種医が全部知るわけじゃないし、知ったとしても忙しければ報告しないだろう。
日本でコロナ死とワクチン死のどっちが多いかね？　
0664１３２人目の素数さん
垢版 |
2021/09/26(日) 10:39:06.97ID:wZBjSjOo
>ワクチン死は半年で1155人報告。
それ、ワクチンとの因果関係は認められてない、単なる時系列的に
接種後にたまたま亡くなっただけかもしれない人数だから。
そういうのを「ワクチン死」などと言ってる時点で、オツムの程度
が知れる。というか、悪質なデマ。
0665１３２人目の素数さん
垢版 |
2021/09/26(日) 11:12:06.20ID:tJw5/+iI
>>664
それを言ったらコロナ死もコロナ感染後の死者だよ。
たまたま死んだ人が病院搬送後に感染がわかったケースが含まれている。
0666１３２人目の素数さん
垢版 |
2021/09/26(日) 11:19:22.87ID:tJw5/+iI
どこのスレでもそうだが、ワクチンに疑問を持つ人に対してワクチン推進派の暴言や罵詈雑言が酷いな。
下品極まりない。どういう連中が書き込んでいるんだろうね？
0667１３２人目の素数さん
垢版 |
2021/09/26(日) 11:27:31.89ID:tJw5/+iI
ワクチン接種後の死亡についてほぼ全てについて因果関係が否定されていない。
0668１３２人目の素数さん
垢版 |
2021/09/26(日) 11:40:48.24ID:UPviXamU
Prob(接種あり│コロナ感染∧死亡）＜Prob(接種なし│コロナ感染∧死亡）
∴
Prob(コロナ感染∧死亡│接種あり）＜Prob(コロナ感染∧死亡│接種なし）
∴
Prob(死亡│接種あり）＜Prob(死亡│接種なし）

というトンデモな条件付き確率の
誘導をしたいのだろう。さてと
ワザワザ誘導するウラ側を察する
つまり、因果関係は不明でないのだ
0669１３２人目の素数さん
垢版 |
2021/09/26(日) 15:08:31.33ID:nA3+2fZa
>>659
接種 0〜5日の死亡数★と
接種 25〜30の死亡数☆の比率、つまり
★÷☆がパッ見10倍位ある感じだ。

然るに
ワクチン接種で死亡確率が10倍以上で
あり、死因（根本的死因）はワクチン
でよるものであるのは自明である。

ワクチン接種で、お無くなりなられた方にお悔やみ申し上げる。

特に接種後二週間は死亡リスクが危険
さらにその後も、ポアソン分布のようにな確率分布で、命が無くなる。

これより、ワクチンの成分は、
生分解されないい異物が血管内でを
グルグルするロシアンルーレットの
状態であろう。

血栓できやすい人が、爆弾抱えて
いつ発症するか分からない
そんな感じである。
0670１３２人目の素数さん
垢版 |
2021/09/26(日) 15:32:12.99ID:/hejQo8M
>>659
ま、不幸にして死ぬ人はいるよ
でも、だから打たなきゃ死なない、というのも誤りだけどね
ワクチン打たないと言い張ったあげく、コロナにかかって死ぬ
哀れな人が沢山いるからね　調べて見　バイアス抜きで
0671１３２人目の素数さん
垢版 |
2021/09/26(日) 15:37:05.42ID:/hejQo8M
>>669
>ワクチンの成分は、
>生分解されないい異物が血管内をグルグルする
>ロシアンルーレットの状態であろう。

科学用語使えば妄想も真実になると思う素人ってイタイタシイね
0672１３２人目の素数さん
垢版 |
2021/09/26(日) 16:14:56.76ID:KWgwjBUI
>>671　噂の酸化グラフェンかもよ
0673>>671の続き
垢版 |
2021/09/26(日) 16:23:50.45ID:KWgwjBUI
グラフェンは、もしかすると
オキシダーゼによる過酸化水素で
分解される可能性もあるが、相当な
炎症を起こすかもね。それ以前に分解
されない。あとステンレスだとしても
いつか頸動脈を突き破ればそこで命は亡くなる。時間の経過と死亡人数が
ポアソン分布になるのも当然と言える
もっとも、ステンレスはワクチン成分
でなく異物（コンタミという）だが
酸化グラフェンはウラ成分である可能性
が濃厚だ。ウラ成分がトッピングされてるかは目視では分からないと思う

それを文系クンにも理解出来る用語に
翻訳したのが「ロシアンルーレット」
というワケ。ルーレットだから
基本的な確率が理解できれば意味が
理解出来ると思うよ
0674１３２人目の素数さん
垢版 |
2021/09/26(日) 16:35:48.93ID:wZBjSjOo
>>665
そういうトンデモ発言をするからバカにされる。

はっきり感染と因果関係のある死者は認定されてるし、その
占める割合も相当高いんだから同列視はできない。
１か０かじゃないんだよ。
0675１３２人目の素数さん
垢版 |
2021/09/26(日) 16:38:46.80ID:wZBjSjOo
>>666
反ワクチン派の暴言もかなりなものだから比較する意味はない。
そもそも、下品とか上品とか、そんなに気になるか？よっぽど
お上品な育ちでいらっしゃるのねw
0676１３２人目の素数さん
垢版 |
2021/09/26(日) 16:41:38.82ID:/hejQo8M
>>673
生物学も化学も分かってなさそうな素人の妄想

精神科で診てもらうことを強くお勧めします
0679１３２人目の素数さん
垢版 |
2021/09/26(日) 16:46:15.19ID:wZBjSjOo
>>669
接種後の時間経過が長いとどんどん網から漏れてくだけなんじゃね？

コロナ以外の原因で死んだ人が、ワクチン接種を受けてたかどうか
いちいち確認するとは思えん。

そもそも、毎日一定割合で人が死ぬんだからワクチン接種とは
無関係な死者が一定数つねに発生してないとおかしい。つまり、
平常時の死亡割合から推定される死者数と比較しないと意味ない。
それより統計的に有意に多いかどうかで議論すべき。
0680１３２人目の素数さん
垢版 |
2021/09/26(日) 16:54:16.66ID:tJw5/+iI
>>674
ワクチン死は医者がワクチンと関係ないと思ったら報告しないし報告する動機がない。
コロナ死は死者がコロナ感染していたら報告することになっている。

もちろんコロナ感染者が肺炎などコロナ感染者に見られる症状があればコロナが原因だろうと医者は考える。
同じようにワクチン接種するまで元気だった人が接種後数日で死ねばワクチンだろうと考える。考えばきゃ報告しない。

両方とも偶然じゃなさそうなのがほとんどで、しかし両方とも死因が科学的に確定されてるわけでもない。　
0681１３２人目の素数さん
垢版 |
2021/09/26(日) 17:18:39.80ID:tJw5/+iI
ワクチンの感染予防効果は確認されていない。厚労省のホームページにもそう書いてある。
実際、ワクチン接種率と全く無関係に実効再生産数が変動してる。
ワクチンの感染抑制効果は統計的に検出はできない。
0682１３２人目の素数さん
垢版 |
2021/09/26(日) 17:40:54.54ID:7v/foW4/
>>677　君のは治らないと思うよ
0683１３２人目の素数さん
垢版 |
2021/09/26(日) 17:42:50.57ID:/hejQo8M
>>681
ワクチンに感染予防効果が全くない、とは証明されていないｗ

むしろ、高齢者の死亡率がロコツに下がってるのは
ワクチンの感性予防効果および症状軽減効果
と考える以外ないｗ

そして、反ワクチンの連中がコロナに感染してくたばるのを見ると
「ワクチンを打って死ぬより、
　ワクチン打たずにコロナにかかって死ぬほうが
　遥かに確率が大きい」
ということがよくわかる

ま、コロナにかかって悶死したいなら止めないがね
マゾヒズムは変態だが、止める理由はないｗ
0684１３２人目の素数さん
垢版 |
2021/09/26(日) 17:43:11.18ID:7v/foW4/
>>676　追伸
長文ご閲覧ありがとう😜
0685１３２人目の素数さん
垢版 |
2021/09/26(日) 17:43:56.28ID:/hejQo8M
>>682
ボクは統合失調症じゃないからな　治るとかいう以前の問題ｗ
0686１３２人目の素数さん
垢版 |
2021/09/26(日) 17:52:00.06ID:tJw5/+iI
>>683
君は重症化予防の話と混同してるようだね。
感染予防効果が確認できないレベルなのは統計データから明らかなんだよ。
実効再生産数の変動にワクチンの影響は検出できない。
検出できない程ワクチンの感染予防効果が低い。あるいは他の影響が大きすぎる。
0687１３２人目の素数さん
垢版 |
2021/09/26(日) 17:54:02.78ID:7v/foW4/
>>683
コロナ　→　肺炎　→　酸欠　→　悶死
と考える人も多いけど、
酸欠でも血中の二酸化炭素の排出が
阻害されない場合には、苦しまない
というワケで、
コロナ　→　肺炎　→　酸欠　→安楽死
という可能性もある。
なお、ワクチン接種により、
コロナで死亡したとの情報はないから
モチロン、厚労省の見解は疑わしい
それも極めて特にだ。
WHO 　米国CD 武田・モデルナそして
厚労省のデマ情報に騙されるのは
いないとは思うが、未知なるテクノロジ
でワープしたワクチンは未知なる
毒性を如何に隠蔽🉑するか楽しみだ

それにしても、レベル低いな
0688１３２人目の素数さん
垢版 |
2021/09/26(日) 18:00:17.46ID:7v/foW4/
それは、それとして、デルタ株の
基本再生産指数R0って2より遥かに
大きいだろ。8とかのハズのようだ
水疱瘡レベルの感染力　なのに
最近減少に転じてキタァァァ

だからさぁ、厚労省もさ、水疱瘡の
ワクチン廃止すればいいと思うよ。
てか、チミ達は厚労省の工作員？
米CDCの工作員？、多分モデルナの
工作員なのかな。お仕事いつもお疲れ様😺
0689１３２人目の素数さん
垢版 |
2021/09/26(日) 18:03:43.24ID:7v/foW4/
>>685　話は通じてるし、統合失調症
ではないと思うよ。
どっちかというと、厚労省とワクチン
の説明の内容が糖質的な文章だね。
専門用語が多いから判別つかないけど
0690１３２人目の素数さん
垢版 |
2021/09/26(日) 18:05:17.40ID:7v/foW4/
あ、糖質で閃いた。政府広報TvCM
「あなたの大切な人の為にワクチン接種検討」というセリフ
これどう考えても、糖質だろ。
0691１３２人目の素数さん
垢版 |
2021/09/26(日) 18:26:32.23ID:tJw5/+iI
まあ関連妄想的な台詞だが、
「ワクチン打てば感染しないしさせない」とデマを言うわけにもいかないから、
視聴者が勝手に誤解するように仕向ける詐欺の手法だね。
0692１３２人目の素数さん
垢版 |
2021/09/27(月) 19:54:04.58ID:hJMWf7Xa
こいつスゲー馬鹿だな。>ID:tJw5/+iI

ワクチン懐疑派には、マトモなやつはいないのか？
0693１３２人目の素数さん
垢版 |
2021/09/27(月) 19:55:43.55ID:hJMWf7Xa
あ、こいつはもっと馬鹿。>ID:7v/foW4/
ってか狂ってる。
0694１３２人目の素数さん
垢版 |
2021/09/27(月) 19:55:47.37ID:bbho3tUq
陰謀論でワクチン打てず次々死んでいくパヨク達
0695１３２人目の素数さん
垢版 |
2021/10/20(水) 18:38:56.51ID:dJBDWcuJ
期日前投票

小選挙区は立憲、比例区は共産　にした

あとは創価学会員の投票を妨害せねば
0696１３２人目の素数さん
垢版 |
2021/10/28(木) 16:11:57.83ID:PXorMPhs
おまえ、共産党や立民が政権とって政治が立ち行くと思ってんのか？
頭おかしいだろ。
0699１３２人目の素数さん
垢版 |
2021/11/10(水) 21:20:02.58ID:81odK9Zx
僕ンとこは実効再生産数が、また１を超えました；;
0700１３２人目の素数さん
垢版 |
2021/11/11(木) 22:07:27.85ID:jQZ/cZe6
そろそろ次の波が始まる頃か。また正月あけにピークだろ。
0701１３２人目の素数さん
垢版 |
2021/11/12(金) 11:05:51.78ID:a4AqLFYV
ドイツは一日5万人越え
死者は200人以上
0702１３２人目の素数さん
垢版 |
2021/11/12(金) 11:37:04.14ID:1m8m55Ln
ワクチン関係無く周期的に増減してる説に一票
0703１３２人目の素数さん
垢版 |
2021/11/12(金) 14:18:54.10ID:b3Drxw85
こういうのって、通常の感染数理モデルでは再現できないよね。
結局、感染拡大期と収束期の粗い近似にしか役に立たないような。
0704１３２人目の素数さん
垢版 |
2021/11/13(土) 06:37:21.56ID:S5vJt4eU
インフルエンザの季節性変動って数理モデルで説明されるんだっけ?
0705１３２人目の素数さん
垢版 |
2021/11/13(土) 11:30:10.86ID:adeg5Qhk
感染拡大から集団免疫獲得して感染収束にいたるという
単純なモデルでは現実と合わんよね。
0706１３２人目の素数さん
垢版 |
2021/11/13(土) 12:57:11.68ID:KQwrHI70
99%以上の確率でこの感染者数を超えないとかでも準備するリソースの目安にはなるんじゃね？
0707１３２人目の素数さん
垢版 |
2021/11/13(土) 13:02:21.26ID:mtRXTHpu
検査が現実認識を歪めてる
0708１３２人目の素数さん
垢版 |
2021/11/13(土) 13:13:53.75ID:KQwrHI70
検査して陽性だとわかった人の中で
入院が必要になったり重症化する人の割合を
基礎疾患とか年齢や性別の属性の条件付で計算したら
準備とか対策の参考になるだろう
0710１３２人目の素数さん
垢版 |
2021/11/14(日) 13:22:46.37ID:s64In/86
>>708
行動パターンも知りたいな
手洗ってないとか
声がでかいとか
0712１３２人目の素数さん
垢版 |
2021/11/16(火) 20:32:35.56ID:0krwOF83
>>701
それ未接種者間の感染が大半なんだよね
特にアルプス以北は厳寒の冬で、通気が悪い家屋・施設なうえに、ハグやキスが習慣化し、しかも土足で住んでる地域だと、ワクチン以外での感染抑制は難しいな
0713１３２人目の素数さん
垢版 |
2021/11/28(日) 06:07:34.48ID:FLFMtEuF
オミクロン株の動向を予測せよ
0714１３２人目の素数さん
垢版 |
2021/11/28(日) 11:19:59.83ID:9vAiZBlq
デルタ株がどうのと騒いでたら、いつの間にかミュー株が
出てきて、、、と思ってたらもうすでにオミクロンかよ！
オメガまで１年ももたんぞ。一周したらどうすんのかね？
0715１３２人目の素数さん
垢版 |
2021/11/28(日) 12:54:43.48ID:9vAiZBlq
なんだ、ニューとグザイは跳ばしたのかよ。
ふざけんなよ！＞WHO
0717１３２人目の素数さん
垢版 |
2021/11/28(日) 22:54:38.16ID:KCt4/hpL
ガセかも知らんがνは英語のnewとかぶり、ξは習Xiとかぶるから
というがオミクロンはオーやゼロに見えてわざとらしい
0721１３２人目の素数さん
垢版 |
2021/12/05(日) 11:14:44.32ID:BTu/4R0O
テドロスはほんとどうしようもないな。
WHOにはほとほと愛想が尽きたわ。ろくでもない団体に成り下がった。
0723１３２人目の素数さん
垢版 |
2021/12/06(月) 11:37:39.79ID:Sh0yO9Al
>>722
それくらい、習菌平が世界中で軍事的・経済的な侵略を続け、生物兵器として放った武漢ウイルスは極めて罪深いって事でしょ？
0724１３２人目の素数さん
垢版 |
2021/12/06(月) 23:05:45.64ID:Z8a9cISR
習(近平)→Xi(Jinping)→ξ→クサイ→臭（菌屁）

見事な変換
0725１３２人目の素数さん
垢版 |
2021/12/07(火) 09:11:14.03ID:27qsa/a1
>>700
自分も同じ見解です
忘年会と新年会で増加しそうですね

再来週から1月中は要注意
0726１３２人目の素数さん
垢版 |
2021/12/07(火) 12:52:43.34ID:kDhERtC0
マスクせずに会食するからなぁ。一番の感染源はそこだろう。
韓国の感染爆発もマスク規制外したからじゃね？
0727１３２人目の素数さん
垢版 |
2021/12/07(火) 13:39:22.94ID:27qsa/a1
企業も会社員のワクチン接種を済ませ、会場も貸し切りで対応するようですが、どこまで効果があるのか微妙ですね
0729１３２人目の素数さん
垢版 |
2021/12/27(月) 20:39:43.23ID:q9B4HDE+
新型コロナの影響で感染拡大の予想が大幅に変わったようです

ワクチン接種により、1月の中旬に少し感染者が増えるだけだった予想が、
新型コロナの感染力を追加した結果、1月下旬から爆発的な感染が予想(AI)

感染数を考慮するなら2~3月頃に緊急事態宣言がされそう
ワクチン接種により死亡率は低いようなので、死者数だけを考慮するなら3月頃に宣言が出されると自分は予想してます
0730１３２人目の素数さん
垢版 |
2021/12/28(火) 00:48:13.92ID:j1//kXKm
すでに感染拡大が一週間で２−３倍のペースになってるから、
そのまま外挿すれば、2月上旬にはステージ３、中旬にはステージ４
を超えて緊急事態宣言になりうる（重症化率にもよるけど）ってのは
AIを使わなくても誰でも予想できる。
0731１３２人目の素数さん
垢版 |
2022/01/08(土) 12:25:39.06ID:cnn/Ayz7
自民党支持してるネトウヨって経済成長は放棄してるの？
あるいは失政し続けてきた政党が突然、パフォーマンス上がると信じてるの？
0732１３２人目の素数さん
垢版 |
2022/01/08(土) 13:11:50.00ID:+V76oJv+
自問自答しろよ、バカｗ
0733１３２人目の素数さん
垢版 |
2022/01/08(土) 13:46:59.23ID:esdRA8zR
Ａ君、B君、C君のテストの平均は
７８点で、A君とC君の平均はB君の点数より３点低く、A君の点数はC君の点数より14点低い。
Ａ君の点数は何点ですか？
0734１３２人目の素数さん
垢版 |
2022/01/08(土) 15:41:34.76ID:+V76oJv+
小学生でも暗算で解けるやろ
0735１３２人目の素数さん
垢版 |
2022/01/08(土) 16:20:17.90ID:33QQIc2w
オミクロンの係数が異常過ぎて恐怖しているいるんだが・・・。

ある日突然日光を奪われる、高層建築で暗闇を創造する企業
【権力】日蝕の積水ハウス【日蝕】
0736１３２人目の素数さん
垢版 |
2022/01/08(土) 16:58:56.26ID:+V76oJv+
倍加時間が１日とかそんなとこだもんな。
さすがにこれが続くとは思えんが。
0737１３２人目の素数さん
垢版 |
2022/01/24(月) 08:46:41.82ID:c9+BUFa1
1ヶ月前のAIでの予想よりも新型コロナの感染が広がるタイミングが早く、尚且つ感染率も高いような感じですかね

ワクチン接種率が8割くらいなのに重症者の数が少なく、重症化率も予想よりだいぶ低いようです

メディア関連でも研究者や政府も困惑しているのか、事前に決めていた対応策を打ち出しても感染拡大に歯止めが効かなく、死亡者数(死亡率)と重症者数(重症化率)が低く、緊急事態宣言を出すタイミングに苦慮しているのがみられます
0738１３２人目の素数さん
垢版 |
2022/01/24(月) 08:55:42.17ID:tuabMKon
重症化率が低いってのは確かなんだろうけど、
桁違いってほどでもないのでは？

年齢によって重症化率が大幅に違うことや、重症化するまで
に１週間以上の時間がかかることを勘案すると、単純に
新規重症数/新規感染者数とやってしまうと、大幅に重症化を
過小評価することになるので、気をつけたい。特に、倍加時間
が極端に短いオミクロン株だとなおさら。
0739１３２人目の素数さん
垢版 |
2022/02/02(水) 10:21:53.30ID:oCb1GOvZ
沖縄県疫学統計・解析委員会は、新型コロナウイルスワクチンの接種回数が多いほど、
どの年代でも入院する割合が低いとの調査結果をまとめた。

入院リスクが高い８０歳以上の感染者を見ると、未接種か１回接種の人は５４・２％（７２人中３９人）が入院していたが、
２回接種では３６・０％（４２０人中１５１人）、３回接種では２２・９％（３５人中８人）だった。

https://www.okinawatimes.co.jp/articles/-/904206#:~:text=%E3%83%AF%E3%82%AF%E3%83%81%E3%83%B3

https://www.okinawatimes.co.jp/mwimgs/0/3/-/img_035e382f25735ba2e904bd12a84c1fdd216866.jpg
0742１３２人目の素数さん
垢版 |
2022/02/23(水) 10:46:25.20ID:e6ewm245
英コロナ規制撤廃のニュースを今の日本で放送するのは不味いです。日本の専門家と連携連絡がとれてなかったのかな

感染拡大を防ぐ為に改めてまん延防止等重点措置対策を呼び掛ける必要があると思います。日本での感染者数が落ちつくまでは、一人ひとり徹底した対応が求められます
0743１３２人目の素数さん
垢版 |
2022/02/23(水) 13:01:05.14ID:R8Dl0MjT
『厚労省さん。多すぎる超過死亡の原因を隠蔽ですか？』
://twitter.com/akahiguma/status/1469462333986861059?s=20
国会議員は、打つと死ぬとわかっているのでワクチン接種をしないのですか？
://twitter.com/eGiZm4HfLw1TZX9/status/1478280884013142020?s=20
ワクチン後の長期副反応、後遺症に対しても対策してくれませんか？
健康被害が出ても分からないから国は何もしないよーって。ほんと詐欺にあった気分
https://twitter.com/chaosvice/status/1487153626527059968?s=20&;t=KHfaWKvYQEdPmth49GP9Kg
母親がワクチン副反応で寝込んでます。もう１ヶ月続いてます。
接種した病院では、うちでは、どうする事も出来ないと言われました。
マスコミ関係の方、どうして、重篤な副反応について報道しないのですか？
亡くなった野球選手についても、軽く報じた程度だし。
://twitter.com/s_hisamax_shi/status/1424357221312303109?s=20&t=epzeOwpC7N7dYA-Id5HF5A
大規模薬害が発生していると疑われる徴候がこれだけ出てるのにすべてデマ扱いのNHK
://twitter.com/Yuta64745343/status/1425525922220756992?s=20
全く報道されないし、おかしいでしょ、せめて報道くらいして欲しい。
://twitter.com/someone5963/status/1420151187098136577?s=20
「コロナで亡くなった場合は繰り返し報道されますが、ワクチン死の場合はなぜかほとんど報道されません。」
://twitter.com/KDystopia/status/1444670342408990721?s=20
こういう大事なのを公表しないで、感染何人だーとか、病床がーとか、止めろよ!!
://twitter.com/2001Kayu/status/1488145997595230208?s=20&t=idsdQFt1H_KhLruikoAa7w

https://twitter.com/5chan_nel (5ch newer account)
0748１３２人目の素数さん
垢版 |
2023/10/07(土) 00:01:10.51ID:TaKrIzpv
(z-z₀)(z-z₀)'=R²
zz'-z₀'z-z₀z'+z₀z₀'=R² 円
zz'+cz+c'z'=R²-|c|²>0
|z|²+cz+c'z'=k>0
0772１３２人目の素数さん
垢版 |
2023/10/12(木) 13:27:51.78ID:KSkJdhHw
|a+b|≤|a|+|b| 三角不等式
√∑(aₖ+bₖ)²≤√∑aₖ²+√∑bₖ²

(a･b)²≤|a|²|b|² シュワルツの不等式
(∑aₖbₖ)²≤(∑aₖ²)(∑bₖ²)
積和²
d²=∑aₖ²=|a|²が基本。
d=√(∑aₖ²)
0833１３２人目の素数さん
垢版 |
2023/10/23(月) 20:04:43.03ID:SW08Krwk
x²≡a modp
解を持てばaはpの平方剰余
解を持たなければaはpの平方非剰余
0871１３２人目の素数さん
垢版 |
2023/12/03(日) 01:12:05.38ID:zGNPPwTX
A∪(B∩C)=(A∪B)∩(A∪C)
a×(b+c)=(a×b)+(a×c)
a+(b×c)=a+bc
a+b\a+c\a
(a∪b)∩(a∪c)
(a+b\a)×(a+c\a)
=a+(b+c)\a=a+b×c
(a+b)×(a+c)=a+bc
(a+b\a)×(a+c\a)=a+(b×c)\a
=a+b×c
0872１３２人目の素数さん
垢版 |
2023/12/03(日) 14:18:22.23ID:zGNPPwTX
A∧(B∨C)=(A∧B)∨(A∧C)
x∈A∧x∈B∪C
⇔x∈A∩B∨x∈A∩C
⇔(A∩B)∪(A∩C)

A∪(B∩C)=(A∪B)∩(A∪C)
x∈A⇒OK
x∈B∩C⇒OK
(A+B\A)∩(A+C\A)
=A+(B\A)(C\A)=A∪(B∩C)
0902１３２人目の素数さん
垢版 |
2023/12/29(金) 00:55:37.99ID:OeEIU6X0
a+b+1
○○、b=-1=e○、b=-a-2○
群をなす
単位元は-1なので逆元も-1に合わせる。
0904１３２人目の素数さん
垢版 |
2023/12/29(金) 01:16:50.91ID:OeEIU6X0
a+b+a×b
○、○
(a+b+ab)￮c=a+b+c+ab+ac+bc+abc
a￮(b+c+bc)=a+b+c+bc+ab+ac+abc
単位元は0、✕
Monoid
0916１３２人目の素数さん
垢版 |
2023/12/29(金) 02:33:57.92ID:OeEIU6X0
n=3の時、、仮定より成り立つ
n≥k≥3の時、成り立つと仮定する
n=k+1の時、最も左にある括弧に着目すると括弧内の元の個数はk個以下なので仮定により括弧は外せる。括弧内の積はGの元なのでそれをAと置いて1個の元とみなせる
するとAを含む全体の元の個数はk個以下となるので他の括弧も全て外せる。
0917１３２人目の素数さん
垢版 |
2023/12/29(金) 03:02:10.98ID:OeEIU6X0
2乗してx軸対称になる
2θ=2π-θ、θ=2π/3、4π/3、
a=-1/2、b=√3/2>0
E、A、A²群をなす
A³=Eである。
0918１３２人目の素数さん
垢版 |
2023/12/29(金) 03:13:24.96ID:OeEIU6X0
a²+b²≠0より
ρR(θ)、0≤θ<2π、ρ>0と置ける
群をなすE=e、A⁻¹=ρ⁻¹R(θ)
0925１３２人目の素数さん
垢版 |
2024/01/16(火) 07:57:39.11ID:JZv7NMsA
rc²s/(r²c⁴+s²)より
s=rc²の時、r²c⁴/(r²c⁴+r²c⁴)=1/2
s=2rc²の時、2r²2c⁴/(c⁴+4c⁴)=2/5
より0以外の値であったとしても連続とは成り得ない。
0926１３２人目の素数さん
垢版 |
2024/01/16(火) 07:58:07.42ID:JZv7NMsA
x²=yの時、x⁴/(x⁴+x⁴)=1/2
2x²=yの時、2x⁴/(x⁴+4x⁴)=2/5より
f(0, 0)をどのように定めても不連続である。
0927１３２人目の素数さん
垢版 |
2024/01/16(火) 08:22:23.07ID:JZv7NMsA
(1-cos(r²))/r²=2sin²(r²/2)/r²
=(sin(r²/2)/(r²/2))²×(r²/2)
→1×0/2=0より連続。
0933１３２人目の素数さん
垢版 |
2024/01/20(土) 02:29:27.34ID:EtTz22xt
1～2001
3→20個
4→15個
12→5個
30個。15、30、45、60、20、40
24×33+9=801個
3、6、９、12、、18、21
4、8、16、
0934１３２人目の素数さん
垢版 |
2024/01/20(土) 02:39:51.51ID:EtTz22xt
17個はmod 5で{0, 1, 2, 3, 4}である。C₀=∅の時、少なくとも1つの類は元を5個含む。C₁～C₄=∅でも同じ。
∅が存在しない時、各類から1個ずつ取れる。
0935１３２人目の素数さん
垢版 |
2024/01/20(土) 03:13:25.59ID:EtTz22xt
(n, 1)、…、(n, n/2)
nは偶数≥2

∑[0, n]=2ⁿ
よってS=偶数-1=奇数。
よって奇数は奇数個ある。
nは奇数≥3
1～(n-1)/2、(n+1)/2～n
同様
0937１３２人目の素数さん
垢版 |
2024/01/20(土) 05:16:35.26ID:EtTz22xt
13→36
14→48
15→12
16→48
17→36
180×2=360

24→48
25→24
26→24
27→48
144×2=288

31→36
35→12
36→24
37→12
84×2=168

41→48
42→48
46→48
47→48
192
1008通り

1010101
1010010
1001010
1001001
0101010
0101001
0100101
7×3!×4!=1008





4、6、8
32-4=28通り
0940１３２人目の素数さん
垢版 |
2024/01/21(日) 00:21:17.64ID:Z5hAZTL9
実行する、回転と鏡映
123321
245542
356653
356653
245542
123321

4個→146、8個→235
4³×8³=2¹⁵=32768通り
0945１３２人目の素数さん
垢版 |
2024/01/21(日) 19:10:04.14ID:Z5hAZTL9
n=3
、9
9-1/2=4
32/4+4/2=10

1176
49-1/2=24
1152/4+24/2=288+12=、300
2k²+2k、(2k+1)²(2k²+2k)
(2k²+2k)(4k²+4k)
2k²(k+1)²+k(k+1)=
k(k+1)(2k²+2k+1)
0948１３２人目の素数さん
垢版 |
2024/02/21(水) 14:37:51.85ID:CS94lxgF
x>0、a∈R、
f(x)=xᵃ…xᵃ
T[1, 0, 0]=x+y+z+w(n-1)!
T[a, a, a, a]=n!xᵃ…xᵃ
T[1, 2]=x¹y²+x²y¹
T[3, 0, 0]=2(x³+y³+z³)
T[2, 1, 0]=x²y+x²z+y²x+y²z+z²x+z²y

階段関数
1000１３２人目の素数さん
垢版 |
2024/03/06(水) 06:58:32.72ID:0HjWUJZx
1∈H∧1∈I
∀x, y∈H: xy∈H
∀x, y∈I: xy∈I
もしxy∉Iとすると
不合理
またH∩I≠∅
なぜならば1∈H∧1∈I
xy∉Hとしても同様
10011001
垢版 |
Over 1000Thread
このスレッドは１０００を超えました。
新しいスレッドを立ててください。
life time: 1309日 20時間 1分 0秒
10021002
垢版 |
Over 1000Thread
5ちゃんねるの運営はUPLIFT会員の皆さまに支えられています。
運営にご協力お願いいたします。


───────────────────
《UPLIFT会員の主な特典》
★ 5ちゃんねる専用ブラウザからの広告除去
★ 5ちゃんねるの過去ログを取得
★ 書き込み規制の緩和
───────────────────

会員登録には個人情報は一切必要ありません。
4 USD/mon. から匿名でご購入いただけます。

▼ UPLIFT会員登録はこちら ▼
https://uplift.5ch.net/

▼ UPLIFTログインはこちら ▼
https://uplift.5ch.net/login
レス数が1000を超えています。これ以上書き込みはできません。
ニューススポーツなんでも実況