数理モデル
読んだ。 [0907.0455] The Peter Principle Revisited: A Computational Studyこんな動画を見た。 www.youtube.com 2010年のイグノーベル賞で、Peterの法則というものをシミュレーションで示した、という話。社員地位構造が階層的な場合、昇進するときに能力…
読んだ。INTENSIVIST Vol.10 No.3 2018 (特集:人工呼吸器)メディカルサイエンスインターナショナルAmazonINTENSIVIST Vol.4 No.4 2012 (特集:呼吸器離脱)メディカルサイエンスインターナショナルAmazon安かったので呼吸器離脱を買っては積んでいて、某新型肺…
結論から言うと、0.うんぬん%以内の増加は40-50歳代の年齢層であるようだが、それよりも80歳以上の高齢者での死亡率増加が大きく、かつ感染者数が増えており重症のままベッドを占拠しておりながらもすぐに死亡してベッドを空ける、というわけではないので、…
新型コロナウイルス 国内感染の状況 こちらの2021年1月13日時点の各年代別感染者数と死亡者数のデータがあったので単純に推定してみる。その1 drcパッケージにあるLL.4関数は、4パラメータモデルのロジスティック曲線に当てはめる。 これで簡単にやるその2 r…
西浦先生が日本の実効再生産数を推定した。 コードはrstanで下記から取れる。 https://nbviewer.jupyter.org/github/contactmodel/COVID19-Japan-Reff/tree/master/解説動画を見逃したのでコードと関連論文からのお勉強になるが、肝としては、 ・知りたいの…
こんなことをした。 mikuhatsune.hatenadiary.com 集中治療学会が、人工呼吸器を要している患者、ECMOをしている患者、ECMOで死亡した患者、ECMOから回復した患者、など日ベースで公開してる。 これに、毎日の感染者や死亡者のデータをくっつけて、SIRを使っ…
こんなツイーヨを観測した。コロナウィルス騒動で全国のオペ室今どうなってるのかアンケート7回目全国の麻酔科の先生におたずねします。今（4/24）現在、オペ室は、— DAJ (@dajhiroki) 2020年4月24日 コロナウィルス騒動で全国のオペ室今どうなってるのかア…
慶応大学病院で、入院中で（新型肺炎COVID-19について）無症状の患者67人をPCR検査すると、6%（=4人）の患者で陽性だった、という話。 某感染症の先生は、めちゃめちゃ単純な推定を使って、一番楽観的な推定で47万人、一番悲観的な推定で430万人、東京に無症…
諦めていた。 mikuhatsune.hatenadiary.com しかし、なんか出来そうになった。 mikuhatsune.hatenadiary.com（型の）コロナウイルスの流行具合から、相互に免疫があるものとして、新型肺炎の今後数年の流行がどうなっていくかをSEIRSモデルで、かつ2種類のウ…
諦めていた。 mikuhatsune.hatenadiary.com というのも、SEIRSモデルでS→Eになるパラメータ が、時間依存的なパラメータとして で定義されるが、これはrstanではできなさそうである。というのが、integrate_ode関数が時間依存のパラメータをどう頑張っても取…
読んだ。 Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. - PubMed - NCBI COVID-19がこのままだと2022年も続いているのではないか、といろいろ話題になった論文。 githubがおいてあるが、ここにあるデータは GitHub…
読んだ。 A Simulation on Potential Secondary Spread of Novel Coronavirus in an Exported Country Using a Stochastic Epidemic SEIR Model. - PubMed - NCBI COI：筆者はこの著者とは直接の関係はないので、純粋に統計解析のツッコミです。こんなツイー…
読んだ。 Eurosurveillance | Estimating the infection and case fatality ratio for coronavirus disease (COVID-19) using age-adjusted data from the outbreak on the Diamond Princess cruise ship, February 2020 ダイヤモンド・プリンセス号のデータ…
読んだ。 Ascertainment rate of novel coronavirus disease (COVID-19) in Japan | medRxiv ascertainment rate という、感染者数（PCR陽性ベース）がどれくらいか、つまり、1だと実際の報告数が潜在的な患者数と同一で、>1だと過剰に報告されている、軽症…
読んだ。 www.ncbi.nlm.nih.gov 巷を賑わせているCOVID-19だが、厚生労働省がダイヤモンド・プリンセス号のPCR陽性者数を逐一ネットに挙げていたので、この論文にもあるようにそこからデータを取ってきて、COVID-19の感染力を推定しようと思った。 Basic rep…
という質問を受けたので1万年と2000年ぶりくらいに遺伝統計やった。 雰囲気としては、こんな感じの値を求めたい。Indian J Endocrinol Metab. 2014 Nov-Dec; 18(6): 850–854. のTable 1 より引用www.ncbi.nlm.nih.gov ここで、データとしては被験者たちのgen…
フェンタニルの持続静注をすることで術後鎮痛を行うことがあるが、フェンタニルを低用量で開始するだけでは目標とする鎮痛濃度に達するにはものすごい時間がかかる。 例えば、身長160cm、体重50kg、年齢40歳の男性に20ug/hr でフェンタニルの持続静注をする…
Shinyというhtml 上でいろいろパラメータをいじったものを即座に計算して描出するインタラクティブでなんかすごいアプリがあったのだが、rmarkdown と同様に手をつけることなく放置していた。 せっかちな人はこちら → https://yfujii08.shinyapps.io/pkpdsim…
高校サッカーを見ていた。2017年度は前橋育英が初優勝で幕を閉じた。 どの試合だったか忘れてしまったが、2点差がついたときに解説が「2点差は危険」ということを言っていた。 調べてみると、やはりよく言われていることのようだが、実際にデータをとってみ…
C92 1日目の「あんかけどうふ」というサークルは、デレマスのみくにゃんなどで有名な杏仁豆腐氏の作品が出るサークルということで、叶姉妹のサークル、けものフレンズのたつき監督と並んで3大徹夜サクチケ転売屋の巣窟、という下馬評があったりなかったり。 …
読んだ。 MAGIC: A diffusion-based imputation method reveals gene-gene interactions in single-cell RNA-sequencing data コードはPython で書かれている。computational flowcytometry のDana Peer ラボ。RNA-seq のデータ行列が取るであろう高次元空間…
親戚が集まったときの会で、ふいにそんなことを言われた。その人としては減量をがんばったということを言いたかったのだろうが、自分自身、糖質制限ダイエットに対してエビデンスを持っていないことと、1日で2kg 減るっていうのは、偶然「2日間連続で体重を…
こんなツイートがあった。 window.twttr = (function(d, s, id) { var js, fjs = d.getElementsByTagName(s)[0], t = window.twttr || {}; if (d.getElementById(id)) return t; js = d.createElement(s); js.id = id; js.src = "https://platform.twitter.c…
読んだ。 PLoS One. 2016 Dec 9;11(12):e0167480. doi: 10.1371/journal.pone.0167480. J(^o^)PAN I/II という顔文字の入った臨床研究なんて初めて見た!! とかそういう話はどうでもいいです。 救急外来では、午前中に受診者がいないと、ふいに誰かが「今日は…
読んだ。 Proc Natl Acad Sci U S A. 2004 Sep 7;101(36):13285-90. antimicrobial cycling という概念がある。細菌感染症に対して抗生物質を使うのだが、うまく治療していても耐性菌は生じることもあるし、適当に使ってしまえばさらに耐性菌のリスクがあが…
multi-armed bandit problem というものがある。 元々はギャンブラーがどのスロットマシーンをやれば、利益を最大化できるかという話から始まったらしいが、最近は教科学習として流行って、webコンテンツのクリック率だとかそういう分野で使われている。 こ…
アレルA とアレルa がある。 アレルA の存在確率がのとき、アレルa の存在確率はである。 1世代交配して得られる遺伝子型はである。 これは、アレルA とアレルA を持つ者同士は強制的に子孫残してね、というような神の見えざる手が働かない、つまり、完全に…
読んだ Behav Ecol Sociobiol. 2007 Dec;62(2):245-253 dimorphic allometry という二形性の成長を見ていて、その変化点らしきところを推定する。 論文ではT. dichotomus septentrionalis と呼ばれる、いわゆるカブトムシの体の大きさ(body) と角(horn) につ…
読んだ。BMC Cell Biology201011:24 細胞の移動追跡を自動化してやってくれるソフトウェアTime Lapse Analyserの話。 マニュアルでやると細胞全体はもちろんチェックできないし、代表的な20細胞を選んでtracking してもらったら軒並みAMD (average mean dist…
Interrupted time-series analysisという解析法をやろうと思ったのでやってみたメモ。 J Clin Pharm Ther. 2002 Aug;27(4):299-309.(PDF)を読めばいい。 データとしては、次の図のような日時と処方数のデータについて、とあるイベントが起きたときのそのイベ…
