打順最適化をして確率計算で出した得点分布と、対戦シミュレーションの結果がおかしい話。 （リンクでは総当たり戦しか書いていない。通りの対戦もすぐできる。） チーム1とチーム2が、1イニングである得点k点（）を取る確率がだとすると （実際にはまで得点…

こちらで打順をいじくって、1試合の得点分布がどうなるかを調べた。 結果、パワーヒッターの打順が(1,9)、(2,5,8)、(3,7)、(4,6)では同一の結果になることがわかった。 (1,9)はおそらく、1番打者が第1イニングで打席に立った後は9番打者の後ろに並ぶことにな…

Tigars<- rbind(c(0,0,0,0,1), c(0.6508972,0.2593801,0.05709625,0.004893964,0.027732463), c(0.6504065,0.2926829,0.04471545,0.010162602,0.002032520), c(0.6991304,0.2034783,0.05391304,0.010434783,0.033043478), c(0.7039007,0.1861702,0.02659574,…

打順を変える続き。 パワーヒッターの打順をチーム内でいろいろ変えたときの得点分布を作成する。 プログラムの概要。 パワーヒッターの効果が出やすいように、パワーヒッターは必ずHRを打つように設定する。 2010年の阪神のデータから、パワーヒッターが1番…

RMKUに混じって、CRANなどへの一般公開を念頭に置いた、関数とパッケージの作成について。

避難所 # パラメータ library(MCMCpack) library(gtools) library(expm) Nbase<- 4 Nout<- 3 Ninning<- 9 people<- 9 hitp<- 1-colMeans(team)[1] Vnum<- 100 cut<- 0.95 # 比較してみたいチーム Tigars<- rbind(c(0,0,0,0,1), c(0.6508972,0.2593801,0.0570…

誰が先頭打者になるかによって、1イニングで取る得点が変わりそうなのでその準備。 1,2,3,…,9という順列をひたすら 2,3,4,…,9,1 2,3,4,…,9,1,2 9,1,2,…,8 とベタに並び変えた仮想打順チームをarray関数を用いて作っておこうか。 よさげな関数がなかったので…

どのイニングに誰が先頭打者になりやすいかの続き。 ：nイニング目で打者その1、その2、…がそれぞれ、そのイニングの先頭打者になる確率ベクトル ：1イニングで打席に立ちうる打者人数の確率推移行列 として、 の形にしたい。 打者9人なら、1イニング目で1番…

野球の続き。チームの人数を複数にしたとき、誰が先頭打者になって、そこから何人が打席に立つかで、1イニングに取る得点が変わるだろう、ということを考えた。 この日は、全イニングの平均を出したので、1番バッターが多いのは当たり前だね。 シミュレーシ…

自分のテーマとして、野球の続き。 ・シミュレーションと確率過程での計算が、ずれているのか、同じ傾向（状況?経過?）を示しているのかどうか。 ・異なる打撃確率を持つpeople人で1チームが構成されている場合に、どう確率過程を計算するか。 ・次のテーマ …

冬コミで買った、嫁の切り抜きシール。 PCの天板に貼ってみた。 かわいすぎて生きるのが楽しい。 けど、3次元世界はつらい。

とりあえず応募してみた。

オレも使ってみよう。 薬物動態学の続き。 薬を投与するのにはいくつか経路がある。 経口投与では、腸管で吸収されるときに初回通過効果を受けて、投与した内のいくらか損失する。 静脈注射では、血管内に薬物を打ち込むので、初回通過効果は受けない。投与…

野球の続き。 今まで推移行列や漸化式を使って、打者にある打撃確率を与えたらk点を確率的に何点取る、というのをやってきた。 ただ、これでは1チーム中に複数の打者がいるときの得点が出せない。いまのところいいやり方が思いついていないし。 例えば9人チ…

風邪をひいたようだ。 看病募集中。 あったかいおかゆをﾌｰﾌｰしてｱｰﾝ。 薬物動態学、というのを学んでいる。 血中濃度をモニターしないといけない薬がいくつかあるとかで そういうとき、というか薬を処方したら問題になるのが、患者がきちんと飲んでいるか。 …

コミケが終わってからその足で鷲宮神社で年越しをしようと。 全く話が見えていないと思われる先生とNKのために解説すると、らき☆すたというアニメのキャラクターの実家がこの神社で、アニメのOPや正月回でこの神社が出たことから聖地として有名。 本当は歴史…