#dplyr #magrittr のお勉強

R

##dplyr##magrittr http://d.hatena.ne.jp/teramonagi/20140814/1408024760http://d.hatena.ne.jp/teramonagi/20140403/1396529027http://rpubs.com/teramonagi/1112http://meme.biology.tohoku.ac.jp/students/iwasaki/rstats/tidyr.htmlhttp://meme.biolog…

2014-08-20

FIZZ BUZZ～デスマコロシアム

R

「第5回デスマコロシアム」問題のトーナメント結果発表です！──優勝者は…！ #デスマコロシアム｜CodeIQ MAGAZINE 標準出力に、余計な記号やスペースなどなしに 12fizz4buzzfizz78fizzbuzz11fizz1314FIZZBUZZ1617fizz19buzzfizz2223fizzbuzz26fizz2829FIZZBUZ…

2014-08-18

比率の差の信頼区間の手法色々

R

単なる覚書 ####比率の差の信頼区間多分Miettinen and Nurminenlibrary(PropCIs)diffscoreci(10, 100, 20, 150, 0.95) ###正確な比率の差の信頼区間時間がかかるので注意#library(ExactCIdiff)#BinomCI(100,150,10,20,0.95,CItype = "Two.sided", precisio…

2014-08-12

男女別の生命表から平均寿命を計算する

R

厚生労働省の所にデータが有ったので、やり直す。 ##男女別生命表#データフレーム作成厚生労働省H22国勢調査による完全生命表dat <- data.frame(numeric( 226 ))dat $ x = c( 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,2…

2014-08-12

生命表から平均寿命を計算する

R

mean survial time; AUC of survival curve r-statstics-fan Tuesday, August 12, 2014 1年生存率と5年生存率 - 驚異のアニヲタ社会復帰への道に触発されて、５年生存率ではなく平均寿命やMSTを計算したくなった。前回の記事を使って、簡単にコピペでデー…

2014-08-12

Rのデータフレームをブログ用にテキスト化してコピペで読み込めるようにする

R

1年生存率と5年生存率 - 驚異のアニヲタ社会復帰への道に触発されて、平均寿命とか計算したくなった。平均寿命といえば生存曲線下面積RMST(Restricted mean survival time)～比例ハザード性が仮定できない生存曲線での代替手法 - r-statistics-fanの日記 …

2014-08-09

Rで多倍長精度計算～とりあえず１００００の階乗を計算してみる

R

多倍長精度計算パッケージ Rmpfr - 裏 RjpWiki #追記2014.9.20 よく見ると自分の関数は小数点以下も出力されていて読みにくい。裏 RjpWikiさんの所で、より良い方法の記事が出ました。 <a href="http://blog.goo.ne.jp/r-de-r/e/0dee8bc64bfb50b39e796e207d50969e" data-mce-href="http://blog.goo.ne.jp/r-de-r/e/0dee8bc64bfb50b39e796e207d50969e">多倍長精度計算結果を見やすく出力 - 裏 RjpWiki</a> 多倍長精度計算…

2014-08-07

子供がやってた算数問題(その３くらい？)

R

□＝１～９（重複しない）ちょ、難しいんだけど。唯一解なのか分からない。総当りで確認する。 library(e1071)f1 <- function(){s <- permutations(9)ans <- s[,1] / (10*s[,2] + s[,3]) + s[,4] / (10*s[,5] + s[,6]) + s[,7] / (10*s[,8] + s[,9])return(…

2014-08-04

生存曲線下面積RMST(Restricted mean survival time)～比例ハザード性が仮定できない生存曲線での代替手法

R

生存曲線下面積RMST(Restricted mean survival time)というのを聞いた。論文の多くは田舎病院では入手できなかったが、下記は読めた。 Royston, P. & Parmar, M.K., 2013. Restricted mean survival time: an alternative to the hazard ratio for the desi…

2014-07-30

アンドロイドのパターンロック~GIGAZINE

R

Android OSのパターンロック認証はアルファベット3文字程度のセキュリティ - GIGAZINE 面白そうな記事があり、自分でも検証したくなった。パターン認証の各点を１２３４５６７８９と定義する。試行錯誤して条件を探ったが、＃４つ以上の点をつなぐ＃…

2014-07-21

子供がやっていた算数問題：シンプソンのパラドックス

R

最近の小学生はこんなこともやるのかと感心した問題。うろ覚えだが、こんな感じの問題であった A塾とB塾どちらに入ろうか検討している。A塾の生徒もB塾の生徒も同時に受けた統一模試の結果について問い合わせた。 A塾に聞いた所、男性も女性もB塾より平均点…

2014-07-14

ネタバレ

// コメントでsentaro0525さんに元論文を当てられてしまった。すごすぎる。ぼかしていたのに。そっから辿ると、【速報告知】渾身のガチ企画：『夏の統計的因果推論祭り』を開催します！（7月10日＠東大本郷）＊5/23登録締切りました＊ - Take a Risk：林…

2014-07-13

表からC-index　傾向スコア　操作変数

R

// C-index // C-index r-statistics-fan Sunday, July 13, 2014 某論文で、傾向スコアと操作変数を使った解析をやっていた。しかし、どちらもデータベースベースドの研究にありがちなこれは調整できてへんやろ感が激しい。傾向スコアの前提が信頼できる…

2014-07-09

７を数える（訂正）

R

0を数える場合の部分を書き直した。 f5 <- function(suuji, z){ if (nchar(z) == 1) { if (suuji == 0){ return(0) }else{ return(as.numeric(as.numeric(z) >= suuji))}} else{ x <- rev(as.numeric(unlist(strsplit(z, "")))) #分解してベクトルに keta <-…

2014-07-09

7を数える。検証

R

CodeIQ 7 を数えるの解答例 - 裏 RjpWiki 合っているか自信がないので、裏 RjpWikiさんのコードと自分のコードを比較して検証してみる。 simn <- 100000res <- matrix(0, nrow = simn, ncol = 2)for (i in 1:simn){ res[i,1] <- func2(i) res[i,2] <- as.in…

2014-07-07

CodeIQ 7を数える

R

// CodeIQ - 裏 RjpWiki で面白そうなことをやっていた。締め切りは過ぎたようなので、自分でもやってみる。かぞえるのではなく、計算する感じにしたので、違反かもしれぬ。一応文字列で入力することにより相当な桁数まで対応。 f5(カウントする数字、"12…

2014-07-05

RでTDM,AUC, C-max,Time above MIC

R

// RでTDM,AUC, C-max,Time above MIC // RでTDM,AUC, C-max,Time above MIC r-statistics-fan Saturday, July 05, 2014 一応rで薬剤の血中濃度を描くことができた。 C-maxやAUCやTime above MICも計算できる。もっと正確に計算することは出来るだろうが、 …