統計

日本版敗血症診療ガイドライン2020 (J-SSCG2020) の血糖管理でネットワークメタアナリシスをやってみる

The Japanese Clinical Practice Guidelines for Management of Sepsis and Septic Shock 2020 (J-SSCG 2020) 「日本版敗血症診療ガイドライン2020 (J-SSCG2020)」正式版 公開のお知らせ|日本集中治療医学会 日本版敗血症診療ガイドライン2020 (J-SSCG2020)…

中心極限定理を確かめる

どんな分布(注)からでもサンプリングされた標本の平均をさらに繰り返して標本の分布を考えたとき、その分布は正規分布になる。 というのでいろいろ分布を変えて試してみる。 Rでは接頭語rdpqに続いて分布名があるので、標準で実装されている関数をapropos…

症例報告・基礎研究・臨床研究の論文の書き方

小児科診療 2020年 07 月号 [雑誌]診断と治療社Amazon小児科診療 (83巻7号) | 医書.jp 症例報告だけでも医者人生が変わるという華々しい話が面白かった。臨床婦人科産科 2020年 11月号 今月の臨床 論文作成の戦略 アクセプトを勝ちとるために医学書院Amazon…

帰無仮説検定が正しいのか正しくないのか

検定とp値がいつまで経ってもわからないので久々にシミュレーションした。 対照群 と治療群 があって、なんらかの指標となる値に変化があるかどうか調べたい。 比較するのは、母集団(真の分布)である(ここでは真の分布は正規分布であると自分は分かってい…

Beyond ER ドクターカー&ヘリと心肺蘇生

読んだ。BeyondER(ビヨンダー) Vol.2 No.1 2023メディカルサイエンスインターナショナルAmazonBeyondER (2巻1号) | 医書.jp ヘリには乗らないのでドクターヘリやドクターカーの運用やエビデンスを勉強した。院外心停止の蘇生も復習がてら読んだ。ショック不…

小児科診療で読んだもの

小児科診療 2020年 09 月号 [雑誌]診断と治療社Amazon小児科診療 (83巻9号) | 医書.jp 新生児の手術もあるのでそういう子たちがその後どうなるか調べた。小児科診療 2021年 08 月号 [雑誌]診断と治療社Amazon小児科診療 (84巻8号) | 医書.jp 心奇形と発生に…

JOHNS 気管切開と医学統計

読んだ。JOHNS Vol.36 No.2(202 特集:気管切開術マニュアル東京医学社AmazonJOHNS (36巻2号) | 医書.jp COI:なし 気管切開について勉強しようと思った。 甲状軟骨切開法というのがあることを知った。気管切開の歴史なども興味深かった。 基本的には手術室も…

NEJMに載るようなRCTでも再解析したら結果が異なるのですが??

面白い話題を見つけたので考えてみようと思ったらデータもあったのでやってみた。 ことの発端はここで 産科編 妊娠中期 妊娠中、虫歯の治療はできますか? (周産期医学 49巻13号) | 医書.jp 妊婦に歯科治療を行うと37週未満の早産が減らせるかどうかというRCT…

多カテゴリへの回帰

二値ならbinaryだが多カテゴリのときにどうするか、ということでnnet::multinomかmlogit::mlogitが使えるがnnet::multinomのほうが使いやすかった。 m <- nnet::multinom(Species ~ ., data=iris) # weights: 18 (10 variable) initial value 164.791843 ite…

COVID-19と無症候性低酸素血症

書いた。 Asymptomatic Hypoxemia as a Characteristic Symptom of Coronavirus Disease: A Narrative Review of Its Pathophysiology. COVID 2022, 2(1), 47-59. COVID-19において、低酸素血症なのにぜんぜん苦しそうにない、いわゆるhappy hypoxia もしく…

ランダムに昇進を決めることが全体の価値を最も高める、というPeterの法則とイグノーベル賞の話

読んだ。 [0907.0455] The Peter Principle Revisited: A Computational Studyこんな動画を見た。 www.youtube.com 2010年のイグノーベル賞で、Peterの法則というものをシミュレーションで示した、という話。社員地位構造が階層的な場合、昇進するときに能力…

legend 内部で数式を混在させながら左揃えで2行にわたってテキストを書きたい

結論から言うと、substitute 関数を使って数式を記述し、atop を分子と分母について2行にわたって記述したいテキストに対してそれぞれ使う。ただし、分数表記でいうと微妙に高さがずれるので、atop は中央の括線を書かないので気にしない、ということにすれ…

懲役刑と罰金刑の関係

※法律素人なので懲役刑もしくは罰金刑のどちらがよい/悪いを推奨しているわけではありません。実際に刑を受ける場合は弁護士に聞いてください。〇〇年の懲役または△△円の罰金、というフレーズをときたま見かけるが、法律的に一定の指標があるのかどうか気に…

3連PCR

新型肺炎で2回PCRしたけど陰性で、でも感染症の専門家がやっぱり疑わしいからダメ押しでもう一回PCRやって、というのでやったら陽性だった と言われた。検査には感度と特異度という検査特性が存在するので、検査が陽性だからといって本当にその疾患が陽性で…

新型肺炎の変異株で若年者は本当に重症化しやすいのか

結論から言うと、0.うんぬん%以内の増加は40-50歳代の年齢層であるようだが、それよりも80歳以上の高齢者での死亡率増加が大きく、かつ感染者数が増えており重症のままベッドを占拠しておりながらもすぐに死亡してベッドを空ける、というわけではないので、…

Interrupted time series analysis でARIMAを使って真面目にやる

読んだ。 Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventionsInterrupted time series analysis は介入の前後でlevel(増減、つまり切片)とtrend…

メタアナリシスで、サンプルサイズが最も大きいわけではないのに、weightが最大になるのはおかしくないですか?

という質問を受けた。 結論から言うとおかしくない。メタアナリシスのweight は各研究内の分散 と研究間の分散により (fixed model の場合)もしくは(random effect model の場合) で決まるから、分散が小さい、すなわち推定精度の高い研究はweight が大…

判別分析の判別線がおかしい

という相談を受けた。 解析したいデータは2変数なので、例として 判別分析モデルの応用 がついてきたが、3変数の場合はirisを用いて データ解析・マイニングとR言語 もついてきて、こんな感じでやりたい、と言われた。0/1のtype は以下の線形モデルで表現さ…

新型肺炎の年齢階級別死亡率

新型コロナウイルス 国内感染の状況 こちらの2021年1月13日時点の各年代別感染者数と死亡者数のデータがあったので単純に推定してみる。その1 drcパッケージにあるLL.4関数は、4パラメータモデルのロジスティック曲線に当てはめる。 これで簡単にやるその2 r…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 12話

侑と歩夢回だったが孤立したノードがたくさん。 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 4話 - 驚異のアニヲタ社会復帰の予備 虹…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 11話

歩夢のレズメンヘラっぷりが炸裂したが果林としずくが相変わらずハブられていた。 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 4話 -…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 9話と10話

9話 彼方ちゃん孤立 10話 エマ孤立 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 4話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スク…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 8話

しずく回だったがかすみんが活躍しすぎ。 冒頭の侑が歩夢、果林、愛さんを呼ぶシーンがなければこの人たちは孤立するところだった。 彼方さんは孤立した。 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 7話

彼方回だったが姉妹のいざこざを同好会でやられても…と思った。 果林とエマはだいぶ出番ないが大丈夫か? 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイドル同好会の名前呼び…

ハザード比

カプランマイヤーによる生存曲線をやっていてCox proportional hazard モデルによるハザード比を出そうとして、シミュレーションデータで遊んでいただけの話。 prodlim パッケージにあるSimSurv 関数は、変数X1とX2があるが、とりあえずこれらは考えず、まっ…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 6話

りなりー回だったが今回はしずくもグラフに参加していた。 けどやっぱり陽キャギャルの愛さんのことが好きになってしまう。 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイド…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 5話

エマさん回だったわけだが、しずくがはぶられているが大丈夫か? 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 4話 - 驚異のアニヲタ社…

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 4話

愛さんがキーパーソンすぎるわけだが、あんな娘現実におったらこっちが勝手に惚れてしまって振られてしまうやつじゃんェ… 虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話 - 驚異のアニヲタ社会復帰の予備

虹ヶ咲学園スクールアイドル同好会の名前呼び合いグラフをかく 1〜3話

虹ヶ咲学園スクールアイドル同好会が始まった。 ラブライブ!虹ヶ咲学園スクールアイドル同好会 以前、ラブライブのμ'sとAqours のメンバーの名前呼びグラフを作ったので、アニメを観ながら名前を呼んだ回数をカウントし、スクリプトをそのまま流用して同じ…

ソシャゲの生存率

こんな記事を観測した。 【悲報】ソシャゲの「10年生存率」、ほぼ0と判明する | やらおん!こんな面白い解析ネタがあったのに乗り遅れた。 しかし、図の生存曲線がどのように作られたかが謎である。カプランマイヤーなのだろうか。 しかも、10年で「ほぼ」0%…