2010-11-03から1日間の記事一覧

行列plot用変換関数

R

こんなん探せば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…

遺伝的浮動をplot3dで描く

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