r-statistics-fanの日記

統計好き人間の覚書のようなもの

2014-07-01から1ヶ月間の記事一覧

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

R

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

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

R

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

ネタバレ

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

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

R

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

7を数える(訂正)

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 <-…

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…

CodeIQ 7を数える

R

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

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も計算できる。 もっと正確に計算することは出来るだろうが、 …

TDM

R

鬼のように忙しい怒涛の2週間であった。 やっと落ち着いてきたので、ちょっとRで遊ぶ。 時間のあるときに遊びでTDMのシミュレーションをしたい。 とりあえず、基本事項として波形を書いてみる。今後続きがあるかどうかは自分でもわからないが、%TAMなどの計…