こちらのつづき 考えている野球のゲームはヒットの種類に制約がある すると可能な状況にも制約が出てくる 以下に可能な状況を表す整数（のベクトル）を発生させるプログラミング 漸化式で作る ヒットの種類が４種類で、5項間の漸化式 a,b,c,dは初期条件 uを…

ファイル名や拡張子を一気に変換したい ファイル 一括変換 で検索 こちら 拡張子を ".txt" から".tex"へ変換 for f in *.txt; do mv $f ${f%.txt}.tex; done マッチング、その後の処理をいろいろ変える 処理法方はこちらを参考 追記2010/10/31 こちらで touc…

こちらの続き xを状況を表す整数とする 状況xの時、何人の打者が出塁したかを表すベクトルVをつくる 並べて書くと一番右のベクトル 再帰的に作る 初期条件をいれておいて V<-c(0) 再帰条件を繰り返す V<-c(V,V+1) 終わりはないが 繰り返す回数は となるまで …

こちらのつづき こちらも参考にする ポアソン配置に時間軸をつけて考える 時間に対してアーラン分布、キョリに関してワイブル分布を出そうという魂胆 配置はN次元空間内に、体積V、時間t、密度h、点の数をxとして点の分布が で与えられる。 １ 時間について …

ポアソン分布と指数分布の関係について調べていてこちらを参考にさせていただいた プロットするときはこちらの記事 1 ポアソン分布 仮定 1.任意の時間と時間でそれぞれに発生する事象は巣顔に独立である 2.時間hに発生する事象P(1|h)について とする 3.時間h…

こちらの記事の野球っぽいゲーム 塁上にいる人のパターンと得点について 第N塁に人がいることを1、いないことを0と書くことにしてのとき、考えられる状況は通り これを整数の２進法表示と同一視して、上から順に整数０から7を対応させる 仮想的にもゆるす ヒ…

２変数鵜関数の3dプロット x,y,zの値を書き出さないといけないのだが 陽関数であっても(x,y,z)の対応が保たれるよう（x,y）の集合を作らないといけない library(rgl) #関数を定義 f <- function(x,y){ x^2-y^2 } n1<- 0:30/10 #x座標の集合 n2<- -20:20/10 #…

JIS形式のファイル testJIS.Rnw をシェル上でSweaveで処理する echo "library(tools); Sweave(\"testJIS.Rnw\")" | R --no-save --no-restore testJIS.tex というファイルがつくられる この testJIS.tes は "日本語(ISO 2022-JP)"で保存されている Rのエンコ…

いろいろな形式で保存されたtest.RnwをRのエディタでひらくときの文字化け シェル上で以下のコマンドを実行した open: test.Rnw #Rが起動し内部エディタで開く cat test.Rnw #shell上で内容を書き出す echo "library(tools); Sweave(\"test.Rnw\")" | R --no…

catでシェルウインドウ上に書き出したときに文字化けしたのはターミナルの環境設定によるものだった 環境設定 詳細 文字エンコーディング から読み込むファイルの保存形式に変更する

test.RnwというスクリプトをRとTeXで処理してdvi形式のファイルを作りたい そのあとdviからpdfに変換する 1. RのSweave()関数でLaTeXのスクリプトに変換する シェル上でRを起動してSweave()を実行するとtest.tex がはき出される R Sweave("test.Rnw") q() あ…

ポアソン分布 （ただし） Rではこちら dpois(k, lambda) rpois(N, lambda) N<-1000 l<-2 x<-rpois(N,l) plot(sort(x)) 指数分布 , Rではこちら N<-1000 lambda<-2 x<- rexp(N, lambda) plot(sort(x))

Rで行列を扱うときに便利な関数にapply()のファミリーがある(R-tips) apply(X,margin,関数) m <- matrix(rbinom(100,1,0.2),nrow=5) m apply(m,1,sum) marginの値 1:行にかんして 2:列にかんして

Coteditorのダウンロード こちらからCotEditor_1.0.1.dmgをダウンロード Thunderbird こちらからダウンロード、インストール TeXとTeXShop こちらに書いてある手順に従って、snow leopardに新規インストール こちらからMacTeX.mpkg.zip をダウンロード、イン…

二項係数 choose(n,r) u<-choose(n,0:n) plot(u) 二項分布 dbinom(r,n,p) n<-100 p<-1/6 v<-dbinom(0:n,n,p) plot(v)