対応の壊れたゲーム結果を復元する 対応の壊れたゲーム結果を復元する - ryamadaの遺伝学・遺伝統計学メモ ここで、面白そうな問題が出ていた。 ボーリングの第一ゲームは、各人のスコアが分かっている。 第二ゲーム以降は、どのスコアが誰のものなのか分か…

レアカード リボライブ! - 驚異のアニヲタ社会復帰への道でレアカードのコンプリートまでの回数をシミュレーションしていた。 勉強のために自分でもやってみる。 ## UR card rare <- 66 #レアカード枚数 rarep <- 0.01 #いずれかのレアカードが出る確率 card …

// 統計クイズ（あなたの直感は正しいだろうか） ＃ｔ検定の問題 各群10人ずつのデータが有る。 a群のアウトカムは平均10.7 b群のアウトカムは平均15.4であった。 アウトカムは数が大きいほど良いとする。 あなたはｂ群の方が良いんじゃないかと思っている。…

// 多人数でババ抜きをすると何ターンかかるか（訂正） ばば抜きのシミュレーション - 裏 RjpWikiさんの記事のお陰で、間違いが見つかった。 訂正版を記載する。 今回も非常に勉強になった。blogをはじめて本当に良かった。 それにしても、一から作りなおし…

// 二種類の生存曲線を混ぜる。分子標的薬と細胞障害性抗癌剤。IPASS試験をRで科学する 二種類の生存曲線が混ざった状況だとどうなるか。分子標的薬と細胞障害性抗癌剤。 生存曲線の比較にはCOX回帰などが行われる。 期間を通じて比例ハザード性が保たれてい…

// がんもどき。予後の異なる二重ガンの一方を放置すると 前回の記事と同じだが、図がわかりにくい気がしてきた。これまでは、パワーポイントに貼ってから、色や凡例をつけてたけど、良い機会なので凡例や見出しの勉強をしよう。というわけで図に凡例をつけ…

// ババ抜きを多人数ですると何ターンかかるか（粘りの珍百景をシミュレートする） ＃＃＃＃＃＃＃ 裏 RjpWikiさんのブログをみて、ざっと見たところ、 当方のコードに間違いを発見しました ＝＞とりあえず、間違いは直しました。 訂正後の記事はこちら ＃＃…

がんもどき。予後の異なる二重ガンの一方を放置すると Rによるシミュレーションを実施臨床に役立てた。 Aという癌があり、２年生存率が90％である。手術すれば100％治癒とする。 Bという癌があり、術後だが２年生存率17.4％である。 Aという癌の治療方針につ…

library(XLConnect) としようとすると、エラーが出る。 指定されたモジュールが見つかりませんとか言われる。 でも、rJAVA.dllはある。おかしい。 仕方がないので、これまではEZRで読み込んでからRの 形式で保存たり、ExelでCSVに変換したりして解析していた…

先日の続き 今度は混合効果モデルで推定する。 すなわち各大学にランダム効果を仮定する。 a <- matrix(scan("clipboard", quiet=TRUE), byrow=TRUE, ncol=3) df <- data.frame(total = a[,1], yaba = a[,2], univ = a[,3]) df$id <- c(1:length(df$yaba)) f…

2014-03-05 - 驚異のアニヲタ社会復帰への道 ここで、国立が苦戦との分析がある。 glmで解析してみた。 やはり多分私立の書きまつがいと思う。 a <- matrix(scan("clipboard", quiet=TRUE), byrow=TRUE, ncol=3) df <- data.frame(total = a[,1], yaba = a[,…

リスペクトしているサイトの一つ裏 RjpWikiで、コードの 改良がなされていた。 なんと、 gtools の permutations() で爆速のコードを実現している。 permutations()でも、問題に特化して探索範囲を狭くすれば 実行時間の桁が変わることが分かった。 勉強にな…

子供がやってた算数問題(改) スキルのある諸先輩方は、きれいな図を書いている。 以前やった算数問題は、テキストだけでむりやり書いたが、イマイチわかりにくい。そこで、分かりやすくRで図を書いて再掲する。早い話Rで図を書く練習である。しかし、余白の…

// 20140301世界一受けたい授業、イヌvs.ネコ 今日は子どもと世界一受けたい授業を見ていた。 ネコvs.犬の対決があった。 スポーツなんかはたまたまその時審査員を買収さえすれば勝ちっていうのもアリの世界みたいですし、矛盾な番組でも色々あったみたいで…