中枢神経に作用する薬剤をマウスに投与して、マウスの行動変化を測定した。 前シナプスから後シナプスへ、神経伝達物質を介して刺激が伝わることでマウスが行動する。ここでは、ドーパミンが行動を規定している（だろう）ものとして、実験を行った。 前シナ…

例えば a<- c(0,4,5,3,20,0,0,2,8) というベクトルがあったとして、この要素に0と等しいものがいくつあるか調べたい。 a == 0 [1] TRUE FALSE FALSE FALSE FALSE TRUE TRUE FALSE FALSE となるが、これを [1] TRUE FALSE 3 6 的に返してくれる関数はないでし…

ヒットを打つ（1種類）確率がpのとき、Nout.num回アウトになるまでに打席に立つバッターの加算確率的な何か。 Vnumは打者人数。 Nout.numは1イニング変わるためのアウトの数。 hitpは0.1、0.2、...0.9とヒットを打つ確率が変わる。ここではシングルヒットの…

NKがlistなるものを多用している。なんじゃこりゃぁぁぁあぁぁぁあぁぁ。 ﾋﾞｸﾝﾋﾞｸﾝｯ。 ということでlistの覚書。 ベクトルを格納するっぽい。スカラーも可。 v1<- 1:3 v2<- 4:6 v<- list(v1,v2) v [[1]] [1] 1 2 3 [[2]] [1] 4 5 6 ここからが慣れないとちょ…

