不思議な生物すたちゅーのX用すたンプです。
こいつ 非同次/確率的な強度λを持つ指数分布に従う確率変数τ(死亡/消滅時刻)の生成 - My Life as a Mock Quant の応用で、 時点1〜5までは無敵(不死) になるような確率分布からのサンプリングをやってみた。 具体的には lambda <- function(t){ifelse(t > 1 & t < 5, 0, 0.75)} というように、デフォルト強度(ポアソン分布の平均値に相当する量)を一時的に0と設定する。こうすると単位時間当たりの死亡率が0になるので、この間は無敵となるのだ。んで、この形からサンプリングすると、以下のように絶対死なない”無敵タイム付”の指数分布からのサンプリングになる。ちゃんと1〜5の間は確率密度が0になるので、この間は死なないよと。 コードは以下、コピペすれば同じ結果を再現できるぞい。 #int_0^t \lambda(u)du - \x
例えば、以下のような内容のR Markdownファイル(hoge.Rmd)を用意していたとする。 げっへっへ、本文だぜー!! ```{r, cache=TRUE} #Rのコード x <- 1:10 ``` また、本文だぜー! ここから、Rのコードだけを抽出するためにはknitrパッケージにあるpurl関数を用いると良い。documentation引数に応じて出力が変わる(チャンク・オプション、本文のコメント化など)ので、以下にその結果をメモる。 documentation=0 コマンド purl("hoge.Rmd", documentation=0) 結果 #Rのコード x <- 1:10 documentation=1(デフォルト) コマンド purl("hoge.Rmd", documentation=1) 結果 ## ----, cache=TRUE----------------
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く