R
// CodeIQ - 裏 RjpWiki で面白そうなことをやっていた。 締め切りは過ぎたようなので、自分でもやってみる。 かぞえるのではなく、計算する感じにしたので、違反かもしれぬ。 一応文字列で入力することにより相当な桁数まで対応。 f5(カウントする数字、"12…
// 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も計算できる。 もっと正確に計算することは出来るだろうが、 …
鬼のように忙しい怒涛の2週間であった。 やっと落ち着いてきたので、ちょっとRで遊ぶ。 時間のあるときに遊びでTDMのシミュレーションをしたい。 とりあえず、基本事項として波形を書いてみる。今後続きがあるかどうかは自分でもわからないが、%TAMなどの計…
STAP細胞が見つかるまでひたすら適当な細胞を生成し続けるプログラム - 驚異のアニヲタ社会復帰への道 理論的な解が存在するかどうかをまず確認 - 裏 RjpWiki いつもお勉強させて頂いているblogで面白そうなことをやっていた。 コードを参考にしながら、自分…
// 女子大生とペアプロするだけの簡単なお仕事です!|paizaオンラインハッカソンVol.2 で、面白そうな企画をやっていたようだ。(過去形) でも、R言語は選択肢に無い。 そもそも、Rがあったとしても、条件に合うようにインプットや アウトプットをどうすり…
// 子供と電車の中で10パズル(make10)をやっていた。 #2014.7.11追記 並び替えを許可する場合に10を作ることが出来る組合せ 並び替えを許可しない場合に10を作ることが出来る組合せ の出力を下の方に追加しました これをみれば、例えば1189の場合、並び…
// 探偵ファイルのロト6の攻略法を分析する 探偵ファイル〜スパイ日記〜/ロト6の攻略法!?/松本企業博士 探偵ファイルで、統計的な記事を見た。 気になったので分析してみた。 #離散値なのでbarplotを使った表に書き直しました #2014.7.11追記 結構こ…
// ピタゴラスの定理の整数解(原始ピタゴラス数) 子供がピタゴラスの定理の話をしていた。 いくつか整数解の組合せを知っているようだ。 aとbが1~1000までとすると、どのくらい整数解があるのか気になった。 wikipediaによると原始ピタゴラス数と言うらし…
// 多変量での検定 - 驚異のアニヲタ社会復帰への道 面白そうな解析をやっていた。 アウトカムが多変量の場合の解析である。 >O群とC群で画像のズレがx軸とy軸でとれるらしく ならば、原点が一番理想的なわけで、原点からの距離を指標にすれば 臨床的に意義…
// 恣意的なカテゴリー分け 連続変数で有意にならなかったからカテゴリ変数にして有意にしてる論文を読んだという話を聞いた これは面白い。つか、よくアクセプトされたな。 黒い統計家になった前提で、どの程度、この技が使えるのかやってみよう。 ### 黒い…
Rで自動的に名義変数らしきものを判別して、自動的にt検定とかfisher検定をする 自動的に、お手軽にt検定やFisher exact testなどの結果を返す関数を作った。 いちいちfactor設定するのも面倒なので、upper.col.no以下の水準数要素数ならば 自動的にfactor…
// 多重比較の調整法で、順序固定ゲートキーピング法と いうのがあると聞いた。 ここらへんも参考になった。 Analysis of clinical trials using sas 勉強用 isseing333 プライマリーの帰無仮説が棄却された場合にだけ、 次の検定に進む感じ。見聞きした2つ…
メタアナリシスじゃない信頼区間をforest plotで書いてみたい 同僚に、信頼区間を出すよう頼まれた。 比率なので、二項分布の正確な信頼区間binom.confintを使うことにする。 折角なので、計算したCIをforest plotで書いてみたいと思った。 しかし、低水準関…
SPSSでの標準化係数 ダミー変数の標準化偏回帰係数 - 裏 RjpWiki さんからたどって 重回帰分析で標準化回帰係数βを出力するR関数 - こにしき(言葉、日本社会、教育) さんの記事を見た。 SPSS20は使おうと思えば使えるので、やってみた。 まず、 #### y =…
// 孤独の7 Rで名作虫食い算:孤独の7を解いてみる。 7だけが明らかになっていて、ほかはすべて伏せられているという 名作虫食い算である。ぜひ、ノーヒントで手計算でチャレンジするのを 勧めたい。 孤独の7 ⇐問題文はここのhtmlを引用しました 問題 *7**…
人員が異動で減って、一時的に仕事量が倍になっている。 あと2週間の辛抱だ。 お世話になっている驚異のアニヲタ社会復帰への道さんが、卒業記念に本を書いたとのこと。サイトが怪しいが、500円とちょっとで安い。おすすめである。 購入はここ=> [Rで始…
対応の壊れたゲーム結果を復元する 対応の壊れたゲーム結果を復元する - ryamadaの遺伝学・遺伝統計学メモ ここで、面白そうな問題が出ていた。 ボーリングの第一ゲームは、各人のスコアが分かっている。 第二ゲーム以降は、どのスコアが誰のものなのか分か…
レアカード リボライブ! - 驚異のアニヲタ社会復帰への道でレアカードのコンプリートまでの回数をシミュレーションしていた。 勉強のために自分でもやってみる。 ## UR card rare <- 66 #レアカード枚数 rarep <- 0.01 #いずれかのレアカードが出る確率 card …
// 統計クイズ(あなたの直感は正しいだろうか) #t検定の問題 各群10人ずつのデータが有る。 a群のアウトカムは平均10.7 b群のアウトカムは平均15.4であった。 アウトカムは数が大きいほど良いとする。 あなたはb群の方が良いんじゃないかと思っている。…
// 多人数でババ抜きをすると何ターンかかるか(訂正) ばば抜きのシミュレーション - 裏 RjpWikiさんの記事のお陰で、間違いが見つかった。 訂正版を記載する。 今回も非常に勉強になった。blogをはじめて本当に良かった。 それにしても、一から作りなおし…
// 二種類の生存曲線を混ぜる。分子標的薬と細胞障害性抗癌剤。IPASS試験をRで科学する 二種類の生存曲線が混ざった状況だとどうなるか。分子標的薬と細胞障害性抗癌剤。 生存曲線の比較にはCOX回帰などが行われる。 期間を通じて比例ハザード性が保たれてい…
// がんもどき。予後の異なる二重ガンの一方を放置すると 前回の記事と同じだが、図がわかりにくい気がしてきた。これまでは、パワーポイントに貼ってから、色や凡例をつけてたけど、良い機会なので凡例や見出しの勉強をしよう。というわけで図に凡例をつけ…
// ババ抜きを多人数ですると何ターンかかるか(粘りの珍百景をシミュレートする) ####### 裏 RjpWikiさんのブログをみて、ざっと見たところ、 当方のコードに間違いを発見しました =>とりあえず、間違いは直しました。 訂正後の記事はこちら ##…
がんもどき。予後の異なる二重ガンの一方を放置すると Rによるシミュレーションを実施臨床に役立てた。 Aという癌があり、2年生存率が90%である。手術すれば100%治癒とする。 Bという癌があり、術後だが2年生存率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.犬の対決があった。 スポーツなんかはたまたまその時審査員を買収さえすれば勝ちっていうのもアリの世界みたいですし、矛盾な番組でも色々あったみたいで…