善子が犬を拾って、梨子に一時預かれないか頼むが、犬は迷い犬だったという話。 よしりこ（りこよし?）がｷﾏｼﾀﾜｰな回になっている。 ダイヤさんは第4話でダイヤちゃんと呼んでほしいと言っておきながら、今回は出番の少なさもあいまって誰の名前も呼んでない…

ダイヤさんとよそよそしく呼ばれるのが嫌なダイヤさんが、特に下級生たちと交流を深める話。 ダイヤさんのことを最も多く呼んでいたのは千歌。さすがの主人公。 曜はダイヤさんの妄想のなかでエコーがかかっていたのでこれで呼んだ回数が増えている。 ルビィ…

読んだ。 Public Health. 2017 Apr;145:113-119. Nat Rev Rheumatol. 2016 Aug;12(8):486-96. Biosocial Surveys. Nat Rev Genet. 2013 Jul;14(7):483-95. Mendelian randomization は何かというと、経済学などの因果推論では操作変数法と呼ばれるものを統計…

関数で定義されている変数を取り出したいという相談を受けた。 deparse でできるようだ。 deparse(plot) [1] "function (x, y, ...) " "UseMethod(\"plot\")" 逆に、テキストになっている変数を関数に使いたいとも聞かれた。 eval とparse の合わせ技ででき…

ラブライブ予選に出場することと、学校説明会でライブをすることで入学希望者を増やそうと計画していたが、学校説明会がずれたことでラブライブ予選とダブルブッキングになってしまった。 みかん畑のみかん輸送機を使うことで、ラブライブからの学校説明会に…

某とくべつけんきゅういんの採用発表があった。 結果はお察しである。 不採択になった場合は、どれだけ評価が低かったかという数値と、応募者(応募領域内での)の順位が分かる。 このとき、評点は適当な変換がなされて、T スコアというスカラーが返ってくるが…

ラブライブサンシャインの2期が始まった。 昔、名前の呼び合いグラフ(これとかこれ)をやったので、今回は放送リアルタイムに各話やっていこうと思う。 1話目は学校説明会を開こうとしたが、結局廃校騒ぎになってしまい、その裏で奮闘する鞠莉にスポットがあ…

SNP とかなんでもいいのだが、行列データがあって、ある列のデータのうち条件をみたすものだけを抽出して、和を取りたい。 それをR とかPython を使わずに、シェルというかawk で完結したいのだが、という悩みを同僚が言ってきたのでやってみた。 一発でsum …

読んだ。 ClusterSignificance: a bioconductor package facilitating statistical analysis of class cluster separations in dimensionality reduced data Bioinformatics. 2017 Oct 1;33(19):3126-3128. オミックス解析などをすると複数パラメータの高次…

R の関数にcolors というものがあって、これで色の名前指定でデフォルトで入っている色を使える。 プロットしているときに単純にcol=番号指定だと1:8 で終わってしまうし、かと言っていい色ないかなと毎回悩むのでとりあえず灰色以外を一括してプロットして…

推定するパラメータに順序制約があるとき、rstan ではordered、正に限定するならpositive_ordered が使える。 例えばYaegashi (J Hum Genet. 1998;43(2):85-90.)らは、出生時の染色体異常が母体年齢に応じてどう推移するかをデータを取って調べている。 ここ…

読んだ。 Cancer Inform. 2014 Jan 16;13:13-20. オミックス研究では、各染色体でのCNVや遺伝子発現状態、各種統計量の定量値や分布などをずらずらっと描きたいが、横に長く描くといけてないと誰かが思ったのだろうか、丸く描くやり方がある。 これをcircos …

こんな問題を見かけた。 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.com…

読んだ。独習 統計学24講: 医療データの見方・使い方作者: 鶴田陽和出版社/メーカー: 朝倉書店発売日: 2013/05/29メディア: 単行本この商品を含むブログ (3件) を見るすべての医療系学生・研究者に贈る 独習統計学応用編24講 ─分割表・回帰分析・ロジスティ…

注意：本解析の結果と、実際に声優がそうであるかはまったく関係がありません。 援助交際してそうなアニメキャラランキング2016をstanで考える - 驚異のアニヲタ社会復帰への道 援助交際してそうなアニメキャラランキング2015をstanで考える - 驚異のアニヲ…

