冴えない彼女の育て方が13巻を持って原作完結した。（FDやGSといったサイドストーリーもあるけど） kindle でいつも1ヶ月遅れで読んでいたけど、今回は待ちきれずにkindleポチーと同時に単行本も買ってしまった。冴えない彼女の育てかた13 (ファンタジア文庫…

某とくべつけんきゅういんの採用発表があった。 結果はお察しである。 不採択になった場合は、どれだけ評価が低かったかという数値と、応募者(応募領域内での)の順位が分かる。 このとき、評点は適当な変換がなされて、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講 ─分割表・回帰分析・ロジスティ…

vim を使っていて、caps lock がESC ならうれしくて、ubuntu 14 まではxmodmap で設定していたが、ubuntu 16 からはこれだとうまくいかないようである。 sudo apt-get install gnome-tweak-tool でtweak から、タイピング > Caps Lock キーの動作にいき、Cap…

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

こちらを参考にインストールする。 これとか。

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

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

昔、windows7 があるPC にubuntu をインストールした。 今回、let's note RZ4 の8GB が安く手に入ったので、windows8 がある環境にubuntu16.04 をデュアルブートできるようにしようとしたら、UEFI で手こずった。 結論からいうとRufus を使ってUFEI対応のUSB…

読んだ。統計的因果推論 (統計解析スタンダード)作者: 岩崎学出版社/メーカー: 朝倉書店発売日: 2015/11/10メディア: 単行本（ソフトカバー）この商品を含むブログ (3件) を見るCOI：研究室にあった。 ランダム化比較試験ができないときの、傾向スコアマッチ…

血液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 時間分 を「声優統計コーパス」として無料公開します - 糞糞糞ネット弁慶 音素バランス文という、音声言語研究では非常になんかいい例…

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

読んだ。 MAGIC: A diffusion-based imputation method reveals gene-gene interactions in single-cell RNA-sequencing data コードはPython で書かれている。computational flowcytometry のDana Peer ラボ。RNA-seq のデータ行列が取るであろう高次元空間…

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

なんでもいいのだがスペースで区切られた、横に長いファイルがあるとする。 seq -s" " 1 20 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 これを行から列に変換して縦に長いファイルが欲しい。というのも、paste コマンドなどを使って結合したいから…

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…