じゃんけんのシミュレーション

R

人員が異動で減って、一時的に仕事量が倍になっている。あと２週間の辛抱だ。お世話になっている驚異のアニヲタ社会復帰への道さんが、卒業記念に本を書いたとのこと。サイトが怪しいが、５００円とちょっとで安い。おすすめである。購入はここ＝> [Rで始…

2014-03-24

対応の壊れたゲーム結果を復元する

R

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

2014-03-13

レアカードガチャ

R

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

2014-03-11

統計クイズ１（あなたの直感は正しいだろうか）　t検定編

クイズ黒い統計家 R

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

2014-03-11

多人数でババ抜きをすると何ターンかかるか（訂正）

テレビ R

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

2014-03-10

二種類の生存曲線が混ざったら。分子標的薬と細胞障害性抗癌剤。IPASS試験のような状況をRで統計解析する

R

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

2014-03-10

がんもどき。予後の異なる二重ガンの一方を放置すると（図だけ改）

時系列概要 R

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

2014-03-08

ババ抜きを多人数ですると何ターンくらいかかるか（粘りの珍百景をシミュレートする）

テレビ R

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

2014-03-07

がんもどき。予後の異なる二重ガンの一方を放置すると

R

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

2014-03-06

RでのrJAVAのエラー

R

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

2014-03-06

MCMCglmmを無理やり使ってみる

R

先日の続き今度は混合効果モデルで推定する。すなわち各大学にランダム効果を仮定する。 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

国立が苦戦

R

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

2014-03-05

算数問題その３

R

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

2014-03-03

子供がやってた算数問題(図だけ改)

R

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

2014-03-02

20140301世界一受けたい授業、イヌvs.ネコを統計解析する

テレビ R

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

2014-02-26

算数の問題

R

子供がやってた算数問題子供が面白い算数の問題をやっていた。 a b ｜a-b｜というように上の2つの差の絶対値が下の数になる。このルールで、１～１０の整数を使ってa~jの表を完成する。 a b c d e f g h i j 問題はe＝６の場合の解であった。とりあえず…

2014-02-25

時系列データの描画

時系列 R

時系列データの描画同僚から、時系列データの解析を頼まれた。線形混合モデルを勧めたが、とりあえず図示することが大事である。どんなに良いモデル選択法を用いても、やはり人間の直感/見た目は重要である。しかーし、良い関数を知らないだけかもしれ…

2014-02-25

PSPPワロタ

ネタ

有名な統計ソフトSPSSのオープンソースのクローンが存在する。それがPSPPである。もう名前からして爆笑したのだが、昔のSPSSのファイルを見るのにイイネ！一応基本的な解析はできるし、無料なのがgood！ただ、Rがあるので自分で使うことはないかな。 SPS…

2014-02-17

ロケミツのサイコロの統計その３（４回連続して６がでるとイカサマ？）

ロケみつテレビ

// ６が連続４回でるとイカサマ？これは、前回の記事で使った関数を流用できる prenzoku <- function(m, q, r) { pren <- numeric(m - r + 1) pren[1] <- q^r pren[2:(r + 1)] <- (1 - q) * q^r for (i in (r + 2):(m - r + 1)) { pren[i] <- (1 - q) * q^r…