C92 1日目の「あんかけどうふ」というサークルは、デレマスのみくにゃんなどで有名な杏仁豆腐氏の作品が出るサークルということで、叶姉妹のサークル、けものフレンズのたつき監督と並んで3大徹夜サクチケ転売屋の巣窟、という下馬評があったりなかったり。 …

こんな話を見かけた。 「昔の夏はもっと涼しかっただろ！」と思い、50年前の日本の気温を調べてみたところ、まさかの結果が…！ | netgeek これについて言及していた人がいて、データソースがおかしい、ということのようだ。 気象庁 | 日本の年平均気温 変化…

血液1滴を検査するだけで、各癌腫特有のmicroRNA を検出して癌であるかどうかを95% 以上の精度で判定するらしい。 https://medical-tribune.co.jp/news/2017/0724509666/ ネットでは称賛の嵐のようである。 よくよく記事を読むと、どういう理屈で癌を判定し…

ガンマ分布は、期間 ごとに1回くらい起こるランダムな事象が 回起こるまでの時間の分布としてモデル化されることがあるらしい。 ガンマ分布自体は、shape とscale (もしくは逆数のrate) のふたつのパラメータで決定される。 悪性腫瘍(いっぱんにがん)は、遺…

賭ケグルイを見ている。 黒髪ロングのプリーツスカート黒タイツのはやみんなので視聴意欲がやばい。 ここで、賭ケグルイの1話で、投票じゃんけんという変則じゃんけんで勝負する。 クラス 人がグー 、チョキ 、パー のいずれかの手を投票する。 蛇喰さんと芽…

ダイエットのデータを昔使って遊んだけど、元データを持っている人が更に長期のデータを公開してくれた。 window.twttr = (function(d, s, id) { var js, fjs = d.getElementsByTagName(s)[0], t = window.twttr || {}; if (d.getElementById(id)) return t;…

声優統計コーパスというものがある。 日本声優統計学会 プロの女性声優 3 名が 3 パターンの感情表現で読み上げた音声 2 時間分 を「声優統計コーパス」として無料公開します - 糞糞糞ネット弁慶 音素バランス文という、音声言語研究では非常になんかいい例…

虚構新聞のネタなので、一般人を試験により認定するという話は存在しません。 こんな記事を見かけた。 一般人認定試験、来年度実施を検討 「共謀罪」成立受け 共謀罪について、一般人であるかどうかの認定試験を行って、一般人か否かが決まる。ここで、「合…

R でいい感じのグラデーションを作ろうと思ったらたいていrainbow を使うのだが、rainbow を使うと先頭が赤、末尾が紫なので色が同系統で混同する。 matlab のjet.colors が赤〜青のグラデーションなのでよさそう、と紹介されたのだが、matlab にしかはいっ…

tmux を勧められた。 以前、vim を導入したときにvim-r-plugin を使うのにtmux が必要で、tmux のpane を使えばRstudio っぽく、コードを書く領域とRが実行できる領域が出現して、コマンドショートカットでコードをRコンソールに送ることができる、という話…

5000兆円ではなく、violinplot の上半分だけほしい。 上も下も同じ形なので、どちらかがあればよい。そうすれば、stan の事後分布とかにもいい感じで使える。 たぶんggplot2 にはよさそうなものがありそうだが、vioplot パッケージを流用して使おうと思った…

Rの確率分布をやっていた。確率分布はたくさんある。 Rでデフォルトではいっているものは何があるのだろうと思って探した。R 内の確率分布に関する関数は、累積密度p、確率点q、確率密度d、乱数r が用意されていて、21個とbirthday 分布が用意してあった。た…

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.com/widgets.js"; fjs.paren…

細胞でもなんでもいいが、視野内の物体の動きを軌跡として取得したい。 最近流行りのディープラーニングでは下調べが足りないがいろいろあるっぽい。しかし、これらの手法を使うにしても自前で開発するにしても、正解データと比較して性能評価が必要である。…

親戚が集まったときの会で、ふいにそんなことを言われた。その人としては減量をがんばったということを言いたかったのだろうが、自分自身、糖質制限ダイエットに対してエビデンスを持っていないことと、1日で2kg 減るっていうのは、偶然「2日間連続で体重を…