動的教材は作るのも楽しい Mathematicaのデモンストレーションプロジェクトにより、動的コンテンツの利用と提供が可能になります 利用する場合 こんなコンテンツをローカルにダウンロードして動かしてみることができます そのために必要なフリーソフトはこち…

以前の記事はこちら 基本、２通り バイナリで表示できる部分 ソースをダウンロードして、ソースファイルを探して読んでみる部分 C.Fortanの部分はソースファイルをダウンロードしないと読めない Ｒのソースはすぐに見られるものと、１回手間をかけないと読め…

不確かな部分が多いです･･･ Ｒの構成関数(言語での分類) Ｒコード Ｃコード FORTRANコード Ｒの構成関数(読み込み時期での分類) 基本骨格部分 ビルド時に読み込まれる関数 アドオンパッケージ部分 必要に応じて、後付けで読み込む部分 自作部分 呼び出し Ｒ…

距離とは次のように定義されている ２つの何か（Ａ，Ｂ）の間に定められ 非負(0以上)の値をとり その値は、ＡのＢに対する値もＢのＡに対する値も同じであり ＡとＢとが同じときには0をとり ３つの何か（Ａ，Ｂ，Ｃ）があるときには、Ａ−Ｂ，Ｂ−Ｃ，Ｃ−Ａの…

リンク先１ リンク先２

基本の使い方はこちらを参照。 texのファイル"test.tex"を次のように作ってみる \documentclass{jbook} \author{ 筆者 } \title{索引を作ってみよう} \usepackage{makeidx} \makeindex \begin{document} \maketitle わたしは、このページの北海道\index{ほっ…

作業環境をリムーバブルディスクに入れて、こっちのＰＣ、あっちのＰＣで作業をするためにもって歩くとする。 その作業には、Windowsのコマンドプロンプトの作業もあるので、リムーバブルディスク(たとえばE:)にcd したい。 C：￥Documents and Settings￥ユ…

完全グラフは、すべての頂点間に辺のあるグラフです。 エッジを行列表現して、segments()関数を使って、こんな風に。 DrawPerfectGraph<-function(d=5){ #d:次数 #正d角形の頂点座標 k<-1:d divk<-k/d divk<-divk*2*pi x<-cos(divk) y<-sin(divk) #プロット…

有限個の状態の離散的世代経過を状態間遷移関数でシミュレートすれば、遺伝的ドリフトを試せる。 #染色体N本中p本が着目アレル #それがk倍になった後に、N本を抜き出し #そこにi本の着目アレルがある確率 probDrift<-function(N,p,k,i){ kN<-k*N kp<-k*p kNp…

ヒトの核染色体(普通の染色体)の長さは、１番から22番、Ｘ、Ｙと塩基対単位で 247249719,242951149,199501827,191273063,180857866,170899992,158821424,146274826,140273252,135374737,134452384,132349534,114142980,106368585,100338915,88827254,7877474…

deFinetti<-function(g=c(0.09,0.42,0.49)){ x<-seq(from=-1/sqrt(3),to=1/sqrt(3),by=0.01) y<--1.5*x^2+1/6 plot(x,y,xlim=c(-1/sqrt(3),1/sqrt(3)),ylim=c(-1/3,2/3),type="l") abline(a=2/3,b=sqrt(3)) abline(a=2/3,b=-sqrt(3)) abline(h=-1/3) AA<-c(-…

初歩からの集団遺伝学作者: 安田徳一出版社/メーカー: 裳華房発売日: 2007/12メディア: 単行本 クリック: 10回この商品を含むブログ (4件) を見る オンリーワン・ゲノム―今こそ『遺伝と多様性』を知ろう作者: 鎌谷直之,早未恵理出版社/メーカー: 星の環会発…

hweExact<-function(g){ n<-sum(g) nA<-2*g[1]+g[2] na<-2*g[3]+g[2] evod<-g[2]%%2 obs<-(g[2]-evod)/2+1 maxAa<-min(nA,na)-evod Aa<-seq(from=evod,to=maxAa,by=2) AA<-(nA-Aa)/2 aa<-(na-Aa)/2 pr<-rep(0,length(Aa)) pr<-exp(n*lgamma(2+1)+lgamma(nA+1…

aからbを選ぶ場合の数は、です。 これに対応するＲの関数はchoose()で > choose(10,3) [1] 120 aからbを選ぶのは、aをb+c,ただしb+c=aに分けること。 Nをn1,n2,...,nk, na+n2+...+nk=Nに分けるのは そんな関数は。ただし、n1,...,nkは、全部与えても、nkを省…

20091222の日記にコメントをいただきました(こちら)。Sweaveを使ったＲの記録作りのことです。 Ｒを立ち上げます(WindowsXPで使っています) 「今日の作業」をします 作業ディレクトリが"C:\Documents and Settings\MyName\My Documents\RDiary\rd20100106"と…

モンティ・ホール問題(Wiki) テレビの問題は次の通り ３つの袋があって、１つはあたり、２つははずれという。 今、１つの袋を選んだあとに、残りの２つの袋のうち、はずれの袋が捨てられた。 はじめに選んだ袋と選ばなかった袋のうち捨てられなかった袋との…