ちょっとメモ
- 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,]
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,]