こちらを実行した後で type<- c("p","l","b","c","o","h","s","S","n") titles<- c("p plot","l plot","b plot","c plot","o plot","h plot","s plot","S plot","n plot") par(mfrow=c(3,3)) for(t in 1:length(type)){ plot(t(B),xlab="",ylab="",type=type…

CardioLabというシミュレーションソフトを用いて、交感神経・副交感神経作用薬の効果を調べ、容量反応曲線およびSchild plotを描いた。 シミュレーションと聞いてwktkだったのに、いざやってみると 教員「じゃ、私の言うとおりに数値入力して」 …教員「じゃ…

峰は一峰性か? こちらでは一峰性っぽいが、アウトカウント4つのとき、微妙に峰が分離されているっぽい。 峰を定義しよう。 ある得点を取った試合数、がある。 その得点+1の得点を取った試合数、もある。 後者から前者を引き、これを増分とする。 増分が、あ…

野球で得点分布の山がどうなるのかをちょっと考えよう。 俺の嫁にどれだけ歌わせても嫌と言わないように、コンピュータもどれだけ計算させても嫌と言わないので、とりあえず、シミュレーションするに当たってこちら側で振っているパラメータをいじくってみま…

こんなん探せばpackageにあるだろたぶん… plot.converter<- # 変換しますよ〜 function(A){ # 変数は行列 x<- t(t(c(rep(1:nrow(A),ncol(A))))) # xを用意 y<- c() # yからの… for(i in 1:(T+1)){ y<- c(y,rep(i,N+1)) } y<- t(t(y)) # yを用意 convert.matr…

こちらでperspを使って描いているけど、plot3dでグリグリしたいので N<-50 #集団のサイズ k<-3 #残せる遺伝子の最大数 #k<-"infinity" #残せる遺伝子が無限大の場合 T<-10 #世代数 initial <- 5 #最初に存在する変異アレル数 #推移確率の行列を作る A<-matri…

オレ「大丈夫じゃない、問題だ。」 最近ブログ中のアニヲタ成分が急激に低下している… というわけで今期アニメヨスガノソラについて語るスレになりました。 ヨスガノソラは今期のエロゲ枠です。ストーリーはwikiを参照するとして、要は主人公ハーレムｳﾊﾊwww…

これを使って9イニングのスコア分布を描いた。 10万回試合したとき、得点は u #0点、1点、…最高で22点でした [1] 15620 17983 17440 14927 11344 8155 5529 3526 2300 1382 783 431 278 [14] 153 77 35 23 7 4 1 0 1 1となる。 今、打撃確率が同じチーム同士…

http://www.sciviews.org/Tinn-R/ Rのテキストエディタ

走者ベクトルを付け加えるやり方で進めることにした。 ここで問題となるのが、超低い確率だが3アウトまで行きつかねぇってときにうわ〜。 ということで制限設けて強制終了することで、終了しない事態、を防ぎましょう。 p<- seq(from=0,to=1,by=0.01) #アウ…

薬理学実習なるものがありました。 実習項目はいくつかありましたが、時間の都合上2つだけしかできませんでした。 交感神経・副交感神経作用薬による血圧調節 まあ何をするかざっくり言うと、アドレナリンとアセチルコリンを投与して血圧変化を追いましょう…

とりあえずこちらの3つめのやり方でかいてみた。 library(MCMCpack) batter <- rbind(c(0.71,0.20,0.06,0.01,0.02), # これをいろいろ振ってみる c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06…

g<- function(x){ # xは条件分岐の信号的なもの pattern<- 0 # 走者パターンとか score<- 0 # 得点 switch(x, # switch関数内のxは1以上の整数しか扱ってくれないらしい (pattern<- 1) && (score<- 0), # x=1なら、こうなる (pattern<- 2) && (score<- 1), #…

プログラムを動かすために、必要な定義。 イニング アウト 打者 打席 塁 打撃パターン 得点 ふたつチームを作って勝負させようと思っていたが、ひとつのチームで9イニング勝手にやったら何点取れるかをヒストグラムにして分布を考えるのも手だ。 やり方ひと…

text \begin{equation}A=\begin{pmatrix}a_{11}&a_{12}\\a_{21}&a_{22}\end{pmatrix}\end{equation}

グリコにかわって新しいネタをやろう。 野球ネタです。 4番松井は勝てるのか 野球では4番にホームランを打てるパワーヒッターをよくもってくるけど、戦術としてこれは最適なのか これをシミュレーションで考えてみよう。 つまり 4番パワーヒッターが一番勝ち…

インなんたらさんにくぎゅぅが出ていたﾊｧﾊｧ。神経科学の学生発表がありました。 どの班も手の込んだスライドを作っていて、試験に向けていい復習になった。 学生からの質問も例年に比べて活発で、教官方は非常に満足しておられたっぽい。 自分の班は注意の神…

『らき☆すた』と学ぶ 化学[理論編]が面白いほどわかる本作者: 松原隆志,美水かがみ出版社/メーカー: 中経出版発売日: 2010/10/07メディア: 単行本（ソフトカバー）購入: 12人 クリック: 1,063回この商品を含むブログ (8件) を見る高校生のときの衝撃は忘れな…

この前の先生からのメール。 本来は癌の病態や治療に興味があったから進学したわけであります（今となっては見る影もない…サーセンｗ）。 有名なのはRBのtwo hit theoryだが、まあぼや〜と思いついたのが癌集団中の癌幹細胞のポピュレーションが閾値を超えた…

ねんどろいど 初音ミク (ノンスケール ABS/PVC塗装済み可動フィギュア)出版社/メーカー: グッドスマイルカンパニー(GOOD SMILE COMPANY)発売日: 2009/05/11メディア: おもちゃ＆ホビー購入: 11人 クリック: 750回この商品を含むブログ (223件) を見る第5次再…

バグ取りして修正してみた。 N<-10 #ジャンケンする回数 library(gtools) #とりあえず全通り出す cmb<- combinations(5,N,repeats.allowed=T) num<- nrow(cmb) #A→1、B→2、C→3、D→4、E→5 と変換 ABCDE.count<- matrix(rep(0,5*nrow(cmb)),nc=5) for(i in 1:n…

とか毎回言ってんだろ？ おれもその一人ですけど何か？ 今期はヨスガノソラに期待。 原作探しちゃったよ… 今週のセミナーでは、バグをとろうということで そういうときどうしたらいいかで、小さいものを代入してみようと。 あとは言葉をよく知らなかったとい…

N 50でやってみて、極端なをとりあえず入れてみた。 #極端にグーを出す > probability [1] 6.156864e-22 > p [,1] [,2] [,3] [1,] 1 1.44129e-34 6.978894e-16 #極端にチョキを打す > probability [1] 1.207472e-21 > p [,1] [,2] [,3] [1,] 5.768531e-06 0.…

で、を固定したとき +6歩進む確率 出る回数A +3歩進む確率 出る回数B -3歩進む確率 出る回数C -6歩進む確率 出る回数D 0歩進む確率 出る回数E 一回のグリコでN回ジャンケンすると A〜EをN個重複を許して並べるやり方は通りあって、このうち ・・・（1） を満…

初音ミク/DVD「ミクの日感謝祭 39's Giving DayProject DIVA presents 初音ミク・ソロコンサート〜こんばんは、初音ミクです。〜」出版社/メーカー: ジェネオン・ユニバーサル発売日: 2010/09/01メディア: DVD購入: 4人 クリック: 84回この商品を含むブログ …

プログラミングセミナーでの内容や思いついたことをぐだ〜と書いておくことにしよう。本日の内容というか急に閲覧しはじめてわからない人へ。 あとでどこかにまとめを書く。一回のジャンケンで相手より前にいる歩数の期待値で考えていたが、無限回ジャンケン…