2013-10-01から1ヶ月間の記事一覧

薬物依存

覚えていなくていつも間違えるからメモ。 依存の型 身体依存 精神依存 耐性 代表的な薬物 アルコール ++ +++ + モルヒネ +++ +++ +++ ヘロイン,コデイン バルビツール酸 +++ ++ ++ フェノバルビタール,睡眠薬 コカイン - +++ - 大麻 - ++ - マリファナ,ハシ…

Multiple imputation

欠損値をどうにかこうにか補おうということでimputationがある。 欠損データ分析には欠損値のメカニズムとして MCAR, MAR, MNAR を挙げている。リンクのサンプルデータとして、欠損値を含む変数, 解析上の変数, 別の変数, 欠損値かどうかの確率変数(欠損値な…

就職活動の結果

就職活動の結果が出たわけだが、やはり都市部では人が多く、東北日本海側で人手不足だった。 愛知で50%のフルマッチ、宮崎県ではフルマッチしたプログラムは0件だった。 前回の中間発表で熊本県が謎の人気だったが、今回もフルマッチ病院が多めだったようだ。…

Breathtaking Journey

NEJM. 2012 Aug 2;367(5):452-7. いろいろな感染症の潜伏期一覧 バックナンバー

機械学習の性能評価

機械学習の評価として、Recall, Precision, F値, Matthews correlation coefficient (MCC)があるので、Rでできる機械学習をひたすらやって性能評価をしてみる。 今回は交差検証は行わず、データセットの全てのデータでモデルを構築して、それを元のデータに…

虫食い算

面白い問題を教えてもらった。 虫食い算をRで総当りで解くらしい。 S E N D +) M O R E M O N E Y 行列構造を保って演算するようにがんばる。 f1 <- function(){ library(gtools) perm <- permutations(10, 8, v=0:9, set=TRUE, repeats.allowed=FALSE) # パ…

2変数プロット

bagplotというのを聞いた。aplpackパッケージでできるらしい。 library(aplpack) dat <- cbind(rnorm(100)+100, rnorm(100)+300) dat <- rbind(dat, c(105, 295)) plot(dat) bagplot(dat,factor=2.5,create.plot=TRUE,approx.limit=300, show.outlier=TRUE,s…

At a loss

NEJM. 2012 Jul 5;367(1):67-72. Gloves and Stocking 型の感覚鈍麻になる疾患 DANG THERAPIST D:DM A:alcohol N:nutrient G:GBS T:tumor H:hereditary CMT E:endcrine R:renal A:amyloidosis P:porphirin I:HIV, spiroheta, Lyme S:salcoidosis T:toxic バ…

ネットワークのコミュニティ

ネットワークのコミュニティ検出をやろうやろうと思って長らく放置していたのでやる。 重複がないコミュニティはigraph, 重複があるコミュニティはlinkcommでできるので、重複があるほうをする。 linkcomm クラスというのができるので扱いにくいかと思ったら…

サンプルサイズ

数学いらずの医科統計学PART7 CHAPTER43で、サンプルサイズについて書いてある。 パッケージならpwr、デフォルトならpower.*.test関数群でできる。 最近はやりのびっぐでーたなら、メモリに乗ってPCで解析できる分だけとにかく持っているデータを使えばいい…

The Wolf at the Door

N Engl J Med. 2012 Jun 7;366(23):2216-21. 血球貪食症候群が悪性腫瘍、感染症、自己免疫疾患に二次性に発症する話。 バックナンバー

時系列解析

セミナーで時系列を学んだのだが、初音ミクの投稿数に適応してみる。 Kalman filterで過去から現在までたどって、平滑化で現在から過去へ遡ることでパラメータを最適化していく。 これを使うと欠損値補完もできる。 そして前回はやらなかった、今後の投稿数…