2011-06-22から1日間の記事一覧

ちょっとメモ

1,2,...,Nの順列の逆対応を作るRソース InvSeq<-function(v){ M<-diag(rep(1,length(v))) M<-M[v,] M<-t(M) as.vector(M%*%(1:length(v))) } N<-5 s<-sample(1:N,N) s invs<-InvSeq(s) invs Im<-diag(rep(1,N)) Im[s,]%*%Im[invs,] こんな風にできるそうだ(…